17,254 research outputs found
First order least squares method with weakly imposed boundary condition for convection dominated diffusion problems
We present and analyze a first order least squares method for convection
dominated diffusion problems, which provides robust L2 a priori error estimate
for the scalar variable even if the given data f in L2 space. The novel
theoretical approach is to rewrite the method in the framework of discontinuous
Petrov - Galerkin (DPG) method, and then show numerical stability by using a
key equation discovered by J. Gopalakrishnan and W. Qiu [Math. Comp. 83(2014),
pp. 537-552]. This new approach gives an alternative way to do numerical
analysis for least squares methods for a large class of differential equations.
We also show that the condition number of the global matrix is independent of
the diffusion coefficient. A key feature of the method is that there is no
stabilization parameter chosen empirically. In addition, Dirichlet boundary
condition is weakly imposed. Numerical experiments verify our theoretical
results and, in particular, show our way of weakly imposing Dirichlet boundary
condition is essential to the design of least squares methods - numerical
solutions on subdomains away from interior layers or boundary layers have
remarkable accuracy even on coarse meshes, which are unstructured
quasi-uniform
Krylov implicit integration factor discontinuous Galerkin methods on sparse grids for high dimensional reaction-diffusion equations
Computational costs of numerically solving multidimensional partial
differential equations (PDEs) increase significantly when the spatial
dimensions of the PDEs are high, due to large number of spatial grid points.
For multidimensional reaction-diffusion equations, stiffness of the system
provides additional challenges for achieving efficient numerical simulations.
In this paper, we propose a class of Krylov implicit integration factor (IIF)
discontinuous Galerkin (DG) methods on sparse grids to solve reaction-diffusion
equations on high spatial dimensions. The key ingredient of spatial DG
discretization is the multiwavelet bases on nested sparse grids, which can
significantly reduce the numbers of degrees of freedom. To deal with the
stiffness of the DG spatial operator in discretizing reaction-diffusion
equations, we apply the efficient IIF time discretization methods, which are a
class of exponential integrators. Krylov subspace approximations are used to
evaluate the large size matrix exponentials resulting from IIF schemes for
solving PDEs on high spatial dimensions. Stability and error analysis for the
semi-discrete scheme are performed. Numerical examples of both scalar equations
and systems in two and three spatial dimensions are provided to demonstrate the
accuracy and efficiency of the methods. The stiffness of the reaction-diffusion
equations is resolved well and large time step size computations are obtained
A Multiscale Thermo-Fluid Computational Model for a Two-Phase Cooling System
In this paper, we describe a mathematical model and a numerical simulation
method for the condenser component of a novel two-phase thermosyphon cooling
system for power electronics applications. The condenser consists of a set of
roll-bonded vertically mounted fins among which air flows by either natural or
forced convection. In order to deepen the understanding of the mechanisms that
determine the performance of the condenser and to facilitate the further
optimization of its industrial design, a multiscale approach is developed to
reduce as much as possible the complexity of the simulation code while
maintaining reasonable predictive accuracy. To this end, heat diffusion in the
fins and its convective transport in air are modeled as 2D processes while the
flow of the two-phase coolant within the fins is modeled as a 1D network of
pipes. For the numerical solution of the resulting equations, a Dual
Mixed-Finite Volume scheme with Exponential Fitting stabilization is used for
2D heat diffusion and convection while a Primal Mixed Finite Element
discretization method with upwind stabilization is used for the 1D coolant
flow. The mathematical model and the numerical method are validated through
extensive simulations of realistic device structures which prove to be in
excellent agreement with available experimental data
- …