An implicit high-order hybridizable discontinuous Galerkin method for the incompressible Navier–Stokes equations

We present an implicit high-order hybridizable discontinuous Galerkin method for the steady-state and time-dependent incompressible Navier–Stokes equations. The method is devised by using the discontinuous Galerkin discretization for a velocity gradient-pressure–velocity formulation of the incompressible Navier–Stokes equations with a special choice of the numerical traces. The method possesses several unique features which distinguish itself from other discontinuous Galerkin methods. First, it reduces the globally coupled unknowns to the approximate trace of the velocity and the mean of the pressure on element boundaries, thereby leading to a significant reduction in the degrees of freedom. Moreover, if the augmented Lagrangian method is used to solve the linearized system, the globally coupled unknowns become the approximate trace of the velocity only. Second, it provides, for smooth viscous-dominated problems, approximations of the velocity, pressure, and velocity gradient which converge with the optimal order of k + 1 in the L²-norm, when polynomials of degree k?0 are used for all components of the approximate solution. And third, it displays superconvergence properties that allow us to use the above-mentioned optimal convergence properties to define an element-by-element postprocessing scheme to compute a new and better approximate velocity. Indeed, this new approximation is exactly divergence-free, H (div)-conforming, and converges with order k + 2 for k ? 1 and with order 1 for k = 0 in the L²-norm. Moreover, a novel and systematic way is proposed for imposing boundary conditions for the stress, viscous stress, vorticity and pressure which are not naturally associated with the weak formulation of the method. This can be done on different parts of the boundary and does not result in the degradation of the optimal order of convergence properties of the method. Extensive numerical results are presented to demonstrate the convergence and accuracy properties of the method for a wide range of Reynolds numbers and for various polynomial degrees.     

Two-level discretizations of nonlinear closure models for proper orthogonal decomposition

Proper orthogonal decomposition has been successfully used in the reduced-order modeling of complex systems. Its original promise of computationally efficient, yet accurate approximation of coherent structures in high Reynolds number turbulent flows, however, still remains to be fulfilled. To balance the low computational cost required by reduced-order modeling and the complexity of the targeted flows, appropriate closure modeling strategies need to be employed. Since modern closure models for turbulent flows are generally nonlinear, their efficient numerical discretization within a proper orthogonal decomposition framework is challenging. This paper proposes a two-level method for an efficient and accurate numerical discretization of general nonlinear closure models for proper orthogonal decomposition reduced-order models. The two-level method computes the nonlinear terms of the reduced-order model on a coarse mesh. Compared with a brute force computational approach in which the nonlinear terms are evaluated on the fine mesh at each time step, the two-level method attains the same level of accuracy while dramatically reducing the computational cost. We numerically illustrate these improvements in the two-level method by using it in three settings: the one-dimensional Burgers equation with a small diffusion parameter ν = 10^?3, the two-dimensional flow past a cylinder at Reynolds number Re = 200, and the three-dimensional flow past a cylinder at Reynolds number Re = 1000.     

A deconvolution enhancement of the Navier–Stokes-αβ model

A deconvolution enhancement of the Navier–Stokes-αβ model for turbulent flow is introduced. The energy and energy-dissipation rate for the enhanced model are derived. It is also shown that the consistency error, relative to the Navier–Stokes equations, and the microscale of the enhanced model are less than those of the Navier–Stokes-αβ model. The proposed model is used to simulate the Taylor–Green vortex problem and results show a qualitatively improved representation of the mean-square vorticity when compared to the Navier–Stokes-αβ model. Numerical studies of the energy spectrum and the alignment between the vorticity and the eigenvectors of the stretching tensor for three-dimensional turbulent flows with Re = 200 are used to explore the utility of the model. A benchmark problem of a two-dimensional channel flow over a step for Re = 600 also indicates that this model can be applied to more general flows than those involving periodic boundary conditions.     

A local discontinuous Galerkin method for directly solving Hamilton–Jacobi equations

In this paper we propose a new local discontinuous Galerkin method to directly solve Hamilton–Jacobi equations. The scheme is a natural extension of the monotone scheme. For the linear case with constant coefficients, the method is equivalent to the discontinuous Galerkin method for conservation laws. Thus, stability and error analysis are obtained under the framework of conservation laws. For both convex and noneconvex Hamiltonian, optimal (k + 1)th order of accuracy for smooth solutions are obtained with piecewise kth order polynomial approximations. The scheme is numerically tested on a variety of one and two dimensional problems. The method works well to capture sharp corners (discontinuous derivatives) and have the solution converges to the viscosity solution.     

Hybridizable discontinuous Galerkin methods for partial differential equations in continuum mechanics

We present hybridizable discontinuous Galerkin methods for solving steady and time-dependent partial differential equations (PDEs) in continuum mechanics. The essential ingredients are a local Galerkin projection of the underlying PDEs at the element level onto spaces of polynomials of degree k to parametrize the numerical solution in terms of the numerical trace; a judicious choice of the numerical flux to provide stability and consistency; and a global jump condition that enforces the continuity of the numerical flux to arrive at a global weak formulation in terms of the numerical trace. The HDG methods are fully implicit, high-order accurate and endowed with several unique features which distinguish themselves from other discontinuous Galerkin methods. First, they reduce the globally coupled unknowns to the approximate trace of the solution on element boundaries, thereby leading to a significant reduction in the degrees of freedom. Second, they provide, for smooth viscous-dominated problems, approximations of all the variables which converge with the optimal order of k + 1 in the L²-norm. Third, they possess some superconvergence properties that allow us to define inexpensive element-by-element postprocessing procedures to compute a new approximate solution which may converge with higher order than the original solution. And fourth, they allow for a novel and systematic way for imposing boundary conditions for the total stress, viscous stress, vorticity and pressure which are not naturally associated with the weak formulation of the methods. In addition, they possess other interesting properties for specific problems. Their approximate solution can be postprocessed to yield an exactly divergence-free and H(div)-conforming velocity field for incompressible flows. They do not exhibit volumetric locking for nearly incompressible solids. We provide extensive numerical results to illustrate their distinct characteristics and compare their performance with that of continuous Galerkin methods.     

Two-phase electrohydrodynamic simulations using a volume-of-fluid approach: A comment

This comment refers to the article of Tomar et al. [1], which presents a numerical methodology of a continuum surface force formulation for simulating two-phase electrohydrodynamic flows. The present work shows, that due to the diffusive character of the Laplacian equation (? · (??₀E) = 0) with discontinuous physical properties (?(x, y, z)), different averaging methods (arithmetic and harmonic) for the fluid property in the transition region have to be applied. The correct choice of the averaging method depends on the orientation of the flux to the interface.An additional improvement is made by calculating the electric displacement D at the cell faces. This leads to a numerical solution independent of the spatial resolution as well as of the interfacial smearing. Simulation results of two different test cases show that the error of the numerical solution is in the order of machine precision.     

A comparison of vortex and pseudo-spectral methods for the simulation of periodic vortical flows at high Reynolds numbers

We present a validation study for the hybrid particle-mesh vortex method against a pseudo-spectral method for the Taylor–Green vortex at Re_Γ = 1600 as well as in the collision of two antiparallel vortex tubes at Re_Γ = 10,000. In this study we present diagnostics such as energy spectra and enstrophy as computed by both methods as well as point-wise comparisons of the vorticity field. Using a fourth order accurate kernel for interpolation between the particles and the mesh, the results of the hybrid vortex method and of the pseudo-spectral method agree well in both flow cases. For the Taylor–Green vortex, the vorticity contours computed by both methods around the time of the energy dissipation peak overlap. The energy spectrum shows that only the smallest length scales in the flow are not captured by the vortex method.In the second flow case, where we compute the collision of two anti-parallel vortex tubes at Reynolds number 10,000, the vortex method results and the pseudo-spectral method results are in very good agreement up to and including the first reconnection of the tubes. The maximum error in the effective viscosity is about 2.5% for the vortex method and about 1% for the pseudo-spectral method. At later times the flows computed with the different methods show the same qualitative features, but the quantitative agreement on vortical structures is lost.     

An inverse Sturm–Liouville problem with a fractional derivative

In this paper, we numerically investigate an inverse problem of recovering the potential term in a fractional Sturm–Liouville problem from one spectrum. The qualitative behaviors of the eigenvalues and eigenfunctions are discussed, and numerical reconstructions of the potential with a Newton method from finite spectral data are presented. Surprisingly, it allows very satisfactory reconstructions for both smooth and discontinuous potentials, provided that the order α ∈ (1, 2) of fractional derivative is sufficiently away from 2.     

Optimal Runge–Kutta schemes for discontinuous Galerkin space discretizations applied to wave propagation problems

We study the performance of methods of lines combining discontinuous Galerkin spatial discretizations and explicit Runge–Kutta time integrators, with the aim of deriving optimal Runge–Kutta schemes for wave propagation applications. We review relevant Runge–Kutta methods from literature, and consider schemes of order q from 3 to 4, and number of stages up to q + 4, for optimization. From a user point of view, the problem of the computational efficiency involves the choice of the best combination of mesh and numerical method; two scenarios are defined. In the first one, the element size is totally free, and a 8-stage, fourth-order Runge–Kutta scheme is found to minimize a cost measure depending on both accuracy and stability. In the second one, the elements are assumed to be constrained to such a small size by geometrical features of the computational domain, that accuracy is disregarded. We then derive one 7-stage, third-order scheme and one 8-stage, fourth-order scheme that maximize the stability limit. The performance of the three new schemes is thoroughly analyzed, and the benefits are illustrated with two examples. For each of these Runge–Kutta methods, we provide the coefficients for a 2N-storage implementation, along with the information needed by the user to employ them optimally.     

A fully discrete Galerkin method for high frequency exterior acoustic scattering in three dimensions

Standard Galerkin discretization techniques (with locally- or globally-supported basis functions) for boundary integral equations are inefficient for high frequency three dimensional exterior scattering simulations because they require a fixed number of unknowns per wavelength in each dimension, leading to large CPU time and memory requirements to set up the dense Galerkin matrix, with each entry requiring evaluation of multi-dimensional highly oscillatory integrals. In this work, using globally-supported basis functions, we describe an efficient fully discrete Galerkin surface integral equation algorithm for simulating high frequency acoustic scattering by three dimensional convex obstacles that includes a powerful integration scheme for evaluation of four dimensional Galerkin integrals with high-order accuracy. Such high-order order accuracy for various practically relevant frequencies (k ∈ [1, 100,000]) substantially improves on approximations based on standard asymptotic techniques. We demonstrate the efficiency of our algorithm for spherical and non-spherical convex scattering for several wavenumbers 1 ? k ? 100,000 for low to high order prescribed tolerance. Our fully discrete algorithm requires only mild growth in the number of unknowns and CPU time as the frequency increases.     

A high order Discontinuous Galerkin – Fourier incompressible 3D Navier–Stokes solver with rotating sliding meshes

We present the development of a sliding mesh capability for an unsteady high order (order ? 3) h/p Discontinuous Galerkin solver for the three-dimensional incompressible Navier–Stokes equations. A high order sliding mesh method is developed and implemented for flow simulation with relative rotational motion of an inner mesh with respect to an outer static mesh, through the use of curved boundary elements and mixed triangular–quadrilateral meshes.A second order stiffly stable method is used to discretise in time the Arbitrary Lagrangian–Eulerian form of the incompressible Navier–Stokes equations. Spatial discretisation is provided by the Symmetric Interior Penalty Galerkin formulation with modal basis functions in the x–y plane, allowing hanging nodes and sliding meshes without the requirement to use mortar type techniques. Spatial discretisation in the z-direction is provided by a purely spectral method that uses Fourier series and allows computation of spanwise periodic three-dimensional flows. The developed solver is shown to provide high order solutions, second order in time convergence rates and spectral convergence when solving the incompressible Navier–Stokes equations on meshes where fixed and rotating elements coexist.In addition, an exact implementation of the no-slip boundary condition is included for curved edges; circular arcs and NACA 4-digit airfoils, where analytic expressions for the geometry are used to compute the required metrics.The solver capabilities are tested for a number of two dimensional problems governed by the incompressible Navier–Stokes equations on static and rotating meshes: the Taylor vortex problem, a static and rotating symmetric NACA0015 airfoil and flows through three bladed cross-flow turbines. In addition, three dimensional flow solutions are demonstrated for a three bladed cross-flow turbine and a circular cylinder shadowed by a pitching NACA0012 airfoil.     

Temperature dependence of dc corona and charge-carrier entrainment in a gas flow channel

Current–voltage (I–V) data for positive and negative polarity point-to-plane geometries are reported for gas flows transverse to the axis of the emitters. Air and nitrogen flows of 0–5 m/s were considered in the experiments and temperatures in ranges of 213–493 K in nitrogen and 283–493 K in air. Entrainment of charge carriers from the individual corona and positive–negative polarity emitter-pairs was considered towards understanding the static elimination process for downstream targets in the gas flows.     

A generalized wall boundary condition for smoothed particle hydrodynamics

In this paper we present a new formulation of the boundary condition at static and moving solid walls in SPH simulations. Our general approach is both applicable to two and three dimensions and is very simple compared to previous wall boundary formulations. Based on a local force balance between wall and fluid particles we apply a pressure boundary condition on the solid particles to prevent wall penetration. This method can handle sharp corners and complex geometries as is demonstrated with several examples. A validation shows that we recover hydrostatic equilibrium conditions in a static tank, and a comparison of the classical dam break simulation with state-of-the-art results in literature shows good agreement. We simulate various problems such as the flow around a cylinder and the backward facing step at Re = 100 to demonstrate the general applicability of this new method.     

Superconductivity in CeO1−xFxFeAs with upper critical field of 94 T

We have successfully synthesized Ce based oxypnictide superconductors with fluorine doping (CeO_1?xF_xFeAs) by a two step solid state reaction method. Detailed XRD and EDX confirm the crystal structure and chemical compositions. We observe that an extremely high H_c2(0) of 94 T can be achieved in the x = 0.1 composition. This increase in H_c2(0) is accompanied by a decrease in transition temperature (38.4 K in x = 0.1 composition) from 42.5 K for the x = 0.2 phase. The in-plane Ginzburg–Landau coherence length is estimated to be ～27 Å at x = 0.2 suggesting a moderate anisotropy in this class of superconductors. The Seebeck coefficient confirms the majority carrier to be electrons and strong dominance of electron–electron correlations in this multiband superconductor.     

Proper general decomposition (PGD) for the resolution of Navier–Stokes equations

In this work, the PGD method will be considered for solving some problems of fluid mechanics by looking for the solution as a sum of tensor product functions. In the first stage, the equations of Stokes and Burgers will be solved. Then, we will solve the Navier–Stokes problem in the case of the lid-driven cavity for different Reynolds numbers (Re = 100, 1000 and 10,000). Finally, the PGD method will be compared to the standard resolution technique, both in terms of CPU time and accuracy.     

Low loss propagation in slow light photonic crystal waveguides at group indices up to 60

We have designed slow light photonic crystal waveguides operating in a low loss and constant dispersion window of Δλ = 2 nm around λ = 1565 nm with a group index of n_g = 60. We experimentally demonstrate a relatively low propagation loss, of 130 dB/cm, for waveguides up to 800 μm in length. This result is particularly remarkable given that the waveguides were written on an electron-beam lithography tool with a writefield of 100 μm that exhibits stitching errors of typically 10–50 nm. We reduced the impact of these stitching errors by introducing “slow–fast–slow” mode conversion interfaces and show that these interfaces reduce the loss from 320 dB/cm to 130 dB/cm at n_g = 60. This significant improvement highlights the importance of the slow–fast–slow method and shows that high performance slow light waveguides can be realised with lengths much longer than the writing field of a given e-beam lithography tool.     

Lattice Boltzmann simulations of 2D laminar flows past two tandem cylinders

We apply the lattice Boltzmann equation (LBE) with multiple-relaxation-time (MRT) collision model to simulate laminar flows in two-dimensions (2D). In order to simulate flows in an unbounded domain with the LBE method, we need to address two issues: stretched non-uniform mesh and inflow and outflow boundary conditions. We use the interpolated grid stretching method to address the need of non-uniform mesh. We demonstrate that various inflow and outflow boundary conditions can be easily and consistently realized with the MRT-LBE. The MRT-LBE with non-uniform stretched grids is first validated with a number of test cases: the Poiseuille flow, the flow past a cylinder asymmetrically placed in a channel, and the flow past a cylinder in an unbounded domain. We use the LBE method to simulate the flow past two tandem cylinders in an unbounded domain with Re = 100. Our results agree well with existing ones. Through this work we demonstrate the effectiveness of the MRT-LBE method with grid stretching.     

Insulation room for aero-acoustic experiments at moderate Reynolds and low Mach numbers

A non-expensive insulation box for aero-acoustic experiments at moderate Reynolds numbers Re < 2 × 10⁴ and low Mach numbers M < 0.2 is presented. Its performance is evaluated with particular attention to unwanted noise sources inherent to the flow facility. Objective acoustic parameters of the insulation box are assessed.     

A linear stability analysis for nonlinear,grey, thermal radiative transfer problems

We present a new linear stability analysis of three time discretizations and Monte Carlo interpretations of the nonlinear, grey thermal radiative transfer (TRT) equations: the widely used “Implicit Monte Carlo” (IMC) equations, the Carter Forest (CF) equations, and the Ahrens–Larsen or “Semi-Analog Monte Carlo” (SMC) equations. Using a spatial Fourier analysis of the 1-D Implicit Monte Carlo (IMC) equations that are linearized about an equilibrium solution, we show that the IMC equations are unconditionally stable (undamped perturbations do not exist) if α, the IMC time-discretization parameter, satisfies 0.5 < α ? 1. This is consistent with conventional wisdom. However, we also show that for sufficiently large time steps, unphysical damped oscillations can exist that correspond to the lowest-frequency Fourier modes. After numerically confirming this result, we develop a method to assess the stability of any time discretization of the 0-D, nonlinear, grey, thermal radiative transfer problem. Subsequent analyses of the CF and SMC methods then demonstrate that the CF method is unconditionally stable and monotonic, but the SMC method is conditionally stable and permits unphysical oscillatory solutions that can prevent it from reaching equilibrium. This stability theory provides new conditions on the time step to guarantee monotonicity of the IMC solution, although they are likely too conservative to be used in practice. Theoretical predictions are tested and confirmed with numerical experiments.     

Point-contact Andreev-reflection spectroscopy in ReFeAsO1−xFx (Re = La,Sm): Possible evidence for two nodeless gaps

A deep understanding of the character of superconductivity in the recently discovered Fe-based oxypnictides ReFeAsO_1?xF_x (Re = rare-earth) necessarily requires the determination of the number of the gaps and their symmetry in k space, which are fundamental ingredients of any model for the pairing mechanism in these new superconductors. In the present paper, we show that point-contact Andreev-reflection spectroscopy experiments performed on LaFeAsO_1?xF_x (La-1111) polycrystals with T_c ～ 27 K and SmFeAsO_0.8F_0.2 (Sm-1111) polycrystals with T_c ～ 53 K gave differential conductance curves exhibiting two peaks at low bias and two additional structures (peaks or shoulders) at higher bias voltages, an experimental situation quite similar to that observed by the same technique in pure and doped MgB₂. The single-band Blonder–Tinkham–Klapwijk model is totally unable to properly fit the conductance curves, while the two-gap one accounts remarkably well for the shape of the whole experimental dI/dV vs. V curves. These results give direct evidence of two nodeless gaps in the superconducting state of ReFeAsO_1?xF_x (Re = La, Sm): a small gap, Δ₁, smaller than the BCS value (2Δ₁/k_BT_c ～ 2.2–3.2) and a much larger gap Δ₂ which gives a ratio 2Δ₂/k_BT_c ～ 6.5–9.0. In Sm-1111 both gaps close at the same temperature, very similar to the bulk T_c, and follow a BCS-like behaviour, while in La-1111 the situation is more complex, the temperature dependence of the gaps showing remarkable deviations from the BCS behaviour at T close to T_c.The normal-state conductance reproducibly shows an unusual, but different, shape in La-1111 and Sm-1111 with a depression or a hump at zero bias, respectively. These structures survive in the normal state up to T¹ ～ 140 K, close to the temperatures at which structural and magnetic transitions occur in the parent, undoped compound.