
We study general parabolic equations of the form $u_t = div A(x,t, u,D u) +
div(F^{p2} F)+ f$ whose principal part depends on the solution itself. The
vector field $A$ is assumed to have small mean oscillation in $x$, measurable
in $t$, Lipschitz continuous in $u$, and its growth in $Du$ is like the
$p$Laplace operator. We establish interior Calder\'onZygmund estimates for
locally bounded weak solutions to the equations when $p>2n/(n+2)$. This is
achieved by employing a perturbation method together with developing a
twoparameter technique and a new compactness argument. We also make crucial
use of the intrinsic geometry method by DiBenedetto \cite{D2} and the maximal
function free approach by Acerbi and Mingione \cite{AM}.

In this paper, we investigate regularity for solutions to the linearized
MongeAmp\`ere equations when the nonhomogeneous term has low integrability. We
establish global
$W^{1,p}$ estimates for all $p<\frac{nq}{nq}$ for solutions to the equations
with right hand side in $L^q$ where $n/2<q\leq n$. These estimates hold under
natural assumptions on the domain, MongeAmp\`ere measures and boundary data.
Our estimates are affine invariant analogues of the global $W^{1,p}$ estimates
of N. Winter for fully nonlinear, uniformly elliptic equations.

We study quasilinear elliptic equations of the form $\text{div}
\mathbf{A}(x,u,\nabla u) = \text{div}\mathbf{F} $ in bounded domains in
$\mathbb{R}^n$, $n\geq 1$. The vector field $\mathbf{A}$ is allowed to be
discontinuous in $x$, Lipschitz continuous in $u$ and its growth in the
gradient variable is like the $p$Laplace operator with $1<p<\infty$. We
establish interior $W^{1,q}$estimates for locally bounded weak solutions to
the equations for every $q>p$, and we show that similar results also hold true
in the setting of {\it Orlicz} spaces. Our regularity estimates extend results
which are only known for the case $\mathbf{A}$ is independent of $u$ and they
complement the wellknown interior $C^{1,\alpha}$ estimates obtained by
DiBenedetto \cite{D} and Tolksdorf \cite{T} for general quasilinear elliptic
equations.

The paper is concerned with sticky weak solutions to the equations of
pressureless gases in two or more space dimensions. Various initial data are
constructed, showing that the Cauchy problem can have (i) two distinct sticky
solutions, or (ii) no sticky solution, not even locally in time. In both cases
the initial density is smooth with compact support, while the initial velocity
field is continuous.

In this paper, we establish global $W^{2,p}$ estimates for solutions to the
linearized MongeAmp\`ere equations under natural assumptions on the domain,
MongeAmp\`ere measures and boundary data. Our estimates are affine invariant
analogues of the global $W^{2,p}$ estimates of Winter for fully nonlinear,
uniformly elliptic equations, and also linearized counterparts of Savin's
global $W^{2,p}$ estimates for the MongeAmp\`ere equations.

Let $\Omega\subset \R^n$ be a bounded convex domain and $\phi\in
C(\bar\Omega)$ be a convex function such that $\phi$ is sufficiently smooth on
$\partial\Omega$ and the MongeAmp\`ere measure $\det D^2\phi$ is bounded away
from zero and infinity in $\Omega$. The corresponding linearized
MongeAmp\`ere equation is \[ \trace(\Phi D^2 u) =f, \] where $\Phi := \det
D^2 \phi ~ (D^2\phi)^{1}$ is the matrix of cofactors of $D^2\phi$. We prove a
conjecture in \cite{GT} about the relationship between $L^p$ estimates for $D^2
u$ and the closeness between $\det D^2\phi$ and one. As a consequence, we
obtain interior $W^{2,p}$ estimates for solutions to such equation whenever the
measure $\det D^2\phi$ is given by a continuous density and the function $f$
belongs to $L^q(\Omega)$ for some $q> \max{\{p,n\}}$.

In this paper, we establish several geometric properties of boundary sections
of convex solutions to the MongeAmp\`ere equations: the engulfing and
separating properties and volume estimates. As applications, we prove a
covering lemma of Besicovitch type, a covering theorem and a strong type $pp$
estimate for the maximal function corresponding to boundary sections. Moreover,
we show that the MongeAmp\`ere setting forms a space of homogeneous type.