A new residual-based stabilized finite element method is analyzed for solving the Stokes equations in terms of velocity and pressure, where the $H^{−1}$ norm is introduced in the measurement of the residuals to obtain a symmetric positive definite (SPD) method. The $H^{−1}$ norm is computable and can be always easily realized offline by the continuous linear finite element solution or the preconditioned counterpart of the Poisson Dirichlet problem. Although the $H^{−1}$ norm is computed in the linear element space, no matter what the finite element spaces for the velocity and the pressure are, optimal error bounds can be established when using continuous finite element pairs $R_l$−$R_m$ for velocity and pressure for any $l, m\ge 1$. Numerical experiments are performed to confirm the theoretical results obtained.

In this paper, we design and analyze third order positivity-preserving discontinuous Galerkin (DG) schemes for solving the time-dependent system of Poisson–Nernst–Planck (PNP) equations, which have found much use in diverse applications. Our DG method with Euler forward time discretization is shown to preserve the positivity of cell averages at all time steps. The positivity of numerical solutions is then restored by a scaling limiter in reference to positive weighted cell averages. The method is also shown to preserve steady states. Numerical examples are presented to demonstrate the third order accuracy and illustrate the positivity-preserving property in both one and two dimensions.

In this paper, we extend the high order finite-difference method with subcell resolution (SR) for two-species stiff one-reaction models to multispecies and multireaction chemical reactive flows, which are significantly more difficult because of the multiple scales generated by different reactions. For reaction problems, when the reaction time scale is very small, the reaction zone scale is also small and the governing equations become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present SR method for reactive Euler system is a fractional step method. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with certain computed flow variables in the shock region modified by the Harten subcell resolution idea. Several numerical examples of multispecies and multireaction reactive flows are performed in both one and two dimensions. Studies demonstrate that the SR method can capture the correct propagation speed of discontinuities in very coarse meshes.

Differential optical absorption spectroscopy (DOAS) is a powerful tool for detecting and quantifying trace gases in atmospheric chemistry\cite {Platt_Stutz08}. DOAS spectra consist of a linear combination of complex multi-peak multi-scale structures. Most DOAS analysis routines in use today are based on least squares techniques, for example, the approach developed in the 1970s uses polynomial fits to remove a slowly varying background, and known reference spectra to retrieve the identity and concentrations of reference gases. An open problem is to identify unknown gases in the fitting residuals for complex atmospheric mixtures.

Atmospheric temperature is one of the most important climate variables. This observational study presents detailed descriptions of the temperature variability imprinted in the 9-year brightness temperature data acquired by the Advanced Microwave Sounding Unit- Instrument A (AMSU-A) aboard Aqua since September 2002 over tropical oceans. A non-linear, adaptive method called the Ensemble Joint Multiple Extraction has been employed to extract the principal modes of variability in the AMSU-A/Aqua data. The semi-annual, annual, quasibiennial oscillation (QBO) modes and QBO–annual beat in the troposphere and the stratosphere have been successfully recovered. The modulation by the El Nin˜o/Southern oscillation (ENSO) in the troposphere was found and correlates well with the Multivariate ENSO Index. The longterm variations during 2002–2011 reveal a cooling trend (-0.5 K/decade at 10 hPa) in the tropical stratosphere; the trend below the tropical tropopause is not statistically significant due to the length of our data. A new tropospheric near-annual mode (period *1.6 years) was also revealed in the troposphere, whose existence was confirmed using National Centers for Environmental Prediction Reanalysis air temperature data. The near-annual mode in the troposphere is found to prevail in the eastern Pacific region and is coherent with a near-annual mode in the observed sea surface temperature over the Warm Pool region that has previously been reported. It remains a challenge for climate models to simulate the trends and principal modes of natural variability reported in this work.