Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow velocity”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Flowing plasma rearrangement in the presence of static perturbing fields

Charged particles interacting with electromagnetic waves have a portion of their energy tied up in wave-driven oscillations. When these waves are localized to the exhaust of linear magnetic confinement systems, this ponderomotive effect can be utilized to enhance particle confinement. The same effect can be derived for particles moving via an E×B drift into a region of a static perturbation to the electromagnetic fields which has a large wave vector component in the direction of the motion. In this work, we use a simplified slab model to self-consistently solve for the electromagnetic fields within the fluid flowing plasma of a static flute-like (k∥=0) perturbation and evaluate the resulting ponderomotive potential. We find that two types of perturbations can exist within the flowing plasma, which are an O wave and an X wave in the frame moving with the fluid. In the case of tenuous plasma, these perturbations are magnetostatic or electrostatic multipole-analog perpendicular to the guiding magnetic field in the lab frame, respectfully. For denser plasmas, the O wave-like perturbation is screened at the electron skin depth scale, and the X wave-like perturbation is a combination of a similar perpendicular electric perturbation and parallel magnetic perturbation. The ponderomotive potential generated in the X wave-like case is gyrofrequency-dependent and can be used as either potential barriers or potential wells, depending on the direction of the flow velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Experimental investigation of laser ablation of stone polycrystalline targets

We report the results of an experimental investigation of ablation of stone polycrystalline targets of complex multicomponent composition, which imitate the substance of asteroids. The targets were irradiated by nanosecond pulses of a neodymium laser at an energy density Φ{sub L} of up to 5 × 10{sup 4} J cm{sup −2}. The experiments demonstrated the existence of two ablation regimes, with the boundary between them lying at Φ{sub L} ≈ 4000 J cm{sup −2}. The regime change is characterised by a change in the form of the dependence of the surface mass density of removed target material on the laser energy density and by the appearance of a minimum in the dependence of specific energy of destruction on Φ{sub L}. This is supposedly related to the passage from a one-dimensional plasma plume expansion to the three-dimensional one and the corresponding decrease in the efficiency of energy transfer from the laser beam to the target due to a lowering of laser-produced plasma density. Our experiments also showed the existence of a maximum in impulse coupling coefficient C{sub m} as a function of laser energy density (C{sub m} ≈ 6.3 × 10{sup −5} N W{sup −1} for Φ{sub L} = Φ{sub opt} ≈ 25 J cm{sup −2}). Maxima were also recorded in the dependences of the ablation efficiency and average ablation flow velocity on Φ{sub L}. For Φ{sub L} > Φ{sub opt}, the decrease in the function C{sub m}(Φ{sub L}) turns out to be much steeper than for metals and polymer materials. The difference is supposedly due to the lower strength and lower plasticity of the polycrystalline stone targets. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Time-resolved turbulent dynamo in a laser plasma

Understanding magnetic-field generation and amplification in turbulent plasma is essential to account for observations of magnetic fields in the universe. A theoretical framework attributing the origin and sustainment of these fields to the so-called fluctuation dynamo was recently validated by experiments on laser facilities in low-magnetic-Prandtl-number plasmas (Pm < 1). However, the same framework proposes that the fluctuation dynamo should operate differently when Pm ≳ 1, the regime relevant to many astrophysical environments such as the intracluster medium of galaxy clusters. This paper reports an experiment that creates a laboratory Pm ≳ 1 plasma dynamo. We provide a time-resolved characterization of the plasma’s evolution, measuring temperatures, densities, flow velocities, and magnetic fields, which allows us to explore various stages of the fluctuation dynamo’s operation on seed magnetic fields generated by the action of the Biermann-battery mechanism during the initial drive-laser target interaction. The magnetic energy in structures with characteristic scales close to the driving scale of the stochastic motions is found to increase by almost three orders of magnitude and saturate dynamically. It is shown that the initial growth of these fields occurs at a much greater rate than the turnover rate of the driving-scale stochastic motions. Here, our results point to the possibility that plasma turbulence produced by strong shear can generate fields more efficiently at the driving scale than anticipated by idealized magnetohydrodynamics (MHD) simulations of the nonhelical fluctuation dynamo; this finding could help explain the large-scale fields inferred from observations of astrophysical systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Doppler Backscattering Data Analysis and Integrated Modeling with OMFIT

One Modeling Framework for Integrated Tasks (OMFIT) is a widely used software tool in the magnetic fusion research community. OMFIT provides magnetic fusion energy researchers with a framework for the development of special-purpose physics modules. This paper describes an OMFIT physics module pertaining to the Doppler Backscattering (DBS) fusion plasma diagnostic. DBS measures density fluctuations and flow velocity through plasma scattering of electromagnetic waves. The OMFIT DBS module was developed to analyze experimental DBS data and facilitate modeling of DBS systems installed on multiple tokamak devices. The OMFIT DBS module is designed to support several analysis workflows: detailed analysis of experimental data, experimental planning, and theory-based synthetic diagnostic modeling. The DBS module uses integrated modeling by leveraging other OMFIT physics modules to perform tasks related to DBS, e.g. ray/beam–tracing simulations, edge-localized mode–synchronized data analysis, magnetic equilibrium reconstruction, and fitting kinetic profile data. Furthermore, this paper describes several supported workflows and serves a reference for the OMFIT DBS module.

Doppler backscattering↗
DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗
DOE OSTI2021

Particle image velocimetry analysis with simultaneous uncertainty quantification using Bayesian neural networks

Particle image velocimetry (PIV) is an effective tool in experimental fluid mechanics for extracting flow fields from images. Recently, convolutional neural networks (CNNs) have been used to perform PIV analysis with accuracy on par with classical methods. Here we extend the use of CNNs to analyze PIV data while providing simultaneous uncertainty quantification on the inferred flow field. The method we apply in this paper is a Bayesian convolutional neural network (BCNN) which learns distributions of the CNN weights through variational Bayes. In order to demonstrate the utility of BCNNs for the PIV task, we compare the performance of three distinct BCNN models with simple architectures. The first network estimates flow velocity from image interrogation regions only. Our second model learns to infer velocity from both the image interrogation regions and interrogation region cross-correlation maps. Finally, our best performing network infers velocities from interrogation region cross-correlation maps only. We find that BCNNs using interrogation region cross-correlation maps as inputs perform better than those using interrogation windows only as inputs and discuss reasons why this may be the case. Additionally, we test the best performing BCNN on a full synthetic test image pair and a real image pair from the 1st International PIV Challenge. We show that ~98% of true particle displacements from the full synthetic image pair can be captured within the BCNN's 95% confidence intervals, and that the BCNN's performance on the real image pair is quantitatively similar to that of algorithms tested in the 1st International PIV Challenge. Finally, we show that BCNNs can be generalized to be used with multi-pass PIV algorithms with a moderate loss in accuracy, which may be overcome by future work on finetuning and training schemes. So to our knowledge, this is the first use of Bayesian neural networks to perform PIV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Solitary zonal structures in subcritical drift waves: a minimum model

Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa–Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa–Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same 'equation of state', which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analysing the effects of heating and gas puffing in Proto-MPEX helicon and auxiliary heated plasmas

The Material Plasma Exposure Experiment (MPEX) is being constructed at Oak Ridge National Laboratory to investigate critical fusion reactor issues, such as plasma–material interactions (PMIs) under reactor-relevant conditions and time scales. The linear device Proto-MPEX was used as a test bed to address anticipated research and development issues associated with heating scenarios and establish the physics basis for MPEX. The SOLPS-ITER code suite has been applied to understand plasma and neutral transport in Proto-MPEX and to increase confidence in predictive simulations for MPEX. Coupling between COMSOL and SOLPS is performed to implement a 2D electron heating profile of the helicon source. The simulations show reasonable agreement with the experimental data for plasma with helicon and auxiliary electron cyclotron heating (ECH). Both Bohm and constant diffusion (D$_\perp$: 0.5 m 2 s –1 and $\chi_\perp$: 1 m 2 s –1 ) simulations show similar levels of agreement with respect to the sparse experimental data available, assuming a few per cent impurity concentration. ECH significantly increases the target electron temperature, however, the target electron density is reduced compared to helicon-only heated plasmas due to an increase in flow velocity and radial losses. The simulations show that further increasing the ECH power results in an increase in the target electron density due to increased recycling flux and ionization. The results indicate that ECH significantly enhances the target heat flux, with ECH power of 50 kW increasing the target heat flux from 0.4 to 17 MW m –2 . It is found that a small amount of gas puffing (GP) near the target plate can further increase the target heat fluxes at the higher ECH power cases, but the target heat flux is reduced at higher GP conditions due to a significant reduction in electron temperature via radiation. ECH and GP scenarios can generate a higher target flux, facilitating improved PMI studies with more reactor-relevant plasma conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 &#x00d7; 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E &#x00d7; B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E &#x00d7; B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The impact of the $\text{WHIM}$ on the $\text{IGM}$ thermal state determined from the low- z Lyman $\alpha$ forest

At z ≲ 1, shock heating caused by large-scale velocity flows and possibly violent feedback from galaxy formation, converts a significant fraction of the cool gas (T ~ 10 4 K) in the intergalactic medium (IGM) into warm–hot phase (WHIM) with T > 10 5 K, resulting in a significant deviation from the previously tight power-law IGM temperature–density relationship, T = T 0 (P/$\bar{p}$) γ-1 ⁠. This study explores the impact of the WHIM on measurements of the low-z IGM thermal state, [T 0 , γ], based on the b–NH 1 distribution of the Ly α forest. Exploiting a machine learning-enabled simulation-based inference method trained on Nyx hydrodynamical simulations, we demonstrate that [T 0 , γ] can still be reliably measured from the b–NH 1 distribution at z = 0.1, notwithstanding the substantial WHIM in the IGM. To investigate the effects of different feedback, we apply this inference methodology to mock spectra derived from the IllustrisTNG and Illustris simulations at z = 0.1. The results suggest that the underlying [T 0 , γ] of both simulations can be recovered with biases as low as |Δlog(T 0 /K)| ≲ 0.05 dex, |Δγ| ≲ 0.1, smaller than the precision of a typical measurement. Given the large differences in the volume-weighted WHIM fractions between the three simulations (Illustris 38 percent, IllustrisTNG 10 percent, and Nyx 4 per cent), we conclude that the b–N H1 distribution is not sensitive to the WHIM under realistic conditions. Finally, we investigate the physical properties of the detectable Ly α absorbers, and discover that although their T and Δ distributions remain mostly unaffected by feedback, they are correlated with the photoionization rate used in the simulation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Revealing the Structure and Dynamics of Self-Generated Electric and Magnetic Fields Near Plasma Stagnation in Laser-Driven Hohlraums

By coupling newly developed triparticle charged particle radiography with radiography reconstruction algorithms and novel reconstruction postprocessing techniques, the spatial structure and time evolution of self-generated electric and magnetic fields in laser-driven vacuum hohlraums have been quantitatively revealed. Through high-fidelity data from a series of experiments, it is shown that late in the hohlraum evolution (after the end of laser drive) these fields are strongly correlated in both space and time, providing evidence that their evolution is primarily dominated by advection with the plasma flow. At these late times, plasma flow velocities inferred from both gross radiography analysis and field reconstructions (and corroborated with Thomson scattering measurements) indicate that the plasma is approaching stagnation near the hohlraum axis. Finally, these experiments provide not only new physical insight into spontaneously generated hohlraum fields, but also provide important spatially and temporally resolved information for future benchmarking of numerical codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Pump with Broadband Probe Experiments for Single-Shot Measurements of Plasma Conditions and Crossed-Beam Energy Transfer

Here, a novel technique for measuring plasma conditions using monochromatic pump-broadband-probe laser interactions has been experimentally demonstrated. Originally proposed in Ludwig et al. [Phys. Plasmas 26, 113108 (2019)], this method utilizes crossed-beam energy transfer between the broadband probe and the pump, mediated by plasma ion acoustic waves. The complete energy transfer spectrum can be captured in a single shot, enabling the inference of plasma parameters such as density, electron and ion temperatures, and flow velocity. Compared to Thomson scattering, this technique offers signal enhancements typically larger than 9 orders of magnitude, significantly reducing the required probe laser intensity and facilitating interactions that are linear and measurements that are nonperturbative of the plasma. Furthermore, it provides a powerful tool for advancing studies of crossed-beam energy transfer under conditions relevant to inertial confinement fusion experiments.

Physics - Plasma physics↗
DOE OSTI2021

Nonequilibrium kinetic freeze-out properties in relativistic heavy ion collisions from energies employed at the RHIC beam energy scan to those available at the LHC

In this paper, we investigate the kinetic freeze-out properties in relativistic heavy ion collisions at different collision energies. We present a study of standard Boltzmann-Gibbs blast-wave (BGBW) fits and Tsallis blast-wave (TBW) fits performed on the transverse momentum spectra of identified hadrons produced in Au + Au collisions at collision energies of $\sqrt{s_{\rm{NN}}}=$7.7 –200 GeV at the Relativistic Heavy Ion Collider (RHIC), and in Pb + Pb collisions at collision energies of $\sqrt{s_{\rm{NN}}}=$ 2.76 and 5.02 TeV at the Large Hadron Collider (LHC). The behavior of strange and multistrange particles is also investigated. We found that the TBW model describes data better than the BGBW one overall, and the contrast is more prominent as the collision energy increases as the degree of nonequilibrium of the produced system is found to increase. From TBW fits, the kinetic freeze-out temperature at the same centrality shows a weak dependence of collision energy between 7.7 and 39 GeV, while it decreases as collision energy continues to increase up to 5.02 TeV. The radial flow is found to be consistent with zero in peripheral collisions at RHIC energies but sizable at LHC energies and central collisions at all RHIC energies. We also observed that the strange hadrons, with higher temperature and similar radial flow, approach equilibrium more quickly from peripheral to central collisions than light hadrons. The dependence of temperature and flow velocity on nonequilibrium parameter (q - 1) is characterized by two second-order polynomials. Both a and dξ from the polynomials fit, related to the influence of the system bulk viscosity, increase toward lower RHIC energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Water Energy Resource Innovation on the Cavitation Characteristics

The main purpose of this study is to numerically correlate the amount of generated vapor over a hydrofoil to the lift and drag coefficients acting on it. As such, cavitation characteristics were investigated of a hydrofoil in the cavitating, sub-cavitating, and non-cavitating flows for different angles of attacks (AoA) with the high upstream flow velocity. The hydrofoil was tested in a square water tunnel with water entering the tunnel at various velocities for each AoA ranges from 9.1 m/s to 12.2 m/s. It was found that lift and drag forces acting on the hydrofoil follow the trend of the experimental data quite closely. While the cavitation can be identified by a unique number (averaged vapor volume fraction), the work done created an inverse correlation between this number and the cavitation number at the same angle of attack. The lift force declines with the increase in the vapor content on the hydrofoil surface, meanwhile the drag force peaks at certain vapor volume fraction, and then, a huge reduction occurs with the considerable decrease in the corresponding cavitation number. A fourth-order correlation generated between the lift to drag (L/D) and the cavitation number (σ). It was found the lift-to-drag ratio decreases by the formation of the cavitation over the hydrofoil, thus causing a drop in the efficiency of the turbomachines.

42 ENGINEERING↗
DOE OSTI2021

The Effect of the Degree of Premixedness on Self-Excited Combustion Instability

The use of lean, premixed fuel and air mixtures is a common strategy to reduce NO x emissions in gas turbine combustors. However, this strategy causes an increased susceptibility to self-excited instability, which manifests as fluctuations in heat release rate, flow velocity, and combustor acoustics that oscillate in-phase in a feedback loop. Here, this study considers the effect of the level of premixedness on the self-excited instability in a single-nozzle combustor. In this system, the fuel can be injected inside the nozzle to create a partially-premixed mixture or far upstream to create a fully-premixed mixture, varying the level of premixedness of the fuel and air entering the combustor. When global equivalence ratio is held constant, the cases with higher levels of premixing have higher instability amplitudes. High-speed CH* chemiluminescence imaging shows that the flame for these cases is more compact and the distribution of the heat release rate oscillations is more concentrated near the corner of the combustor in the outer recirculation zones. Rayleigh index images, which are a metric for the relative phase of pressure and heat release rate oscillations, suggest that vortex rollup in the corner region is primarily responsible for determining instability characteristics when premixedness is varied. This finding is further supported through analysis of local flame edge dynamics.

42 ENGINEERING↗
DOE OSTI2023

High-Fidelity Energy Deposition Ignition Model Coupled with Flame Propagation Models at Engine-like Flow Conditions

With the heightened pressure on car manufacturers to increase the efficiency and reduce the carbon emissions of their fleets, more challenging engine operation has become a viable option. Highly dilute, boosted, and stratified charge, among others, promise engine efficiency gains and emissions reductions. At such demanding engine conditions, the spark-ignition process is a key factor for the flame initiation propagation and the combustion event. From a computational standpoint, there exist multiple spark-ignition models that perform well under conventional conditions but are not truly predictive under strenuous engine operation modes, where the underlying physics needs to be expanded. In this paper, a hybrid Lagrangian-Eulerian spark-ignition (LESI) model is coupled with different turbulence models, grid sizes, and combustion models. The ignition model, previously developed, relies on coupling Eulerian energy deposition with a Lagrangian particle evolution of the spark channel, at every time-step. The spark channel is attached to the electrodes and allowed to elongate at a speed derived from the flow velocity. The LESI model is used to simulate spark ignition in a nonquiescent crossflow environment at engine-like conditions, using converge commercial computational fluid dynamics (CFD) solver. The results highlight the consistency, robustness, and versatility of the model in a range of engine-like setups, from typical with Reynolds-averaged Navier-Stokes (RANS) and a larger grid size to high fidelity with large-eddy simulation (LES) and a finer grid size. The flame kernel growth is then evaluated against Schlieren images from an optical constant volume ignition chamber with a focus on the performance of flame propagation models, such as G-equation and thickened flame model, versus the baseline well-stirred reactor model. Finally, future development details are discussed.

Advanced ignition modeling↗
DOE OSTI2023

Optimization of Dimensions of Smooth and Twisted-Tape-Inserted Tubes for Heat Transfer with NaCl/KCl/MgCl 2 Molten Salts by Principle of Entropy Generation Minimization

The entropy generation minimization principle is used as the criterion to optimize the flow and heat transfer of solar collectors and heat exchangers that use molten salts NaCl–KCl–MgCl 2 and KCl–MgCl 2 . The Gnielinski correlation for the Nusselt number versus Reynolds number, as well as the Moody friction factor given by Petukhov, was used for the calculation of the convective heat transfer coefficient and pressure loss due to friction in smooth tubes. For twisted-tap-inserted tube, equations of Nu and friction factor provided by Manglik and Bergles were used. The objective function, the entropy generation rate of the heat transfer system, was expressed as the function of Reynolds number, Prandtl number, heating flux, tube diameter, etc. As a result of the analysis, the optimum Reynolds number was determined and thereby to determine the optimum Nusselt number, convective heat transfer coefficient, friction factor, and tube diameter, which also allows the calculation of optimum flow velocity. The analysis was conducted in the fluid temperature range of 500–700 °C, which covers the operation temperature for supercritical CO 2 power cycles in concentrated solar power (CSP) system. Here, optimized results from the smooth tube and twisted-tap-inserted tube are compared, which is important to the design of solar receivers for CSP systems.

14 SOLAR ENERGY↗
DOE OSTI2023

Granular Flow in Novel Octet Shape–Based Lattice Frame Material

Moving packed-bed heat exchangers in concentrated solar power (CSP) plants involves heat transfer between heated falling particles and supercritical carbon dioxide. The overall effective thermal conductivity of the moving packed bed and particle-side channel contact resistances are still the bottlenecks in achieving the desirable thermal transport levels. To this end, a novel moving packed bed heat exchanger consisting of an Octet lattice packed between the walls of the particle-side channel is proposed in this study. Granular flow analysis in Octet lattice moving packed bed heat exchanger (OLHX) was conducted through experiments and discrete element method (DEM)-based numerical simulations. The experimental images clearly demonstrated stagnation regions upstream of lattice fibers, void regions downstream of the fiber junctions, and wavy-type unobstructed flow on the lateral sides of the fibers. DEM simulations were successful in capturing all these critical flow phenomena. Larger flow velocities were observed on the lateral sides of the fibers in the simulations. Also, when the particles in the silo were emptied, the final images showed an accumulation of particles on the inter-fiber as well as fiber–channel wall junctions. Moreover, the fiber connections resulted in some regions devoid of particle contact on the channel endwall, which means that these regions would suffer from poor thermal exchange. Finally, the overall mass flowrate increased with increasing porosity for a fixed particle diameter.

14 SOLAR ENERGY↗