Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUID AMPLIFIER”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Efficiency of Water Vapor on Top‐of‐Atmosphere Radiation

Earth's climate sensitivity is greatly affected by the compensation between temperature feedback and water vapor (WV) feedback. Using abrupt 4xCO 2 experiments, we show that the global-mean WV feedback is nearly a linear function of the temperature feedback, the slope of which is explained by the longwave radiative efficiency of WV (ϵ). Although ϵ remains constant across models in the global mean, it exhibits substantial spatial variations and is particularly weak in Antarctica, where near-surface inversions decouple the surface from the free troposphere. We introduce a surface–free troposphere temperature difference (SFTD) metric, showing that positive SFTD (e.g., high lifting condensation level) amplifies ϵ, while negative SFTD (e.g., strong surface inversion) suppresses it. These findings provide a clear explanation of how local climate conditions modulate the radiative compensation between temperature and WV feedbacks.

climate feedback↗
DOE OSTI2019

Analog simulator of integro-differential equations with classical memristors

An analog computer makes use of continuously changeable quantities of a system, such as its electrical, mechanical, or hydraulic properties, to solve a given problem. While these devices are usually computationally more powerful than their digital counterparts, they suffer from analog noise which does not allow for error control. We will focus on analog computers based on active electrical networks comprised of resistors, capacitors, and operational amplifiers which are capable of simulating any linear ordinary differential equation. However, the class of nonlinear dynamics they can solve is limited. In this work, by adding memristors to the electrical network, we show that the analog computer can simulate a large variety of linear and nonlinear integro-differential equations by carefully choosing the conductance and the dynamics of the memristor state variable. We study the performance of these analog computers by simulating integro-diferential models related to fluid dynamics, nonlinear Volterra equations for population growth, and quantum models describing non-Markovian memory effects, among others. Finally, we perform stability tests by considering imperfect analog components, obtaining robust solutions with up to 13% relative error for relevant timescales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experimental studies of high-confinement mode plasma response to non-axisymmetric magnetic perturbations in ASDEX Upgrade

Here, the interaction of externally applied small non-axisymmetric magnetic perturbations (MP) with tokamak high-confinement mode (H-mode) plasmas is reviewed and illustrated by recent experiments in ASDEX Upgrade. The plasma response to the vacuum MP field is amplified by stable ideal kink modes with low toroidal mode number n driven by the H-mode edge pressure gradient (and associated bootstrap current) which is experimentally evidenced by an observable shift of the poloidal mode number m away from field alignment (m = qn, with q being the safety factor) at the response maximum. A torque scan experiment demonstrates the importance of the perpendicular electron flow for shielding of the resonant magnetic perturbation, as expected from a two-fluid MHD picture. Two significant effects of MP occur in H-mode plasmas at low pedestal collisionality, $\nu _{\text{ped}}^{\ast}\leqslant 0.4$ : (a) a reduction of the global plasma density by up to $61 \% $ and (b) a reduction of the energy loss associated with edge localised modes (ELMs) by a factor of up to 9. A comprehensive database of ELM mitigation pulses at low ${{\nu}^{\ast}}$ in ASDEX Upgrade shows that the degree of ELM mitigation correlates with the reduction of pedestal pressure which in turn is limited and defined by the onset of ELMs, i. e. a modification of the ELM stability limit by the magnetic perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

INDUCED MAXIMUM MAGNETIC FIELD IN COSMIC OUTFLOW SYSTEM BY A RELATIVISTIC CURRENT FILAMENTATION INSTABILITY: EXACT ANALYTICAL MODEL

We present an analytical study of current filamentation instability (CFI) in a fully relativistic cold plasma system, including arbitrary currents. For our purposes, we employ the cold fluid equations, together with Maxwell's equations as well as the plasma shell concept and boost frame method, to obtain an exact solution of the instability growth rate. A simple relation is found for the maximum growth rate of the CFI (for any arbitrary current system), which remarkably is used to calculate the large magnitude of an induced magnetic field in astrophysical environments such as active galactic nuclei (AGNs), microquasars, supernova remnants (SNRs), and stellar winds. We find that the magnetic field is amplified in the SNR up to the level required to justify the recent discovery of the year-scale variability in the X-ray emission of SNRs. Also, the maximum magnetic field of two and three orders higher (or two orders lower) than that of the SNR has been derived for microquasars and AGNs (or stellar winds), respectively. Moreover, making use of the exact analytical solution of the CFI, it is shown that the maximum magnetic field up to around 10{sup 8} G can be detected from a classical cold counterstreaming system after a saturation time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Buoyancy instability of homologous implosions

With this study, I consider the hydrodynamic stability of imploding ideal gases as an idealized model for inertial confinement fusion capsules, sonoluminescent bubbles and the gravitational collapse of astrophysical gases. For oblate modes (short-wavelength incompressive modes elongated in the direction of the mean flow), a second-order ordinary differential equation is derived that can be used to assess the stability of any time-dependent flow with planar, cylindrical or spherical symmetry. Upon further restricting the analysis to homologous flows, it is shown that a monatomic gas is governed by the Schwarzschild criterion for buoyant stability. Under buoyantly unstable conditions, both entropy and vorticity fluctuations experience power-law growth in time, with a growth rate that depends upon mean flow gradients and, in the absence of dissipative effects, is independent of mode number. If the flow accelerates throughout the implosion, oblate modes amplify by a factor (2C) |N 0 |t i , where C is the convergence ratio of the implosion, N 0 is the initial buoyancy frequency and t i is the implosion time scale. If, instead, the implosion consists of a coasting phase followed by stagnation, oblate modes amplify by a factor exp(π|N 0 |t s ), where N 0 is the buoyancy frequency at stagnation and t s is the stagnation time scale. Even under stable conditions, vorticity fluctuations grow due to the conservation of angular momentum as the gas is compressed. For non-monatomic gases, this additional growth due to compression results in weak oscillatory growth under conditions that would otherwise be buoyantly stable; this over-stability is consistent with the conservation of wave action in the fluid frame. The above analytical results are verified by evolving the complete set of linear equations as an initial value problem, and it is demonstrated that oblate modes are the fastest-growing modes and that high mode numbers are required to reach this limit (Legendre mode ℓ ≳ 100 for spherical flows). Finally, comparisons are made with a Lagrangian hydrodynamics code, and it is found that a numerical resolution of ~30 zones per wavelength is required to capture these solutions accurately. This translates to an angular resolution of ~(12/ℓ)°, or ≲ 0.1° to resolve the fastest-growing modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Generation of Large-Volume High-Pressure Plasma by Spatio-Temporal Control of Space Charge

Due to the fundamental limitations of scaling, very little progress has been made towards achieving a large-volume dense non-equilibrium plasma near atmospheric pressures. Commercially available state-of-the art glow-like atmospheric plasma devices for industrial applications are either narrow tubes or wide slits with narrow openings. Often multiple sources are put together in various configurations to process larger surfaces. The traditional methods of exciting electrodes create spatially fixed electric fields. As a result, the space charge at atmospheric pressure tends to be spatially confined resulting in non-uniformity which eventually leads to instability as the discharge is scaled. Theoretical work done under this project showed that it is possible to generate a spatially rotating electric field by exciting a set of electrodes with phase staggered sinusoidal waveforms. The modeling and simulations were done using plasma fluid models. It was shown that such a field can produce a uniform plasma. At the conclusion of the project, experimental proof of concept with an eight-electrode system in various gases (Air, Helium and Argon) was demonstrated. Power measurements and spectral investigation show that the concept can be used to generate a uniform stable plasma. This plasma source has the potential of opening new applications of nonthermal plasma including combustion of carbon-free fuel. The current limitation of the proposed method in scaling to higher volume and pressure is the need for multiple high voltage amplifiers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Vorticity screening by dense layers

We investigate variable density layers in the presence of vortices, by analogy to line charges in the presence of a dielectric slab. We adapt solutions for dielectric layers from the literature to the variable density fluid case, obtaining exact hypergeometric expressions for flow velocities induced everywhere by a vortex in the vicinity of a finite-thickness slab of different densities. The dimensionless Atwood number, which appears naturally in other variable density phenomena, such as the Rayleigh–Taylor instability, reappears here as the natural ordering parameter for the influence of the slab, which acts primarily to screen the influence of vorticity across the slab to a lower effective circulation. The solution takes the form of scaling correction factors that might be applied to computational models or evaluated directly. Extensions, where vorticity is localized within the dense layer itself and amplified outside of it, are considered, as well as some speculative applications such as inertial confinement fusion capsule designs.

42 ENGINEERING↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Hydrodynamic instability at impact interfaces and planetary implications

Abstract Impact-induced mixing between bolide and target is fundamental to the geochemical evolution of a growing planet, yet aside from local mixing due to jetting – associated with large angles of incidence between impacting surfaces – mixing during planetary impacts is poorly understood. Here we describe a dynamic instability of the surface between impacting materials, showing that a region of mixing grows between two media having even minimal initial topography. This additional cause of impact-induced mixing is related to Richtmyer-Meshkov instability (RMI), and results from pressure perturbations amplified by shock-wave refraction through the corrugated interface between impactor and target. However, unlike RMI, this new impact-induced instability appears even if the bodies are made of the same material. Hydrocode simulations illustrate the growth of this mixing zone for planetary impacts, and predict results suitable for experimental validation in the laboratory. This form of impact mixing may be relevant to the formation of stony-iron and other meteorites.

58 GEOSCIENCES↗
DOE OSTI2022

Nanoconfinement matters in humidified CO 2 interaction with metal silicates

With enigmatic observations of enhanced reactivity of wet CO 2 -rich fluids with metal silicates, the mechanistic understanding of molecular processes governing carbonation proves critical in designing secure geological carbon sequestration and economical carbonated concrete technologies. Here, we use the first principle and classical molecular simulations to probe the impact of nanoconfinement on physicochemical processes at the rock–water–CO 2 interface. We choose nanoporous calcium–silicate–hydrate (C–S–H) and forsterite (Mg 2 SiO 4 ) as model metal silicate surfaces that are of significance in cement chemistry and geochemistry communities, respectively. We show that while a nanometer-thick interfacial water film persists at unsaturated conditions consistent with in situ infrared spectroscopy, the phase behavior of the water–CO 2 mixture changes from its bulk counterpart depending on the surface chemistry and nanoconfinement. We also observe enhanced solubility at the interface of water and CO 2 phases, which could amplify the CO 2 speciation rate. Additionally, through free energy calculations, we show that CO 2 could be found in a metastable state near the C–S–H surface, which can potentially react with surface water and hydroxyl groups to form carbonic acid and bicarbonate. These findings support the explicit consideration of nanoconfinement effects in reactive and non-reactive pore-scale processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

On the acoustic signature of tandem airfoils: The sound of an elastic airfoil in the wake of a vortex generator

The acoustic signature of an acoustically compact tandem airfoil setup in uniform high-Reynolds number flow is investigated. The upstream airfoil is considered rigid and is actuated at its leading edge with small-amplitude harmonic pitching motion. The downstream airfoil is taken passive and elastic, with its motion forced by the vortex-street excitation of the upstream airfoil. The non-linear near-field description is obtained via potential thin-airfoil theory. It is then applied as a source term into the Powell-Howe acoustic analogy to yield the far-field dipole radiation of the system. To assess the effect of downstream-airfoil elasticity, results are compared with counterpart calculations for a non-elastic setup, where the downstream airfoil is rigid and stationary. Depending on the separation distance between airfoils, airfoil-motion and airfoil-wake dynamics shift between in-phase (synchronized) and counter-phase behaviors. Consequently, downstream airfoil elasticity may act to amplify or suppress sound through the direct contribution of elastic-airfoil motion to the total signal. Resonance-type motion of the elastic airfoil is found when the upstream airfoil is actuated at the least stable eigenfrequency of the downstream structure. This, again, results in system sound amplification or suppression, depending on the separation distance between airfoils. With increasing actuation frequency, the acoustic signal becomes dominated by the direct contribution of the upstream airfoil motion, whereas the relative contribution of the elastic airfoil to the total signature turns negligible.

42 ENGINEERING↗
DOE OSTI2021

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Coupled geomechanical investigation of depletion-induced fault reactivation

Fault reactivation during subsurface fluid production pose significant challenges to safe and sustainable resource extraction. Here, this study presents a three-dimensional coupled geomechanical framework to investigate the processes driving fault reactivation, capturing the interactions between reservoir dynamics and geomechanical responses. Verification against theoretical estimations based on linear poroelasticity confirms the model's capacity in representing reservoir background stress responses. However, the study reveals that relying solely on background stress states can underestimate or overestimate fault reactivation potential, emphasizing the importance of including localized stress perturbations such as differential compaction and stress redistribution. Applied to a fault (M1) inspired by the geological characteristics of the Groningen field, the model shows slip initiation at 2965 m depth with 16.0 MPa depletion, aligning with field observations where seismicity occurred at approximately 3 km depth after 15.8 MPa depletion. Parametric studies reveal: (1) inelastic reservoir compaction delays fault reactivation and mitigates fault slip by reducing stress concentration, (2) higher intermediate in-situ stress magnitudes decrease the Coulomb Failure Stress (CFS) increase rate and reduce fault slip, (3) larger fault offsets amplify shear stress near the offset zone, promoting earlier reactivation and longer rupture propagation, and (4) fault permeability significantly influences pressure diffusion, with low-permeability faults leading to sharper stress changes and earlier fault destabilization. These insights highlight the critical role of geological and mechanical parameters in fault reactivation and provide a predictive framework for mitigating induced seismicity risks.

Coupled geomechanical modeling↗
DOE OSTI2014

Simulating Blade-Strike on Fish passing through Marine Hydrokinetic Turbines

The study reported here evaluated the occurrence, frequency, and intensity of blade strike of fish on an axial-flow marine hydrokinetic turbine by using two modeling approaches: a conventional kinematic formulation and a proposed Lagrangian particle- based scheme. The kinematic model included simplifying assumptions of fish trajectories such as distribution and velocity. The proposed method overcame the need for such simplifications by integrating the following components into a computational fluid dynamics (CFD) model: (i) advanced eddy-resolving flow simulation, (ii) generation of ambient turbulence based on field data, (iii) moving turbine blades in highly transient flows, and (iv) Lagrangian particles to mimic the potential fish pathways. The test conditions to evaluate the blade-strike probability and fish survival rate were: (i) the turbulent environment, (ii) the fish size, and (iii) the approaching flow velocity. The proposed method offered the ability to produce potential fish trajectories and their interaction with the rotating turbine. Depending upon the scenario, the percentile of particles that registered a collision event ranged from 6% to 19% of the released sample size. Next, by using a set of experimental correlations of the exposure-response of living fish colliding with moving blades, the simulated collision data were used as input variables to estimate the survival rate of fish passing through the operating turbine. The resulting survival rates were greater than 96% in all scenarios, which is comparable to or better than known survival rates for conventional hydropower turbines. The figures of strike probability and mortality rate were amplified by the kinematic model. The proposed method offered the advantage of expanding the evaluation of other mechanisms of stress and injury on fish derived from hydrokinetic turbines and related devices.

CFD↗
DOE OSTI2017

Impact of numerical choices on water conservation in the E3SM Atmosphere Model Version 1 (EAM V1)

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations for sea level rise projection. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1° horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model is negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in the new model results in a very thin model layer at the Earth’s surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for this model.

Zhang, Kai↗
DOE OSTI2018

Impact of numerical choices on water conservation in the E3SM Atmosphere Model version 1 (EAMv1)

Abstract. The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1∘ horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model becomes negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors in early V1 versions decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in V1 results in a very thin model layer at the Earth's surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for V1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

P-wave velocity anisotropy related to sealed fractures reactivation tracing the structural diagenesis in carbonates

Fracture properties are important in carbonate reservoir characterization, as they are responsible for a large part of the fluid transfer properties at all scales. It is especially true in tight rocks where the matrix transfer properties only slightly contribute to the fluid flow. Open fractures are known to strongly affect seismic velocities, amplitudes and anisotropy. Here, we explore the impact of fracture evolution on the geophysical signature and directional V p anisotropy of fractured carbonates through diagenesis. For that purpose, we studied a meter-scale, parallelepiped quarry block of limestone using a detailed structural and diagenetic characterization, and numerous V p measurements. The block is affected by two en-échelon fracture clusters, both being formed in opening mode (mode 1) and cemented, but only one being reactivated in shear. We compared the diagenetic evolution of the fractures, which are almost all 100% filled with successive calcite cements, with the P-wave velocities measured across this meter-scale block of carbonate, which recorded the tectonic and diagenetic changes of a South Provence sedimentary basin. We found that a directional V p anisotropy magnitude as high as 8–16% correlates with the reactivated fractures’ cluster dip angle, which is explained by the complex filling sequence and softer material present inside the fractures that have been reactivated during the basin's tectonic inversion. We show that although a late karstification phase preferentially affected these reactivated fractures, it only amplified the pre-existing anisotropy due to tectonic shear. We conclude that V p anisotropy measurements may help to identify the fracture sealing/opening processes associated with polyphased tectonic history, the anisotropy being independent of the current stress-state. Finally, this case shows that velocity anisotropies induced by fractures resulted here from a cause that is different from how these features have often been interpreted: selective reactivation of sealed fractures clusters rather than direction of currently open ones.

58 GEOSCIENCES↗
DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗