Source author record

Nadeem A. Malik

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

physics.flu-dyn physics.comp-ph math-ph math.MP math.GM physics.chem-ph physics.geo-ph

Catalog footprint

Research overview

8works

7topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2020arXiv

A Comparison of Turbulence Generated by 3DS Sparse Grids With Different Blockage Ratios and Different Co-Frame Arrangements

A new type of grid turbulence generator, the 3D sparse grid (3DS), is a co-planar arrangement of co-frames each containing a different length scale of grid elements [Malik, N. A. US Patent No. US 9,599,269 B2 (2017)] and possessing a much bigger parameter space than the flat 2D fractal square grid (2DF). Using DNS we compare the characteristics of the turbulence (mean flow, turbulence intensity, energy spectrum) generated by different types of 3DS grids. The peak intensities generated by 3DS can exceed the peaks generated by the 2DF by 80\%; we observe that a 3DS with blockage ratio 24\% produces turbulence similar to the 2DF with blockage ratio 32\% implying lower energy input for the same turbulence.

preprint2020arXiv

Objective quantification of Particle Pair Diffusion in Homogeneous Isotropic Turbulence

Turbulence consists of interacting flow structures covering a wide range of length and time scales. A long-standing question looms over pair diffusion of particles in close proximity i.e. particle pair diffusion at small separations: what range of turbulence length scales governs pair diffusion? Here, we attempt to answer this question by addressing pair diffusion by both fine scales and larger scale coherent structures in which the fine scales are embedded - we unavoidably encounter a combination of both local and non-local interactions associated with the small and large length scales. The local structures possess length scales of the same order of magnitude as the pair separation $l$ , and they induce strong relative motion between the particle pair. However, the non-local structures, possessing length scales much larger than $l$ , also induce (via Biot-Savart) significant relative motion, an effect ignored in prior studies (based on Richardson-Obukhov R-O theory). This fundamentally changes the interpretation of the pair diffusion process, giving the pair diffusivity $K$ scaling as $K\sim l^{1.556}$ -- agreeing within $1\%$ of data. The `R-O constant' $g_l$ is shown to be not a constant, although widely assumed to be a constant. However, new constants (representing pair diffusivity $G_K$ and pair separation $G_l$ ) are identified, which we show to asymptote to $G_K\approx 0.73$ and $G_l\sim 0.01$ at high Reynolds numbers. As an application, we show that the radius of a cloud of droplets in a spray is smaller by an order of magnitude as compared to R-O theory.

preprint2016arXiv

A numerical study of the the response of transient inhomogeneous flames to pressure fluctuations and negative stretch in contracting hydrogen/air flames

Transient premixed hydrogen/air flames contracting through inhomogeneous fuel distributions and subjected to stretch and pressure oscillations are investigated numerically using an implicit method which couples the fully compressible flow to the realistic chemistry and multicomponent transport properties. The impact of increasing {\em negative} stretch is investigated through the use of planar, cylindrical and spherical geometries, and a comparison with the results from {\em positively} stretched expanding H2/air flames (MALIK2010) and CH4/air flames (MALIK2012a) is made. A flame relaxation number $n_R=τ_R/τ_L$ ( $τ_R$ is the time that the flame takes to return to the mean equilibrium conditions after initial disturbance; $τ_L$ is a flame time scale) decreases by 10\% with increasing {\em negative} stretch, in contrast to the two expanding flames where $n_R$ decreased by 40\% with increasing {\em positive} stretch. $n_R$ appears to much more sensitive to variations in positive/negative curvature than to the thermo-chemistry of different flame types. $n_R$ may thus be a useful indicator of the strength of flame-curvature coupling. The spectra of pressure fluctuations $E_p(ω)$ scale close to $\simω^{-3}$ , which is steeper than in the expanding H2/air flames where $E_p(ω)\simω^{-2}$ . Rapid transport ('flapping') of the flame front by the largest convective velocities induced by the random pressure fluctuations is prominent because of the lack of gas velocity ahead of the flame, $u_g=0$ . 'Memory effects' between fuel consumption and the rate of heat release is obscured by the flapping, although the contracting flames display short time-lags in the mean.

preprint2016arXiv

A realistic transport model with pressure dependent parameters for gas flow in tight porous media with application to determining shale rock properties

Shale gas recovery has seen a major boom in recent years due to the increasing global energy demands; but the extraction technologies are very expensive. It is therefore important to develop realistic transport modelling and simulation methods, for porous rocks and porous media, that can compliment the field work. Here, a new nonlinear transport model for single phase gas flow in tight porous media is derived, incorporating many important physical processes that occur in such porous systems: continuous flow, transition flow, slip flow, Knudsen diffusion, adsorption and desorption in to and out of the rock material, and a correction for high flow rates (turbulence). This produces a nonlinear advection-diffusion type of partial differential equation (PDE) with pressure dependent model parameters and associated compressibility coefficients, and highly nonlinear apparent convective flux (velocity) and apparent diffusivity. An important application is to the determination of shale rock properties, such as porosity and permeability, by history matching of the the simulation results to data from pressure-pulse decay tests in a shale rock core sample [Pong K., Ho C., Liu J., Tai Y. Non-linear pressure distribution in uniform microchannels. ASME Fluids Eng. Div. (FED) Vol. 197, 51--56, (1994)]. The estimates of the rock porosity and the permeability from our model simulations are realistic of shale rocks, more realistic than obtained from previous models, and illustrates the potential of the modelling strategy presented here in producing accurate simulations of shale gas flow in tight reservoirs.

preprint2016arXiv

On Turbulent Particle Pair Diffusion

Richardson's theory of turbulent particle pair diffusion [Richardson, L. F. Proc. Roy. Soc. Lond. A 100, 709--737, 1926], based upon observational data, is equivalent to a locality hypothesis in which the turbulent pair diffusivity $(K)$ scales with the pair separation $(σ_l)$ with a 4/3-power law, $K\sim σ_l^{4/3}$ . Here, a reappraisal of the 1926 dataset reveals that one of the data-points is from a molecular diffusion context; the remaining data from geophysical turbulence display an unequivocal non-local scaling, $K \sim σ_l^{1.564}$ . Consequently, the foundations of pair diffusion theory have been re-examined, leading to a new theory based upon the principle that both local and non-local diffusional processes govern pair diffusion in homogeneous turbulence. Through a novel mathematical approach the theory is developed in the context of generalised power law energy spectra, $E(k)\sim k^{-p}$ for $1<p\le 3$ , over extended inertial subranges. The theory predicts the scaling, $K(p)\sim σ_l^{γ_p}$ , with $γ_p$ intermediate between the purely local and the purely non-local scalings, i.e. $(1+p)/2<γ_p\le 2$ . A Lagrangian diffusion model, Kinematic Simulations [Kraichnan, R. H., Phys. Fluids 13, 22-31, 1970; Fung et al., J. Fluid Mech. 236, 281-318, 1992], is used to examine the predictions of the new theory all of which are confirmed. The simulations produce the scalings, $K\sim σ_l^{1.545}$ to $\sim σ_l^{1.570}$ , in the accepted range of intermittent turbulence spectra, $E(k)\sim k^{-1.72}$ to $\sim k^{-1.74}$ , in close agreement with the revised 1926 dataset.

preprint2016arXiv

Turbulent Particle Pair Diffusion, Locality Versus Non-locality: Numerical Simulations

Particle pair (relative) diffusion in a field of homogeneous turbulence with generalised power-law energy spectra, $E(k)\sim k^{-p}$ for $1< p\le 3$ and $k_1\le k\le k_η$ with $k_η/k_1=10^6$ , is investigated numerically using Kinematic Simulation (Kraichnan 1970, Fung et al 1992). If $Δ=|{\bf x_2}(t)-{\bf x}_1(t)|$ and $σ_Δ=\sqrt{\langleΔ^2\rangle}$ , (the angled brackets is the ensemble average over particle pairs), then we find that: (1) The pair diffusivity scales like $D_p\sim σ_Δ^{γ_p}$ with $γ_p$ obtained from the simulations such that $γ^l_p<γ_p<2$ , where $γ^l_p =(1+p)/2$ is the Richardson locality scaling. The range of $Δ$ over which these scalings are observed diminishes from above and from below as $p$ increases towards $3$ .\\ (2) $M(p)=γ_p/γ^l_p>1$ in the range $1<p<3$ , and $M$ has a peak at $p_m\approx 1.8$ . This suggests that for spectra close to this in the range $1.5<p<2$ , which includes Kolmogorov turbulence, local and non-local correlations play comparable roles in the pair diffusion process.\\ (3) The mean square separation scales like $\langleΔ^2\rangle_p \sim τ_p^{χ_p}$ where $χ_p=1/(1-γ_p/2)$ and $τ_p$ is an adjusted travel time.\\ (4) For Kolmogorov turbulence $p=5/3$ , we observe $D_{5/3}\sim σ_Δ^{1.53}$ , and $\langleΔ^2\rangle_{5/3}\simτ^{4.2}$ . \\ (5) At $p_*\approx 1.4$ ( $E(k)\sim k^{-1.4}$ ) we observe $D_{p_*}\simσ_Δ^{4/3}$ , and $\langleΔ^2\rangle_{p_*}\simτ^3$ ; these are different to the Richardson $4/3$ -law and Richardson-Obukov $t^3$ -regime which occur for $p=5/3$ . These results are consistent with Malik's (2014) theory, and supports the principle upon which the theory is based that both local and non-local correlations are effective in the pair diffusion process inside the inertial subrange.

preprint2015arXiv

Residual sweeping effect in turbulent particle pair diffusion in a Lagrangian diffusion model

Thomson, D. J. \& Devenish, B. J. [{\em J.~Fluid Mech.} 526, 277 (2005)] and others have suggested that sweeping effects make Lagrangian properties in Kinematic Simulations (KS), Fung et al [Fung J. C. H., Hunt J. C. R., Malik N. A. \& Perkins R. J. {\em J.~Fluid Mech.} 236, 281 (1992)], unreliable. Here it is shown through a novel analysis based upon analysing pairs of particle trajectories in a frame of reference moving with the large energy containing scales of motion that the normalized integrated error $e^I_K$ in the turbulent pair diffusivity ( $K$ ) due to the sweeping effect decreases with increasing pair separation ( $σ_l$ ), such that $e^I_K\to 0$ as $σ_l/η\to \infty$ ; and $e^I_K\to \infty$ as $σ_l/η\to 0$ . $η$ is the Kolmogorov turbulence microscale. There is an intermediate range of separations $1<σ_l/η< \infty$ in which the error $e^I_K$ remains negligible. Simulations using KS shows that in the swept frame of reference, this intermediate range is large covering almost the entire inertial subrange simulated, $1<σ_l/η< 10^5$ , implying that the deviation from locality observed in KS therefore cannot be atributed to sweeping errors and could be real. This is important for pair diffusion theory and modeling.

preprint2014arXiv

Approximate solutions of a time-fractional diffusion equation with a source term using the variational iteration method

We consider a time fractional differential equation of order $α$ , $0<α<1$ , $\frac{\partial c(x,t)}{\partial t}={}^C_0\mathcal{D}_t^α[(Ac)(x,t)]+q(x,t) ,\quad x > 0, t > 0, \quad c(x,0)=f(x).$ where ${}^C_0\mathcal{D}_t^α$ is the Caputo fractional derivative of order $α$ , $A$ is a linear differential operator, $q(x,t)$ is a source term, and $f(x)$ is the inital condition. Approximate (truncated) series solutions are obtained by means of the Variational Iteration Method (VIM). We find the series solutions for different cases of the source term, in a form that is readily implementable on the computer where symbolic computation platform is available. The error in truncated solution $c_n$ diminishes exponentially fast for a given $α$ as the number of terms in the series increases. VIM has several advantages over other methods that produce solutions in the series form. The truncated VIM solutions often converge rapidly requiring only a few terms for fast and accurate approximations.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Nadeem A. Malik

Research overview

Is this your work?

8 publications

A Comparison of Turbulence Generated by 3DS Sparse Grids With Different Blockage Ratios and Different Co-Frame Arrangements

Objective quantification of Particle Pair Diffusion in Homogeneous Isotropic Turbulence

A numerical study of the the response of transient inhomogeneous flames to pressure fluctuations and negative stretch in contracting hydrogen/air flames

A realistic transport model with pressure dependent parameters for gas flow in tight porous media with application to determining shale rock properties

On Turbulent Particle Pair Diffusion

Turbulent Particle Pair Diffusion, Locality Versus Non-locality: Numerical Simulations

Residual sweeping effect in turbulent particle pair diffusion in a Lagrangian diffusion model

Approximate solutions of a time-fractional diffusion equation with a source term using the variational iteration method

See the researcher in context

Building this map preview