We study a nonlinear weighted least-squares finite element method for the Navier- Stokes equations governing non-Newtonian fluid flows by using the Carreau- Yasuda model. The Carreau-Yasuda model is used to describe the shear-thinning behavior of blood. We prove that the least-squares approximation converges to linearized solutions of the non-Newtonian model at the optimal rate. By using continuous piecewise linear finite element spaces for all variables and by appropriately adjusting the nonlinear weighting function, we obtain optimal L2-norm error convergence rates in all variables. Numerical results are given for a Carreau fluid in the 4-to-1 contraction problem, revealing the shear-thinning behavior. The physical parameter effects are also investigated.

We implemented an adaptively refined least-squares finite element approach for the Navier-Stokes equations that govern generalized Newtonian fluid flows using the Carreau model. To capture the flow region, we developed an adaptive mesh refinement approach based on the least-squares method. The generated refined grids agree well with the physical attributes of the flows. We also proved that the least-squares approximation converges to the linearized versions solutions of the Carreau model at the best possible rate. Model problems considered in the study are the flow past a planar channel and 4-to-1 contraction problems. We presented the numerical results of the model problems, revealing the efficiency of the proposed scheme, and investigated the physical parameter effects.

We consider the numerical solution of the time-fractional diffusion-wave equation on two-dimensional and three-dimensional unbounded spatial domains. Introduce an artificial boundary and find the exact and approximate artificial boundary conditions for the given problem, which lead to a bounded computational domain. Using the exact or approximating boundary conditions on the artificial boundary, the original problem is reduced to an initial-boundary-value problem on the bounded computational domain which is respectively equivalent to or approximates the original problem. Finite difference methods are used to solve the reduced problems on the bounded computational domain and the stability of these finite difference methods are proved. The numerical results demonstrate that the method given in this paper is effective and feasible.

A time-dependent Hermite-Galerkin spectral method (THGSM) is investigated in this paper for the nonlinear convection-diffusion equations in the unbounded domains. The time-dependent scaling factor and translating factor are introduced in the definition of the generalized Hermite functions (GHF). As a consequence, the THGSM based on these GHF has many advantages, not only in theoretical proofs, but also in numerical implementations. The stability and spectral convergence of our proposed method have been established in this paper. The Korteweg-de Vries-Burgers (KdVB) equation and its special cases, including the heat equation and the Burgers' equation, as the examples, have been numerically solved by our method. The numerical results are presented, and it surpasses the existing methods in accuracy. Our theoretical proof of the spectral convergence has been supported by the numerical results.

Global geometrical optics method is a new semi-classical approach for the high frequency linear waves proposed by the author in [13]. In this paper, we rederive it in a more concise way. It is shown that the right candidate of solution ansatz for the high frequency wave equations is the extended WKB function, other than the WKB function used in the classical geometrical optics approximation. A new and main contribution of this paper is an interface analysis for the Helmholtz equation when the incident wave is of extended WKB-type. We derive asymptotic expressions for the reflected and/or transmitted propagating waves in the general case. These expressions are valid even when the incident rays include caustic points.