The problem of compressive-sensing (CS) L2-L1-TV reconstruc- tion of magnetic resonance (MR) scans from undersampled k-space data has been addressed in numerous studies. However, the regularization parameters in models of CS L2-L1-TV reconstruction are rarely studied. Once the regu- larization parameters are given, the solution for an MR reconstruction model is xed and is less eective in the case of strong noise. To overcome this shortcoming, we present a new alternating formulation to replace the standard L2-L1-TV reconstruction model. A weighted-average alternating minimization method is proposed based on this new formulation and a convergence analysis of the method is carried out. The advantages of and the motivation for the pro- posed alternating formulation are explained. Experimental results demonstrate that the proposed formulation yields better reconstruction results in the case of strong noise and can improve image reconstruction via exible parameter slection.

In this paper, we consider local multiscale model reduction for problems with multiple scales in space and time. We developed our approaches within the framework of the Generalized Multiscale Finite Element Method (GMsFEM) using spacetime coarse cells. The main idea of GMsFEM is to construct a local snapshot space and a local spectral decomposition in the snapshot space. Previous research in developing multiscale spaces within GMsFEM focused on constructing multiscale spaces and relevant ingredients in space only. In this paper, our main objective is to develop a multiscale model reduction framework within GMsFEM that uses spacetime coarse cells. We construct spacetime snapshot and offline spaces. We compute these snapshot solutions by solving local problems. A complete snapshot space will use all possible boundary conditions; however, this can be very expensive. We propose using

nthispaperwewillpresentthelocalstabilityanalysisandlocalerrorestimatefor the local discontinuous Galerkin (LDG) method, when solving the time-dependent singularly perturbed problems in one dimensional space with a stationary outflow boundary layer. Based on a general framework on the local stability, we obtain the optimal error estimate out of the local subdomain, which is nearby the outflow boundary point and has the width of O(h log(1/h)), for the semi-discrete LDG scheme and the fully-discrete LDG scheme with the second order explicit Runge–Kutta time-marching. Here h is the maximum mesh length. The numerical experiments are given also.

In this paper, a novel penalty method based on the immersed boundary formulation is proposed for simulating the transient Stokes flow with an inextensible interface enclosing a suspended solid particle. The main idea of this approach relies on the penalty techniques by modifying the constitutive equation of Stokes flow to weaken the incompressibility condition, relating the surface divergence to the elastic tension $\sigma$ to relax the interface's inextensibility, and connecting the particle surface-velocity with the particle surface force ${\mbf F}$ to regularize the particle's rigid motion. The advantage of these regularized governing equations is that when they are discretized by the standard centered difference scheme on a staggered grid, the resulting linear system can easily be reduced by eliminating the unknowns $p_h$, $\sigma_h$ and ${\mbf F}_h$ directly, so that we just need to solve a smaller linear system of the velocity approximation ${\mbf u}_h$. This advantage is preserved and even enhanced when such approach is applied to the transient Stokes flow with multiple compound vesicles. Moreover, this smaller linear system is symmetric and negative-definite, which enables us to use efficient linear solvers. Another important feature of the proposed method is that the discretization scheme is unconditionally stable in the sense that an appropriately defined energy functional associated with the discrete system is decreasing and hence bounded in time. We numerically test the accuracy and stability of the immersed boundary discretization scheme. The tank-treading and tumbling motions of inextensible interface with a suspended solid particle in the simple shear flow will be studied extensively. The simulation of the motion of multiple compound vesicles will be performed as well. Numerical results illustrate the superior performance of the proposed penalty method.

Numerical simulations of waves in highly heterogeneous media have important applications, but direct computations are prohibitively expensive. In this paper, we develop a new generalized multiscale finite element method with the aim of simulating waves at a much lower cost. Our method is based on a mixed Galerkin type method with carefully designed basis functions that can capture various scales in the solution. The basis functions are constructed based on some local snapshot spaces and local spectral problems defined on them. The spectral problems give a natural ordering of the basis functions in the snapshot space and allow systematically enrichment of basis functions. In addition, by using a staggered coarse mesh, our method is energy conserving and has block diagonal mass matrix, which are desirable properties for wave propagation. We will prove that our method has spectral convergence, and