A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model. © 2011 Elsevier Ltd.

Immersed Boundary Method for generalised finite volume and finite difference Navier-Stokes solvers

In Immersed Boundary Methods (IBM) the effect of complex geometries is introduced through the forces added in the Navier-Stokes solver at the grid points in the vicinity of the immersed boundaries. Most of the methods in the literature have been used with Cartesian grids. Moreover many of the methods developed in the literature do not satisfy some basic conservation properties (the conservation of torque, for instance) on non-uniform meshes. In this paper we will follow the RKPM method originated by Liu et al. [1] to build locally regularized functions that verify a number of integral conditions. These local approximants will be used both for interpolating the velocity field and for spreading the singular force field in the framework of a pressure correction scheme for the incompressible Navier-Stokes equations. We will also demonstrate the robustness and effectiveness of the scheme through various examples. Copyright © 2010 by ASME.

Finite element solution of Navier Stokes equations adapted to a priori error estimates

The paper is based on qualitative properties of the solution of the Navier-Stokes equations for incompressible fluid, and on properties of their finite element solution. In problems with corner-like singularities (e.g. on the well-known L-shaped domain) usually some adaptive strategy is used. In this paper we present an alternative approach. For flow problems on domains with corner singularities we use the a priori error estimates and asymptotic expansion of the solution to derive an algorithm for refining the mesh near the corner. It gives very precise solution in a cheap way. We present some numerical results.

Adjoint algorithms for the Navier-Stokes equations in the low Mach number limit

This paper describes a derivation of the adjoint low Mach number equations and their implementation and validation within a global mode solver. The advantage of using the low Mach number equations and their adjoints is that they are appropriate for flows with variable density, such as flames, but do not require resolution of acoustic waves. Two versions of the adjoint are implemented and assessed: a discrete-adjoint and a continuous-adjoint. The most unstable global mode calculated with the discrete-adjoint has exactly the same eigenvalue as the corresponding direct global mode but contains numerical artifacts near the inlet. The most unstable global mode calculated with the continuous-adjoint has no numerical artifacts but a slightly different eigenvalue. The eigenvalues converge, however, as the timestep reduces. Apart from the numerical artifacts, the mode shapes are very similar, which supports the expectation that they are otherwise equivalent. The continuous-adjoint requires less resolution and usually converges more quickly than the discrete-adjoint but is more challenging to implement. Finally, the direct and adjoint global modes are combined in order to calculate the wavemaker region of a low density jet. © 2011 Elsevier Inc.

Energy stable and momentum conserving hybrid finite element method for the incompressible Navier-Stokes equations

A hybrid method for the incompressible Navier-Stokes equations is presented. The method inherits the attractive stabilizing mechanism of upwinded discontinuous Galerkin methods when momentum advection becomes significant, equal-order interpolations can be used for the velocity and pressure fields, and mass can be conserved locally. Using continuous Lagrange multiplier spaces to enforce flux continuity across cell facets, the number of global degrees of freedom is the same as for a continuous Galerkin method on the same mesh. Different from our earlier investigations on the approach for the Navier-Stokes equations, the pressure field in this work is discontinuous across cell boundaries. It is shown that this leads to very good local mass conservation and, for an appropriate choice of finite element spaces, momentum conservation. Also, a new form of the momentum transport terms for the method is constructed such that global energy stability is guaranteed, even in the absence of a pointwise solenoidal velocity field. Mass conservation, momentum conservation, and global energy stability are proved for the time-continuous case and for a fully discrete scheme. The presented analysis results are supported by a range of numerical simulations. © 2012 Society for Industrial and Applied Mathematics.

Anti-Stokes photoluminescence in ZnO microcrystal

Low temperature (10 K) strong anti-Stokes photoluminescence (ASPL) of ZnO microcrystal excited by low power cw 532 nm laser is reported here. Energy upconversion of 1.1 eV is obtained in our experiment with no conventional nonlinear effect. Through the study of the normal photoluminescence and temperature dependence of ASPL we conclude that the green band luminescence in ZnO is related to deep donor to valance band transition. Using the two-step two-photon absorption model, we provide a plausible mechanism leading to the ASPL phenomenon in our experiment. (c) 2006 American Institute of Physics.

The relationship between chemical bond properties and Stokes shift of Eu2+ in some silicate host lattices

The covalency of each bond in divalent europium doped hosts CaSiO3, SrSiO3, BaSiO3, Sr2LiSiO4F, Ba5SiO4Cl6 and Ba5SiO4Br6 were calculated by using the complicate crystal chemical bond theory. The relationship between the Stokes shift and the bond properties of Eu2+ in these crystals was discussed. The result demonstrates that, in the isostructural crystals that being doped with Eu2+, there is a more precise connection between the magnitude of Stokes shift and the mean covalency of the dopant site.

Third-order Stokes wave solutions for interfacial internal waves in a three-layer density-stratified fluid

Interfacial internal waves in a three-layer density-stratified fluid are investigated using a singular method, and third-order asymptotic solutions of the velocity potentials and third-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory. as expected, the third-order solutions describe the third-order nonlinear modification and the third-order nonlinear interactions between the interfacial waves. The wave velocity depends on not only the wave number and the depth of each layer but also on the wave amplitude.

Second-order Stokes solutions for internal waves in three-layer density-stratified fluid

In this paper, internal waves in three-layer stratified fluid are investigated by using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory. As expected, the first-order solutions are consistent with ordinary linear theoretical results, and the second-order solutions describe the second-order modification on the linear theory and the interactions between the two interfacial waves. Both the first-order and second-order solutions derived depend on the depths and densities of the three-layer fluid. It is also noted that the solutions obtained from the present work include the theoretical results derived by Umeyama as special cases.