A novel discontinuous Galerkin (DG) method is developed to solve timedependent bi-harmonic type equations involving fourth derivatives in one and multiple space dimensions. We present the spatial DG discretization based on a mixed formulation and central interface numerical fluxes so that the resulting semi-discrete schemes are L^2 stable even without interior penalty. For time discretization, we use Crank–Nicolson so that the resulting scheme is unconditionally stable and second order in time. We present the optimal L^2 error estimate of O(h^{k+1}) for polynomials of degree k for semi-discrete DG schemes, and the L2 error of O(h^{k+1} + (\Delta t)^2) for fully discrete DG schemes. Extensions to more general fourth order partial differential equations and cases with non-homogeneous boundary conditions are provided. Numerical results are presented to verify the stability and accuracy of the schemes. Finally, an application to the one-dimensional Swift–Hohenberg equation endowed with a decay free energy is presented.

In this paper, we investigate the boundary layer behavior of solutions to the one dimensional Broadwell model of the nonlinear Boltzmann equation for small mean free path. We consider the analogue of Maxwell's diffusive and the reflexive boundary conditions. It is found that even for such a simple model, there are boundary layers due to purely kinetic effects which cannot be detected by the corresponding Navier-Stokes system. It is also found numerically that a compressive boundary layer is not always stable in the sense that it may detach from the boundary and move into the interior of the gas as a shock layer.

In this paper, we study existence of global entropy weak solutions to a critical-case unstable thin lm equation in one-dimensional case ht + @x(h n @xxxh) + @x(h n+2@xh) = 0; where n  1. There exists a critical mass Mc = 2 p 6 3 found by Witelski etal. (2004 Euro. J. of Appl. Math. 15, 223{256) for n = 1. We obtain global existence of a non-negative entropy weak solution if initial mass is less than Mc. For n  4, entropy weak solutions are positive and unique. For n = 1,a nite time blow-up occurs for solutions with initial mass larger than Mc.For the Cauchy problem with n = 1 and initial mass less than Mc, we show that at least one of the following long-time behavior holds: the second moment goes to innity as the time goes to innity or h(
; tk) * 0 in L1 (R) for some subsequence tk ! 1.

The global well-posedness of the Boltzmann equation with initial data of large amplitude has remained a long-standing open problem. In this paper, by developing a new L x L v 1 L x , v approach, we prove the global existence and uniqueness of mild solutions to the Boltzmann equation in the whole space or torus for a class of initial data with bounded velocity-weighted L x L v 1 L x , v norm under some smallness condition on the L x L v 1 L x , v norm as well as defect mass, energy and entropy so that the initial data allow large amplitude oscillations. Both the hard and soft potentials with angular cut-off are considered, and the large time behavior of solutions in the L x L v 1 L x , v norm with explicit rates of convergence are also studied.

The focus of this work is on rigorous mathematical analysis of the topological derivative based detection algorithms for the localization of an elastic inclusion of vanishing characteristic size. A filtered quadratic misfit is considered, and the performance of the topological derivative imaging functional resulting therefrom is analyzed. Our analysis reveals that the imaging functional may not attain its maximum at the location of the inclusion. Moreover, the resolution of the image is below the diffraction limit. Both phenomena are due to the coupling of pressure and shear waves propagating with different wave speeds and polarization directions. A novel imaging functional based on the weighted Helmholtz decomposition of the topological derivative is, therefore, introduced. It is thereby substantiated that the maximum of the imaging functional is attained at the location of the inclusion and the resolution is enhanced and proves to be the diffraction limit. Finally, we investigate the stability of the proposed imaging functionals with respect to measurement and medium noises.