High-frequency response of a calculation methodology for gas dynamics based on Independent Time Discretisation

A calculation methodology to solve the one-dimensional governing equations system is presented. This calculation methodology is based on the Independent Time Discretisation (ITD) of the ducts composing the system. The purpose is the improvement of the trade-off between the accuracy and the computational cost that the current 1D gas dynamic models can yield. The ITD methodology is applied to the specific problem of noise prediction in internal combustion engines in order to evaluate its performance in the frequency domain. The application of the ITD methodology to the well-known acoustic configurations which are representative of the main attenuation mechanisms in commercial mufflers shows its ability. The potential is evaluated in terms of reduction of the computational cost and the accuracy and robustness provided by the results as a function of the spatial mesh size and the family of finite difference numerical method applied.     

An adaptive mesh method for 1D hyperbolic conservation laws

An adaptive method is developed for solving one-dimensional systems of hyperbolic conservation laws, which combines the rezoning approach with the finite volume weighted essentially non-oscillatory (WENO) scheme. An a posteriori error estimate, used to equidistribute the mesh, is obtained from the differences between respective numerical solutions of 5th-order WENO (WENO5) and 3rd-order ENO (ENO3) schemes. The number of grids can be adaptively readjusted based on the solution structure. For higher efficiency, mesh readjustment is performed every few time steps rather than every time step. In addition, a high order conservative interpolation is used to compute the physical solutions on the new mesh from old mesh based on the finite volume ENO reconstruction. Extensive examples suggest that this adaptive method exhibits more accurate resolution of discontinuities for a similar level of computational time comparing with that on a uniform mesh.     

Construction of locally conservative fluxes for the SUPG method

We consider the construction of locally conservative fluxes by means of a simple postprocessing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields nonconservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element method for solving nondominating advection diffusion equations and the streamline upwind/Petrov‐Galerkin method for solving advection dominated problems. We then describe the postprocessing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The postprocessing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the postprocessing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift‐diffusion equations. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1971–1994, 2015     

ALTERING CONNECTIVITY WITH LARGE DEFORMATION MESH FOR LAGRANGIAN METHOD AND ITS APPLICATION IN MULTIPLE MATERIAL SIMULATION

A new approach for treating the mesh with Lagrangian scheme of finite volume method is presented. It has been proved that classical Lagrangian method is difficult to cope with large deformation in tracking material particles due to severe distortion of cells, and the changing connectivity of the mesh seems especially attractive for solving such issues. The mesh with large deformation based on computational geometry is optimized by using new method. This paper develops a processing system for arbitrary polygonal unstructured grid, the intelligent variable grid neighborhood technologies is utilized to improve the quality of mesh in calculation process, and arbitrary polygonal mesh is used in the Lagrangian finite volume scheme. The performance of the new method is demonstrated through series of numerical examples, and the simulation capability is efficiently presented in coping with the systems with large deformations.     

不可压缩流动的并行数值方法

不可压缩流动的数值模拟是计算流体力学的重要组成部分. 基于有限元离散方法, 本文设计了不可压缩Navier-Stokes (N-S)方程支配流的若干并行数值算法. 这些并行算法可归为两大类: 一类是基于两重网格离散方法, 首先在粗网格上求解非线性的N-S方程, 然后在细网格的子区域上并行求解线性化的残差方程, 以校正粗网格的解; 另一类是基于新型完全重叠型区域分解技巧, 每台处理器用一局部加密的全局多尺度网格计算所负责子区域的局部有限元解. 这些并行算法实现简单, 通信需求少, 具有良好的并行性能, 能获得与标准有限元方法相同收敛阶的有限元解. 理论分析和数值试验验证了并行算法的高效性     

基于PDE和几何曲率流驱动扩散的图像分析与处理

本文介绍由变分优化模型导出的偏微分方程(PDEs)模型与几何曲率流驱动扩散在图像恢复方面的应用，以及多种非线性异质扩散模型，讨论了PDEs模型在图像分析与处理方面的优点，理论与实验结果表明，要恢复得到商质量的图像，PDEs模型的利用是极为必要的．文中还介绍了求解PDEs模型的数值方案．其中，曲率计算是一个关键问题，其结果直接参与自适应扩散的控制．详细总结了基于有限差分和水平集方法，求解藕合非线性异质扩散模型方程的数值方案，追求高质量图像、高精度计算方法、降低计算复杂性是本文处理方法不断进步的发展动力。     

A Lyapunov function for a two-chemical species version of the chemotaxis model

Since the accuracy of finite element solutions of partial differential equations is generally mesh dependent, especially when solutions have singularities and discontinuities, a proper mesh generation is often important and sometimes crucial for an accurate numerical approximation of such problems. In this paper, the mesh transformation method is applied to the boundary value problems of elliptic partial differential equations, and it is proved that the method leads to the optimal finite element solutions. AMS subject classification (2000) 73C50, 65K10, 65N12, 65N30     

The numerical simulation of a typical fluid-structure interaction problem

To determine the dynamic response of a structure under the influence of the fluid flow one must solve a coupled computational fluid dynamics (CFD) and computational structural dynamics (CSD) mathematical problem. This paper presents the comparison of two methods for the calculation of the fluid-structure interaction. The first one is of explicit-implicit type and uses a staggered time advancement of the fluid and structure problems. The second uses a fully implicit discretization in the physical time of the fluid-structure equations and an explicit advancement in the dual-time. The physical fluid-structure problem is accompanied by the equations of the mesh motion, which are written as for a pseudo-structural system with its own dynamics. Representative numerical results are presented for the two degrees of freedom tipical section in unsteady transonic flow. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Uniformly convergent novel finite difference methods for singularly perturbed reaction–diffusion equations

In this paper, we construct a kind of novel finite difference (NFD) method for solving singularly perturbed reaction–diffusion problems. Different from directly truncating the high‐order derivative terms of the Taylor's series in the traditional finite difference method, we rearrange the Taylor's expansion in a more elaborate way based on the original equation to develop the NFD scheme for 1D problems. It is proved that this approach not only can highly improve the calculation accuracy but also is uniformly convergent. Then, applying alternating direction implicit technique, the newly deduced schemes are extended to 2D equations, and the uniform error estimation based on Shishkin mesh is derived, too. Finally, numerical experiments are presented to verify the high computational accuracy and theoretical prediction.     

一种健壮的低耗散通量分裂格式

随着计算流体力学的快速发展,设计精确、高效并且健壮的数值格式变得尤为重要.通过对3种流行的通量分裂方法(AUSM、Zha-Bilgen和Toro-Vázquez)的对流通量和压力通量进行特征分析,构造了一种简单、低耗散并且健壮的通量分裂格式(命名为R-ZB格式).采用Zha-Bilgen分裂方法将Euler方程的通量分裂成对流通量和压力通量,其中对流通量采用迎风方法来计算,压力通量采用低耗散的HLL格式来计算,从而克服了原始的HLL格式不能精确分辨接触间断的缺点.数值实验表明,该文给出的R-ZB格式不仅保留了原始Zha-Bilgen格式简单高效、能够精确分辨接触间断等优点,而且具有更好的健壮性,在计算二维问题时不会出现数值激波不稳定现象.     

An $h$-Adaptive Runge-Kutta Discontinuous Galerkin Method for Hamilton-Jacobi Equations

In [35, 36], we presented an $h$-adaptive Runge-Kutta discontinuous Galerkin method using troubled-cell indicators for solving hyperbolic conservation laws. A tree data structure (binary tree in one dimension and quadtree in two dimensions) is used to aid storage and neighbor finding. Mesh adaptation is achieved by refining the troubled cells and coarsening the untroubled "children". Extensive numerical tests indicate that the proposed $h$-adaptive method is capable of saving the computational cost and enhancing the resolution near the discontinuities. In this paper, we apply this $h$-adaptive method to solve Hamilton-Jacobi equations, with an objective of enhancing the resolution near the discontinuities of the solution derivatives. One- and two-dimensional numerical examples are shown to illustrate the capability of the method.     

Resolución numérica de las ecuaciones del movimiento de edificios de varias plantas con no linealidades severas

This paper presents a new algorithm for solving the equations of motion of multi-storey buildings that incorporate frictional energy dissipators as seismic protection. The behavior of the dissipators is represented by Coulomb dry friction models; they introduce severe nonlinearities in the dynamic behavior of the structure every time that the contact conditions (stick or slip) change in the dissipators. These nonlinearities complicate the resolution of the equations of motion as it usually is described by lumped masses models whose degrees of freedom are the displacements of the floors and, as the stick or slip conditions change, the degrees of freedom must be modified: for blocking conditions they are only the displacements of the storeys while under sliding conditions the displacements of the dissipators have to be also considered. In previous articles the accuracy of the proposed algorithm has been verified by comparison with experimental results; as well, the computational efficiency of the algorithm has been confirmed by comparing the required resources (in terms of computation time and of memory allocation) with those of other algorithms. The objectives of this paper are to describe in detail the numerical solution of the equations of motion and present representative examples confirming the ability of the algorithm to reproduce the dynamic behavior of buildings with friction dissipators and reporting preliminarily about the usefulness of such devices to reduce the oscillations of the structure to be protected.     

Three-and four-step implicit absolutely stable fourth-order Runge-Kutta schemes

Two types of implicit fourth-order Runge-Kutta schemes are constructed for first-order ordinary differential equations, multidimensional transfer equations, and compressible gas equations. The absolute stability of the schemes is proved by applying the principle of frozen coefficients. Adaptive artificial viscosity ensuring good time convergence and oscillations damping near discontinuities is used in solving gas dynamics equations. The comparative efficiency of the schemes is illustrated by numerical results obtained for compressible gas flows.     

A two‐level finite element method for the stationary Navier‐Stokes equations based on a stabilized local projection

In this article, we propose a two‐level finite element method to analyze the approximate solutions of the stationary Navier‐Stokes equations based on a stabilized local projection. The local projection allows to circumvent the Babuska‐Brezzi condition by using equal‐order finite element pairs. The local projection can be used to stabilize high equal‐order finite element pairs. The proposed method combines the local projection stabilization method and the two‐level method under the assumption of the uniqueness condition. The two‐level method consists of solving a nonlinear equation on the coarse mesh and solving a linear equation on fine mesh. The nonlinear equation is solved by the one‐step Newtonian iteration method. In the rest of this article, we show the error analysis of the lowest equal‐order finite element pair and provide convergence rate of approximate solutions. Furthermore, the numerical illustrations coincide with the theoretical analysis expectations. From the view of computational time, the results show that the two‐level method is effective to solve the stationary Navier‐Stokes equations. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Numeric solution of advection–diffusion equations by a discrete time random walk scheme

Explicit numerical finite difference schemes for partial differential equations are well known to be easy to implement but they are particularly problematic for solving equations whose solutions admit shocks, blowups, and discontinuities. Here we present an explicit numerical scheme for solving nonlinear advection–diffusion equations admitting shock solutions that is both easy to implement and stable. The numerical scheme is obtained by considering the continuum limit of a discrete time and space stochastic process for nonlinear advection–diffusion. The stochastic process is well posed and this guarantees the stability of the scheme. Several examples are provided to highlight the importance of the formulation of the stochastic process in obtaining a stable and accurate numerical scheme.     

A finite element method for contact/impact

Ideas from the analysis of differential-algebraic equations are applied to the numerical solution of frictionless contact/impact problems in solid mechanics. Index-one and two formulations for dynamic contact–impact within the context of the finite element method are derived. The resulting equations are shown to stabilize the kinematic fields at the contact interface, at the expense of a small energy loss, which is shown to decrease consistently with mesh refinement. This energy dissipation is shown to be necessary for the establishment of persistent contact. A Newmark-type time integration scheme is derived from the proposed formulation, and shown to yield excellent results in modeling the transition to contact/impact.     

Modeling and simulation of a gas distribution pipeline network

This research study focuses on the modeling and simulation of a gas distribution pipeline network with a special emphasis on gas ducts. Gas ducts are the most important components of such kind of systems since they define the major dynamic characteristics. Isothermal, unidirectional flow is usually assumed when modeling the gas flow through a gas duct. This paper presents two simplified models derived from the set of partial differential equations governing the dynamics of the process. These models include the inclination term, neglected in most related papers. Moreover, two numerical schemes are presented for the integration of such models. Also, it is shown how the pressure drop along the pipe has a strong dependency with the inclination term. To solve the system dynamics through the proposed numerical schemes a based MATLAB-Simulink library was built. With this library it is possible to simulate the behavior of a gas distribution network from the individual simulation of each component. Finally, the library is tested through three application examples, and results are compared with the existing ones in the literature.     

A TWO-GRID METHOD FOR THE C0 INTERIOR PENALTY DISCRETIZATION OF THE MONGE-AMP(E)RE EQUATION

The purpose of this paper is to analyze an efficient method for the solution of the nonlinear system resulting from the discretization of the elliptic Monge-Ampère equation by a $C^0$ interior penalty method with Lagrange finite elements. We consider the two-grid method for nonlinear equations which consists in solving the discrete nonlinear system on a coarse mesh and using that solution as initial guess for one iteration of Newton's method on a finer mesh. Thus both steps are inexpensive. We give quasi-optimal $W^{1,\infty}$ error estimates for the discretization and estimate the difference between the interior penalty solution and the two-grid numerical solution. Numerical experiments confirm the computational efficiency of the approach compared to Newton's method on the fine mesh.     

A hybrid linking approach for solving the conservation equations with an adaptive mesh refinement method

Accurate numerical computation of complex flows on a single grid requires very fine meshes to capture phenomena occurring at both large and small scales. The use of adaptive mesh refinement (AMR) methods, significantly reduces the involved computational time and memory. In the present article, a hybrid linking approach for solving the conservation equations with an AMR method is proposed. This method is essentially a coupling between the h-AMR and the multigrid methods. The efficiency of the present approach has been demonstrated by solving species and Navier–Stokes equations.     

Numerical solution of the fully non-linear weakly dispersive serre equations for steep gradient flows

We demonstrate a numerical approach for solving the one-dimensional non-linear weakly dispersive Serre equations. By introducing a new conserved quantity the Serre equations can be written in conservation law form, where the velocity is recovered from the conserved quantities at each time step by solving an auxiliary elliptic equation. Numerical techniques for solving equations in conservative law form can then be applied to solve the Serre equations. We demonstrate how this is achieved. The system of conservation equations are solved using the finite volume method and the associated elliptic equation for the velocity is solved using a finite difference method. This robust approach allows us to accurately solve problems with steep gradients in the flow, such as those generated by discontinuities in the initial conditions.The method is shown to be accurate, simple to implement and stable for a range of problems including flows with steep gradients and variable bathymetry.