
We show that weak solutions to conormal derivative problem for elliptic
equations in divergence form are continuously differentiable up to the boundary
provided that the mean oscillations of the leading coefficients satisfy the
Dini condition, the lower order coefficients satisfy certain suitable
conditions, and the boundary is locally represented by a $C^1$ function whose
derivatives are Dini continuous. We also prove that strong solutions to oblique
derivative problem for elliptic equations in nondivergence form are twice
continuously differentiable up to the boundary if the mean oscillations of
coefficients satisfy the Dini condition and the boundary is locally represented
by a $C^1$ function whose derivatives are double Dini continuous. This in
particular extends a result of M. V. Safonov (Comm. Partial Differential
Equations 20:13491367, 1995)

We prove the uniqueness of H\"older continuous weak solutions via duality
argument and vanishing viscosity method for the KellerSegel system of porous
medium type equations coupled to the Stokes system in dimensions three. An
important step is the estimate of the Green function of parabolic equations
with lower order terms of variable coefficients, which seems to be of
independent interest.

We construct fundamental solutions of secondorder parabolic systems of
divergence form with bounded and measurable leading coefficients and divergence
free firstorder coefficients in the class of $BMO^{1}_x$, under the
assumption that weak solutions of the system satisfy a certain local
boundedness estimate. We also establish Gaussian upper bound for such
fundamental solutions under the same conditions.

We extend and improve the results in \cite{DK16}: showing that weak solutions
to full elliptic equations in divergence form with zero Dirichlet boundary
conditions are continuously differentiable up to the boundary when the leading
coefficients have Dini mean oscillation and the lower order coefficients verify
certain conditions. Similar results are obtained for nondivergence form
equations. We extend the weak type(1, 1) estimates in \cite{DK16} and
\cite{Es94} up to the boundary and derive a Harnack inequality for nonnegative
adjoint solutions to nondivergence form elliptic equations, when the leading
coefficients have Dini mean oscillation.

Under various conditions, we establish Schauder estimates for both divergence
and nondivergence form secondorder elliptic and parabolic equations involving
H\"older seminorms not with respect to all, but only with respect to some of
the independent variables. A novelty of our results is that the coefficients
are allowed to be merely measurable with respect to the other independent
variables.

We show that any weak solution to elliptic equations in divergence form is
continuously differentiable provided that the modulus of continuity of
coefficients in the $L^1$mean sense satisfies the Dini condition. This in
particular answers a question recently raised by Yanyan Li and allows us to
improve a result of Brezis. We also prove a weak type$(1,1)$ estimate under a
stronger assumption on the modulus of continuity. The corresponding results for
nondivergence form equations are also established.

Fluidsolid interaction has been a challenging subject due to their strong
nonlinearity and multidisciplinary nature. Many of the numerical methods for
solving FSI problems have struggled with nonconvergence and numerical
instability. In spite of comprehensive studies, it has been still a challenge
to develop a method that guarantees both convergence and stability.
Our discussion in this work is restricted to the interaction of viscous
incompressible fluid flow and a rigid body. We take the monolithic approach by
Gibou and Min that results in an extended Hodge projection. The projection
updates not only the fluid vector field but also the solid velocities. We
derive the equivalence of the extended Hodge projection to the Poisson equation
with nonlocal Robin boundary condition. We prove the existence, uniqueness,
and regularity for the weak solution of the Poisson equation, through which the
Hodge projection is shown to be unique and orthogonal. Also, we show the
stability of the projection in a sense that the projection does not increase
the total kinetic energy of fluid and solid. Also, we discusse a numerical
method as a discrete analogue to the Hodge projection, then we show that the
unique decomposition and orthogonality also hold in the discrete setting. As
one of our main results, we prove that the numerical solution is convergent
with at least the first order accuracy. We carry out numerical experiments in
two and three dimensions, which validate our analysis and arguments.

We study an elliptic equation with measurable coefficients arising from
photoacoustic imaging in inhomogeneous media. We establish Holder continuity
of weak solutions and obtain pointwise bounds for Green's functions subject to
Dirichlet or Neumann condition.

Vecer derived a degenerate parabolic equation with a boundary condition
characterizing the price of Asian options with generally sampled average. It is
well understood that there exists a unique probabilistic solution to such a
problem but it remained unclear whether the probabilistic solution is a
classical solution. We prove that the probabilistic solutions to Vecer's PDE
are regular.

The aim of this paper is to investigate Green's function for parabolic and
elliptic systems satisfying a possibly nonlocal Robintype boundary condition.
We construct Green's function for parabolic systems with timedependent
coefficients satisfying a possibly nonlocal Robintype boundary condition
assuming that weak solutions of the system are locally H\"{o}lder continuous in
the interior of the domain, and as a corollary we construct Green's function
for elliptic system with a Robintype condition. Also, we obtain Gaussian bound
for Robin Green's function under an additional assumption that weak solutions
of Robin problem are locally bounded up to the boundary. We provide some
examples satisfying such a local boundedness property, and thus have Gaussian
bounds for their Green's functions.

We study the $3 \times 3$ elliptic systems $\nabla \times (a(x) \nabla\times
u)\nabla (b(x) \nabla \cdot u)=f$, where the coefficients $a(x)$ and $b(x)$
are positive scalar functions that are measurable and bounded away from zero
and infinity. We prove that weak solutions of the above system are H\"older
continuous under some minimal conditions on the inhomogeneous term $f$. We also
present some applications and discuss several related topics including
estimates of the Green's functions and the heat kernels of the above systems.

We consider the elliptic system of linear elasticity with bounded measurable
coefficients in a domain where the second Korn inequality holds. We construct
heat kernel of the system subject to Dirichlet, Neumann, or mixed boundary
condition under the assumption that weak solutions of the elliptic system are
H\"older continuous in the interior. Moreover, we show that if weak solutions
of the mixed problem are H\"older continuous up to the boundary, then the
corresponding heat kernel has a Gaussian bound. In particular, if the domain is
a two dimensional Lipschitz domain satisfying a corkscrew or nontangential
accessibility condition on the set where we specify Dirichlet boundary
condition, then we show that the heat kernel has a Gaussian bound. As an
application, we construct Green's function for elliptic mixed problem in such a
domain.

We study the Neumann Green's function for second order parabolic systems in
divergence form with timedependent measurable coefficients in a cylindrical
domain $\mathcal{Q}=\Omega\times (\infty,\infty)$, where $\Omega\subset
\mathbb{R}^n$ is an open connected set such that a multiplicative Sobolev
embedding inequality holds there. Such a domain includes, for example, a
bounded Sobolev extension domain, a special Lipschitz domain, and an unbounded
domain with compact Lipschitz boundary. We construct the Neumann Green's
function in $\mathcal{Q}$ under the assumption that weak solutions of the
systems satisfy an interior H\"older continuity estimate. We also establish
global Gaussian bounds for Neumann Green's function under an additional
assumption that weak solutions with zero Neumann data satisfy a local
boundedness estimate. In the scalar case, such a local boundedness estimate is
a consequence of De GiorgiMoserNash theory holds for equations with bounded
measurable coefficients in Sobolev extension domains, while in the vectorial
case, one may need to impose further regularity assumptions on the coefficients
of the system as well as on the domain to obtain such an estimate. We present a
unified approach valid for both the scalar and vectorial cases and discuss some
applications of our results including the construction of Neumann functions for
second order elliptic systems with measurable coefficients in two dimensional
domains.

We consider secondorder linear parabolic operators in nondivergence form
that are intrinsically defined on Riemannian manifolds. In the elliptic case,
Cabr\'e proved a global KrylovSafonov Harnack inequality under the assumption
that the sectional curvature of the underlying manifold is nonnegative. Later,
Kim improved Cabr\'e's result by replacing the curvature condition by a certain
condition on the distance function. Assuming essentially the same condition
introduced by Kim, we establish KrylovSafonov Harnack inequality for
nonnegative solutions of the nondivergent parabolic equation. This, in
particular, gives a new proof for LiYau Harnack inequality for positive
solutions to the heat equation in a manifold with nonnegative Ricci curvature.

We study Neumann functions for divergence form, second order elliptic systems
with bounded measurable coefficients in a bounded Lipschitz domain or a
Lipschitz graph domain. We establish existence, uniqueness, and various
estimates for the Neumann functions under the assumption that weak solutions of
the system enjoy interior H\"older continuity. Also, we establish global
pointwise bounds for the Neumann functions under the assumption that weak
solutions of the system satisfy a certain natural local boundedness estimate.
Moreover, we prove that such a local boundedness estimate for weak solutions of
the system is in fact equivalent to the global pointwise bound for the Neumann
function. We present a unified approach valid for both the scalar and the
vectorial cases.

We obtain sharp $L^p$ and H\"older estimates for the Neumann function of the
operator $\nabla \cdot \gamma \nabla  ik$ on a bounded domain. We also obtain
quantitative description of its singularity. We then apply these estimates to
quantitative photoacoustic imaging in inhomogeneous media. The problem is to
reconstruct the optical absorption coefficient of a diametrically small anomaly
from the absorbed energy density.

We construct Green's functions for divergence form, second order parabolic
systems in nonsmooth timevarying domains whose boundaries are locally
represented as graph of functions that are Lipschitz continuous in the spatial
variables and 1/2H\"older continuous in the time variable, under the
assumption that weak solutions of the system satisfy an interior H\"older
continuity estimate. We also derive global pointwise estimates for Green's
function in such timevarying domains under the assumption that weak solutions
of the system vanishing on a portion of the boundary satisfy a certain local
boundedness estimate and a local H\"older continuity estimate. In particular,
our results apply to complex perturbations of a single real equation.

We establish global Gaussian estimates for the Green's matrix of divergence
form, second order parabolic systems in a cylindrical domain under the
assumption that weak solutions of the system vanishing on a portion of the
boundary satisfy a certain local boundedness estimate and a local H\"older
estimate. From these estimates, we also derive global estimates for the Green's
matrix for elliptic systems with bounded measurable coefficients in two
dimensional domains. We present a unified approach valid for both the scalar
and vectorial cases and discuss several applications of our result.

We establish global pointwise bounds for the Green's matrix for divergence
form, second order elliptic systems in a domain under the assumption that weak
solutions of the system vanishing on a portion of the boundary satisfy a
certain local boundedness estimate. Moreover, we prove that such a local
boundedness estimate for weak solutions of the system is equivalent to the
usual global pointwise bound for the Green's matrix. In the scalar case, such
an estimate is a consequence of De GiorgiMoserNash theory and holds for
equations with bounded measurable coefficients in arbitrary domains. In the
vectorial case, one need to impose certain assumptions on the coefficients of
the system as well as on domains to obtain such an estimate. We present a
unified approach valid for both the scalar and vectorial cases and discuss
several applications of our result.

We establish Schauder estimates for both divergence and nondivergence form
secondorder elliptic and parabolic equations involving H\"older seminorms not
with respect to all, but only with respect to some of the independent
variables.

We study a certain one dimensional, degenerate parabolic partial differential
equation with a boundary condition which arises in pricing of Asian options.
Due to degeneracy of the partial differential operator and the nonsmooth
boundary condition, regularity of the generalized solution of such a problem
remained unclear. We prove that the generalized solution of the problem is
indeed a classical solution.

We establish existence and pointwise estimates of fundamental solutions and
Green's matrices for divergence form, second order strongly elliptic systems in
a domain $\Omega \subseteq \mathbb{R}^n$, $n \geq 3$, under the assumption that
solutions of the system satisfy De GiorgiNash type local H\"{o}lder continuity
estimates. In particular, our results apply to perturbations of diagonal
systems, and thus especially to complex perturbations of a single real
equation.

We study Green's matrices for divergence form, second order strongly elliptic
systems with bounded measurable coefficients in two dimensional domains. We
establish existence, uniqueness, and pointwise estimates of the Green's
matrices.

We study boundary blowup solutions of semilinear elliptic equations
$Lu=u_+^p$ with $p>1$, or $Lu=e^{au}$ with $a>0$, where $L$ is a second order
elliptic operator with measurable coefficients. Several uniqueness theorems and
an existence theorem are obtained.

We establish existence and various estimates of fundamental matrices and
Green's matrices for divergence form, second order strongly parabolic systems
in arbitrary cylindrical domains under the assumption that solutions of the
systems satisfy an interior H\"{o}lder continuity estimate. We present a
unified approach valid for both the scalar and the vectorial cases.