Fully HOC Scheme for Mixed Convection Flow in a Lid-Driven Cavity Filled with a Nanofluid

A fully higher-order compact (HOC) finite difference scheme on the 9-point two-dimensional (2D) stencil is formulated for solving the steady-state laminar mixed convection flow in a lid-driven inclined square enclosure filled with water-$Al_2O_3$ nanofluid. Two cases are considered depending on the direction of temperature gradient imposed (Case I, top and bottom; Case II, left and right). The developed equations are given in terms of the stream function-vorticity formulation and are non-dimensionalized and then solved numerically by a fourth-order accurate compact finite difference method. Unlike other compact solution procedure in literature for this physical configuration, the present method is fully compact and fully higher-order accurate. The fluid flow, heat transfer and heat transport characteristics were illustrated by streamlines, isotherms and averaged Nusselt number. Comparisons with previously published work are performed and found to be in excellent agreement. A parametric study is conducted and a set of graphical results is presented and discussed to elucidate that significant heat transfer enhancement can be obtained due to the presence of nanoparticles and that this is accentuated by inclination of the enclosure at moderate and large Richardson numbers.     

A Fast High Order Iterative Solver for the Electromagnetic Scattering by Open Cavities Filled with the Inhomogeneous Media

The scattering of the open cavity filled with the inhomogeneous media is studied. The problem is discretized with a fourth order finite difference scheme and the immersed interface method, resulting in a linear system of equations with the high order accurate solutions in the whole computational domain. To solve the system of equations, we design an efficient iterative solver, which is based on the fast Fourier transformation, and provides an ideal preconditioner for Krylov subspace method. Numerical experiments demonstrate the capability of the proposed fast high order iterative solver.     

Solitary wave solution to Aw-Rascle viscous model of traffic flow

A traveling wave solution to the Aw-Rascle traffic flow model that includes the relaxation and diffusion terms is investigated. The model can be approximated by the well-known Kortweg-de Vries (KdV) equation. A numerical simulation is conducted by the first-order accurate Lax-Friedrichs scheme, which is known for its ability to capture the entropy solution to hyperbolic conservation laws. Periodic boundary conditions are applied to simulate a lengthy propagation, where the profile of the derived KdV solution is taken as the initial condition to observe the change of the profile. The simulation shows good agreement between the approximated KdV solution and the numerical solution.     

Numerical simulation of free surface flows with the level set method using an extremely high-order accuracy WENO advection scheme

Three-dimensional dynamic gas–liquid flow simulations that accurately track the phase interface are numerically challenging. This article presents a numerical study of the performance of the level-set phase–interface tracking method when combined with extremely high order (7th to 11th) weighted essentially non-oscillatory (WENO) advection schemes for gas–liquid free surface flows. Comparisons between simulation results and prior benchmark results suggest that such a combination of methods can be satisfactorily applied to the level-set and Navier-Stokes equations for free surface flow simulations when volume conservation is enforced at every time step, and minor numerical oscillations are suppressed through use of an artificial viscosity term. In particular, simulations of solid body rotation, the unsteady flow following an ideal dam break, tank sloshing, and the rise of a single bubble all agree with analytical or experimental results to within ± 3.12% when the level-set method is combined with an 11th order WENO scheme. Furthermore, use of an 11th order WENO advection scheme actually has a computational cost advantage because, for the same accuracy, it can be used on a coarser grid when compared with a more-common second-order advection scheme; computational savings of up to 87% are possible.     

应用于计算气动声学的优化有限紧致格式

对于含间断的计算气动声学问题,数值计算的格式不仅要求低耗散低色散的设计,对短波具有较高的分辨率,还要求能捕捉激波.中心紧致格式具有高精度,具有无耗散和低色散特征,但不能捕捉间断和激波;WENO格式处理间断较为成功,而耗散和色散误差相对较大.有限紧致格式可以将紧致格式与WENO格式相结合构造成混合格式,利用光滑因子之间的关系对激波区域进行自动判断,将传统的全域求解的紧致格式划分为有限的局部紧致求解,间断点上的激波捕捉铜梁自动作为局部紧致求解的边界通量,在在光滑区域具有紧致格式的高精度低耗散性质,在激波附近不产生非物理振荡.本文利用有限紧致格式思想,构造了新的适合于气动声学问题的优化有限紧致格式,将其应用于计算气动声学一维标准测试问题,对相关格式的模拟性能进行了评估,显示该格式在宽频声波传播和含有间断的声波传播模拟方面具有优势.     

An Adaptive Grid Eulerian Method for the Computation of Free Surface Flows

A numerical scheme for the prediction of free surface flows is presented and investigated. The method is based on an adaptive grid Eulerian finite-volume method, where non-orthogonal boundary-fitted moving grids are employed to follow the free surface. The underlying flow solver consists in a pressure-correction scheme of SIMPLE type with multigrid acceleration, which is iteratively combined with the moving grid technique. Several numerical examples are considered to illustrate the capabilities of the approach.     

Two-stage fourth-order gas-kinetic scheme for three-dimensional Euler and Navier-Stokes solutions

ABSTRACT

For the one-stage third-order gas-kinetic scheme (GKS), successful applications have been achieved for the three-dimensional compressible flows [Pan, L., K. Xu, Q. Li, and J. Li. 2016. “An Efficient and Accurate Two-stage Fourth-order Gas-kinetic Scheme for the Navier-Stokes Equations.” Journal of Computational Physics 326: 197–221]. The high-order accuracy of the scheme is obtained by integrating a multidimensional time-accurate gas distribution function over the cell interface within a time step without using Gaussian quadrature points and Runge-Kutta time-stepping technique. However, to the further increase of the order of the scheme, such as the fourth-order one, the one step formulation becomes very complicated for the multidimensional flow. Recently, a two-stage fourth-order GKS with high efficiency has been constructed for two-dimensional inviscid and viscous flow computations ([Li, J., and Z. Du. 2016. “A Two-stage Fourth Order Time-accurate Discretization for Lax-Wendroff Type Flow Solvers I. Hyperbolic Conservation Laws.” SIAM Journal on Scientific Computing 38: 3046–3069]; Pan et al. 2016), and the scheme uses the time accurate flux function and its time derivatives. In this paper, a fourth-order GKS is developed for the three-dimensional flows under the two-stage framework. Based on the three-dimensional WENO reconstruction and flux evaluation at Gaussian quadrature points on a cell interface, the high-order accuracy in space is achieved first. Then, the two-stage time stepping method provides the high accuracy in time. In comparison with the formal third-order GKS [Pan, L., and K. Xu. 2015. “A Third-order Gas-kinetic Scheme for Three-dimensional Inviscid and Viscous Flow Computations.” Computers & Fluids 119: 250–260], the current fourth-order method not only improves the accuracy of the scheme, but also reduces the complexity of the gas-kinetic flux solver greatly. More importantly, the fourth-order GKS has the same robustness as the second-order shock capturing scheme [Xu, K. 2001. “A Gas-kinetic BGK Scheme for the Navier-Stokes Equations and its Connection with Artificial Dissipation and Godunov Method.” Journal of Computational Physics 171: 289–335]. Numerical results validate the outstanding reliability and applicability of the scheme for three-dimensional flows, such as the cases related to turbulent simulations.     

On the Accuracy of Upwind and Symmetric TVD Schemes in Simulating Low Mach Number Flow

The accuracy of MUSCL upwind and Yee-Roe-Davis symmetric TVD schemes for simulating low Mach number flow is studied through a numerical experiment of the 2-D lid driven cavity problem. The steady slate solution is reached by using a marching approach based on the pseudocompressibilty method in conjunction with implicit approximate factorization. A finite volume discretization of the conservation equations is used with a four level multigrid method to accelerate the convergence. The tests performed which were in the range of 100 ≤ Re ≤ 5000, show that the Yee-Roe-Davis symmetric scheme generates results in very good agreement with the benchmark results over this range of Re. The MUSCL upwind scheme accuracy deteriorates with the increasing Re.     

基于自编码网络流形学习的毛竹笋不溶性膳食纤维含量红外光谱建模

提出了一种结合自编码网络(AN)流形学习和偏最小二乘(PLS)法的红外光谱建模方法AN-PLS。AN-PLS方法首先用AN算法对红外光谱数据进行非线性降维,再结合PLS建立回归模型。利用该方法建立了毛竹笋中不溶性膳食纤维含量的近红外光谱和中红外光谱回归模型。结果表明,用AN-PLS方法建立的回归模型,比用其他常用光谱数据预处理方法结合PLS及用单独PLS算法建立的模型具有更小的预测均方根误差RMSEP和更高的决定系数R2,因此,AN-PLS具有较优的建模与预测能力,利用近红外光谱和中红外光谱技术结合AN-PLS建模,可实现毛竹笋中不溶性膳食纤维含量的准确测量。     

一维守恒的Lagrangian ADER格式

基于欧拉框架下ADER格式,构造一维守恒只有一个时间步的、高精度中心型拉格朗日ADER(LADER)格式.构造r阶LADER格式包括:从欧拉方程出发推导拉格朗日框架下积分形式的方程、采用WENO方法高精度重构节点处守恒量和从1阶到r-1阶的空间导数、求拉氏框架下这些变量的Godunov值,并计算1阶到r-1阶的时间全导数,最后高精度离散积分形式的流通量函数.对光滑流场的模拟表明,LADER格式达到设计的精度;对含强间断的流场模拟表明,数值解在间断附近基本无振荡.