An official website of the United States government
We propose a novel framework for simulating ink as a particle-laden flow using particle flow maps. Our method addresses the limitations of existing flow-map techniques, which struggle with dissipative forces like viscosity and drag, thereby extending the application scope from solving the Euler equations to solving the Navier-Stokes equations with accurate viscosity and laden-particle treatment. Our key contribution lies in a coupling mechanism for two particle systems, coupling physical sediment particles and virtual flow-map particles on a background grid by solving a Poisson system. We implemented a novel path integral formula to incorporate viscosity and drag forces into the particle flow map process. Our approach enables state-of-the-art simulation of various particle-laden flow phenomena, exemplified by the bulging and breakup of suspension drop tails, torus formation, torus disintegration, and the coalescence of sedimenting drops. In particular, our method delivered high-fidelity ink diffusion simulations by accurately capturing vortex bulbs, viscous tails, fractal branching, and hierarchical structures.
more »
« less
Virtual Experiments of Particle Mixing Process with the SPH-DEM Model
Particle mixing process is critical for the design and quality control of concrete and composite production. This paper develops an algorithm to simulate the high-shear mixing process of a granular flow containing a high proportion of solid particles mixed in a liquid. DEM is employed to simulate solid particle interactions; whereas SPH is implemented to simulate the liquid particles. The two-way coupling force between SPH and DEM particles is used to evaluate the solid-liquid interaction of a multi-phase flow. Using Darcy’s Law, this paper evaluates the coupling force as a function of local mixture porosity. After the model is verified by two benchmark case studies, i.e., a solid particle moving in a liquid and fluid flowing through a porous medium, this method is applied to a high shear mixing problem of two types of solid particles mixed in a viscous liquid by a four-bladed mixer. A homogeneity metric is introduced to characterize the mixing quality of the particulate mixture. The virtual experiments with the present algorithm show that adding more liquid or increasing liquid viscosity slows down the mixing process for a high solid load mix. Although the solid particles can be mixed well eventually, the liquid distribution is not homogeneous, especially when the viscosity of liquid is low. The present SPH-DEM model is versatile and suitable for virtual experiments of particle mixing process with different blades, solid particle densities and sizes, and liquid binders, and thus can expedite the design and development of concrete materials and particulate composites.
more »
« less
- PAR ID:
- 10276150
- Date Published:
- Journal Name:
- Materials
- Volume:
- 14
- Issue:
- 9
- ISSN:
- 1996-1944
- Page Range / eLocation ID:
- 2199
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract. Secondary organic aerosols (SOAs) can exist in liquid, semi-solid, or amorphous solid states. Chemical transport models (CTMs), however, usually assume that SOA particles are homogeneous and well-mixed liquids, with rapid establishment of gas–particle equilibrium for simulations of SOA formation and partitioning. Missing the information of SOA phase state and viscosity in CTMs impedes accurate representation of SOA formation and evolution, affecting the predictions of aerosol effects on air quality and climate. We have previously developed a parameterization to estimate the glass transition temperature (Tg) of an organic compound based on volatility and to predict viscosity of SOA. In this study, we apply this method to predict the phase state of SOA particles over China in summer of 2018 using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). The simulated Tg of dry SOA (Tg,org) agrees well with the value estimated from ambient volatility measurements at an urban site in Beijing. For the spatial distributions of Tg,org, simulations show that at the surface the values of Tg,org range from ∼287 to 305 K, with higher values in northwestern China, where SOA particles have larger mass fractions of low-volatility compounds. Considering water uptake by SOA particles, the SOA viscosity shows a prominent geospatial gradient in which highly viscous or solid SOA particles are mainly predicted in northwestern China. The lowest and highest SOA viscosity values both occur over the Qinghai–Tibet Plateau, where the solid phase state is predicted over dry and high-altitude areas and the liquid phase state is predicted mainly in the south of the plateau with high relative humidity during the summer monsoon season. Sensitivity simulations show that, including the formation of extremely low-volatility organic compounds, the percent time that a SOA particle is in the liquid phase state decreases by up to 12 % in southeastern China during the simulated period. With an assumption that the organic and inorganic compounds are internally mixed in one phase, we show that the water absorbed by inorganic species can significantly lower the simulated viscosity over southeastern China. This indicates that constraining the uncertainties in simulated SOA volatility distributions and the mixing state of the organic and inorganic compounds would improve prediction of viscosity in multicomponent particles in southeastern China. We also calculate the characteristic mixing timescale of organic molecules in 200 m SOA particles to evaluate kinetic limitations in SOA partitioning. Calculations show that during the simulated period the percent time of the mixing timescale longer than 1 h is >70 % at the surface and at 500 hPa in most areas of northern China, indicating that kinetic partitioning considering the bulk diffusion in viscous particles may be required for more accurate prediction of SOA mass concentrations and size distributions over these areas.more » « less
-
null (Ed.)Particle size and settling speed determine the rate of particulate mass transfer from the ocean surface to the sea bed. Turbulent shear in the ocean can act on large, faster-settling flocculated particles to break them into slower-settling primary particles or sub-aggregates. However, it is difficult to understand the disruption behavior of aggregates and their response to varying shear forces due to the complex ocean environment. A study was conducted to simulate the disruption behavior of marine aggregates in the mixed layer of the ocean. The breakup process was investigated by aggregating and disrupting flocs of bentonite clay particles in a rotating and oscillating cylindrical tank 10 cm in diameter filled with salt water. This laboratory tank, which operated based on an extension of Stokes’ second problem inside a cylinder, created laminar oscillating flow superimposed on a constant rotation. This motion allowed the bentonite particles to aggregate near the center of the tank but also exposed large aggregates to high shear forces near the wall. A high-speed camera system was used, along with particle tracking measurements and image processing techniques, to capture the breakup of the large particle aggregates and locate their radial position. The breakup response of large aggregates and the sizes of their daughter particles after breakup were quantified using the facility. The disruption strength of the aggregated particles is presented and discussed relative to their exposure to varying amounts of laminar shear.more » « less
-
Vat photopolymerization (VP) additive manufacturing involves selectively curing low‐viscosity photopolymers via exposure to ultraviolet light in a layer‐wise fashion. Dispersing filler materials in the photopolymer enables tailored end‐use properties, but also increases the viscosity and the timescale associated with interparticle network structural recovery postshear. These rheological properties influence self‐leveling and recoating of the liquid photopolymer mixture during VP. Herein, viscosity of photopolymer and rigid spherical glass microparticles (filler) is experimentally determined as a function of filler fraction, filler size distribution (mono‐ and polydisperse), shear rate, and temperature, which are important VP process parameters. Employing existing viscosity models for mono‐ and polydisperse polymer mixtures demonstrates that particle–particle interactions and the formation of nonspherical clusters of particles strongly affect the viscosity of both monodisperse and polydisperse mixtures with particle volume fractions > 0.05 due to agglomeration/deagglomeration of clusters at elevated shear rates. Consequently, unmodified viscosity models, which assume uniformly dispersed, rigid, spherical particles, are applicable only for mixtures with particle volume fractions < 0.05. It is shown that modifying model parameters such as the fluidity limit and intrinsic viscosity, which explicitly account for nonspherical clusters of particles, improves agreement between viscosity models and experiments, in particular when using a fractal approach.more » « less
-
A Eulerian—Lagrangian model has been developed to simulate particle attachment to surfaces with arc-shaped ribs in a two-dimensional channel flow at low Reynolds numbers. Numerical simulation has been performed to improve the quantitative understanding of how rib geometries enhance shear rates and particle-surface interact for various particle sizes and flow velocities. The enhanced shear rate is attributed to the wavy flows that develop over the ribbed surface and the weak vortices that form between adjacent ribs. Varying pitch-to-height ratio can alter the amplitude of the wavy flow and the angle of attack of the fluid on the ribs. In the presence of these two competing factors, the rib geometry with a pitch-to-height ratio of two demonstrates the greatest shear rate and the lowest fraction of particle attachment. However, the ribbed surfaces have negligible effects on small particles at low velocities. A force analysis identifies a threshold shear rate to reduce particle attachment. The simulated particle distributions over the ribbed surfaces are highly non-uniform for larger particles at higher velocities. The understanding of the effect of surface topography on particle attachment will benefit the design of surface textures for mitigating particulate fouling in a wide range of applications.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3390/ma14092199
