addtoresetequationsection addtoresetlemmasection addtoresetremarksection addtoresettheoremsection

The shock formation and optimal regularities of the resulting shock curves for 1-D scalar conservation laws

Yin Huicheng^1,*,   Zhu Lu^2,
[Yin Huicheng<EMAIL_ADDRESS><EMAIL_ADDRESS>and Zhu Lu<EMAIL_ADDRESS>are
supported by the NSFC (No.11571177, No.11731007, No.12001162).]

1.  School of Mathematical Sciences and Institute of Mathematical Sciences,

Nanjing Normal University, Nanjing, 210023, China.

2. College of Science, Hohai University, Nanjing, 210098, China.

The study on the shock formation and the regularities of the resulting shock surfaces
for hyperbolic conservation laws is a basic problem in the nonlinear partial differential equations.
In this paper, we are concerned with the shock formation and the optimal regularities of the resulting shock curves
for the 1-D conservation law $\p_tu+\p_xf(u)=0$ with the smooth initial data $u(0,x)=u_0(x)$.
If $u_0(x)\in C^{1}(\Bbb R)$ and $f(u)\in C^2(\Bbb R)$,
it is well-known that the solution $u$ will blow up on the time $T^*=-\f{1}{\min{g'(x)}}$
when $\min{g'(x)}<0$ holds for $g(x)=f'(u_0(x))$. Let $x_0$ be a local minimum point of $g'(x)$ such that $g'(x_0)=\min{g'(x)}<0$ and
$g''(x_0)=0$, $g^{(3)}(x_0)>0$ (which is called the generic nondegenerate condition), then by Theorem 2 of [11],
a weak entropy  solution $u$ together with the shock curve $x=\vp(t)\in C^2[T^*, T^*+\ve)$ starting from
the blowup point $(T^*, x^*=x_0+g(x_0)T^*)$ can be locally constructed. When the generic nondegenerate condition is violated,
namely, when $x_0$ is a local minimum point of $g'(x)$  such that $g''(x_0)=g^{(3)}(x_0)=...=g^{(2k_0)}(x_0)=0$ but $g^{(2k_0+1)}(x_0)>0$
for some $k_0\in\Bbb N$ with $k_0\ge 2$; or $g^{(k)}(x_0)=0$ for any $k\in\Bbb N$ and $k\ge 2$, we will
study the shock formation and the optimal regularity of the shock curve $x=\vp(t)$, meanwhile,
some precise descriptions on the behaviors of $u$ near the  blowup point $(T^*, x^*)$
are given. Our main aims are to show that: around the blowup point, the shock really appears whether the initial data are
degenerate with finite orders or with infinite orders;  the optimal regularities
of the shock solution and the resulting shock curve  have the explicit relations with the degenerate degrees of the initial data.

Keywords: Shock formation, shock curve, entropy condition, Rankine-Hugoniot condition, hyperbolic conservation law

0.2 true cm

Mathematical Subject Classification 2000: 35L05, 35L72

§ MAIN RESULT

The study on the blowup and shock formation of smooth solutions to the
hyperbolic conservation laws is a basic problem in the nonlinear partial differential equations,
which has made much progress for the multi-dimensional cases in recent years (see [2]-[3],
[6]-[9], [12]-[15]). In the present paper, we are concerned with the
shock formation and the optimal regularities of the resulting shock curves
for the 1-D conservation law

\begin{equation}\label{0-1}
\left\{
\begin{aligned}
&\p_tu+\p_xf(u)=0,\ (t,x)\in\bR_+\times\bR,\\
&u(0,x)=u_0(x),\ x\in\bR,
\end{aligned}
\right.
\end{equation}

where $f(u)\in C^2(\Bbb R)$ and $u_0(x)\in C^1(\Bbb R)$.
It is well-known that the $C^1$ solution $u$ of (<ref>) will blow up at the time $T^*=-\f{1}{\min{g'(x)}}$
with $g(x)=f'(u_0(x))$ and $\min_{x\in\Bbb R}{g'(x)}<0$. If we further assume $g(x)\in L^{\infty}(\Bbb R)\cap C^p(\Bbb R)$
with $p\ge 4$, and let $x_0$ be a local minimum point of $g'(x)$ such that

\begin{align}\label{0-2}
g'(x_0)=\min_{x\in\bR}{g'(x)}<0,\quad g''(x_0)=0, \quad g^{(3)}(x_0)>0,
\end{align}

which is called the generic nondegenerate condition in [1], then by Theorem 2 of [11],
a weak entropy  solution $u$  of (<ref>) together with the shock curve $x=\vp(t)$ starting from
the blowup point $(T^*, x^*=x_0+g(x_0)T^*)$ can be locally obtained as follows:

\begin{align}\label{0-3}
&\vp(t)\in C^p(T^*, T^*+\ve)\cap C^{\f{p}{2}}[T^*, T^*+\ve).
\end{align}

(ii) In some part of the neighbourhood of $(T^*,x^*)$ near $x=\vp(t)$, for $t\ge T^*$ and $x\not=\vp(t)$,

\begin{equation}\label{0-4}
\left\{
\begin{aligned}
&|u(t,x)-u(T^*,x^*)|\le C((t-T^*)^3+(x-x^*)^2)^{\f16},\\
\end{aligned}
\right.
\end{equation}

When the generic nondegenerate condition (<ref>) is violated,
namely, if $x_0$ is a local minimum point of $g'(x)$ such that

\begin{equation}\label{0-5}
\left\{
\begin{aligned}
&g(x)\in L^{\infty}(\Bbb R)\cap C^{2k+2}(\Bbb R)\quad\text{for
$k\in\Bbb N$ with $k\ge 2$,}\\
&g'(x_0)=\min_{x\in\bR}{g'(x)}<0,\quad g''(x_0)=g^{(3)}(x_0)=...=g^{(2k)}(x_0)=0,\quad g^{(2k+1)}(x_0)>0,
\end{aligned}
\right.
\end{equation}

\begin{equation}\label{0-6}
\left\{
\begin{aligned}
&g(x)\in L^\infty(\Bbb R)\cap C^{\infty}(\Bbb R),\\
&g'(x_0)=\min_{x\in\bR}{g'(x)}<0,\quad g^{(k)}(x_0)=0\quad \text{for any $k\in\Bbb N$ and $k\ge 2$},
\end{aligned}
\right.
\end{equation}

we will
study the shock formation and the optimal regularity of the resulting shock $x=\vp(t)$ from the blowup point $(T^*, x^*)$,
meanwhile, some precise descriptions on the behaviors of  the solution  $u$ around the blowup point $(T^*, x^*)$
(rather than only some part near the shock curve) will be given.

Without loss of generality and for convenience, we set $x_0=0$ in (<ref>) and (<ref>). In addition, under condition (<ref>), near $x_0=0$ we assume

\begin{align}\label{0-7}
\end{align}

where $r(x)\in C^{2k+2}$ satisfies that $r^{(j)}(x)=O(x^{2k-j+2})$ for $0\le j\le 2k+2$;
under condition (<ref>), we choose a class of initial data

\begin{align}\label{0-8}
\end{align}

where $p>0$ and $r_0(x)\in C^\i\cap L^\i$ with

O(x^2-j),     j=0,1,2,

O(1),     j≥3

for $x$ near $0$.

Starting from the blowup point $(1,0)$ of (<ref>), let the formed shock curve  $\G$ be denoted by $x=\vp(t)$
if the shock really appears.
On the left hand side and right hand side of $\G$ for $t\ge 1$, the weak entropy solution  $u$ is represented
by $u_-$ and $u_+$ respectively (see Figure 1 below).

Shock formation

It follows from the Rankine-Hugoniot condition and entropy condition
on $\G$ that

\begin{align}\label{Rankine-Hugoniot Condition}
\end{align}

where $[u](t,\vp(t))=u_+(t,\vp(t))-u_-(t,\vp(t))$ is the jump of $u$ across $\G$, and

\begin{align}\label{entropy condition}
\end{align}

Our main results are

Under assumption (<ref>), there exists a unique solution  $u\in C^1((0,1)\times\bR)\cap C([0,1]\times\bR)$ to
problem (<ref>) together with (<ref>)-(<ref>) for $t\ge 1$.

(1) $\vp(t)\in C^\f{k+1}{k}[1,1+\ve)$ and
$u\in C^1((1,1+\ve)\times\bR)\setminus \{x=\vp(t)\})$ for some $\ve>0$.

(2) near the blowup point $(1,0)$, the behaviors of $u$ and its derivatives are as follows
|u(t,x)-u(1,0)|   =   O(|t-1|^12k+|x|^12k+1),

|_t u(t,x)|   =   O((|t-1|^12k+|x|^12k+1)^-(2k-1)),

|_x u(t,x)|   =   O((|t-1|^12k+|x|^12k+1)^-2k).

Under assumption (<ref>), there exists a unique solution  $u\in C^1((0,1)\times\bR)\cap C([0,1]\times\bR)$ to
problem (<ref>) together with (<ref>)-(<ref>) for $t\ge 1$.

(i) $\vp(t)\in C^1[1,1+\ve]$ and
$u\in C^1(((1,1+\ve)\times\bR)\setminus \{x=\vp(t)\})$ for some $\ve>0$. In addition, $\vp(t)=O((t-1)|\ln(t-1)|^{-\f{2}{p}})$
near $t=1$ and for $t>1$.

(ii) near the blowup point $(1,0)$, the behaviors of $u$ and its derivatives are as follows
|u(t,x)-u(1,0)|   =   O(|ln|t-1||^-1p+|ln|x||^-1p),

|_t u(t,x)|   =   O(|t-1|^-1|ln|t-1||^-1-1p+|x|^-1|ln|x||^-1-1p),

|_x u(t,x)|   =   O(|t-1|^-1|ln|t-1||^-1+|x|^-1|ln|x||^-1).

If we take $k=1$, then Theorem 1.1 coincides with the result in Theorem 2 of [11].
In addition, the author in  [11] only shows the behaviors of $u$ in some part of the
neighbourhood of the blowup point  $(1,0)$, which corresponds to
the smallness of the variable $|\la|=\f{|x|}{(t-1)^{\f32}}$ for $t>1$. This can be referred to
the proof of Lemma 2.1 in [11], where $|\la|$ is required to be small.

The regularities of $\vp(t)$ in Theorem 1.1 and  Theorem 1.2 are optimal.
One can see Remark 3.1 and Remark 4.1 below.

Under the generic nondegenerate assumption of the initial data, for the 1-D $2\times 2$ $p-$ system of
polytropic gases, the authors in [10]-[11] and [5] obtain the formation and  construction
of the shock wave starting from the blowup point under some variant conditions; for the 1-D $3\times 3$
strictly hyperbolic conservation laws with the small initial data or the 3-D full compressible Euler equations with
symmetric structure and small perturbation,
the authors in [4], [16] and [7] also get the formation and  construction
of the resulting shock waves, respectively.

In order to prove Theorem 1.1-1.2,
our focus is to solve the singular and nonlinear ordinary differential equation (<ref>)
as in [11]. Note that the equation (<ref>) is equivalent to
$\vp'(t)=G(t,\vp(t))\triangleq \int_0^1f'(\th u_+(t,\vp(t))+(1-\th)u_-(t,\vp(t)))d\th$, where
the function $G(t,\vp)$ is not Lipschtzian with respect to variable $\vp$ since the first order derivative of
$u_{\pm}(t,x)$ with respect to the variable $x$ admits the strong singularities (see (1.13) and (1.17)). To get the uniqueness and
regularity of $(\vp(t), u_{\pm}(t,x))$, we require to carefully analyze the behavior and regularity
of solution $u$ near the blowup point $(1,0)$. Due to the more degenerate conditions (<ref>)
and (<ref>), we shall introduce some different transformations of $(t,x)$ from that in [11]
(for examples, see (<ref>),  (<ref>), (<ref>) and so on).
By involved computation, the behaviors of solution $u$ around the point $(1,0)$ are
derived and then the optimal regularities of $\vp(t)$ are also established. From our results,
we have known two basic facts for problem (<ref>): (1) Around the blowup point, the shock really appears whether the initial data are
degenerate with finite orders or with infinite orders. (2)  The optimal regularities
of the shock solution and the resulting shock curve  have explicit relations with the degenerate degrees of the initial data.

Our paper is organized as follows. In Section 2, we give some basic analysis on the characteristics envelope
of equation (<ref>) near $(1,0)$, meanwhile, the detailed behaviors of the characteristics near  $(1,0)$
are established. The proofs of Theorem 1.1 and Theorem 1.2 are given in Section 3 and Section 4 respectively.

§ SOME PRELIMINARY

For problem (<ref>), we define the characteristics $x=x(t,y)$ starting from the point
$(0,y)$ as follows

d x(t,y)d t=f'(u(t,x(t,y))),

Then along this characteristics we have

This means that the characteristics $x(t,y)$ is straight and

For any fixed $t>0$, in order to solve $y=y(t,x)$ in (<ref>) such that the solution $u$ in (<ref>)
can be obtained, it is necessary to let

xy(t,y)=1+tg'(y)>0.By assumption (<ref>) or (<ref>), we have that near $x=0$,

(i) for $0\le t< 1$, $\f{\p x}{\p y}(t,y)>0$;

(ii) $\f{\p x}{\p y}(1,y)\ge 0$, and only at $y=0$, $\f{\p x}{\p y}(1,y)=0$.

Thus for $0\le t\le 1$, one can get a function
$y=y(t,x)$ satisfying (<ref>) such that the solution to (<ref>) is

On the other hand, one can compute that for $0\le t<1$,

This means that as $(t,x)$ tends to $(1-,0)$, then $y(t,x)\rightarrow 0$ and
$|\p_x y(t,x)|\rightarrow +\i$.

Let $\ve>0$ be a sufficiently small constant.
Under assumption (<ref>), it is easy to check that for $1<t<1+\ve$ and $y$ near $0$, there exist two roots of $\p_y x(t,y)=0$
with respect to the variable $y$, which are denoted by
$\eta_-(t)$ and $\eta_+(t)$ with $\eta_-(t)<\eta_+(t)$. Set $x_\pm(t)=x(t,\eta_\pm(t))$, we then have

* for $x<x_+(t)$ ($x>x_-(t)$ resp.) and equality (<ref>), there exists a unique root denoted by $y_-(t,x)$ ($y_+(t,x)$ resp.).

* for $x=x_+(t)$ ($x=x_-(t)$ resp.)  and equality (<ref>), there exist two roots denoted by $y_-(t,x)<\eta_+(t)$ ($\eta_-(t)<y_+(t,x)$ resp.).

* for $x_+(t)<x<x_-(t)$  and equality (<ref>), there exist three roots denoted by $y_-(t,x)<y_0(t,x)<y_+(t,x)$.

Ø_-   =   {(t,x):1< t<1+, x<x_-(t)}

Ø_+   =   {(t,x):1< t<1+, x>x_+(t)}

Ø_0   =   {(t,x):1< t<1+, x_+(t)<x<x_-(t)}.

Under (<ref>), we derive some properties of $\eta_\pm(t)$ and $x_\pm(t)$ near the blowup point $(1, 0)$.

There exists an $\ve>0$ sufficiently small such that

(1) $\eta_\pm(t)\in C^{2k+1}(1,1+\ve)\cap C^\f{1}{2k}[1,1+\ve)$ admit the following expansion

(2) $x_\pm(t)=x(t,\eta_\pm(t))\in C^{2k+1}(1,1+\ve)\cap C^\f{2k+1}{2k}[1,1+\ve)$ are the envelopes of the characteristic lines
(<ref>) which form a cusp at $(1,0)$, meanwhile,

(1) Note that $\eta_\pm(t)$ are the solutions of

1+tg'(y)=-(t-1)+(2k+1)ty^2k+t r'(y)=0.
This immediately yields $\eta_\pm\in C^{2k+1}(1,1+\ve]$ by the implicit function theorem.
For $t\to 1+$, set

s=(t-1)^12k, z=ys.
Then (<ref>) becomes

Since $r'(sz)=O(s^{2k+1})$ for $s$ near $0$,
$F(0,z^0_\pm)=0$ holds for $z^0_\pm=\pm(2k+1)^{-\f{1}{2k}}$. By direct computation, we have that
_s F(s,z)   =   2k(2k+1)s^2k-1z^2k-2k s^-2k-1r'(sz)+(s^-2k+1)z r”(sz),

_z F(s,z)   =   (1+s^2k)[2k(2k+1)z^2k-1+s^-2k+1r”(sz)].
Together with $r(sz)=\f{g^{(2k+2)}(0)}{(2k+2)!}(sz)^{2k+2}+o(s^{2k+2})$, this yields
_s F(0,z^0_±)   =   ±g^(2k+2)(0)(2k+1)!(2k+1)^-2k+12k,

_z F(0,z^0_±)   =   ±2k(2k+1)^12k≠0.

By the implicit function theorem, for small $\ve>0$ there exist

such that $F(s,z_\pm(s))=0$ and

Therefore, (<ref>) is shown and then $\eta_\pm(t)\in C^{\f{1}{2k}}[1,1+\ve]$.

(2) By (<ref>) and (<ref>), we have

x_\pm(t)=x(t,\eta_\pm(t))=-(t-1)\eta_\pm(t)+t\eta^{2k+1}_\pm(t)+t r(\eta_\pm(t)).

Together with (<ref>), this yields $x_\pm(t)\in C^{2k+1}(1,1+\ve]\cap C^\f{2k+1}{2k}[1,1+\ve]$
and the expansion (<ref>). In addition, due to $\f{\p}{\p y}x(t,\eta_\pm(t))=0$, then for $t\in[1,1+\ve]$,

This means that the tangent direction of $x=x_\pm(t)$ coincides with the characteristic speed of (<ref>) at $(t,x_\pm(t))$.

Consequently, the proof of (2) is finished.

Under (<ref>)-(<ref>), we have

For $\eta_\pm(t),\ x_\pm(t),\ y_\pm(t,x)$ and $y_0(t,x)$, we can deduce the following properties
for small $\ve>0$:

(1) $\eta_\pm(t)\in C^\i(1,1+\ve]\cap C[1,1+\ve]$ with $\eta_\pm(\tau)=\pm|\ln(t-1)|^{-\f{1}{p}}+O(\f{\ln|\ln(t-1))|}{|\ln(t-1)|})$.

(2) $x_\pm(t)=x(t,\eta_\pm(t))\in C^\i(1,1+\ve]\cap C^{1}[1,1+\ve]$ are the envelopes of the characteristic lines
and form a cusp at $(1,0)$. Moreover we have the expansion $x_\pm(t)=\mp(t-1)|\ln(t-1)|^{-\f{1}{p}}+O(\f{(t-1)\ln|\ln(t-1)|}{|\ln(t-1)|})$.

(3) For any $t\in (1,1+\ve]$, $y_-(t,\cdot)$ is an increasing function from $(-\i,x_-(t)]$ onto $(-\i,\eta_-(t)]$; $y_0(t,\cdot)$
is a decreasing function from $[x_+(t),x_-(t)]$ onto $[\eta_-(t),\eta_+(t)]$; $y_+(t,\cdot)$ is an increasing function from
$[x_+(t),+\i)$ onto $[\eta_+(t),+\i)$. Moreover, $y_m(t,x)\in C^\i(\O_m)\cap C(\bar\O_m)$, where $m=-,+,0$.

(1) Set $\tau=t-1$ for $t\geq1$. Note that $\eta_\pm(t)$ are defined for small $\tau>0$ and  are the solutions of the following equation

1+tg'(y)=(-τ+1|y|^p e^-|y|^-p)+τ|y|^p e^-|y|^-p+(τ+1)e^-|y|^-pr_1(y)=0,
where $r_1(y)=\f{1}{p}+\f{y r_0(y)}{|y|^{p+2}}+r_0'(y)=O(|y|^{
\min\{-p+1,0\}})$. Denote $\o=|\ln\tau|^{-1}$ and $z=|y|^{-p}-(\o^{-1}-\ln \o)$, then (<ref>) becomes

Obviously, $F(0,0)=0$. In addition, by direct computation, we have that for small $|y|$,

\f{\p y}{\p z}=-\f{y|y|^{p}}{p}.

Thus one can check that
yz)).This yields $\f{\p F}{\p z}(0,0)=-1$ by $|y|\lesssim \o^\f{1}{p}$. Since $F(\tau,z)$ is continuous and has the continuous partial derivative $\f{\p F}{\p z}$ near $(0,0)$, it follows from the implicit function theorem that there exists a continuous function $z=z(\tau)$ near $z=0$ to satisfy $F(\tau,z)=0$. This deduces  $\eta_\pm(t)=\pm(|\ln\tau|+\ln|\ln\tau|+z(\tau))^{-\f{1}{p}}=\pm|\ln\tau|^{-\f{1}{p}}+O(\f{\ln|\ln\tau|}{|\ln\tau|})$. On the other hand, for $\tau>0$,
we have $\f{\p x}{\p y}(t,y)\neq 0$ and $g'\in C^\i$. Then by the implicit function theorem, $\eta_\pm(t)\in C^\i(1,1+\ve]$ hold.

(2) By (<ref>) and (<ref>), we have that for small $\tau>0$,

which derives $x_\pm(t)\in C^\i(1,1+\ve]\cap C[1,1+\ve]$. On the other hand, it holds that for $t\in (1,1+\ve]$,
ddtx_±(t)   =   tx(t,η_±(t))+yx(t,η_±(t))ddtη_±(t)

=   g(η_±(t)),which means that the tangent direction of $x=x_\pm(t)$ is same as the characteristic speed of (<ref>) at
$(t,x_\pm(t))$. In addition, at the point $(1,0)$, one has
x'_±(1)   =   lim_t→1_+x_±(t,η_±(t))-0t-1

=   lim_t→1_+[-η_±(t)+t(η_±(t)p+o(η_±(t)))e^O(ln|lnτ||lnτ|)]

=   0

=   g(0).Hence we finish the proof of (2).

(3) For any fixed $t\in (1, 1+\ve]$, due to

\f{\p}{\p y}x(t,y)\left\{
\begin{array}{ll}
>0, & \text{ for }y\in(-\i,\eta_-(t))\cup(\eta_+(t),+\i),\\
=0, & \text{ for }y=\eta_\pm(t),\\
<0, & \text{ for }y\in(\eta_-(t),\eta_+(t)),
\end{array}
\right.

then by the inverse function theorem, $y_m(t,\cdot)$ with $m=-,+,0$ are well defined and satisfy the corresponding monotonicity.
Moreover, $y_+(t,\cdot)\in C^\i(x_+(t),+\i)\cap C[x_+(t),+\i)$, $y_0(t,\cdot)\in C^\i(x_+(t),x_-(t))\cap C[x_+(t),x_-(t)]$
and $y_-(t,\cdot)\in C^\i(-\i,x_-(t))\cap C(-\i,x_-(t)]$.

On the other hand, because of $\p_y x(t,x)\neq 0$ for $(t,x)\notin \{x=x_\pm(t)\}$, thus
it follows from the implicit function theorem
that $y_m(t,x)\in C^\i(\O_m)$, $m=-,+,0$. For the continuity of $y_m(t,x)$ in $\bar \O_m$, we take $y_+(t,x)$ as an example.
By $x_+(t)\in C^1([1,1+\ve])$, we then get

|y_+(\bar t,x_+(\bar t))-y_+(t,x)|=|y_+(\bar t,x_+(\bar t))-y_+(t,x_+(t))|+|y_+(t,x_+(t))-y_+(t,x)|\rightarrow 0

as $(t,x)\rightarrow(\bar t,x_+(\bar t))$ for $\bar t\in[1,1+\ve]$ and $(t,x)\in \O_+$. Thus $y_+(t,x)\in C^\i(\O_+)\cap C(\bar\O_+)$

To study the formation of shock wave and the regularity of the resulting shock $x=\vp(t)$ to equation (<ref>),
it is required to study the properties of $y_\pm(t,x)$ for $(t,x)$ in the cusp domain $\O_0$. Under assumption
(<ref>), motivated by [11], we take the following change of the variables

τ=t-1, s=τ^12k, μ=ys, =xs^2k+1,
and will establish the behavior of $y_\pm(t,x)$ near $(1,0)$ in some sub-domain of $\O_0$.

For small $\ve>0$, under  assumption
(<ref>), there exists some constant $\dl>0$ such that for $(s,\la)\in\{0\le s\le \ve,\ |\la|\le \dl\}$,
$(s,\la)\rightarrow s^jy_\pm(t,x)$ are of $C^{j+2}$ for $j=-1,0,1,\ldots,2k$ and $y_{\pm}(t,x)$ admit the
following expansions
y_+(t,x)   =   s(1+2k-g^(2k+2)(0)2k(2k+2)!s)+O(s^3+s^2),

y_-(t,x)   =   s(-1+2k-g^(2k+2)(0)2k(2k+2)!s)+O(s^3+s^2).

h(y)≜r(y)y^2k+1=∫_0^1 (1-þ)^2k(2k)!g^(2k+1)(þy)dþ-1.
Then $h(y)\in C^{p-2k-1}$ and $h(0)=0$. Furthermore,

y h'(y)=y(2k)!∫_0^1(1-þ)^2kþg^(2k+2)(þy)dþ=-1(2k)!∫_0^1g^(2k+1)(þy)(1-þ)^2k-1[1-(2k+1)þ] dþ.
This derives $y h'(y)\in C^1$. Similarly, $y^j h^{(j)}(y)\in C^{1}$ holds for $j=2,\ldots,2k+1$.
Divided by $s^{2k}$, (<ref>) becomes

For $s=\la=0$, by $G(0,0,\mu)=-\mu+\mu^{2k+1}=-\mu(1-\mu^{2k})=0$
we get the roots $\mu^0_\pm=\pm1$ and $\mu^0_c=0$. Note that

By the implicit function theorem, there exist functions $\mu=\mu_\pm(s,\la)\in C^{1}$ near $(s,\la)=(0,0)$ such that

G(s,,μ_±(s,))=0, μ_±(0,0)=±1,
and then $s^{-1}y_\pm\in C^{1}$. On the other hand, due to
_s G(s,,μ)   =   2ks^2k-1μ^2k+2ks^2k-1μ^2k+1h(sμ)+(1+s^2k)μ^2k+2h'(sμ),

_G(s,,μ)   =   -1,
_s G(0,0,±1)   =   h'(0)

_G(0,0,±1)   =   -1.
It follows from (<ref>), (<ref>) and (<ref>) that

_s μ_±(0,0)=-g^(2k+2)(0)2k(2k+2)!,   _μ_±(0,0)=12k.
Consequently, the expansions (<ref>) and (<ref>) are shown.

We next prove $(s,\la)\rightarrow y_\pm\in C^{2}$. By

s_s G(s,,μ_±)+_μG(s,,μ_±)(s_sμ_±(s,))=0

s_s G(s,,μ_±(s,))
we then have $s\p_s\mu_\pm(s,\la)\in C^{1}$.
Thus from $\p_s y_\pm=\mu_\pm+s\p_s\mu_\pm$, one can see that
$(s,\la)\to y_\pm$ is of $C^{1}$.
In addition, by $\p_\la y_\pm=s\p_\la\mu_\pm(s,\la)$ and similar computation,
we can get $y_\pm\in C^{2}$ with respect to $s$ and $\la$.

Note that

s^j^j_s G(s,,μ_±(s,))   =   (s,,s_sμ_±,…,s^j-1^j-1_sμ_±,

+_μG(s,,μ_±)(s^j^j_sμ_±(s,)), j=2,3,…,2k+1,where $\vG$ is a polynomial with respect to its arguments. Then
$s^j\p^j_s\mu_\pm(s,\la)\in C^{1}$ for $j=2,3,\ldots,2k+1$ by induction.
Similarly, $s^j\p^m_
\la\p^{j-m}_s \mu_\pm(s,\la)\in C^1$ for $1\le m\le j\le 2k+1$. Consequently,  the proof
of $(s,\la)\rightarrow s^jy_\pm(t,x)\in C^{j+2}$ for $j=-1,0,1,\ldots,2k$ is completed.

Under assumption (<ref>), we now study the asymptotic behavior of $y_\pm(t,x)$ near $(1,0)$.
In this case, we take the following change of the variables

τ=t-1, s=|lnτ|^-1p, =xsτ, μ=ys.
Then we obtain

Under assumption (<ref>), there exist some small constants $\ve$, $\dl>0$ such that for $(t,x)\in\O_0\triangleq\{1\le t\le 1+\ve,\ -\dl s\tau< x< \dl s\tau\}$, $y_\pm\in C^{1+p}$ hold for the variables $s$ and $\la$. Furthermore,
y_+(t,x)   =   s(1+lnpps^p+s^pp)+O(s^min{p+2,2p+1}+s^p+1||^2),

y_-(t,x)   =   s(-1-lnpps^p+s^pp)+O(s^min{p+2,2p+1}+s^p+1||^2).

By divided by $s\tau$, (<ref>) can be written as follows

Below we assume $\mu\neq0$. Without loss of generality, one can assume
$\mu>0$ (corresponding to the case of $y_+(t,x)$). We divide the proof of Lemma <ref>
into the two cases of $p\ge1$ and $p\in(0,1)$.

Case 1. $p\geq1$.

Set $\zeta=s^{-p}(1-\mu^{-p})-\ln p$ and $\mu=\left(1-s^p(\zeta+\ln p)\right)^{-\f{1}{p}}$. Then (<ref>) becomes

Obviously $F_1(0,0,0)=0$ and $\ds\lim_{s\to0+,\zeta\to0}\mu=1$. Note that
_sμ=s^p-1(ζ+lnp)(1-s^p(ζ+lnp))^-1p-1,  _ζμ=s^pp(1-s^p(ζ+lnp))^-1p-1.Thanks to $p\ge1$, we have that $\p_s\mu$ is bounded. On the other hand,
_s F_1   =   ps^-p-1μe^ζ-s^-p+(-1+e^ζ+e^ζ-s^-p)_sμ

_F_1   =   -1,

_ζF_1   =   μ(e^ζ+e^ζ-s^-p)
+(-1+e^ζ+e^ζ-s^-p)_ζμ+pe^ζ(1+e^-s^-p)s(r_0(sμ)+s r'_0(sμ)_ζμ).This derives $F_1\in C^1$ and
_s F_1(0,0,0)=p2r”_0(0),  _F_1(0,0,0)=-1,  _ζF_1(0,0,0)=1.

Thus by the implicit function theorem, one can obtain that there exists a unique function $\zeta(s,\la)\in C^1$ satisfying $F_1(s,\la,\zeta(s,\la))=0$ and admitting the following expansion

At this time, we get

μ(s,)   =   (1-s^p(ζ+lnp))^-1p

=   (1-s^p(lnp+p2r”_0(0)s++O(s^2+^2))^-1p

=   1+lnpps^p+s^pp+O(s^p+1+s^p||^2)
and $s^{l}\mu\in C^{l+p}$, $l=0,1$.

If we consider the case of $\mu<0$, then by the same method, one can obtain the expansion

Case 2. $0<p<1$.

Set $\o=s^{p}$, $\zeta=\o^{-1}(1-\mu^{-p})-\ln p$ and $\mu=\left(1-\o(\zeta+\ln p)\right)^{-\f{1}{p}}$. Then (<ref>) becomes

Obviously $F_2(0,0,0)=0$ and $\ds\lim_{\o\to0+,\ \zeta\to0}\mu=1$. Note that
_øμ=ζ+lnpp(1-ø(ζ+lnp))^-1p-1,  _ζμ=øp(1-ø(ζ+lnp))^-1p-1.On the other hand,
_øF_2   =   ø^-2μe^ζ-ø^-1+(-1+e^ζ+e^ζ-ø^-1)_øμ

_F_2   =   -1,

_ζF_2   =   μ(e^ζ+e^ζ-ø^-1)
+(-1+e^ζ+e^ζ-ø^-1)_ζμ+pe^ζ(1+e^-ø^-1)ø^1p(r_0(ø^1pμ)+ø^1p r'_0(ø^1pμ)_ζμ).Then $F_2\in C^1$ and
_øF_2(0,0,0)=0,  _F_2(0,0,0)=-1,  _ζF_2(0,0,0)=1.

Thus by the implicit function theorem, one can obtain that there exists a unique function $\zeta(\o,\la)\in C^1$ satisfying $F_2(\o,\la,\zeta(\o,\la))=0$ and admitting such an expansion

This yields

μ(s,)   =   (1-s^p(ζ(s^p,)+lnp))^-1p

=   (1-s^p(lnp++O(s^2p+^2))^-1p

=   1+lnpps^p+s^pp+O(s^2p+s^p||^2)
and $s^{l}\mu\in C^{l+p}$ for $l=0,1$.

If we consider the case of $\mu<0$, then by the same method, one can obtain the expansion

Consequently, we complete the proof of (<ref>) and (<ref>).

§ PROOF OF THEOREM <REF>

We now construct the shock curve $x=\vp(t)$ of (<ref>) in $\O_0$.
It follows from the Rankine-Hugoniot condition that

This, together with the mean-value theorem, yields

which means that the entropy condition on $x=\vp(t)$ holds.

f(u_0(x))-f(u_0(y))u_0(x)-u_0(y),      if x≠y,

g(x),      if x=y.

Under assumption (<ref>), it is easy to verify $a(x,y)\in C^{2k+2}(\bR^2)$ with $a(x,y)=-\f{1}{2}(x+y)+b(x,y)$,
where $b(x,y)=O(x^2+y^2)\in C^{2k+2}$.

Under assumption (<ref>), for (<ref>) and small $\ve>0$, there exists a solution $x=\vp(t)$ on $[1,1+\ve)$
to satisfy

(1) $\vp(t)$ is a $C^{2}$ function on $s\in[0,\ve)$, where $s=(t-1)^{\f{1}{2k}}$;

(2) $\vp(t)\in C^{\f{k+1}{k}}[1,1+\ve)$.

(1) Set $\la(s)=\f{\vp(t)}{s^{2k+1}}$. Then (<ref>) becomes

s'(s)+(2k+2)(s)=s d(s,(s)),

By the proof procedure of Lemma <ref>, we have $s^jd(s,\la)\in C^{j}$ for $j=0,1,2$.

Then by the same analysis in [11], there exists a unique solution $\la(s)\in C^{1}[0,\ve)$
to (<ref>) and further $s\la'(s)\in C^{1}$. Due to $(s\la(s))'=s\la'(s)+\la(s)$,
then $s\la(s)\in C^{2}$, and $s^2\la'(s)=s^2d-(2k+1)s\la\in C^2$.
Therefore, $s\rightarrow\vp'(t)=\f{1}{2k}[s^2\la'(s)+(2k+1)s\la(s)]$ is
a $C^2$ function.

(2) Let $s\rightarrow0+$ in (<ref>), we have

\la'(0)=\f{d(0,0)}{2k+3}=\f{g^{(2k+2)}(0)}{(2k+3)!}+\lim_{s\rightarrow0_+}\f{b(s\mu_+,s\mu_-)}{s^2}=O(1).

Therefore $\la(s)=O(s)$ and then $\vp(t)=O(s^{2k+2})=O((t-1)^\f{k+1}{k})\in C^{\f{k+1}{k}}[1,1+\ve)$.

Remark 3.1. The regularity of $\vp(t)$ in Lemma <ref> is optimal. Indeed, we consider
the following problem

\left\{
\begin{array}{ll}
&\p_t u+\p_x(\f{1}{2}u^2)=0, \\
&\ds u(0,x)=-x+x^{2k+1}+|x|^{2k+2+\ve},\ \ve>0.
\end{array}
\right.

In this case, we have $g(x)=-x+x^{2k+1}+x^{2k+2}$ and
y_+(t,(t))   =   (t-1)^12k(1+(t)2k(t-1)^2k+12k-(t-1)^12k2k)+o((t-1)^1k)

y_-(t,(t))   =   (t-1)^12k(-1+(t)2k(t-1)^2k+12k-(t-1)^12k2k)+o((t-1)^1k).
It follows from Rankine-Hugoniot condition that
'(t)   =   -y_-(t,(t))+y_+(t,(t))2

=   (t)k(t-1)-(t-1)^1kk+o((t-1)^1k).This derives $\vp(t)\in C^{\f{k+1}{k}}[1,1+\ve)$ which is optimal.

Under assumption (<ref>), for any $c\in(-\f{2k}{(2k+1)^{\f{1}{2k}+1}},+\i)$, there exist $\ve=\ve(c)>0$
and $\dl=\dl(c)>0$
such that for $(s,\la)\in\{0< s<\ve,c-\dl<\la<c+\dl\}$,
$(s,\la)\rightarrow y_+(t,x)$ has the expansion
and for $(s,\la)\in\{0<s<\ve,-c-\dl<\la<-c+\dl\}$,
$(s,\la)\rightarrow y_-(t,x)$ has the expansion
where $\mu_c$ is the unique solution in $(\f{1}{(2k+1)^\f{1}{2k}},+\i)$ of the equation

By (<ref>), (<ref>) and (<ref>),
we have that
_μG(0,±c,±μ_c)   =   -1+(2k+1)μ_c^2k>0,

_s   G(0,±c,±μ_c)   =   μ_c^2k+2h'(0)
=μ_c^2k+2g^(2k+2)(0)(2k+2)! ,

_G(0,±c,±μ_c)   =   -1.
Then by the implicit function theorem one has that there is a unique function
$\mu_\pm(s,\la)$ near $(0,\pm c)$ satisfying $G(s,\la,\mu_\pm(s,\la))\equiv0$ and
admitting the expansion

Thus (<ref>) and (<ref>) are proved.

To study the asymptotic behavior of $y_\pm$ near the $x-$axis,
we now take the following transform

ξ=x^12k+1, η=t-1ξ^2k, ν=yξ.
Under assumption (<ref>), by divided $\xi^{2k+1}$, (<ref>) then becomes

We now have

Under assumption (<ref>), for small $\dl>0$, there exists some small constant $\ve>0$ such that for $(\eta,\xi)\in\{ |\eta|\le \ve,\ 0<\xi<\dl\}$, we can get the expansion of $y_+(t,x)$ on $(\xi,\eta)$
and for $(\eta,\xi)\in\{|\eta|\le \ve,\ -\dl<\xi<0\}$, we can get

It follows from direct computation that $H(0,0,1)=0$ and
_νH   =   -η+(2k+1)(1+ξ^2kη)ν^2k+(2k+1)(1+ξ^2kη)ν^2kh(ξν)

_ηH   =   -ν+ξ^2kν^2k+1+ξ^2kν^2k+1h(ξν),

_ξH   =   2kξ^2k-1ην^2k+1+2kξ^2k-1ην^2k+1h(ξη)+(1+ξ^2kη)ν^2k+1νh'(ξν).
_νH(0,0,1)   =   2k+1,

_ηH(0,0,1)   =   -1,

_ξH(0,0,1)   =   h'(0)=g^(2k+2)(0)(2k+2)!.
By the implicit function theorem, there exists a unique solution $\nu=\nu(\eta,\xi)$
of (<ref>) near $(\eta,\xi)=(0,0)$ such that
Therefore, we obtain (<ref>) for $\xi>0$. Analogously,  (<ref>) holds for $\xi<0$.

Next we consider the asymptotic behavior of $y(t,x)$ near the blowup point $(1, 0)$ in the domain $\{(t,x):\ t<1\}$.
Without confusions, we still use the same notation as for $t>1$. Note that through each point in $\{t<1\}$,
there exists a unique characteristic line. By taking the following transform similar to (<ref>)

τ=1-t, s=τ^12k, μ=ys, =xs^2k+1,
and then by divided $s^{2k+1}$ on two sides of  (<ref>), then (<ref>) becomes

For each $c\in\bR$, there exist $\ve=\ve(c),\dl=\dl(c)>0$
such that for $(s,\la)\in\{0< s<\ve,c-\dl<\la<c+\dl\}$,
$(s,\la)\rightarrow y(t,x)$ has the expansion
where $\mu_c$ is the unique solution of the equation

It follows from direct computation that
_μG(s,,μ)   =   1+(2k+1)(1 -s^2k)μ^2k+(2k+1)(1-s^2k)μ^2k+1sh'(sμ),

_s   G(s,,μ)   =   -2ks^2k-1μ^2k+1-2ks^2k-1μ^2k+1h(sμ)

_G(s,,μ)   =   -1.
_μG(0,c,μ_c)   =   1+(2k+1)μ_c^2k>0,

_s   G(0,c,μ_c)   =   μ_c^2k+2h'(0)=μ_c^2k+2g^(2k+2)(0)(2k+2)!,

_G(0,c,μ_c)   =   -1.
By the implicit function theorem, there exists a $\mu(s,\la)$ near $(c,\mu_c)$ satisfying

from which we can deduce (<ref>).

We start to prove Theorem <ref>.

Proof of Theorem <ref>:

(1) By Lemma <ref>, $\vp(t)\in C^\f{k+1}{k}[1,1+\ve)$ and
$u\in C^1((1,1+\ve)\times\bR)\setminus \{x=\vp(t)\})$ for small $\ve>0$ have been shown.

(2)  Let $\dl,\ \ve>0$ be the constants obtained in Lemma <ref> and denote
Ø_x,+   =   B∩{(t,x): 0<x<^2k+1, |t-1|<x^2k2k+1},

Ø_x,-   =   B∩{(t,x): -^2k+1<x<0, |t-1|<(-x)^2k2k+1},

Ø_0   =   B∩{(t,x): t<1, |x|^12k+1<2^12k(1-t)^12k}.
Let $c_0\in(0,\f{2k}{(2k+1)^{1+\f{1}{2k}}})$ be some fixed constant and denote
Ø_t,+   =   B∩{(t,x): -c_0(t-1)^12k<x^12k+1<2^12k(t-1)^12k},

Ø_t,-   =   B∩{(t,x): -2^12k(t-1)^12k<x^12k+1<c_0(t-1)^12k}.
It is easy to see that for $(t,x)\in\O_{x,+}\cup\O_{t,+}$, $u(t,x)=u_0(y_+(t,x))$;
for $(t,x)\in\O_{x,-}\cup\O_{t,-}$, $u(t,x)=u_0(y_-(t,x))$; for $(t,x)\in\O_{0}$, $u(t,x)=u_0(y(t,x))$. By Heine-Borel property of compactness, there exist $\{c_{j,\pm},\dl_{j,\pm}=\dl_{j,\pm}(c_{j,\pm}),\ve_{j,\pm}=\ve_{j,\pm}(c_{j,\pm})\}_{j=1}^n$ and $\{c_{j,0},\dl_{j,0}=\dl_{j,0}(c_{j,0}),\ve_{j,0}=\ve_{j,0}(c_{j,0})\}_{j=1}^n$ such that

Ø_t,+⊂∪_j=1^nØ_t,+^j, Ø_t,-⊂∪_j=1^nØ_t,-^j, Ø_0⊂∪_j=1^nØ_0^j,
Ø_t,+^j   =   {(t,x): 0<(t-1)^12k<_j,+, c_j,+-_j,+<x(t-1)^2k+12k<c_j,++_j,+},

Ø_t,-^j   =   {(t,x): 0<(t-1)^12k<_j,-, c_j,--_j,-<x(t-1)^2k+12k<c_j,-+_j,-},

Ø_0^j   =   {(t,x): 0<(1-t)^12k<_j,0, c_j,+-_j,+<x(1-t)^2k+12k<c_j,++_j,+},
and these domains satisfy the corresponding properties in Lemma <ref> and <ref>.

Set $B=\{(t,x):\ 0<\sqrt{(t-1)^2+x^2}< \rho\}$, and choose $\rho>0$ sufficiently small such that

We now establish the behaviors of $u$ and its derivatives near $(1,0)$.
It suffices to only consider this in the domains $\O_{x,+}$, $\O_{t,+}^j$ and  $\O_{0}^j$
since the other cases can be treated analogously.

For $(t,x)\in\O_{x,+}$, we have

for $(t,x)\in\O_{t,+}^j$,

and for $(t,x)\in\O_{0}^j$,

Therefore (<ref>) is obtained.

Let's turn to prove the estimates (<ref>) and (<ref>).
Note that

For $(t,x)\in\O_{t,+}^j$, by (<ref>) in Lemma <ref> we have

1+tg'(y_+(t,x))   =   -s^2k+(2k+1)(1+s^2k)y_+^2k(t,x)+O(y_+^2k+1(t,x))

=   (-1+(2k+1)μ_c_j,+^2k)s^2k+O(s^2k+1+s^2k|-c_j,+|)

≳    s^2k

=   (t-1)

≳   |t-1|+|x|^2k2k+1,
where $s=t-1$, and the fact of $-1+(2k+1)\mu_{c_{j,+}}^{2k}>0$ has been used.

For $(t,x)\in\O_{x,+}$, by (<ref>) in Lemma <ref> we have

1+tg'(y(t,x))   =   -ηξ^2k+(2k+1)(1+ηξ^2k)y^2k(t,x)+O(y^2k+1(t,x))

=   ξ^2k+O(ξ^2k+1+ηξ^2k)

≳    ξ^2k

=   x^2k2k+1

≳   |t-1|+|x|^2k2k+1,
where $\xi=x^\f{1}{2k+1}$.

For $(t,x)\in\O_{0,+}^j$, by (<ref>) in Lemma <ref> we arrive at

1+tg'(y(t,x))   =   s^2k+(2k+1)(1+s^2k)y^2k(t,x)+O(y^2k+1(t,x))

=   (1+(2k+1)μ_c_j,0^2k)s^2k+O(s^2k+1+s^2k|-c_j,0|)

≳    s^2k

=   (1-t)

≳   |t-1|+|x|^2k2k+1,
where $s=1-t$, and the fact of $1+(2k+1)\mu_{c_{j,0}}^{2k}>0$ has been used.

Therefore, $1+tg'(y(t,x))\gtrsim|t-1|+|x|^\f{2k}{2k+1}$ holds for $(t,x)\in B$.
In light of (<ref>) and the fact of $g(y(t,x))\sim y(t,x)$ in $B$,
(<ref>)-(<ref>) hold and thus the proof of Theorem <ref>
is completed.

§ PROOF OF THEOREM <REF>

By the characteristics method, we can define $u_\pm(t,x)=u_0(y_\pm(t,x))$. By (<ref>),
the shock curve $x=\vp(t)$ satisfies

a(x,y)≜   f(u_0(x))-f(u_0(y))u_0(x)-u_0(y), if x≠y,

g(x), if x=y.
By (<ref>), we have $a(x,y)$ in $C^\i(\bR^2)$ and

where $b(x,y)=b(y,x)$ and $b(x,y)=O(x^2+y^2)\in C^\i$.

We now study the regularity of the shock wave $x=\vp(t)$ as a function of $s=|\ln(t-1)|^{-\f{1}{p}}$.

Under assumption (<ref>), for (<ref>) and small $\ve>0$, there exists a solution $x=\vp(t)$ on $[1,1+\ve)$
such that

(1) $s\rightarrow\vp(t)$ is of $C^{1}$ on $[0,\ve)$;

(2) $x=\vp(t)$ is of $C^1$ on $[1,1+\ve]$ with the behavior $\vp(t)=O(s^2\tau)$.

(1) Set $\la(s)=\f{\vp(t)}{s\tau}$. Then

\f{d\vp}{dt}(t)=s\left(\f{s^{1+p}}{p}\f{d\la(s)}{ds}+(\f{s^p}{p}+1)\la(s)\right).

Substituting this into (<ref>) yields

=   1sa(sμ_+(s,(s)),sμ_-(s,(s)))

=   -12(μ_+(s,(s))+μ_-(s,(s)))+1sb(sμ_+(s,(s)),sμ_-(s,(s)))

=   -s^pp+O(s+s^2p+1||+s^p+1||^2).
By (<ref>) and (<ref>), we have

Moreover, $s^ld(s,\la)\in C^{l+p}$ holds for $l=0,1$, which is derived by Lemma <ref> and

In addition, (<ref>) in $(s,\la)$ can be written as

s^1+ppd(s)d s+(2ps^p+1)(s)=d(s,(s)),

This yields

(s)=ps^-2∫_0^s ø^1-pe^s^-p-ø^-pd(ø,(ø))dø.
It follows from direct computation that

|(s)|   ≤   s^-2∫_0^ss^2e^s^-pø^-1-pe^-ø^-p|d(ø,(ø))|dø)

≲    (s+s^2p+1_L^ı[0,s]+s^p+1^2_L^ı[0,s])∫_0^se^s^-pø^-1-pe^-ø^-pdø

≲    s+s^2p+1_L^ı[0,s]+s^p+1^2_L^ı[0,s].Thus $\|\la\|_{L^\i[0,s]}\leq Cs$ for $s\in(0,\ve]$ and small $\ve>0$. By the analogous
computation, we can apply the contraction mapping theorem to prove that there exists a continuous solution
$\la$ to the integral equation (<ref>). From (<ref>), we have

'(s)   =   ps^-1-pd(s,(s))-2ps^-3∫_0^sø^1-pe^s^-p-ø^-pd(ø,(ø))dø-ps^-p-3e^s^-p∫_0^sø^2d(ø,(ø))de^-ø^-p

=   -2ps^-3∫_0^sø^1-pe^s^-p-ø^-pd(ø,(ø))dø+ps^-p-3e^s^-p∫_0^s
This derives

|'(s)|   ≲   s^-3∫_0^ss^2ø^-1-pe^s^-p-ø^-p|d(ø,(ø))|dø+s^-p-3∫_0^ss^2+pø^-1-pe^s^-p-ø^-p|d(ø,(ø))|dø

≲   1+_L^ı[0,s]s+s^-1e^s^-p∫_0^sø^-1-pe^-ø^-p(s+s^2'_L^ı[0,s]+^2_L^ı[0,s])dø

≲   1+_L^ı[0,s]s+s'_L^ı[0,s],and then $\la'(s)\in C[0,\ve]$ can be shown. In addition, by $\vp(t)=s\tau\la(s)$ and $s=|\ln\tau|^{-\f{1}{p}}$,
then $\vp(t)=O(s^2\tau)$ holds.

Remark 4.1. The regularity of $\vp(t)$ in Lemma <ref> is optimal. Indeed, we consider Burgers' equation

\left\{
\begin{array}{ll}
&\ds\f{\p u}{\p t}+\f{\p}{\p x}(\f{1}{2}u^2)=0, \\
&\ds u(0,x)=-x+\f{1}{p}e^{-|x|^{-p}}\left(x+x^2\right),\ p>0.
\end{array}
\right.

In this case, $g(x)=u_0(x)=-x+\f{1}{p}e^{-|x|^{-p}}\left(x+x^2\right)$. On the other hand, (<ref>) can be written as

F(s,\la,\zeta)\triangleq G(s,\la,\mu)=\mu(-1+e^\zeta)+\mu e^{\zeta-s^{-p}}+s\mu^{2}e^{\zeta}(1+e^{-s^{-p}})-\la,

which derives $\f{\p F}{\p s}|_{s=\la=\zeta=0}=1$. So we have that
y_+(t,x)   =   |ln(t-1)|^-1p(1+lnpp|ln(t-1)|^-1)+(t)p(t-1)|ln(t-1)|+1p|ln(t-1)|^-1-2p+o(|ln(t-1)|^-1-2p),

y_-(t,x)   =   -|ln(t-1)|^-1p(1+lnpp|ln(t-1)|^-1)+(t)p(t-1)|ln(t-1)|+1p|ln(t-1)|^-1-2p+o(|ln(t-1)|^-1-2p).It follows from  Rankine-Hugoniot condition that
'(t)   =   -y_-(t,(t))+y_+(t,(t))2

=   -(t)p(t-1)|ln(t-1)|-1p|ln(t-1)|^-1-2p+o(|ln(t-1)|^-1-2p).This means that $\vp(t)=O((t-1)|\ln(t-1)|^{-\f{2}{p}})$ is optimal.

Under assumption (<ref>), for any $c\in(-1,+\i)$, there exist $\ve=\ve(c),\dl=\dl(c)>0$
such that for $(s,\la)\in\{0<s<\ve,c-\dl<\la<c+\dl\}$,
$(s,\la)\rightarrow y_+(t,x)$ has the expansion
and for $(s,\la)\in\{0<s<\ve,-c-\dl<\la<-c+\dl\}$,
$(s,\la)\rightarrow y_-(t,x)$ has the expansion

Similarly to Lemma <ref>, we first consider the case of $p\geq1$. By taking $\la=c>-1$ and $s=0$ in (<ref>),
we have the solution $\zeta_{c}=\ln(c+1)$ and $\mu_{c}=1$. Furthermore, direct computation yields
_s F_1(0,ζ_c,c)   =   c(ln(c+1)+lnp)_p^1+p(c+1)2r”_0(0),

_ζF_1(0,ζ_c,c)   =   c+1,

_F_1(0,ζ_c,c)   =   -1,
where $\dl_p^1=\left\{\begin{array}{cc}1,&\ p=1,\\0,&\ p>1.\end{array}\right.$ By the implicit function theorem, we have
and then

from which (<ref>) follows.

For $p\in(0,1]$, recalling $\o=s^p$ and then taking $\o=0$, $\la=c$ and $\zeta=\zeta_c=\ln(c+1)$ in (<ref>),
one can arrive at
_øF_2(0,ζ_c,c)   =   c(ln(c+1)+lnp)p,

_ζF_2(0,ζ_c,c)   =   c+1,

_F_2(0,ζ_c,c)   =   -1,
and then by the implicit function theorem we have

and (<ref>) can be obtained.

On the other hand for $\la=-c$ and $\mu_{-c}=-1$, (<ref>) can be proved by the same method.

Next we consider the behavior of $y(t,x)$ near $x-$axis. Note that for $y>0$,

is a monotonically increasing function of $\xi$ from $[0,+\i)$ to $[0,+\i)$. Then there exists a unique inverse function $h(x)$ of (<ref>) satisfying that for $x>0$ sufficiently small,

h(|x|),     x>0,

-h(|x|),     x<0,

ν=yξ, η=(t-1)h(|x|)|x|.

Under assumption (<ref>), there exist some constants $\ve$, $\dl>0$ small enough such that for $(\eta,\xi)\in\{0<\xi<\dl,-\ve<\eta<\ve\}$,
$(\eta,\xi)\rightarrow y_+(t,x)$ has the expansion
and for $(\eta,\xi)\in\{-\dl<\xi<0,-\ve<\eta<\ve\}$,
$(\eta,\xi)\rightarrow y_+(t,x)$ has the expansion

We only consider the case of $x>0$ and then $y=y_+(t,x)>0$ since the other case can be treated analogously.
By $x=\xi e^{-\xi^{-p}}$, $y=\xi\nu$, $t-1=\eta e^{-\xi^{-p}}$ and (<ref>), (<ref>) becomes

Similarly to Lemma <ref>, we divide the proof procedure into two cases of $p\ge1$ and $0<p<1$.
Firstly we consider $p\ge1$. Set $\th=\xi^{-p}(1-\nu^{-p})-\ln p$ and $\nu=\left(1-\xi^{p}(\th+\ln p)\right)^{-\f{1}{p}}$. Then (<ref>) becomes

+ηνe^þ-ξ^-p+p e^þ(ηe^-ξ^-p+1)ξr_0(ξν).
Note $J_1(0,0,0)=0$, and

_ξν=ξ^p-1(þ+lnp)(1-ξ^p(þ+lnp))^-1p-1, _þν=1pξ^p(1-ξ^p(þ+lnp))^-1p-1
are bounded near $\xi=0$ and $\th=0$ by $p\ge1$. Due to

_ξJ_1   =   (e^þ-η)_ξν+ηe^þ-ξ^-p(_ξν+pνξ^-p-1)

-pe^þξ^2(ηe^-ξ^-p(1-pξ^-p)+1)r_0(ξν)+p e^þ(ηe^-ξ^-p+1)ξr'_0(ξν)(ν+ξ_ξν),

_þJ_1   =   e^þν+(e^þ-η)_þν+ηe^þ-ξ^-p(ν+_þν)+p e^þ(ηe^-ξ^-p+1)ξ(r_0(ξν)+ξr'_0(ξν)_þν),

_ηJ_1   =   -ν+e^þ-ξ^-p(ν+pξr_0(ξν)),
we then obtain

_ξJ_1(0,0,0)=_1^plnp+p2r”_0(0), _þJ_1(0,0,0)=1, _ηJ_1(0,0,0)=-1.
Thus by the implicit function theorem, one can deduce that there exists a unique function $\th=\th(\eta,\xi)$ near $(\eta,\xi)=(0,0)$ satisfying

Recalling $\nu=\left(1-\xi^p(\th+\ln p)\right)^{-\f{1}{p}}$, we then have

For $p\in(0,1)$, set $\vs=\xi^p$ and then
$\nu=(1-\vs(\th+\ln p))^{-\f{1}{p}}$.

In this case, (<ref>) becomes

+ηνe^þ-^-1+p e^þ(ηe^-^-1+1)^1pr_0(^1pν).

By $J_2(0,0,0)=0$ and

_ν=1p(þ+lnp)(1-(þ+lnp))^-1p-1, _þν=p(1-(þ+lnp))^-1p-1,
under  assumption (<ref>), we have that $J_2(\eta,\xi,\th)\in C^\f{1}{p}$ and

_ξJ_2   =   (e^þ-η)_ν+ηe^þ-^-1(_ν+ν^-2)

-e^þ^1p+1(ηe^-^-1(1-p^-1)+1)r_0(^1pν)+p e^þ(ηe^-^-1+1)^1pr'_0(^1pν)(1p^1p-1ν+^1p_ν),

_þJ_2   =   e^þν+(e^þ-η)_þν+ηe^þ-^-1(ν+_þν)+p e^þ(ηe^-^-1+1)^1p(r_0(^1pν)+^1p r'_0(^1pν)_þν),

_ηJ_2   =   -ν+e^þ-^-1(ν+p^1pr_0(^1pν)).
This yields

_J_2(0,0,0)=lnpp, _þJ_2(0,0,0)=1, _ηJ_2(0,0,0)=-1.
By the implicit function theorem, we know that there exists a unique solution $\th=\th(\eta,\vs)$
for $(\eta,\vs)$ near $(0,0)$, which satisfies $\th(0,0)=0$ and

By $\xi=\vs^\f{1}{p}$, then

ν(η,ξ)   =   (1-(þ+lnp))^-1p

=   (1-lnp-η+O(^2+η^2))^-1p

=   1+lnppξ^p+1pξ^pη+O(ξ^2p+ξ^pη^2).
Therefore we finish the proof of (<ref>).

For $x<0$, we can transform (<ref>) to $H(\eta,-\xi,-\nu)=0$. Analogously, we can obtain (<ref>) about $y_-(t,x)$
and $x<0$.

Next we consider the behavior of $y(t,x)$ for $t<1$ near $(1,0)$. Without of confusion, we still denote

τ=1-t, s=|lnτ|^-1p, =xsτ, μ=ys,
as in the case of $t>1$. By divided $s\tau$, (<ref>) then becomes

Under assumption (<ref>), we have that

(1) for any $c>1$, there exist $\ve=\ve(c),\dl=\dl(c)>0$
such that for $(s,\la)\in\{0<s<\ve,1<c-\dl<\la<c+\dl\}$,
$(s,\la)\rightarrow y(t,x)$ has the expansion

(2) for any $0\le c<1$, there exist $\ve=\ve(c),\dl=\dl(c)>0$
such that for $(s,\la)\in\{0<s<\ve,c-\dl<\la<c+\dl<1\}$,
$(s,\la)\rightarrow y(t,x)$ has the expansion

(3) for any $c<-1$, there exist $\ve=\ve(c),\dl=\dl(c)>0$
such that for $(s,\la)\in\{0<s<\ve,c-\dl<\la<c+\dl<-1\}$,
$(s,\la)\rightarrow y(t,x)$ has the expansion

(4) for any $-1<c<0$, there exist $\ve=\ve(c),\dl=\dl(c)>0$
such that for $(s,\la)\in\{0<s<\ve,-1<c-\dl<\la<c+\dl<1\}$,
$(s,\la)\rightarrow y(t,x)$ has the expansion

We only prove the cases (1), (2) for $c\ge0$ since  the cases (3), (4) can be obtained by the same way.

(1) If $c>1$, it is similar to the proof of Lemma <ref> that we adopt the variable transformation

and then $\mu=(1-s^p(\zeta+\ln p))^{-\f{1}{p}}$.

At first, we assume $p\ge 1$. Then (<ref>) becomes

It is easy to see that for $c> 1$, $F_1(0,c,\ln(c-1))=0$ and

_s F_1(s,,ζ)   =   (1+e^ζ-e^ζ-s^-p)_sμ+ps^-p-1μe^ζ-s^-p

+p e^ζs^2(e^-s^-p(1-ps^-p)-1)r_0(sμ)

_F_1(s,,ζ)   =   -1,

_ζF_1(s,,ζ)   =   (1+e^ζ-e^ζ-s^-p)_ζμ+μe^ζ(1-e^-s^-p)+pe^ζ(1-e^-s^-p)s(r_0(sμ)+r'_0(sμ)_ζμ).This yields
_s F_1(0,c,ln(c-1))=c(lnp+ln(c-1))_1^p+p(c-1)2r”_0(0), _F_1(0,c,ln(c-1))=-1,

By the implicit function theorem, there exists a unique solution $\zeta(s,\la)$ satisfying $\zeta(0,c)=\ln(c-1)$ and $\zeta(s,\la)$ is $C^\i$ near $(0,c)$ with

μ   =   (1-s^p(ζ+lnp))^-1p

=   1+ln(c-1)+lnpps^p+s^p(-c)p(c-1)+O(s^p+1+s^p(-c)^2).

Secondly, we consider the case of $p\in(0,1)$. Let $\o=s^{p}$
and then $\mu=\left(1-\o(\zeta+\ln p)\right)^{-\f{1}{p}}$. In this case, (<ref>) becomes

It is obvious that for $c> 1$, $F_2(0,c,\ln(c-1))=0$. By direct computation, we have that

_øF_2(ø,,ζ)   =   (1+e^ζ-e^ζ-ø^-1)_øμ+s^-2μe^ζ-ø^-1

_F_2(ø,,ζ)   =   -1,

_ζF_2(ø,,ζ)   =   (1+e^ζ-e^ζ-ø^-1)_ζμ+μe^ζ(1-e^-ø^-1)+
pe^ζ(1-e^-ø^-1)ø^1p(r_0(ø^1pμ)+r'_0(ø^1pμ)_ζμ).This yields
_s F_1(0,c,ln(c-1))=c(lnp+ln(c-1))p, _F_1(0,c,ln(c-1))=-1, _ζF_1(0,c,ln(c-1))=c-1.
By the implicit function theorem, there exists a unique solution $\zeta(s,\la)$ satisfying $\zeta(0,c)=\ln(c-1)$ and $F_2(s,\la,\zeta(s,\la))=0$ with

μ   =   (1-s^p(ζ+lnp))^-1p

=   1+ln(c-1)+lnpps^p+s^p(-c)p(c-1)+O(s^p+1+s^p(-c)^2).
Together with (<ref>) and (<ref>), this derives (<ref>).

(2) For $0<c<1$, define

It is clear that $G(0,c,c)=0$ and

_s G(s,,μ)   =   s^-p-1(μ+pr_0(sμ)s)

_G(s,,μ)   =   -1,

_μG(s,,μ)   =   1+s^-pμ^-p-1(μ+pr_0(sμ)s)

(e^s^-p(1-μ^-p)-e^-s^-pμ^-p).Then in light of $c<1$, it follows that

_s G(0,c,c)   =   0,

_G(0,c,c)   =   -1,

_μG(0,c,c)   =   1.Thus by the implicit function theorem, there exists a unique solution $\mu=\mu(s,\la)$ satisfying that $\mu(0,c)=c$ and

Therefore, by $y=s\mu$, we finish the proof of (<ref>).

Proof of Theorem <ref>:

(1) It can be obtained by Lemma <ref>.

(2) Since we don't get the behavior of $y(t,x)$ for $c=\pm1$ in Lemma <ref>, we have to choose the domain $\O_{-,t,+}^0$ and $\O_{-,t,-}^0$ as follows
Ø_t,-,+^0   =   {(t,x): 0<s<_0, 1-_0<xsτ<1+_0},

Ø_t,-,-^0   =   {(t,x): 0<s<_0, -1-_0<xsτ<-1+_0},
where $\tau=1-t$, $s=|\ln \tau|^{-\f{1}{p}}$ and $\ve_0$, $\dl_0>0$ sufficiently small.

We only consider the behavior in $\O_{-,t,+}^0$ since the treatment
in $\O_{-,t,-}^0$ is similar. By monotonicity of $y(t,\cdot)$ for each fixed $t\in[0,1]$, we know that
for $(t,x)\in\O_{-,t,+}^0$,

Let's firstly turn to consider $y(t,(1+\dl_0)s\tau)$. Note that $\mu=\f{y(t,(1+\dl_0)s\tau)}{s}$ satisfies

For $p\ge1$, we let

and then $\mu=\left(1-s^p(\zeta+\ln p)\right)^{-\f{1}{p}}$.
In this case, (<ref>) becomes

Note that if $s\to 0+$, then $\zeta\to(\ln{\dl_0})+$ and $\mu\to1+$.
In addition,

and then by taking $s\to0+$, we have

_sμ(0,ln_0)=_1^p(ln_0+lnp), _ζμ(0,ln_0)=0.
On the other hand, it follows from direct computation that
_s J_1(s,ζ)   =   (1+e^ζ-e^ζ-s^-p)_sμ-ps^-p-1μe^ζ-s^-p

_ζJ_1(s,ζ)   =   μe^ζ(1-e^-s^-p)+(1+e^ζ-e^ζ-s^-p)_ζμ+pe^ζ(e^-s^-p-1)s(r_0(sμ)+s_ζμr'_0(sμ)).Together with (<ref>), this yields

_s J_1(0,ln_0)=(1+_0)(ln_0+lnp)+p_02r”_0(0), _ζJ_1(0,ln_0)=1.

Thus by the implicit function theorem, there exists a function $\zeta(s)$ satisfying $J(s,\zeta(s))=0$ such that

ζ(s)=ln_0-((1+_0)(ln_0+lnp)+p_02r”_0(0))s+O__0(s^2),  s∈(0,_0],
where $\ve_0=\ve_0(\dl_0)>0$ is small enough. Therefore

μ(s)=(1-s^p(ζ+lnp))^-1p=1+ln_0+lnpps^p+O__0(s^p+1),  s∈(0,_0].

For $p\in(0,1)$, we take $\o=s^p$ and then (<ref>) becomes

where $\mu=\left(1-\o(\zeta+\ln p)\right)^{-\f{1}{p}}$. By direct computation, one has
_øJ_2(ø,ζ)   =   (1+e^ζ-e^ζ-ø^-1)_øμ-ø^-2μe^ζ-ø^-1

_ζJ_2(ø,ζ)   =   μe^ζ(1-e^-ø^-1)+(1+e^ζ-e^ζ-ø^-1)_ζμ+pe^ζ(e^-ø^-1-1)ø^1p(r_0(ø^1pμ)+ø^1p_ζμr'_0(ø^1pμ)),and

This yields

and then

So we can obtain

and then

By (<ref>) and (<ref>), we have that for $p>0$ and $0<s<\ve_0$,

where $\ve_{0}>0$ is a small constant depending on $\dl_0$. It follows that for $s\in[0,\ve_0)$,

Secondly, we consider the behavior of $y(t,(1-\dl_0)s\tau)$. For $x=(1-\dl_0)s\tau$, $\mu=\f{y(t,(1-\dl_0)s\tau)}{s}$ satisfies

It is easy to know $\mu=1-\dl_0$ as $s\to0_+$. Furthermore, by direct computation, one has
_s L(s,μ)   =    s^-p-1(μ+psr_0(sμ))

+e^s^-p(1-μ^-p)-e^-s^-pμ^-ps^2(-r_0(sμ)+s μr'_0(sμ)),

_μL(s,μ)   =   (1+1pe^s^-p(1-μ^-p)

(ps^-p-1μ^-p-1r_0(sμ)+r'_0(sμ)).By assumption (<ref>), we have

_s L(0,1-_0)=0,  _μL(0,1-_0)=1.
Then there exists a function

such that $L(s,\mu(s))=0$ for $s>0$ sufficiently small and dependent on $\dl_0$. It follows that
for $s\in[0,\ve_0)$,

Thus for small fixed $\dl_0>0$, we have that for $s\in[0,\ve_0)$ and $x\in((1-\dl_0)s\tau,(1+\dl_0)s\tau)$,

Recalling $u(t,x)=u_0(y(t,x))$, we have that

and by

we have that for $(t,x)\in\O_{t,-,+}^0$,
ux(t,x)   =   u'_0(y(t,x))1+tg'(y(t,x))≲|ln(1-t)|^-1(1-t)^-1,

ut(t,x)   =   -u'_0(y(t,x))g(y(t,x))1+tg'(y(t,x))≲|ln(1-t)|^-1-1p(1-t)^-1.

We decompose the neighbourhood $B$ of $(1,0)$ into $\O_{x,+},\ \O_{x,-},\ \O^{j}_{t,+},\ \O^{j}_{t,-,m}$ for $j=1,2,\ldots,N$ and $\O^{j}_{t,-,+},\ \O^{j}_{t,-,-}$ for $j=0,1,2,\ldots,N$ as follows (see Figure 2 below).

Decomposition of $B$

Ø_x,+   =   {(t,x): x>0, 0<ξ<, -<η<},

Ø_x,-   =   {(t,x): x<0, 0<ξ<, -<η<},
where $(\xi,\eta)$ are defined in (<ref>). Taking the suitable constants $c_{1,+}<c_{2,+}<\ldots<c_{N-1,+}<c_{N,+}$, $\{\ve_{j,+}\}_{j=1}^N$ and $\{\dl_{j,+}\}_{j=1}^N$, and setting

Ø^j_t,+={(t,x): 0<s<_j,+, c_j,+-_j,+<<c_j,++_j,+}, s=|ln(t-1)|^-1 p, =xs(t-1),
where $B\cap\{t\ge0\}\subset \left(\cup_{j=1}^N\O^{j}_{t,+}\right)\cup\O_{x,+}\cup\O_{x,+}$.
By choosing small $\dl_0>0$, we can define $\O^{0}_{t,-,+}$ and $\O^{0}_{t,-,-}$ as in (<ref>) and (<ref>).
Meanwhile, taking some suitable constants $\{c_{1,-,m}\}_{j=1}^N$, $\{\ve_{j,-,m}\}_{j=1}^N$ and $\{\dl_{j,-,m}\}_{j=1}^N$ where $m=+,0,-$,
and setting

Ø^j_t,-,m={(t,x): 0<s<_j,-,m, c_j,-,m-_j,-,m<<c_j,-,m+_j,-,m}, s=|ln(1-t)|^-1 p, =xs(1-t),
where $c_{N,-,-}<\ldots<c_{1,-,-}<-1<c_{1,-,0}<c_{1,-,0}<\ldots<c_{N,-,0}<1<c_{1,-,+}<\ldots<c_{N,-,+}$, such
that $B\cap\{t\le0\}\subset \left(\cup_{j=0}^N\O^{j}_{t,-,+}\right)\cup
\left(\cup_{j=1}^N\O^{j}_{t,-,0}\right)\cup\left(\cup_{j=0}^N\O^{j}_{t,-,-}\right)\cup
\O_{x,+}\cup\O_{x,+}$ holds. Note that in $\O_{x,\pm}$, $y(t,x)\sim\xi$; and in others, $y(t,x)\sim s$.
Then similarly to (<ref>), we can obtain (<ref>). On the other hand,
in $\O_{x,\pm}$, $1+tg'(y(t,x))\gtrsim|x||\ln|x||$; and in others, $1+tg'(y(t,x))\gtrsim |t-1||\ln|t-1||$ as in (<ref>).
Then similarly to the proof for (<ref>) and (<ref>), we can establish (<ref>) and (<ref>).

[1]  S. Alinhac, Blowup of small data solutions for a class of quasilinear wave equations in two space dimensions.
Ann. of Math. 149 (1999), no. 1, 97-127.

[2] T. Buckmaster, S. Shkoller, V. Vicol, Formation of  shocks for 2D isentropic compressible Euler,
arXiv: 1907.03784v2, 8 Jul 2019.

[3] T. Buckmaster, S. Shkoller, V. Vicol, Formation of point shocks for 3D compressible Euler,
arXiv: 1912.04429v2, 22 Jun 2020.

[4] Chen Shuxing, Xin Zhouping, Yin Huicheng, Formation and construction of
shock wave for quasilinear hyperbolic system and its application to inviscid compressible flow.
The Institute of Mathematical Sciences at CUHK, 2010, Research Reports: 2000-10(069)

[5] Chen Shuxing, Dong Liming, Formation of shock for the p-system with general smooth initial
data. Sci. in China, Ser. A, 44 (2001), no. 9, 1139-1147.

[6] D. Christodoulou, The formation of shocks in 3-dimensional fluids.
EMS Monographs in Mathematics. European Mathematical Society (EMS), Zürich, 2007.

[7] D. Christodoulou, A. Lisibach, Shock development in spherical symmetry.
Ann. PDE 2 (2016), no. 1, Art. 3, 246 pp.

[8] D. Christodoulou, Miao Shuang, Compressible flow and Euler's equations,
Surveys of Modern Mathematics, 9. International Press, Somerville, MA; Higher
Education Press, Beijing, 2014.

[9] G. Holzegel, S. Klainerman, J. Speck, Wong Willie Wai-Yeung, Small-data shock formation in solutions
to 3D quasilinear wave equations: an overview. J. Hyperbolic Differ. Equ. 13 (2016), no. 1, 1-105.

[10]   Kong Dexing, Formation and propagation of singularities for
$2\times 2$ quasilinear hyperbolic systems. Trans. Amer. Math. Soc. 354 (2002), no. 8, 3155-3179.

[11]    M.P. Lebaud, Description de la formation d'un choc dans le
$p-$système. J. Math. Pures Appl. (9) 73 (1994), no. 6, 523-565.

[12] J. Luk, J. Speck, Shock formation in solutions to the 2D compressible Euler equations in the presence of non-zero vorticity. Invent. Math. 214 (2018), no. 1, 1-169.

[13] F. Merle, P. Raphaël, I. Rodniaski, J. Szefel, On the
implosion of a three dimensional compressible fluid, arXiv: 1912.11009v2, 13 Jun 2020.

[14] Miao Shuang, Yu Pin, On the formation of shocks for quasilinear wave equations.
Invent. Math. 207 (2017), no. 2, 697-831.

[15] J. Speck, Shock formation for 2D quasilinear wave systems featuring multiple speeds: blowup
for the fastest wave, with non-trivial interactions up to the singularity.
Ann. PDE 4 (2018), no. 1, Art. 6, 131 pp.

[16] Yin Huicheng,  Formation and construction of a shock wave for 3-D compressible Euler equations
with the spherical initial data. Nagoya Math. J. 175 (2004), 125-164.

[17]  Xin Zhouping, Some current topics in nonlinear conservation laws. Some current topics on nonlinear conservation laws, xiii-xxxi, AMS/IP Stud. Adv. Math., 15, Amer. Math. Soc., Providence, RI, 2000.