A collocated finite volume scheme to solve free convection for general non-conforming grids

We present a new collocated numerical scheme for the approximation of the Navier–Stokes and energy equations under the Boussinesq assumption for general grids, using the velocity–pressure unknowns. This scheme is based on a recent scheme for the diffusion terms. Stability properties are drawn from particular choices for the pressure gradient and the non-linear terms. Convergence of the approximate solutions may be proven mathematically. Numerical results show the accuracy of the scheme on irregular grids.     

Second-order sign-preserving conservative interpolation (remapping) on general grids

An accurate conservative interpolation (remapping) algorithm is an essential component of most arbitrary Lagrangian–Eulerian (ALE) methods. In this paper we describe a local remapping algorithm for a positive scalar function. This algorithm is second-order accurate, conservative, and sign preserving. The algorithm is based on estimating the mass exchanged between cells at their common interface, and so is equally applicable to structured and unstructured grids. We construct the algorithm in a series of steps, clearly delineating the assumptions and errors made at each step. We validate our theory with a suite of numerical examples, analyzing the results from the viewpoint of accuracy and order of convergence.     

A general scheme for microscopic theories

We discuss from an operational point of view some fundamental concepts of the micro-scopic physical theories, with the aim of providing a background for a successive investigation of the microscopic space-time structure. We consistently develop the remark that a frame of reference is determined by physical objects which may interact with the objects under investigation. As it is not clear that a state can be prepared by means of physical operations, we do not use the concept of physical state for the foundation of the theory. In our approach, the primitive concepts are the measurement procedures, following which one gets a numerical result, and the transformation procedures, which have the aim of building a frame of reference. We discuss several rules which allow us to define new procedures in terms of known procedures. The statistical laws of physics are formulated in terms of an order relation between measurement procedures, which defines also an equivalence relation. The equivalence classes of measurement procedures are called measurements. We define also equivalence classes of transformation procedures, called transformations. The mathematical structure of the set of measurements and of the set of transformations is discussed in detail. We consider measurements with an arbitrary finite number of possible results, as this enables us to give a rigorous definition of compatibility. Finally, we point out that all the physical theories necessarily contain ideal measurements and transformations which do not correspond to any known physical procedure. The introduction of these ideal objects permits a considerable simplification of the mathematical structure of the theory, but reduces its physical content.     

A compatible and conservative spectral element method on unstructured grids

We derive a formulation of the spectral element method which is compatible on very general unstructured three-dimensional grids. Here compatible means that the method retains discrete analogs of several key properties of the divergence, gradient and curl operators: the divergence and gradient are anti-adjoints (the negative transpose) of each other, the curl is self-adjoint and annihilates the gradient operator, and the divergence annihilates the curl. The adjoint relations hold globally, and at the element level with the inclusion of a natural discrete element boundary flux term.     

A dynamic finite volume scheme for large-eddy simulation on unstructured grids

In recent years there has been considerable progress in the application of large-eddy simulation (LES) to increasingly complex flow configurations. Nevertheless a lot of fundamental problems still need to be solved in order to apply LES in a reliable way to real engineering problems, where typically finite-volume codes on unstructured meshes are used. A self-adaptive discretisation scheme, in the context of an unstructured finite-volume flow solver, is investigated in the case of isotropic turbulence at infinite Reynolds number. The Smagorinsky and dynamic Smagorinsky sub-grid scale models are considered. A discrete interpolation filter is used for the dynamic model. It is one of the first applications of a filter based on the approach presented by Marsden et al. In this work, an original procedure to impose the filter shape through a specific selection process of the basic filters is also proposed. Satisfactory results are obtained using the self-adaptive scheme for implicit LES. When the scheme is coupled with the sub-grid scale models, the numerical dissipation is shown to be dominant over the sub-grid scale component. Nevertheless the effect of the sub-grid scale models appears to be important and beneficial, improving in particular the energy spectra. A test on fully developed channel flow at Re_τ = 395 is also performed, comparing the non-limited scheme with the self-adaptive scheme for implicit LES. Once again the introduction of the limiter proves to be beneficial.     

A conservative numerical scheme for solving an autonomous Hamiltonian system

A new numerical scheme is proposed for solving Hamilton’s equations that possesses the properties of symplecticity. Just as in all symplectic schemes known to date, in this scheme the conservation laws of momentum and angular momentum are satisfied exactly. A property that distinguishes this scheme from known schemes is proved: in the new scheme, the energy conservation law is satisfied for a system of linear oscillators. The new numerical scheme is implicit and has the third order of accuracy with respect to the integration step. An algorithm is presented by which the accuracy of the scheme can be increased up to the fifth and higher orders. Exact and numerical solutions to the two-body problem, calculated by known schemes and by the scheme proposed here, are compared.     

An efficient cell-centered diffusion scheme for quadrilateral grids

A new algorithm for solution of diffusion equations in two dimensions on structured quadrilateral grids is proposed. The algorithm is based on a semi-implicit method for the time discretization and has a nine-point local stencil in space. Our scheme is fast, quite accurate and demonstrates good spatial convergence. The presented numerical tests show that it is well suited for hydrocodes with cell-centered principal variables.     

Enriched multi-point flux approximation for general grids

It is well known that the two-point flux approximation, a numerical scheme used in most commercial reservoir simulators, has O(1) error when grids are not K-orthogonal. In the last decade, the multi-point flux approximations have been developed as a remedy. However, non-physical oscillations can appear when the anisotropy is really strong. We found out the oscillations are closely related to the poor approximation of pressure gradient in the flux computation.In this paper, we propose the control volume enriched multi-point flux approximation (EMPFA) for general diffusion problems on polygonal and polyhedral meshes. Non-physical oscillations are not observed for realistic and strongly anisotropic heterogeneous material properties described by a full tensor. Exact linear solutions are recovered for grids with non-planar interfaces, and a first and second order convergence are achieved for the flux and scalar unknowns, respectively.     

A conservative level-set based method for compressible solid/fluid problems on fixed grids

A three-dimensional Eulerian method is presented for simulating dynamic systems comprising multiple compressible solid and fluid components where internal boundaries are tracked using level-set functions. Aside from the interface interaction calculation within mixed cells, each material is treated independently and the governing constitutive laws solved using a conservative finite volume discretisation based upon the solution of Riemann problems to determine the numerical fluxes. The required reconstruction of mixed cell volume fractions and cut cell geometries is presented in detail using the level-set fields. High-order accuracy is achieved by incorporating the weighted-essentially non-oscillatory (WENO) method and Runge–Kutta time integration. A model for elastoplastic solid dynamics is employed formulated using the tensor of elastic deformation gradients permitting the equations to be written in divergence form. The scheme is demonstrated using selected one-dimensional initial value problems for which exact solutions are derived, a two-dimensional void collapse, and a three-dimensional simulation of a confined explosion.     

A general scheme for computation of instanton effects

A general approach to the computation of instanton effects without invoking an ad hoc cut-off on their size is illustrated here by the computation of instanton contributions to the static $\text{Q}$Q potential and vacuum polarization. Our results suggest that the instanton effects are much smaller than perturbative effects in the region where both may be reliably computed.     

Method for efficient prevention of gravity wave decoupling on rectangular semi-staggered grids

Generation of short gravity wave noise often occurs on semi-staggered rectangular grids as a result of sub-grid decoupling when there is a strong forcing in the mass field. In this study a numerical scheme has been proposed to prevent the generation of the gravity wave decoupling. The proposed numerical method provides efficient communication between decoupled gravity wave finite-difference solutions by a simple averaging of a term in the height tendency in the continuity equation. The proposed method is tested using a shallow water sink model developed for the purpose of this study. It has been demonstrated that this method outperforms other existing approaches. The new scheme is time efficient, based on explicit time integration and can be easily implemented. The proposed method is applicable in hydrodynamic models specified on semi-staggered B or E grids.     

A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part I: On a rectangular collocated grid system

A consistent, conservative and accurate scheme has been designed to calculate the current density and the Lorentz force by solving the electrical potential equation for magnetohydrodynamics (MHD) at low magnetic Reynolds numbers and high Hartmann numbers on a finite-volume structured collocated grid. In this collocated grid, velocity (u), pressure (p), and electrical potential (φ) are located in the grid center, while current fluxes are located on the cell faces. The calculation of current fluxes on the cell faces is conducted using a conservative scheme, which is consistent with the discretization scheme for the solution of electrical potential Poisson equation. A conservative interpolation is used to get the current density at the cell center, which is used to conduct the calculation of Lorentz force at the cell center for momentum equations. We will show that both “conservative” and “consistent” are important properties of the scheme to get an accurate result for high Hartmann number MHD flows with a strongly non-uniform mesh employed to resolve the Hartmann layers and side layers of Hunt’s conductive walls and Shercliff’s insulated walls. A general second-order projection method has been developed for the incompressible Navier–Stokes equations with the Lorentz force included. This projection method can accurately balance the pressure term and the Lorentz force for a fully developed core flow. This method can also simplify the pressure boundary conditions for MHD flows.     

Quasi-optical grids with thin rectangular patch/aperture elements

Theoretical analysis is presented for an efficient and accurate performance evaluation of quasi-optical grids comprised of thin rectangular patch/aperture elements with/without a dielectric substrate/superstrate. The convergence rate of this efficient technique is improved by an order of magnitude with the approximate edge conditions incorporated in the basis functions of the integral equation solution. Also presented are the interesting applications of this efficient analytical technique to the design and performance evaluation of the coupling grids and beam splitters in the optical systems as well as thermal protection sunshields used in the communication systems of satellites and spacecrafts.     

A third-order finite-volume residual-based scheme for the 2D Euler equations on unstructured grids

A residual-based (RB) scheme relies on the vanishing of residual at the steady-state to design a transient first-order dissipation, which becomes high-order at steady-state. Initially designed within a finite-difference framework for computations of compressible flows on structured grids, the RB schemes displayed good convergence, accuracy and shock-capturing properties which motivated their extension to unstructured grids using a finite volume (FV) method. A second-order formulation of the FV–RB scheme for compressible flows on general unstructured grids was presented in a previous paper. The present paper describes the derivation of a third-order FV–RB scheme and its application to hyperbolic model problems as well as subsonic, transonic and supersonic internal and external inviscid flows.     

A general multi-secret visual cryptography scheme

In a (k, n) visual cryptography scheme (VCS), a secret image is encoded into n shadow images that are distributed to n participants. Any k participants can reveal the secret image by stacking their shadow images, and less than k participants have no information about the secret image. In this paper we consider the case when the secret image is more than one, and this is a so-called multi-secret VCS (MVCS). The previous works on MVCS are all the simple 2-out-of-2 cases. In this paper, we discuss a general (k, n)-MVCS for any k and n. This paper has three main contributions: (1) our scheme is the first general (k, n)-MVCS, which can be applied on any k and n, (2) we give the formal security and contrast conditions of (k, n)-MVCS and (3) we theoretically prove that the proposed (k, n)-MVCS satisfies the security and contrast conditions.     

A general central limit theorem for FKG systems

A central limit theorem is given which is applicable to (not necessarily monotonic) functions of random variables satisfying the FKG inequalities. One consequence is convergence of the block spin scaling limit for the magnetization and energy densities (jointly) to the infinite temperature fixed point of independent Gaussian blocks for a large class of Ising ferromagnets whenever the susceptibility is finite. Another consequence is a central limit theorem for the density of thesurface of the infinite cluster in percolation models.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A nominally second-order accurate finite volume cell-centered scheme for anisotropic diffusion on two-dimensional unstructured grids

In this paper, we describe a second-order accurate cell-centered finite volume method for solving anisotropic diffusion on two-dimensional unstructured grids. The resulting numerical scheme, named CCLAD (Cell-Centered LAgrangian Diffusion), is characterized by a local stencil and cell-centered unknowns. It is devoted to the resolution of diffusion equation on distorted grids in the context of Lagrangian hydrodynamics wherein a strong coupling occurs between gas dynamics and diffusion. The space discretization relies on the introduction of two half-edge normal fluxes and two half-edge temperatures per cell interface using the partition of each cell into sub-cells. For each cell, the two half-edge normal fluxes attached to a node are expressed in terms of the half-edge temperatures impinging at this node and the cell-centered temperature. This local flux approximation can be derived through the use of either a sub-cell variational formulation or a finite difference approximation, leading to the two variants CCLADS and CCLADNS. The elimination of the half-edge temperatures is performed locally at each node by solving a small linear system which is obtained by enforcing the continuity condition of the normal heat flux across sub-cell interface impinging at the node. The accuracy and the robustness of the present scheme is assessed by means of various numerical test cases.