• 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:1349--1367, 1995)
  • We prove the uniqueness of H\"older continuous weak solutions via duality argument and vanishing viscosity method for the Keller-Segel 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 second-order parabolic systems of divergence form with bounded and measurable leading coefficients and divergence free first-order 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 non-divergence 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 non-negative adjoint solutions to non-divergence form elliptic equations, when the leading coefficients have Dini mean oscillation.
  • Under various conditions, we establish Schauder estimates for both divergence and non-divergence form second-order elliptic and parabolic equations involving H\"older semi-norms 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 non-divergence form equations are also established.
  • Fluid-solid 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 non-convergence 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 non-local 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 photo-acoustic 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 Robin-type boundary condition. We construct Green's function for parabolic systems with time-dependent coefficients satisfying a possibly nonlocal Robin-type 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 Robin-type 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 non-tangential 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 time-dependent 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 Giorgi-Moser-Nash 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 second-order linear parabolic operators in non-divergence form that are intrinsically defined on Riemannian manifolds. In the elliptic case, Cabr\'e proved a global Krylov-Safonov 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 Krylov-Safonov Harnack inequality for nonnegative solutions of the non-divergent parabolic equation. This, in particular, gives a new proof for Li-Yau 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 photo-acoustic 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 non-smooth time-varying domains whose boundaries are locally represented as graph of functions that are Lipschitz continuous in the spatial variables and 1/2-H\"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 time-varying 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 Giorgi-Moser-Nash 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 non-divergence form second-order elliptic and parabolic equations involving H\"older semi-norms 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 non-smooth 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 Giorgi-Nash 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 blow-up 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.