Finite element methods developed for unfitted meshes have been widely applied to various interface problems. However, many of them resort to non-conforming spaces for approximation, which is a critical obstacle for the extension to $\textbf{H}(\text{curl})$ equations. This essential issue stems from the underlying Sobolev space $\textbf{H}^s(\text{curl};\,\Omega)$, and even the widely used penalty methodology may not yield the optimal convergence rate. One promising approach to circumvent this issue is to use a conforming test function space, which motivates us to develop a Petrov–Galerkin immersed finite element (PG-IFE) method for $\textbf{H}(\text{curl})$-elliptic interface problems. We establish the Nédélec-type IFE spaces and develop some important properties including their edge degrees of freedom, an exact sequence relating to the $H^1$ IFE space and optimal approximation capabilities. We analyse the inf-sup condition under certain assumptions and show the optimal convergence rate, which is also validated by numerical experiments.

We consider weak solutions $u:\Omega _T\to {\open R}^N$ to parabolic systems of the type

$$u_t-{\rm div}\;a(x,t,Du) = 0\quad {\rm in}\;\Omega _T = \Omega \times (0,T),$$

$x\mapsto a(x,t,\xi )$

$D_xa(\cdot ,\cdot ,\xi )$

$L^\alpha (0,T;L^\beta (\Omega ))$

$\alpha ,\beta $

$$\displaystyle{{p(n + 2)-2n} \over {2\alpha }} + \displaystyle{n \over \beta } = 1-\kappa $$

$$2 \les p < q \les p + \displaystyle{{2\kappa } \over {n + 2}}.$$

$u_t\in L^{p/(q-1)}_{{\rm loc}}(\Omega _T)$

Pinning effect on current-induced magnetic transverse domain wall dynamics in nanostrip is studied for its potential application to new magnetic memory devices. In this study, we carry out a series of calculations by solving generalized Landau-Lifshitz equation involving a current spin transfer torque in one and two dimensional models. The critical current for the transverse wall depinning in nanostrip depends on the size of artificial rectangular defects on the edges of nanostrip. We show that there is intrinsic pinning potential for a defect such that the transverse wall oscillates damply around the pinning site with an intrinsic frequency if the applied current is below critical value. The amplification of the transverse wall oscillation for both displacement and maximum value of m3 is significant by applying AC current and current pulses with appropriate frequency. We show that for given pinning potential, the oscillation amplitude as a function of the frequency of the AC current behaves like a Gaussian distribution in our numerical study, which is helpful to reduce strength of current to drive the transverse wall motion.

We prove continuity in generalized parabolic Morrey spaces of sublinear operators generated by the parabolic Calderón—Zygmund operator and by the commutator of this operator with bounded mean oscillation (BMO) functions. As a consequence, we obtain a global -regularity result for the Cauchy—Dirichlet problem for linear uniformly parabolic equations with vanishing mean oscillation (VMO) coefficients.

Weighted interior penalty discontinuous Galerkin method is developed to solve the two-dimensional non-equilibrium radiation diffusion equation on unstructured mesh. There are three weights including the arithmetic, the harmonic, and the geometric weight in the weighted discontinuous Galerkin scheme. For the time discretization, we treat the nonlinear diffusion coefficients explicitly, and apply the semi-implicit integration factor method to the nonlinear ordinary differential equations arising from discontinuous Galerkin spatial discretization. The semi-implicit integration factor method can not only avoid severe timestep limits, but also takes advantage of the local property of DG methods by which small sized nonlinear algebraic systems are solved element by element with the exact Newton iteration method. Numerical results are presented to demonstrate the validity of discontinuous Galerkin method for high nonlinear and tightly coupled radiation diffusion equation.

This article extends the finite element method of lines to a parabolic initial boundary value problem whose diffusion coefficient is discontinuous across an interface that changes with respect to time. The method presented here uses immersed finite element (IFE) functions for the discretization in spatial variables that can be carried out over a fixed mesh (such as a Cartesian mesh if desired), and this feature makes it possible to reduce the parabolic equation to a system of ordinary differential equations (ODE) through the usual semi-discretization procedure. Therefore, with a suitable choice of the ODE solver, this method can reliably and efficiently solve a parabolic moving interface problem over a fixed structured (Cartesian) mesh. Numerical examples are presented to demonstrate features of this new method.

In this paper, an immersed interface method is presented to simulate the dynamics of inextensible interfaces in an incompressible flow. The tension is introduced as an augmented variable to satisfy the constraint of interface inextensibility and the resulting augmented system is solved by the GMRES method. In this work, the arclength of the interface is locally and globally conserved as the enclosed region undergoes deformation. The forces at the interface are calculated from the configuration of the interface and the computed augmented variable, and then applied to the fluid through the related jump conditions. The governing equations are discretized on a MAC grid via a second-order finite difference scheme which incorporates jump contributions and solved by the conjugate gradient Uzawa-type method. The proposed method is applied to several examples including the deformation of a liquid capsule with inextensible interfaces in a shear flow. Numerical results reveal that both the area enclosed by interface and arclength of interface are conserved well simultaneously. These provide further evidence on the capability of the present method to simulate incompressible flows involving inextensible interfaces.

In this paper, a novel implementation of immersed interface method combined with Stokes solver on a MAC staggered grid for solving the steady two-fluid Stokes equations with interfaces. The velocity components along the interface are introduced as two augmented variables and the resulting augmented equation is then solved by the GMRES method. The augmented variables and /or the forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines and are then applied to the fluid through the jump conditions. The Stokes equations are discretized on a staggered Cartesian grid via a second order finite difference method and solved by the conjugate gradient Uzawa-type method. The numerical results show that the overall scheme is second order accurate. The major advantages of the present IIM-Stokes solver are the efficiency and flexibility in terms of types of fluid flow and different boundary conditions. The proposed method avoids solution of the pressure Poisson equation, and comparisons are made to show the advantages of time savings by the present method. The generalized two-phase Stokes solver with correction terms has also been applied to incompressible two-phase Navier-Stokes flow.

The immersed interface method is modified to compute Schrödinger equation with discontinuous potential. By building the jump conditions of the solution into the finite difference approximation near the interface, this method can give at least second order convergence rate for the numerical solution on uniform cartesian grids. The accuracy of this algorithm is tested via several numerical examples.

In this paper we establish a priori${{h}^{1}}$-estimates in a bounded domain for parabolic equations with vanishing $\text{LMO}$ coefficients.