This paper studies the magneto-heat coupling model which describes iron loss of conductors and energy exchange between magnetic field and Ohmic heat. The temperature influences Maxwell's equations through the variation of electric conductivity, while electric eddy current density provides the heat equation with Ohmic heat source. It is in this way that Maxwell's equations and the heat equation are coupled together. The system also incorporates the heat exchange between conductors and cooling oil which is poured into and out of the transformer. We propose a weak formulation for the coupling model and establish the well-posedness of the problem. The model is more realistic than the traditional eddy current model in numerical simulations for large power transformers. The theoretical analysis of this paper paves a way for us to design efficient numerical computation of the transformer in the future.

We introduce a third order adaptive mesh method to arbitrary high order Godunov approach. Our adaptive mesh method consists of two parts, i.e., mesh-redistribution algorithm and solution algorithm. The mesh-redistribution algorithm is derived based on variational approach, while a new solution algorithm is developed to preserve high order numerical accuracy well. The feature of proposed Adaptive ADER scheme includes that 1). all simulations in this paper are stable for large CFL number, 2). third order convergence of the numerical solutions is successfully observed with adaptive mesh method, and 3). high resolution and non-oscillatory numerical solutions are obtained successfully when there are shocks in the solution. A variety of numerical examples show the feature well.

Propagation at a finite speed is established for non-negative weak solutions to a thin-film approximation of the two-phase Muskat problem. The expansion rate of the support matches the scale invariance of the system. Moreover, we determine sufficient conditions on the initial data for the occurrence of waiting time phenomena.

Ill-posedness for the compressible Navier–Stokes equations has been proved by Chen et al. [On the ill-posedness of the compressible Navier–Stokes equations in the critical Besov spaces, Revista Mat. Iberoam.31 (2015), 1375–1402] in critical Besov space $L^{p}$$(p>6)$ framework. In this paper, we prove ill-posedness with the initial data satisfying

$$\begin{eqnarray}\displaystyle \Vert \unicode[STIX]{x1D70C}_{0}-\bar{\unicode[STIX]{x1D70C}}\Vert _{{\dot{B}}_{p,1}^{\frac{3}{p}}}\leqslant \unicode[STIX]{x1D6FF},\quad \Vert u_{0}\Vert _{{\dot{B}}_{6,1}^{-\frac{1}{2}}}\leqslant \unicode[STIX]{x1D6FF}. & & \displaystyle \nonumber\end{eqnarray}$$

$L(a)\unicode[STIX]{x1D6E5}u$

$u\cdot \unicode[STIX]{x1D6FB}u$

In this work, we investigate the numerical approximation of the compressible Navier-Stokes equations under the framework of discontinuous Galerkin methods. For discretization of the viscous and heat fluxes, we extend and apply the symmetric direct discontinuous Galerkin (SDDG) method which is originally introduced for scalar diffusion problems. The original compressible Navier-Stokes equations are rewritten into an equivalent form via homogeneity tensors. Then, the numerical diffusive fluxes are constructed from the weak formulation of primal equations directly without converting the second-order equations to a first-order system. Additional numerical flux functions involving the jump of second order derivative of test functions are added to the original direct discontinuous Galerkin (DDG) discretization. A number of numerical tests are carried out to assess the practical performance of the SDDG method for the two dimensional compressible Navier-Stokes equations. These numerical results obtained demonstrate that the SDDG method can achieve the optimal order of accuracy. Especially, compared with the well-established symmetric interior penalty (SIP) method [18], the SDDG method can maintain the expected optimal order of convergence with a smaller penalty coefficient.

This paper describes an adaptive preconditioner for numerical continuation of incompressible Navier–Stokes flows based on Stokes preconditioning [42] which has been used successfully in studies of pattern formation in convection. The preconditioner takes the form of the Helmholtz operator I–ΔtL which maps the identity (no preconditioner) for Δt≪1 to Laplacian preconditioning for Δt≫1. It is built on a first order Euler time-discretization scheme and is part of the family of matrix-free methods. The preconditioner is tested on two fluid configurations: three-dimensional doubly diffusive convection and a two-dimensional projection of a shear flow. In the former case, it is found that Stokes preconditioning is more efficient for , away from the values used in the literature. In the latter case, the simple use of the preconditioner is not sufficient and it is necessary to split the system of equations into two subsystems which are solved simultaneously using two different preconditioners, one of which is parameter dependent. Due to the nature of these applications and the flexibility of the approach described, this preconditioner is expected to help in a wide range of applications.

A full multigrid method with coarsening by a factor-of-three to distributed control problems constrained by Stokes equations is presented. An optimal control problem with cost functional of velocity and/or pressure tracking-type is considered with Dirichlet boundary conditions. The optimality system that results from a Lagrange multiplier framework, form a linear system connecting the state, adjoint, and control variables. We investigate multigrid methods with finite difference discretization on staggered grids. A coarsening by a factor-of-three is used on staggered grids that results nested hierarchy of staggered grids and simplified the inter-grid transfer operators. A distributive-Gauss-Seidel smoothing scheme is employed to update the state- and adjoint-variables and a gradient update step is used to update the control variables. Numerical experiments are presented to demonstrate the effectiveness and efficiency of the proposed multigrid framework to tracking-type optimal control problems.

This paper is concerned with the invisibility cloaking in acoustic wave scattering from a new perspective. We are especially interested in achieving the invisibility cloaking by completely regular and isotropic mediums. It is shown that an interior transmission eigenvalue problem arises in our study, which is the one considered theoretically in Cakoni et al. (Transmission eigenvalues for inhomogeneous media containing obstacles, Inverse Problems and Imaging, 6 (2012), 373–398). Based on such an observation, we propose a cloaking scheme that takes a three-layer structure including a cloaked region, a lossy layer and a cloaking shell. The target medium in the cloaked region can be arbitrary but regular, whereas the mediums in the lossy layer and the cloaking shell are both regular and isotropic. We establish that if a certain non-transparency condition is satisfied, then there exists an infinite set of incident waves such that the cloaking device is nearly invisible under the corresponding wave interrogation. The set of waves is generated from the Herglotz approximation of the associated interior transmission eigenfunctions. We provide both theoretical and numerical justifications.

Incorrect propagation speed of discontinuities may occur by straightforward application of standard dissipative schemes for problems that contain stiff source terms for underresolved grids even for time steps within the CFL condition. By examining the dissipative discretized counterpart of the Euler equations for a detonation problem that consists of a single reaction, detailed analysis on the spurious wave pattern is presented employing the fractional step method, which utilizes the Strang splitting. With the help of physical arguments, a threshold values method (TVM), which can be extended to more complicated stiff problems, is developed to eliminate the wrong shock speed phenomena. Several single reaction detonations as well as multispecies and multi-reaction detonation test cases with strong stiffness are examined to illustrate the performance of the TVM approach.

The lobe dynamics andmass transport between separation bubble and main flow in flow over airfoil are studied in detail, using Lagrangian coherent structures (LCSs), in order to understand the nature of evolution of the separation bubble. For this problem, the transient flow over NACA0012 airfoil with low Reynolds number is simulated numerically by characteristic based split (CBS) scheme, in combination with dual time stepping. Then, LCSs and lobe dynamics are introduced and developed to investigate themass transport between separation bubble and main flow, from viewpoint of nonlinear dynamics. The results show that stable manifolds and unstable manifolds could be tangled with each other as time evolution, and the lobes are formed periodically to induce mass transport between main flow and separation bubble, with dynamic behaviors. Moreover, the evolution of the separation bubble depends essentially on the mass transport which is induced by lobes, ensuing energy and momentum transfers. As the results, it can be drawn that the dynamics of flow separation could be studied using LCSs and lobe dynamics, and could be controlled feasibly if an appropriate control is applied to the upstream boundary layer with high momentum.

We propose a stochastic Galerkin method using sparse wavelet bases for the Boltzmann equation with multi-dimensional random inputs. Themethod uses locally supported piecewise polynomials as an orthonormal basis of the random space. By a sparse approach, only a moderate number of basis functions is required to achieve good accuracy in multi-dimensional random spaces. We discover a sparse structure of a set of basis-related coefficients, which allows us to accelerate the computation of the collision operator. Regularity of the solution of the Boltzmann equation in the random space and an accuracy result of the stochastic Galerkin method are proved in multi-dimensional cases. The efficiency of the method is illustrated by numerical examples with uncertainties from the initial data, boundary data and collision kernel.

We analyze in this paper the pressure splitting scheme of a partitioned semi-implicit coupling algorithm for fluid-structure interaction (FSI) simulation. The semi-implicit coupling algorithm is developed on the ground of the artificial compressibility characteristic-based split (AC-CBS) scheme that serves not only for the fluid subsystem but also for the global FSI system. As the dual-time stepping procedure recommended for quasi-incompressible flows is incorporated into the implicit coupling stage, the fluctuating pressure may be unusually susceptible to the AC coefficient. Moreover, it is not trivial to devise an optimal AC formulation for pressure estimation. Instead, we consider a stabilized second-order pressure splitting scheme in the AC-CBS-based partitioned semi-implicit coupling algorithm. Computer simulation of a benchmark FSI experiment demonstrates that good agreement is exposed between the available and present data.

In this work we discuss the numerical discretization of the time-dependent Maxwell's equations using a fully explicit leap-frog type discontinuous Galerkin method. We present a sufficient condition for the stability and error estimates, for cases of typical boundary conditions, either perfect electric, perfect magnetic or first order Silver-Müller. The bounds of the stability region point out the influence of not only the mesh size but also the dependence on the choice of the numerical flux and the degree of the polynomials used in the construction of the finite element space, making possible to balance accuracy and computational efficiency. In the model we consider heterogeneous anisotropic permittivity tensors which arise naturally in many applications of interest. Numerical results supporting the analysis are provided.

This paper is concerned with the construction of global, non-vacuum, strong, large amplitude solutions to initial–boundary-value problems for the one-dimensional compressible Navier–Stokes equations with degenerate transport coefficients. Our analysis derives the positive lower and upper bounds on the specific volume and the absolute temperature.

In this study an explicit Finite Difference Method (FDM) based scheme is developed to solve the Maxwell's equations in time domain for a lossless medium. This manuscript focuses on two unique aspects – the three dimensional time-accurate discretization of the hyperbolic system of Maxwell equations in three-point non-staggered grid stencil and it's application to parallel computing through the use of Graphics Processing Units (GPU). The proposed temporal scheme is symplectic, thus permitting conservation of all Hamiltonians in the Maxwell equation. Moreover, to enable accurate predictions over large time frames, a phase velocity preserving scheme is developed for treatment of the spatial derivative terms. As a result, the chosen time increment and grid spacing can be optimally coupled. An additional theoretical investigation into this pairing is also shown. Finally, the application of the proposed scheme to parallel computing using one Nvidia K20 Tesla GPU card is demonstrated. For the benchmarks performed, the parallel speedup when compared to a single core of an Intel i7-4820K CPU is approximately 190x.

Pressure-correction projection finite element methods (FEMs) are proposed to solve nonstationary natural convection problems in this paper. The first-order and second-order backward difference formulas are applied for time derivative, the stability analysis and error estimates of the semi-discrete schemes are presented using energy method. Compared with characteristic variational multiscale FEM, pressure-correction projection FEMs are more efficient and unconditionally energy stable. Ample numerical results are presented to demonstrate the effectiveness of the pressure-correction projection FEMs for solving these problems.

Moment models are often used for the solution of kinetic equations such as the Boltzmann equation. Unfortunately, standard models like Grad's equations are not hyperbolic and can lead to nonphysical solutions. Newly derived moment models like the Hyperbolic Moment Equations and the Quadrature-Based Moment Equations yield globally hyperbolic equations but are given in partially conservative form that cannot be written as a conservative system.

In this paper we investigate the applicability of different dedicated numerical schemes to solve the partially conservative model equations. Caused by the non-conservative type of equation we obtain differences in the numerical solutions, but due to the structure of the moment systems we show that these effects are very small for standard simulation cases. After successful identification of useful numerical settings we show a convergence study for a shock tube problem and compare the results to a discrete velocity solution. The results are in good agreement with the reference solution and we see convergence considering an increasing number of moments.

In this paper, a novel second-order two-scale (SOTS) computational method is developed for nonlinear dynamic thermo-mechanical problems of composites with cylindrical periodicity. The non-linearities of these multi-scale problems were caused by the temperature-dependent properties of the composites. Firstly, the formal SOTS solutions for these problems are constructed by the multiscale asymptotic analysis. Then we theoretically explain the importance of the SOTS solutions by the error analysis in the pointwise sense. In addition, a SOTS numerical algorithm is proposed in detail to effectively solve these problems. Finally, some numerical examples verify the feasibility and effectiveness of the SOTS numerical algorithm we proposed.

We present a numerical method to solve the Vlasov-Poisson-Fokker-Planck (VPFP) system using the NRxx method proposed in [4, 7, 9]. A globally hyperbolic moment system similar to that in [23] is derived. In this system, the Fokker-Planck (FP) operator term is reduced into the linear combination of the moment coefficients, which can be solved analytically under proper truncation. The non-splitting method, which can keep mass conservation and the balance law of the total momentum, is used to solve the whole system. A numerical problem for the VPFP system with an analytic solution is presented to indicate the spectral convergence with the moment number and the linear convergence with the grid size. Two more numerical experiments are tested to demonstrate the stability and accuracy of the NRxx method when applied to the VPFP system.

The computational work and storage of numerically solving the time fractional PDEs are generally huge for the traditional direct methods since they require total memory and work, where NT and NS represent the total number of time steps and grid points in space, respectively. To overcome this difficulty, we present an efficient algorithm for the evaluation of the Caputo fractional derivative of order α∈(0,1). The algorithm is based on an efficient sum-of-exponentials (SOE) approximation for the kernel t–1–α on the interval [Δt, T] with a uniform absolute error ε. We give the theoretical analysis to show that the number of exponentials Nexp needed is of order for T≫1 or for TH1 for fixed accuracy ε. The resulting algorithm requires only storage and work when numerically solving the time fractional PDEs. Furthermore, we also give the stability and error analysis of the new scheme, and present several numerical examples to demonstrate the performance of our scheme.