Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid-particle interactions”

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.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗
DOE OSTI2021

Coupling particle scale model and SuperDEM–CFD for multiscale simulation of biomass pyrolysis in a packed bed pyrolyzer

An efficient biomass pyrolysis process requires a comprehensive understanding of the chemical and physical phenomena that occur at multi-length and time scales. In this study, a multiscale computational approach was developed and validated for biomass pyrolysis in a packed-bed reactor by integrating pyrolysis kinetics, a particle scale model, and SuperDEM-CFD in open-source code MFiX. A one-dimensional particle-scale model that discretizes the characteristic length of biomass particle into layers was developed to predict the intraparticle phenomena inside a single particle. The 1D model was validated by comparing it with a single biomass particle pyrolysis experiment. Here, a recently developed SuperDEM-CFD model was employed to simulate the non-spherical particle-particle contact and fluid-particle interaction. The coupled model was applied to simulate the pyrolysis of cubic biomass particles in a packed bed and validated by comparing with experimental data. Simulation with and without particle-scale model was compared, and the effect of the gas-solid heat transfer models was also investigated.

09 BIOMASS FUELS↗
DOE OSTI2023

Revisiting the empirical particle-fluid coupling model used in DEM-CFD by high-resolution DEM-LBM-IMB simulations: A 2D perspective

The work investigates the applicability of the unresolved Computational Fluid Dynamics and Discrete Element Method (CFDDEM) technique based on empirical equations for fluid-particle coupling. We first carry out a series of representative volume element simulations using the high-resolution particle-resolved Lattice Boltzmann method and Discrete Element Method (LBMDEM) coupled by an Immersed Moving Boundary (IMB) scheme. Then, we compare the results obtained by both LBMDEM and empirical equations used in unresolved CFDDEM with analytical solutions. It is found that the existing empirical equations used in solving fluid-particle interactions in 2D CFDDEM fail to accurately calculate the hydrodynamic force applied to solid particles. The underlying reason is that the existing empirical models are obtained based on 3D experimental results and thus are not applicable to 2D problems. Based on the simulation results, a new drag coefficient model is then proposed. The estimated drag forces using the new model are compared favourably with the simulated ones, indicating the good performance of the proposed model.

42 ENGINEERING↗
DOE OSTI2021

In-situ magnetic alignment model for additive manufacturing of anisotropic bonded magnets

In this work, we report a mathematical framework which predicts the degree of alignment (DoA) in an in-situ aligned additively manufactured 3D printed bonded magnets. A multiphysics model is developed which couples the harmonious interactions of magnetic particles in a viscous flowing polymer under the presence of an external magnetic field. The hydrodynamic fluid-particle interaction is paired with the magnetophoretic force to predict the particle trajectory and distribution during extrusion through a 3D printer nozzle. Succeeding the force balance, a magnetohydrodynamic torque equilibrium analysis is performed to predict the net-orientation of the magnetic particles as a function of the applied field strength, viscous forces, and particle-to-particle interactions (P2P). Experimental validation of the DoA predictions is performed using 65 vol% Nd-Fe-B+Sm-Fe-N in Nylon12 (DoA exp = 0.620 and DoA theory = 0.686) and 15 vol% Sm-Co in PLA (DoA exp = 0.830 and DoA theory = 0.863). A parametric analysis is performed to analyze the effect of operating and design parameters like alignment field strength, magnetic loading fraction, extrusion load, and particle size. The model predicts a competing behavior between particle-fluid and particle-particle interactions under the presence of an applied field. Lastly, the model provides a framework to efficiently predict the DoA in tandem with a functionalized-magnetic 3D printer and allows the user to adjust the operating parameters according to the desired DoA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2021

Evaporation and propagation of liquid drop streams at vacuum pressures: Experiments and modeling

Evaporation of streams of liquid droplets in environments at vacuum pressures below the vapor pressure has not been widely studied. Here, experiments and simulations are reported that quantify the change in droplet diameter when a steady stream of ≈100 μm diameter drops is injected into a chamber initially evacuated to <10 –8 bar. In experiments, droplets fall through the center of a 0.8 m long liquid nitrogen cooled shroud, simulating infinity radiation and vapor mass flux boundary conditions. Furthermore, experimentally measured changes in drop diameters vary from ≈0 to 6 μm when the initial vapor pressure is increased from 10 –6 to 10 –3 bar by heating the liquid. Measured diameter changes are predicted by a model based on the Hertz-Knudsen equation. One uncertainty in the calculation is the “sticking coefficient” β. Assuming a constant β for all conditions studied here, predicted diameter changes best match measurements with β≈0.3. This value falls within the range of β reported in the literature for organic liquids. Finally, at the higher vapor pressure conditions considered here, the drop stream disperses transverse to the main flow direction. This spread is attributed to forces imparted by an absolute pressure gradient produced by the evaporating stream.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗
DOE OSTI2020

Cellulose nanofibrils and nanocrystals in confined flow: Single-particle dynamics to collective alignment revealed through scanning small-angle x-ray scattering and numerical simulations

Nanostructured materials made through flow-assisted assembly of proteinaceous or polymeric nanosized fibrillar building blocks are promising contenders for a family of high-performance biocompatible materials in a wide variety of applications. Optimization of these processes relies on improving our knowledge of the physical mechanisms from nano- to macroscale and especially understanding the alignment of elongated nanoparticles in flows. Here, we study the full projected orientation distributions of cellulose nanocrystals (CNCs) and nanofibrils (CNFs) in confined flow using scanning microbeam SAXS. For CNCs, we further compare with a simulated system of dilute Brownian ellipsoids, which agrees well at dilute concentrations. However, increasing CNC concentration to a semidilute regime results in locally arranged domains called tactoids, which aid in aligning the CNC at low shear rates, but limit alignment at higher rates. Similarly, shear alignment of CNF at semidilute conditions is also limited owing to probable bundle or flock formation of the highly entangled nanofibrils. Finally, this work provides a quantitative comparison of full projected orientation distributions of elongated nanoparticles in confined flow and provides an important stepping stone towards predicting and controlling processes to create nanostructured materials on an industrial scale.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Ring patterns generated by an expanding colloidal meniscus

The drop-and-dry is a common technique allowing for creation of periodic nanoparticle (NP) structures for sensing, photonics, catalysis, etc. However, the reproducibility and scalability of this approach for fabrication of NP-based structures faces serious challenges due to the complexity of the simple, at first glance, evaporation process. In this work we study the effect of the spatial confinement on the NP self-assembly under slow solvent evaporation, when the air-liquid-substrate contact line (CL) expands from the center towards the walls of a cylindrical cell, forming a toroid. Using in situ video monitoring of the stick-slip CL motion, we find regular hydrodynamic perturbations in the meniscus, and reveal fine details of the formation of quasiperiodic rings of close packed NP layers. We report that drying of the toroidal NP droplet has a number of important differences from drying of the classical hemispherical colloidal drops. In toroidal drops we observe linear-in-time average meniscus motion, in contrast to the hemispherical drops where the meniscus moves as a square root of time. While both droplet geometries produce NP ring patterns, the ring width for the toroidal drop decreases with increasing ring radius, while it decreases with decreasing the radius of the hemispherical drop. We suggest that free ligands are the main cause of the Marangoni instabilities driving the periodic vorticity in the meniscus. Finally, we show that the usually ignored contact line tension may yield a considerable contribution to the CL pinning causing the CL slip-stick motion and the ring formation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerical calculation of the particle–fluid–particle stress in random arrays of fixed particles

Based on the nearest particle statistics [Zhang, J. Fluid Mech. 910, A16 (2021)], the phase interaction force in a multiphase flow is decomposed into a particle–mean-field force and the divergence of the particle–fluid–particle (PFP) stress. The PFP stress is proportional to the correlation product of the distance from a particle to its nearest neighbor and the force on the particle conditionally averaged on the nearest-neighbor location. In this work, a functional form of the stress is obtained corrected to the first order of the ratio between the interparticle distance and the macroscopic length scale of the flow. Particle-resolved numerical simulations are used to calculate the PFP stress in random arrays of fixed particles and to explore the physics represented by the stress. The numerical results show that the PFP stress is attractive along the direction of the flow and is repulsive in the directions perpendicular to the flow. In the flow regime simulated, this PFP stress can be considered as a macroscopic representation of the drafting–kissing–tumbling mechanism. Finally, the Reynolds stress for the fluid phase is also calculated and compared with the PFP stress.

42 ENGINEERING↗
DOE OSTI2022

Shock interacting with a random array of stationary particles underwater

Accurately predicting the force experienced by particles within a random bed during shock traversal is a challenging problem. In this paper, particle-resolved inviscid simulations of shock propagation over randomly distributed beds are performed, with the goal of quantifying the force on individual particles. Water was considered the continuum material, and then results were compared with previous research completed in air and the differences resulting from these different mediums were examined. Simulations were conducted for four different combinations of incident Mach number and bed volume fraction with the particles remaining stationary. Additionally, time-resolved streamwise and transverse force coefficients were calculated for each particle in the bed. It was observed that while the average force coefficient was similar to that of an isolated particle, the force on individual particles was substantially different. An important observation was that the inviscid drag force on some particles remained consistently positive, but others consistently negative even long after the passage of the shock. This persistent force on the particle was attributed to the potential flow fields that result from neighbors' specific locations upstream and downstream. Similar behavior was observed in the transverse force, with particles experiencing sustained force in the same lateral direction. These persistent forces are inviscid quasisteady contributions which with the viscous counterpart will play an essential role in the particle's long-term dispersion if allowed to move in response to the force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Stokesian dynamics with odd viscosity

Stokesian dynamics is a well-established computational method for simulating dynamics of many particles suspended in a conventional passive fluid medium. Active fluids composed of self-driving particles with broken time-reversal symmetry permit the emergence of a so-called odd viscosity. Here, in this work, we extended the conventional Stokesian dynamics formalism to incorporate the additional hydrodynamic effects due to odd viscosity, which enables simulating collective behavior of many particles suspended in an active fluid medium with both even viscosity and odd viscosity.

42 ENGINEERING↗
DOE OSTI2022

Experimental Demonstration of the Thin-Film Liquid-Metal Jet as a Charge Stripper

For high-power heavy ion accelerators, the development of a suitable charge stripper, which can handle intense beams, is essential. This Letter describes the first experimental demonstration of a heavy ion liquid lithium charge stripper. A 10–20 μm thick liquid lithium jet flowing at >50 m/s was formed and confirmed stable when bombarded by various heavy ion beams, while increasing the charge state of the incoming beams to the desired charge state range. Finally, this demonstration proved the existing power limitation with the conventional strippers can be overcome by the liquid-metal stripper, opening completely new possibilities in high-power accelerator development.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Rectified Rotational Dynamics of Mobile Inclusions in Two-Dimensional Active Nematics

We investigate the dynamics of mobile inclusions embedded in 2D active nematics. The interplay between the inclusion shape, boundary-induced nematic order, and autonomous flows powers the inclusion motion. Disks and achiral gears exhibit unbiased rotational motion, but with distinct dynamics. In comparison, chiral gear-shaped inclusions exhibit long-term rectified rotation, which is correlated with dynamics and polarization of nearby +1/2 topological defects. Furthermore, the chirality of defect polarities and the active nematic texture around the inclusion correlate with the inclusion’s instantaneous rotation rate. Inclusions provide a promising tool for probing the rheological properties of active nematics and extracting ordered motion from their inherently chaotic motion.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bayesian estimations of orientation distribution functions from small-angle scattering enable direct prediction of mechanical stress in anisotropic materials

Properties of soft materials are influenced by their anisotropic structuring under nonequilibrium fields. Although anisotropic structure-property relationships have been extensively explored theoretically, comparison to experiments requires determination of the microstructural orientation probability distribution function (OPDF) of microstructural elements. Small angle scattering (SAS) measurements encode information about the OPDF, but tools to navigate this connection are incomplete. Here, we develop and validate an explicit framework to link arbitrary OPDFs to SAS measurements. Specifically, we propose, validate, and apply a method, maximum a posteriori scattering inference (MAPSI), whereby the OPDF may be obtained from SAS measurements using a Bayesian estimation method. Using this method, we obtain estimates of the full 3D OPDF for two model semidilute fd-virus (rodlike) dispersions at concentrations that are approximately equal to and twice the overlap concentration. From the OPDF, we calculate its second and fourth moments and compare these to predictions for a dilute suspension of rigid rods and to a recent theory for semidilute suspensions. Finally, we use both the theoretical and measured moments to calculate the stress, both for dilute and semidilute suspensions. These predictions are not only compared to each other, but also to measured values of the shear stress, and point to new insights into the behavior of suspensions of highly elongated particles in the transition between dilute and semidilute behavior. We also use this new framework to provide perspective on the connection between scalar parameterizations of scattering and the OPDF that have frequently been used in the past. Here, the new tools developed in this work provide an unprecedented path toward experimental validation of dynamical theories of rodlike colloids and polymers, and for measurement of nonequilibrium structures and stresses of other complex fluids and soft materials with SAS.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A coupled polygonal DEM-LBM technique based on an immersed boundary method and energy-conserving contact algorithm

Here we present a framework of coupling polygonal discrete elements and the lattice Boltzmann method using a direct forcing immersed boundary scheme. In this technique, an energy-conserving contact algorithm is utilized to handle the interactions between convex and concave polygonal particles. The surface of a polygon is represented by discrete boundary points which includes vertices of polygonal particles and/or points interpolated from vertices. The fluid-particle coupling is obtained through the interactions of the boundary points and the imaginary fluid particles using a direct-forcing immersed boundary method. Validations of the proposed technique are made by single particle and multiple arbitrarily-shaped particle sedimentation tests, and the effect of particle shape is illustrated using a drafting-kissing-tumbling benchmark.

42 ENGINEERING↗