Re.ned Modeling of the pressure redistribution/scrambling under consistent consideration of scalar turbulence effects in turbulent combustion

In practical turbulent flow problems of engineering importance the coupling between velocity and scalar turbulence along with the variable density plays a non negligible role. For computations using second moment closure approach, the pressure redistribution/scrambling is the most critical term to be modeled as well known. Almost all existing models consist in rescating models derived on a constant density basis in a density weighted form. With regard to turbulent premixed combustion in fact, the application of such models to a range of transient one‐dimensional and two‐dimensional premixed flames in the flamelet regime has been found to yield unsatisfactory results, see [1]. As pointed out by Sadiki [2], the use of the Favre method must be consistently considered as far as open thermodynamic systems are concerned. Furthermore, the need for maintaining certain invariance properties, physical and mathematical realizability conditions in formulating turbulence models is well accepted. Because turbulent processes are irreversible, these efforts demand a carefull consideration of thermodynamic concepts. Based on the results in [1] and following [2], this work aims to derive a physically consistent formulation of the pressure redistribution/scrambling term under consideration of the variable density. Considering the case of premixed flames, the thermochemistry is included by means of a single reactive scalar ‐ the reaction progress variable. The accuracy of the model extensions proposed is demonstrated by comparing the numerical results with experimental data in opposed jet premixed flame configuration.     

Betti numbers in multidimensional persistent homology are stable functions

Multidimensional persistence mostly studies topological features of shapes by analyzing the lower level sets of vector‐valued functions, called filtering functions. As is well known, in the case of scalar‐valued filtering functions, persistent homology groups can be studied through their persistent Betti numbers, that is, the dimensions of the images of the homomorphisms induced by the inclusions of lower level sets into each other. Whenever such inclusions exist for lower level sets of vector‐valued filtering functions, we can consider the multidimensional analog of persistent Betti numbers. Varying the lower level sets, we obtain that persistent Betti numbers can be seen as functions taking pairs of vectors to the set of non‐negative integers. In this paper, we prove stability of multidimensional persistent Betti numbers. More precisely, we prove that small changes of the vector‐valued filtering functions imply only small changes of persistent Betti numbers functions. This result can be obtained by assuming the filtering functions to be just continuous. Multidimensional stability opens the way to a stable shape comparison methodology based on multidimensional persistence. In order to obtain our stability theorem, some other new results are proved for continuous filtering functions. They concern the finiteness of persistent Betti numbers for vector‐valued filtering functions and the representation via persistence diagrams of persistent Betti numbers, as well as their stability, in the case of scalar‐valued filtering functions. Finally, from the stability of multidimensional persistent Betti numbers, we obtain a lower bound for the natural pseudo‐distance. Copyright © 2013 John Wiley & Sons, Ltd.     

Advected turbulent line thermal driven by concentration difference

The spatial evolution of an advected line thermal driven by concentration difference––a turbulent buoyant body of fluid, for which small density difference is caused by a proportional variation in scalar concentration, horizontally introduced at no excess momentum into a horizontal ambient current, is studied using the standard two-equation k– model with a buoyancy expansion. The numerical results show that the advected line thermal is characterized longitudinally by a flat trajectory with scalar dilution taking place essentially near the jet exit, and transversely by a vortex-pair flow and a kidney-shaped concentration structure with double peak maxima corresponding to stronger buoyancy effect; the computed flow details and scalar mixing characteristics can be described by self-similar relations beyond a dimensionless distance of around 10; the aspect ratio for the kidney-shaped sectional thermal is found to be around 1.2–1.4; the predicted flow feature and mixing rate are well supported by asymptotic dimensional analysis and related experimental data. The analogy between a steady advected line thermal and corresponding time-dependent line thermal is also found reasonable by a special exploration into the horizontal velocity distribution and the significance of horizontal diffusion effect of the advected line thermal. Only about half of the vertical momentum resulting from the buoyancy effect is found contained in the advected line thermal, corresponding to an added virtual mass coefficient of approximately 1 for the sectional thermal.     

Large-eddy-simulation of 3-dimensional Rayleigh-Taylor instability in incompressible fluids

The 3-dimensional incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation ( LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinusoidal perturbation and random perturbation are simulated. A full treatment of the whole evolution process of the instability is addressed. To verify the reliability of the LES code, the averaged turbulent energy as well as the flux of passive scalar are calculated at both the resolved scale and the subgrid scale. Our results show good agreement with the experimental and other numerical work. The LES method has proved to be an effective approach to the Rayleigh-Taylor instability.     

Non‐Gaussian invariant measures for the Majda model of decaying turbulent transport

The problem of turbulent transport of a scalar field by a random velocity field is considered. The scalar field amplitude exhibits rare but very large fluctuations whose typical signature are fatter than Gaussian tails for the probability distribution of the scalar. The existence of such large fluctuations is related to clustering phenomena of the Lagrangian paths within the flow. This suggests an approach to turn the large deviation problem for the scalar field into a small deviation, or small ball, problem for some appropriately defined process measuring the spreading with time of the Lagrangian paths. Here, such a methodology is applied to a model proposed by Majda consisting of a white‐in‐time linear shear flow and some generalizations of it where the velocity field has finite, or even infinite, correlation time. The non‐Gaussian invariant measure for the (reduced) scalar field is derived and, in particular, it is shown that the one‐point distribution of the scalar has stretched exponential tails, with a stretching exponent depending of the parameters in the model. Different universality classes for the scalar behavior are identified which, all other parameters being kept fixed, display a one‐to‐one correspondence with a exponent measuring time persistence effects in the velocity field. © 2001 John Wiley & Sons, Inc.     

Quaternionic spherical wave functions

The scalar spherical wave functions (SWFs) are solutions to the scalar Helmholtz equation obtained by the method of separation of variables in spherical polar coordinates. These functions are complete and orthogonal over a sphere, and they can, therefore, be used as a set of basis functions in solving boundary value problems by spherical wave expansions. In this work, we show that there exists a theory of functions with quaternionic values and of three real variables, which is determined by the Moisil–Theodorescu‐type operator with quaternionic variable coefficients, and which is intimately related to the radial, angular and azimuthal wave equations. As a result, we explain the connections between the null solutions of these equations, on one hand, and the quaternionic hyperholomorphic and anti‐hyperholomorphic functions, on the other. We further introduce the quaternionic spherical wave functions (QSWFs), which refine and extend the SWFs. Each function is a linear combination of SWFs and products of ‐hyperholomorphic functions by regular spherical Bessel functions. We prove that the QSWFs are orthogonal in the unit ball with respect to a particular bilinear form. Also, we perform a detailed analysis of the related properties of QSWFs. We conclude the paper establishing analogues of the basic integral formulae of complex analysis such as Borel–Pompeiu's and Cauchy's, for this version of quaternionic function theory. As an application, we present some plot simulations that illustrate the results of this work. Copyright © 2016 John Wiley & Sons, Ltd.     

Revisiting the Concentration Problem of Vector Fields within a Spherical Cap: A Commuting Differential Operator Solution

We propose a novel basis of vector functions, the mixed vector spherical harmonics that are closely related to the functions of Sheppard and Török and help us reduce the concentration problem of tangential vector fields within a spherical cap to an equivalent scalar problem. Exploiting an analogy with previous results published by Grünbaum, Longhi and Perlstadt, we construct a differential operator that commutes with the concentration operator of this scalar problem and propose a stable and convenient method to obtain its eigenfunctions. Having obtained the scalar eigenfunctions, the calculation of tangential vector Slepian functions is straightforward.     

Asymptotic behaviour of commutation errors and the divergence of the Reynolds stress tensor near the wall in the turbulent channel flow

The derivation of the space averaged Navier–Stokes equations for the large eddy simulation (LES) of turbulent incompressible flows introduces two groups of terms which do not depend only on the space averaged flow field variables: the divergence of the Reynolds stress tensor and commutation errors. Whereas the former is studied intensively in the literature, the latter terms are usually neglected. This note studies the asymptotic behaviour of these terms for the turbulent channel flow at a wall in the case that the commutation errors arise from the application of a non‐uniform box filter. To perform analytical calculations, the unknown flow field is modelled by a wall law (Reichardt law and 1/αth power law) for the mean velocity profile and highly oscillating functions model the turbulent fluctuations. The asymptotics show that near the wall, the commutation errors are at least as important as the divergence of the Reynolds stress tensor. Copyright © 2006 John Wiley & Sons, Ltd.     

Domain decomposition preconditioners for multiscale problems in linear elasticity

We analyze two‐level overlapping Schwarz domain decomposition methods for vector‐valued piecewise linear finite element discretizations of the PDE system of linear elasticity. The focus of our study lies in the application to compressible, particle‐reinforced composites in 3D with large jumps in their material coefficients. We present coefficient‐explicit bounds for the condition number of the two‐level additive Schwarz preconditioned linear system. Thereby, we do not require that the coefficients are resolved by the coarse mesh. The bounds show a dependence of the condition number on the energy of the coarse basis functions, the coarse mesh, and the overlap parameters, as well as the coefficient variation. Similar estimates have been developed for scalar elliptic PDEs by Graham et al. 1 The coarse spaces to which they apply here are assumed to contain the rigid body modes and can be considered as generalizations of the space of piecewise linear vector‐valued functions on a coarse triangulation. The developed estimates provide a concept for the construction of coarse spaces, which can lead to preconditioners that are robust with respect to high contrasts in Young's modulus and the Poisson ratio of the underlying composite. To confirm the sharpness of the theoretical findings, we present numerical results in 3D using vector‐valued linear, multiscale finite element and energy‐minimizing coarse spaces. The theory is not restricted to the isotropic (Lamé) case, extends to the full‐tensor case, and allows applications to multiphase materials with anisotropic constituents in two and three spatial dimensions. However, the bounds will depend on the ratio of largest to smallest eigenvalue of the elasticity tensor.     

A joint SSPDF–CMC method for simulating turbulent methane–air jet combustion

A joint single scalar probability density function and conditional moment closure (SSPDF–CMC) method is proposed for modeling a turbulent methane–air jet flame. In general, the probability density function (PDF) of passive scalar (such as mixture fraction) is non-Gaussian and not fully determined by the advecting velocity field, therefore the presumed shape of PDF of mixture fraction assumed as clipped Gaussian distribution or beta function in normal conditional moment closure (CMC) method is incorrect. In SSPDF–CMC method, the PDF of mixture fraction is obtained using a Monte-Carlo method to solve a PDF transport equation. An assumption that the averaged scalar advection is approximately equal to the averaged scalar dissipation in the wake of a grid-generated turbulence flow is adopted to model the averaged scalar dissipation. The predictions using the proposed method are compared with those using the conventional CMC method and the experimental data. It is seen that the predicted Favre conditional averaged statistics and Favre unconditional averaged statistics using the proposed method are in better agreement with the measurement data than those using the conventional CMC method. The predicted conditional or unconditional mean NO even using the SSPDF model is only in fair agreement with the experiments. It shows that the first-order closure for the conditional reaction rate of NO should be improved.     

Spatiospectral concentration of vector fields on a sphere

We construct spherical vector bases that are bandlimited and spatially concentrated, or, alternatively, spacelimited and spectrally concentrated, suitable for the analysis and representation of real-valued vector fields on the surface of the unit sphere, as arises in the natural and biomedical sciences, and engineering. Building on the original approach of Slepian, Landau, and Pollak we concentrate the energy of our function bases into arbitrarily shaped regions of interest on the sphere, and within certain bandlimits in the vector spherical-harmonic domain. As with the concentration problem for scalar functions on the sphere, which has been treated in detail elsewhere, a Slepian vector basis can be constructed by solving a finite-dimensional algebraic eigenvalue problem. The eigenvalue problem decouples into separate problems for the radial and tangential components. For regions with advanced symmetry such as polar caps, the spectral concentration kernel matrix is very easily calculated and block-diagonal, lending itself to efficient diagonalization. The number of spatiospectrally well-concentrated vector fields is well estimated by a Shannon number that only depends on the area of the target region and the maximal spherical-harmonic degree or bandwidth. The spherical Slepian vector basis is doubly orthogonal, both over the entire sphere and over the geographic target region. Like its scalar counterparts it should be a powerful tool in the inversion, approximation and extension of bandlimited fields on the sphere: vector fields such as gravity and magnetism in the earth and planetary sciences, or electromagnetic fields in optics, antenna theory and medical imaging.     

Matrix‐valued Szegő polynomials and quantum random walks

We consider quantum random walks (QRW) on the integers, a subject that has been considered in the last few years in the framework of quantum computation. We show how the theory of CMV matrices gives a natural tool to study these processes and to give results that are analogous to those that Karlin and McGregor developed to study (classical) birth‐and‐death processes using orthogonal polynomials on the real line. In perfect analogy with the classical case, the study of QRWs on the set of nonnegative integers can be handled using scalar‐valued (Laurent) polynomials and a scalar‐valued measure on the circle. In the case of classical or quantum random walks on the integers, one needs to allow for matrix‐valued versions of these notions. We show how our tools yield results in the well‐known case of the Hadamard walk, but we go beyond this translation‐invariant model to analyze examples that are hard to analyze using other methods. More precisely, we consider QRWs on the set of nonnegative integers. The analysis of these cases leads to phenomena that are absent in the case of QRWs on the integers even if one restricts oneself to a constant coin. This is illustrated here by studying recurrence properties of the walk, but the same method can be used for other purposes. The presentation here aims at being self‐contained, but we refrain from trying to give an introduction to quantum random walks, a subject well surveyed in the literature we quote. For two excellent reviews, see [1, 9]. See also the recent notes [20]. © 2009 Wiley Periodicals, Inc.     

A nonlinear inequality and application to global asymptotic stability of perturbed systems

The goal of this work is to present a new nonlinear inequality which is used in a study of the Lyapunov uniform stability and uniform asymptotic stability of solutions to time‐varying perturbed differential equations. New sufficient conditions for global uniform asymptotic stability and/or practical stability in terms of Lyapunov‐like functions for nonlinear time‐varying systems is obtained. Our conditions are expressed as relation between the Lyapunov function and the existence of specific function which appear in our analysis through the solution of a scalar differential equation. Moreover, an example in dimensional two is given to illustrate the applicability of the main result. Copyright © 2014 John Wiley & Sons, Ltd.     

Characterizations for Optimality Conditions of General Robust Optimization Problems

In this paper, by virtue of the image space analysis, the general scalar robust optimization problems under the strictly robust counterpart are considered, among which, the uncertainties are included in the objective as well as the constraints. Besides, on the strength of a corrected image in a new type, an equivalent relation between the uncertain optimization problem and its image problem is also established, which provides an idea to tackle with minimax problems. Furthermore, theorems of the robust weak alternative as well as sufficient characterizations of robust optimality conditions are achieved on the frame of the linear and nonlinear (regular) weak separation functions. Moreover, several necessary and sufficient optimality conditions, especially saddle point sufficient optimality conditions for scalar robust optimization problems, are obtained. Finally, a simple example for finding a shortest path is included to show the effectiveness of the results derived in this paper.     

On quaternionic analysis for the Schrödinger operator with a particular potential and its relation with the Mathieu functions

It has been found recently that there exists a theory of functions with quaternionic values and in two real variables, which is determined by a Cauchy–Riemann‐type operator with quaternionic variable coefficients, and that is intimately related to the so‐called Mathieu equations. In this work, it is all explained as well as some basic facts of the arising quaternionic function theory. We establish analogues of the basic integral formulas of complex analysis such as Borel–Pompeiu's, Cauchy's, and so on, for this version of quaternionic function theory. This theory turns out to be in the same relation with the Schrödinger operator with special potential as the usual holomorphic functions in one complex variable, or quaternionic hyperholomorphic functions, or functions of Clifford analysis, are with the corresponding Laplace operator. Moreover, it is similar to that of α‐hyperholomorphic functions and the Helmholtz operator. Copyright © 2012 John Wiley & Sons, Ltd.     

Local approximation of functions with several variables

This paper presents the problem of local approximation of scalar functions with several variables, including points of non-differentiability. The procedure considers that the mapping displays rates of change of power type with respect to linear changes in the coordinate domain, and the exponents are not necessarily integer. The approach provides a formula describing the local variability of scalar fields which contains and generalizes Taylor’s formula of first order. The functions giving the contact are Müntz polynomials. The knowledge of their coefficients and exponents enables the finding of local extremes including cases of non-smoothness. Sufficient conditions for the existence of global maxima and minima of concave-convex functions are obtained as well.     

Modeling Sediment Deposition Patterns Arising From Suddenly Released Fixed-Volume Turbulent Suspensions

Models presented in several recent papers [1–3] dealing with particle transport by, and deposition from, bottom gravity currents produced by the sudden release of dilute, well‐mixed fixed‐volume suspensions have been relatively successful in duplicating the experimentally observed long‐time, distal, areal density of the deposit on a rigid horizontal bottom. These models, however, fail in their ability to capture the experimentally observed proximal pattern of the areal density with its pronounced dip in the region initially occupied by the well‐mixed suspension and its equally pronounced local maximum at roughly the one‐third point of the total reach of the deposit. The central feature of the models employed in [1–3] is that the particles are always assumed to be vertically well‐mixed by fluid turbulence and to settle out through the bottom viscous sublayer with the Stokes settling velocity for a fluid at rest with no re‐entrainment of particles from the floor of the tank. Because this process is assumed from the outset in the models of [1–3], the numerical simulations for a fixed‐volume release will not take into account the actual experimental conditions that prevail at the time of release of a well‐mixed fixed‐volume suspension. That is, owing to the vigorous stirring that produces the well‐mixed suspension, the release volume will initially possess greater turbulent energy than does an unstirred release volume, which may only acquire turbulent energy as a result of its motion after release through various instability mechanisms. The eddy motion in the imposed fluid turbulence reduces the particle settling rates from the values that would be observed in an unstirred release volume possessing zero initial turbulent energy. We here develop a model for particle bearing gravity flows initiated by the sudden release of a fixed‐volume suspension that takes into account the initial turbulent energy of mixing in the release volume by means of a modified settling velocity that, over a time scale characteristic of turbulent energy decay, approaches the full Stokes settling velocity. Thereafter, in the flow regime, we assume that the turbulence persists and, in accord with current understanding concerning the mechanics of dense underflows, that this turbulence is most intense in the wall region at the bottom of the flow and relatively coarse and on the verge of collapse (see [22]) at the top of the flow where the density contrast is compositionally maintained. We capture this behavior by specifying a “shape function” that is based upon experimental observations and provides for vertical structure in the volume fraction of particles present in the flow. The assumption of vertically well‐mixed particle suspensions employed in [1–5] corresponds to a constant shape function equal to unity. Combining these two refinements concerning the settling velocity and vertical structure of the volume fraction of particles into the conservation law for particles and coupling this with the fluid equations for a two‐layer system, we find that our results for areal density of deposits from sudden releases of fixed‐volume suspensions are in excellent qualitative agreement with the experimentally determined areal densities of deposit as reported in [1, 3, 6]. In particular, our model does what none of the other models do in that it captures and explains the proximal depression in the areal density of deposit.     

AN AUTOMATED METHOD FOR EXTENDING SENSITIVITY ANALYSIS TO MODEL FUNCTIONS

Abstract Parameter sensitivity analysis is a relatively well‐developed field compared with function sensitivity analysis. How sensitive are model conclusions to the choice of functions used in the right‐hand side of difference or differential equation models? Most work in this area has been scenario based where alternative functions are tested. This requires knowledge of plausible alternatives and is usually restricted to well‐known functions. In this work, a method is proposed for function sensitivity analysis, which it is hoped will provide a starting point for discussion and stimulating further research on this important but neglected topic. The method provides information on the sensitivity of a model's results to changes in the shape of the functions. This is done in an automated way without the need to specify alternative functional forms. The proposed method is illustrated on a small ecosystem model.     

Error Estimates for the Adaptive Computation of a Scalar Three Well Problem

We investigate the numerical approximation of Young measure solutions appearing as generalised solutions in scalar non‐convex variational problems. A priori and a posteriori error estimates for a macroscopic quantity, i.e., the stress, are given. Numerical experiments for a scalar three well problem, occurring as a subproblem in the theory of phase transitions in crystalline solids, show that the computational effort can be significantly reduced using an adaptive mesh‐refinement strategy combined with an active set technique by Carstensen and Roubíček.