
We consider wellposedness of the boundary value problem in presence of an
inclusion with complex conductivity $k$. We first consider the transmission
problem in $\mathbb{R}^d$ and characterize solvability of the problem in terms
of the spectrum of the NeumannPoincar\'e operator. We then deal with the
boundary value problem and show that the solution is bounded in its $H^1$norm
uniformly in $k$ as long as $k$ is at some distance from a closed interval in
the negative real axis. We then show with an estimate that the solution depends
on $k$ in its $H^1$norm Lipschitz continuously. We finally show that the
boundary perturbation formula in presence of a diametrically small inclusion is
valid uniformly in $k$ away from the closed interval mentioned before. The
results for the single inclusion case are extended to the case when there are
multiple inclusions with different complex conductivities: We first obtain a
complete characterization of solvability when inclusions consist of two
disjoint disks and then prove solvability and uniform estimates when imaginary
parts of conductivities have the same signs. The results are obtained using the
spectral property of the associated NeumannPoincar\'e operator and the
spectral resolution.

The AllenCahn equation is solved numerically by operator splitting Fourier
spectral methods. The basic idea of the operator splitting method is to
decompose the original problem into subequations and compose the approximate
solution of the original equation using the solutions of the subproblems.
Unlike the first and the second order methods, each of the heat and the
freeenergy evolution operators has at least one backward evaluation in higher
order methods. We investigate the effect of negative time steps on a general
form of third order schemes and suggest three third order methods for better
stability and accuracy. Two fourth order methods are also presented. The
traveling wave solution and a spinodal decomposition problem are used to
demonstrate numerical properties and the order of convergence of the proposed
methods.

The FENE dumbbell model consists of the incompressible NavierStokes equation
and the FokkerPlanck equation for the polymer distribution. In such a model,
the polymer elongation cannot exceed a limit $\sqrt{b}$, yielding all
interesting features near the boundary. In this paper we establish the local
wellposedness for the FENE dumbbell model under a class of Dirichlettype
boundary conditions dictated by the parameter $b$. As a result, for each $b>0$
we identify a sharp boundary requirement for the underlying density
distribution, while the sharpness follows from the existence result for each
specification of the boundary behavior. It is shown that the probability
density governed by the FokkerPlanck equation approaches zero near boundary,
necessarily faster than the distance function $d$ for $b>2$, faster than $dln
d$ for $b=2$, and as fast as $d^{b/2}$ for $0<b<2$. Moreover, the sharp
boundary requirement for $b\geq 2$ is also sufficient for the distribution to
remain a probability density.

We prove global wellposedness for the microscopic FENE model under a sharp
boundary requirement. The wellposedness of the FENE model that consists of the
incompressible NavierStokes equation and the FokkerPlanck equation has been
studied intensively, mostly with the zero flux boundary condition. Recently it
was illustrated by C. Liu and H. Liu [2008, SIAM J. Appl. Math.,
68(5):13041315] that any preassigned boundary value of a weighted
distribution will become redundant once the nondimensional parameter $b>2$. In
this article, we show that for the wellposedness of the microscopic FENE model
($b>2$) the least boundary requirement is that the distribution near boundary
needs to approach zero faster than the distance function. Under this condition,
it is shown that there exists a unique weak solution in a weighted Sobolev
space. Moreover, such a condition still ensures that the distribution is a
probability density. The sharpness of this boundary requirement is shown by a
construction of infinitely many solutions when the distribution approaches zero
as fast as the distance function.