Space–time integrated least squares: a time‐marching approach

A time‐marching formulation is derived from the space–time integrated least squares (STILS) method for solving a pure hyperbolic convection equation and is numerically compared to various known methods. Copyright © 2004 John Wiley & Sons, Ltd.     

A numerical study of the target system of an ADSS with different flow guides

The mechanical design of the target module of an accelerator driven subcritical nuclear reactor system (ADSS) calls for an analysis of the related thermal-hydraulic issues because of large amount of heat deposition in the spallation region during the course of nuclear interactions with the molten lead bismuth eutectic (LBE) target. The LBE also should carry the entire heat generated as a consequence of the spallation reaction. The problem of heat removal by the LBE is a challenging thermal-hydraulic issue. For this, one has to examine the flows of low Prandtl number fluids (LBE) in a complex ADSS geometry. In this study, the equations governing the laminar flow and thermal energy are solved numerically using the streamline upwind Petrov-Galerkin (SUPG) finite element (FE) method. The target systems with a straight and a nozzle guide have been considered. The principal purpose of the analysis is to trace the flow and temperature distribution and thereby to check the suitability of the flow guide in avoiding the recirculation or stagnation zones in the flow space that may lead to hot spots.      

类升力体外形高超声速粘性绕流拓扑结构分析

采用交替方向隐式分解的隐式NND格式求解全N-S方程模拟“类升力体”外形在高超声速下的大攻角流动,给出了“类升力体”外形表面极限流线随攻角变化的拓扑结构及40°攻角下垂直于体轴的横截面流线拓扑结构。结果表明:类升力体外形三维流场结构十分复杂,攻角从0°～40°变化时,背风面表面极限流线依次由不分离、开式分离向起始于鞍、结点组合的高阶奇点的分离方式转化,翼面横向分离亦随攻角增大而增大;垂直于体轴的横截面流动拓扑结构与文献[2]给出的理论分析一致。     

The topology of the sectional streamlines

The unsteady, three-dimensional full Navier-Stokes equations are solved using a Beam-Warming implicit algorithm in this paper. Computations of the flow over a 76° sweep delta wing at 36.5° angle of attack is presented. The sectional streamlines are depicted and the evolution of the instantaneous crossflow topology of the leading-edge vortex is analyzed. It is found that, along the axis, the topology of the primary vortex alters several times starting from stable focus near the apex to unstable focus, and lasts back to stable focus near wake edge; The stable limit cycle and unstable limit cycle are shown in this evolution. These various altering topologies stem from the stretching and compression of the vortex core.     

非线性对流扩散问题的交替方向差分流线扩散法

1引言Peaceman,Douglas等人于1955年提出了差分格式的交替方向法。随后,Douglas,Dupont于1972年又提出了有限元格式的交替方向法[1]。其基本思想是:对两个或三个空间变量的二阶抛物型和双曲型问题,将交替方向法与Galerkin方法相结合,通过算子分裂技术,把高维问题转化为一系列低维问题,交替地沿各空间变量的方向求解。[2]、[3]和[4]给出了对更一般扩散问题(带对流项的抛物方程)的数值求解和误差分析。     

On streamline diffusion arising in Galerkin FEM with predictor/multi‐corrector time integration

In the present paper, the author shows that the predictor/multi‐corrector (PMC) time integration for the advection–diffusion equations induces numerical diffusivity acting only in the streamline direction, even though the equations are spatially discretized by the conventional Galerkin finite element method (GFEM). The transient 2‐D and 3‐D advection problems are solved with the PMC scheme using both the GFEM and the streamline upwind/Petrov Galerkin (SUPG) as the spatial discretization methods for comparison. The solutions of the SUPG‐PMC turned out to be overly diffusive due to the additional PMC streamline diffusion, while the solutions of the GFEM‐PMC were comparatively accurate without significant damping and phase error. A similar tendency was seen also in the quasi‐steady solutions to the incompressible viscous flow problems: 2‐D driven cavity flow and natural convection in a square cavity. Copyright © 2002 John Wiley & Sons, Ltd.     

三维热传导型半导体器件瞬态模拟问题的Crank-Nicolson差分-流线扩散有限元法及其数值分析

本文研究三维热传导型半导体器件瞬态模拟问题的数值方法.针对数学模型中各方程不同的特点,分别提出不同的有限元格式.特别针对浓度方程组是对流为主扩散问题的特点,使用Crank-Nicolson差分-流线扩散计算格式,提高了数值解的稳定性.得到的L2误差估计关于空间剖分步长是拟最优的,关于时间步长具有二阶精度.     

A note on the calculation of strain histories in orthogonal streamline coordinate systems

An extension of a previous work concerning the calculation of strain histories along streamlines is made to get more complete and useful expressions of Finger's strain tensor in a cylindrical (or Cartesian) coordinate system as well as in an orthogonal streamline coordinate system. One of the results shows that Winter's tracking model is correct.Relations among the recent three results of Winter, Adachi and Crochet et al. are presented clearly. Moreover useful applications of Frenet-Serret's formula to the study of the deformation and flow kinematics along streamlines are shown in comparison with the ordinary tensor approach.     

An adaptive wavelet space‐time SUPG method for hyperbolic conservation laws

This article concerns with incorporating wavelet bases into existing streamline upwind Petrov‐Galerkin (SUPG) methods for the numerical solution of nonlinear hyperbolic conservation laws which are known to develop shock solutions. Here, we utilize an SUPG formulation using continuous Galerkin in space and discontinuous Galerkin in time. The main motivation for such a combination is that these methods have good stability properties thanks to adding diffusion in the direction of streamlines. But they are more expensive than explicit semidiscrete methods as they have to use space‐time formulations. Using wavelet bases we maintain the stability properties of SUPG methods while we reduce the cost of these methods significantly through natural adaptivity of wavelet expansions. In addition, wavelet bases have a hierarchical structure. We use this property to numerically investigate the hierarchical addition of an artificial diffusion for further stabilization in spirit of spectral diffusion. Furthermore, we add the hierarchical diffusion only in the vicinity of discontinuities using the feature of wavelet bases in detection of location of discontinuities. Also, we again use the last feature of the wavelet bases to perform a postprocessing using a denosing technique based on a minimization formulation to reduce Gibbs oscillations near discontinuities while keeping other regions intact. Finally, we show the performance of the proposed combination through some numerical examples including Burgers’, transport, and wave equations as well as systems of shallow water equations.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2062–2089, 2017