Non-uniform Particle Distribution Research Articles

Experimental simulation evaluation and control measures of paving segregation for the stone matrix asphalt mixtures

The aggregate segregation for asphalt mixtures, characterized by non-uniform particle distribution, significantly undermines the quality and performance of pavements. This investigation aims to explore and reveal the mechanisms underlying segregation, including raw material properties and construction procedures. Novel evaluation methods, predictive metrics, and control measures were introduced using the developed apparatus. The results indicate that physical properties, particle morphology, and gradation are fundamental factors driving segregation. Spatial hindrance, interfacial friction, collision, and self-organization collectively coordinate particle spatial migration and distribution. Excessive coarse particles induce interstitial displacement, while that of fine particles trigger third-body lubrication. The lubricating-bonding factor η was proposed to characterize asphalt regulation. Segregation prediction based on composite angularity facilitates informed raw material selection. In addition, measures to improve paving quality and efficiency during loading, unloading, and auger mixing were proposed. This research contributes to the understanding of segregation mechanisms, enabling early prevention and effective management strategies.

An Equalized Flow Velocity Strategy for Perovskite Colloidal Particles in Flexible Perovskite Solar Cells.

The non-uniform distribution of colloidal particles in perovskite precursor results in an imbalanced response to the shear force during flexible printing process. Herein, it is observed that the continuous disordered migration occurring in perovskite inks significantly contributes to the enlargement of colloidal particles size and diminishes the crystallization activity of the inks. Therefore, a molecular encapsulation architecture by glycerol monostearate to mitigate colloidal particles collisions in the precursor ink, while simultaneously homogenizing the size distribution of perovskite colloids to minimize their diffusion disparities, is devised. The utilization of colloidal particles with a molecular encapsulation structure enables the achievement of uniform deposition during the printing process, thereby effectively balancing the crystallization rate and phase transition in the film and facilitating homogeneous crystallization of perovskite films. The large-area flexible perovskite device (1.01 cm2 and 100 cm2) fabricated through printing processes, achieves an efficiency of 24.45% and 15.87%, respectively, and manifests superior environmental stability, maintaining an initial efficiency of 91% after being stored in atmospheric ambiences for 150 days (unencapsulated). This work demonstrates that the dynamic evolution process of colloidal particles in both the precursor ink and printing process represents a crucial stride toward achieving uniform crystallization of perovskite films.

Lane formation in gravitationally driven colloid mixtures consisting of up to three different particle sizes.

Brownian dynamics simulations are utilized to study segregation phenomena far from thermodynamic equilibrium. In the present study, we expand upon the analysis of binary colloid mixtures and introduce a third particle species to further our understanding of colloidal systems. Gravitationally driven, spherical colloids immersed in an implicit solvent are confined in two-dimensional linear microchannels. The interaction between the colloids is modeled by the Weeks-Chandler-Andersen potential, and the confinement of the colloids is realized by hard walls based on the solution of the Smoluchowski equationin half space. In binary and ternary colloidal systems, a difference in the driving force is achieved by differing colloid sizes but fixed mass density. We observe for both the binary and ternary systems that a driving force difference induces a nonequilibrium phase transition to lanes. For ternary systems, we study the tendency of lane formation to depend on the diameter of the medium-sized colloids. Here we find a sweet spot for lane formation in ternary systems. Furthermore, we study the interaction of two differently sized colloids at the channel walls. Recently we observed that driven large colloids push smaller colloids to the walls. This results in small particle lanes at the walls at early simulation times. In this work we additionally find that thin lanes are unstable and dissolve over very long time frames. Furthermore, we observe a connection between lane formation and the nonuniform distribution of particles along the channel length. This nonuniform distribution occurs either alongside lane formation or in shared lanes (i.e., lanes consisting of two colloid types).

Nanoimprinted TiO2 Metasurfaces with Reduced Meta-Atom Aspect Ratio and Enhanced Performance for Holographic Imaging.

Metasurface holograms, with the capability to manipulate spatial light amplitudes and phases, are considered next-generation solutions for holographic imaging. However, conventional fabrication approaches for meta-atoms are heavily dependent on electron-beam lithography (EBL), a technique known for its expensive and time-consuming nature. In this paper, a polarization-insensitive metasurface hologram is proposed using a cost-effective and rapid nanoimprinting method with titanium dioxide (TiO2) nanoparticle loaded polymer (NLP). Based on a simulation, it has been found that, despite a reduction in the aspect ratio of meta-atoms of nearly 20%, which is beneficial to silicon master etching, NLP filling, and the mold release processes, imaging efficiency can go up to 54% at wavelength of 532 nm. In addition, it demonstrates acceptable imaging quality at wavelengths of 473 and 671 nm. Moreover, the influence of fabrication errors and nanoimprinting material degradation in terms of residual layer thickness, meta-atom loss or fracture, thermal-induced dimensional variation, non-uniform distribution of TiO2 particles, etc., on the performance is investigated. The simulation results indicate that the proposed device exhibits a high tolerance to these defects, proving its applicability and robustness in practice.

Capture Efficiency Analysis in the Circular Restricted Three-body Problem

Temporary capture efficiency is studied in the framework of the circular restricted three-body problem in two steps. First, a non-uniform distribution of test particles around the secondary’s orbit is obtained by fully accounting the secondary’s gravitational influence. Second, the capture efficiency is computed based on the non-uniform distribution. Several factors influencing the result are discussed. By studying the capture efficiency in the circular restricted three-body problem of different mass ratios, a power-law relation between the capture efficiency (p) and the mass ratio (μ) is established, which is given by p ≈ 0.27 × μ 0.53, within the range of 3.0035 ×10−6 ≤ μ ≤ 3.0034 × 10−5. Taking the Sun–Earth system as an example, the influence from the orbit eccentricity of the secondary on the non-uniform distribution and the capture efficiency is studied. Our studies find that the secondary’s orbit eccentricity has a negative influence on the capture efficiency.

Particle image based simultaneous velocity and particle concentration measurement

The aim of this study is the expansion of the application of particle image velocimetry (PIV) to include the determination of particle concentration within the visualized area, in addition to velocity analysis. The assessment of particle concentration is valuable in various lab-scale experiments involving particle dispersion. Additionally, it plays a crucial role in evaluating the quality of PIV images. The research investigates two particle image-based concentration techniques: the exponential averaging-based sliding method and the Voronoi cell-based method on the particle images. The exponential averaging method provides a straightforward approach, utilizing a constant length scale for sliding average application to particle images. However, this method may result in broadened interfaces or a ‘marker-shot’ effect at low concentrations, making it less suitable for scenarios involving highly non-uniform particle distributions, such as concentrated jet emissions into ambient environments. Consequently, detecting interfaces in such cases requires additional effort for reliable results. In contrast, the Voronoi cell-based technique offers the advantage of spatially adaptive resolution, making it well-suited for variable concentration distributions and situations where interface detection is crucial. To comprehensively evaluate the performance of these techniques, a synthetic test case was generated to simulate a diffusion problem featuring an initial step in concentration distribution. Both the exponential averaging and Voronoi cell-based methods were applied and compared using this synthetic test case. Additionally, the effect of particle–particle overlap is analyzed theoretically and experimentally with uniform concentration and comparison with particle counter measurements. A modified Voronoi method is introduced, providing flexibility in capturing a wide range of concentration regions and features. An example experimental scenario involving a turbulent puff was considered demonstrating the application of the developed methods. The results demonstrate that the Voronoi method effectively captures small structures with high concentrations while providing reliable results in regions with low concentrations.

On the Preferential Concentration of Particles in Turbulent Channel Flow: The Effect of the Added-Mass Factor

Preferential concentration, observed in turbulent flows when particle response times are of the same order of the flow’s characteristic timescales, manifests as non-uniform particle distributions in space. Unraveling its governing mechanisms holds crucial implications for both natural and industrial processes reliant on particle-laden flows. Focusing on particles with small inertia, this study employs Direct Numerical Simulations coupled with Lagrangian particle tracking to investigate the influence of the added-mass factor on the preferential concentration of particles denser than the fluid in the one-way coupling regime. It is shown how the added-mass factor β affects particle distribution within the channel through the statistical correlations between particle concentration and typical descriptors of the flow topology. The results suggest that increasing values of β (corresponding to lighter particles) significantly reduce the effectiveness of turbophoresis in producing particle accumulation in the near-wall region. Resulting in a gradual decorrelation between particle concentration and both the strain-rate and the vorticity tensors, higher values of β lead to a more uniform particle distribution, regardless of the Stokes number.

Photocatalytic degradation of 2,4-dichlorophenol using nanomaterials silver halide catalysts

In this study, the photocatalytic activity of nanomaterials Ag/AgX (X = Cl, Br, I) is reported. Highly efficient silver halide (Ag/AgX where X = Cl, Br, I) photocatalysts were synthesized through a hydrothermal method. The samples were characterized using a range of techniques such as X-ray diffraction (XRD), scanning electron microscopy (SEM), and Brunauer–Emmett–Teller (BET) to check their structural, morphology, textural and optical properties. In addition, the photocatalytic activity of photocatalysts was evaluated through the degradation of 2,4-dichlorophenol (2,4-DCP) under UV and visible light irradiation. XRD analysis confirmed the presence of a single-phase structure (pure phase) in the synthesized photocatalysts. SEM micrographs showed agglomeration with a non-uniform distribution of particles, which is a characteristic of surfactant-free precipitation reactions in aqueous media. The Ag/AgBr photocatalyst exhibited the best degradation efficiency, resulting in 83.37% and 89.39% photodegradation after 5 h of UV and visible light irradiation, respectively. The effect of catalyst loading, initial solution pH, and 2,4-DCP concentration was investigated for the best-performing Ag/AgBr photocatalyst. The degradation kinetics were best described by the pseudo-first-order Langmuir–Hinshelwood model. The photocatalytic capacity of Ag/AgBr decreased by 50% after five reuse cycles. SEM images revealed heightened levels of photodegradation on the catalyst surface. The study proved the feasibility of using simple synthesis methods to produce visible light active photocatalysts capable of degrading refractory phenolic pollutants in aqueous systems.

Study on motion behaviour of coal water slurry particles under vibration conditions

This study uses computational fluid dynamics (CFD) to simulate the motion behavior of coal-water slurry (CWS) particles under vibration conditions. Results show that without vibration, the coal particles settle to the bottom of the container due to gravity, forming distinct regions of low, transition, and high concentration. Under low vibration intensity, the slight energy activates particle motion and promotes settlement. Under high vibration intensity, violent reciprocating motion causes severe shaking and uncertainty, resulting in non-uniform particle distribution and back mixing. Particle velocity distribution is significantly influenced by vibration, with higher frequencies and amplitudes resulting in greater velocities. Additionally, particle pseudo-temperature is higher in the near-wall area due to collisions with particles and wall surfaces. Low frequency and amplitude promote particle settlement, increasing compactness in the bottom area and reducing concentration in the top area. The findings provide valuable insights into CWS behavior under vibration conditions, which can be used to optimize the design and operation of CWS handling and transportation systems.

Numerical evaluation of space form design for outpatient waiting space to reduce airborne infection risk in hospitals

The outpatient waiting space (OWS), with high density and frequent variations of occupants, represents a high-risk space for the spread of respiratory infectious diseases (RIDs) in hospitals. While previous studies have proposed strategies related to occupant behavior and environmental quality, there is a limited exploration of the impact of space form design on infection risk in OWSs. This study developed a numerical method for assessing the airborne infection risk, considering the non-uniform distribution of infectious particles and variations in occupant distributions. Validated computational fluid dynamics (CFD) simulations were employed to calculate the particle distributions. The intake fraction, accumulated throughout the waiting process, was then integrated into the Wells-Riley model for risk evaluation. Using the typical ventilation form of ceiling supply and ceiling return, with an air change per hour (ACH) set to 6 h−1 as the baseline, 12 typical space forms, including four space shapes (L shape, T shape, U shape, and O shape) and three entrance layouts (single-side entrance, middle entrance, and double-side entrances), were selected for comparisons. The results showed significant variations in infectious particle distributions across the 12 space forms, with average infection probabilities (IPs) ranging from 3.47 % to 19.75 %, suggesting that appropriate space form design has the potential to reduce IP by 16.28 %. These findings offer a novel perspective on mitigating airborne infection risk in hospitals by emphasizing the significance of architectural design. Moreover, precise space design guidance is proposed for OWS to safeguard human health against the spread of RIDs.

A novel counter-rotating twin-tool friction stir processing technology for fabrication of Al/SiC surface composite with improved particle distribution

The fabrication of surface composite through friction stir processing (FSP) introduces certain challenges, including insufficient heat generation and non-uniform distribution of reinforcement particles. The implementation of multi-pass FSP can partially reduce these issues; however, it comes at the cost of increased manufacturing time and labour costs. In order to address these challenges comprehensively, a novel counter-rotating twin-tool friction stir processing (CRTTFSP) technology was utilized for fabrication of Al/SiC surface composite. The presence of hard SiC particles notably improved the microhardness of the composite. High heat generation and counter-rotation of the twin-tool, contributed to better particle distribution, which led to consistent microhardness values. Hence, this technology offers uniform mechanical properties of the fabricated composite and also eliminates the need of repetitive processing over the same zone, resulting in considerable reductions in both manufacturing time and associated costs making it more suitable for industrial applications.

A numerical analysis of particle encapsulation in a flow-focusing droplet generation device

In this paper, the process of encapsulating solid particle(s) into liquid droplets in a high-throughput flow-focusing microchannel is investigated numerically. Open source software is used, which computes fluid flow in an Eulerian framework and particle dynamics with a Lagrangian approach. Previous studies have demonstrated that if no action is taken, particles suspended in a liquid passing through a flow-focusing microchannel will be encapsulated at random. This is perhaps unsurprising, but in one such study, less than 35% of droplets were found to contain exactly one particle. The two aims of this study are (i) to explore the flow patterns arising in a microfluidic channel and (ii) to elucidate the effect of salient governing parameters on encapsulation efficiency (i.e., the fraction of droplets encapsulating one particle) by focusing on ordering the particles before reaching the droplet generation section. Following validation against experimental reference data, the capillary number is varied across the three droplet generation regimes: squeezing, dripping, and jetting. We demonstrate that under certain conditions, an encapsulation frequency of 100% can be achieved with ordered particles, but in most cases, this is significantly lower. We examine the flow field to help understand how this non-uniform distribution of particles occurs. Notably, we find the dripping to be the best option for particle encapsulation and in this case extend the study to explore the effect of junction angle, finding that an angle of 60° is the most favorable. Improved understanding of the encapsulation process derived from this study can help to improve design of high-throughput droplet generation microfluidic systems.

Assessing the conveyance of fluidized natural gas hydrate in a horizontal tubing

This study evaluates the transportability of clayey-silt natural gas hydrate (NGH) solid particles during solid fluidization in horizontal tubing, a crucial step in recovering marine NGH. Non-uniform solid particle distribution in horizontal flows can cause unstable operations and financial losses. Thus, assessing transportability is necessary to avoid these issues and harness the potential of NGH as an alternative energy source. This study introduces a model based on each phase's two-layer approach and momentum balance equation, considers potential hydrates (re)formation, and overcomes previous models’ limitations. Moreover, it assumes that the critical settling velocity marks the moment when the lower bed of solid particles stagnates, thereby leading to the disappearance of frictional pressure in the lower bed. The proposed model is evaluated and compared with existing models using a laboratory dataset and the South China Sea offshore field dataset and found to provide the highest prediction accuracy and precision. Moreover, it reveals that the tubing diameter, solid particle diameter, and solid particle concentration are the major influential parameters affecting the solid particle conveyance. Consequently, the proposed model is suitable for evaluating the transportability of those solid particles in horizontal tubing, making it ideal for assessing their transportability in pipelines.

A model for proppant dynamics in a perforated wellbore

This paper presents a model to simulate behavior of particle-laden slurry in a horizontal perforated wellbore with the goal of quantifying fluid and particle distribution between the perforations. There are two primary phenomena that influence the result. The first one is the non-uniform particle distribution within the wellbore’s cross-section and how it changes along the flow. The second phenomenon is related to the ability of particles to turn from the wellbore to a perforation. Consequently, the paper considers both of these phenomena independently at first, and then they are combined to address the whole problem of flow in a perforated wellbore. A mathematical model for calculating the particle and velocity profiles within the wellbore is developed. The model is calibrated against available laboratory data for various flow velocities, particle diameters, pipe diameters, and particle volume fractions. It predicts a steady-state solution for the particle and velocity profiles, as well as it captures the transition in time from a given state to the steady-state solution. The key dimensionless parameter that quantifies the latter solution is identified and is called dimensionless gravity. When it is small, the particles are fully suspended and the solution is uniform. At the same time, when the aforementioned parameter is large, then the solution is strongly non-uniform and resembles a flowing bed state. A mathematical model for the problem of particle turning is developed and is calibrated against available experimental and computational data. The key parameter affecting the result is called turning efficiency. When the efficiency is close to one, then most of the particles that follow the fluid streamlines going into the perforation are able enter the hole. At the same time, zero efficiency corresponds to the case of no particles entering the perforation. Solutions for the both sub-problems are combined to develop a model for the perforated wellbore. Results are compared (not calibrated) to a series of laboratory and field scale experiments for perforated wellbores. Comparison with the available computational results is presented as well. In addition, the comparison is presented in view of the parametric space defined by the dimensionless gravity and turning efficiency. Such a description allows to explain seemingly contradictory results observed in different tests and also allows to highlight parameters for which perforation orientation plays a significant role.

A multiscale computational framework using active learning to model complex suspension flows

We devise a novel computational framework to achieve efficient multiscale modeling for suspension flows. The modeling framework comprises a particle-resolving direct numerical simulation for microscopic computation, a single-fluid continuum model to capture the bulk flow behavior at the macroscale, and a Gaussian process regression that connects the two modeling components. The microscopic calculation simulates detailed flow fields around individual suspended particles. Along with a coarse-grained method, we are able to obtain the mean volume fractions and stresses for discrete particles, which are then further used to calculate the rheological properties of the bulk flow fields at the macroscale under specific forcing conditions. Using Gaussian process regression, a few data points from the microscopic calculations can be further interpolated/extrapolated to form complete constitutive relationships that cover the whole forcing range in macroscopic calculations (i.e., continuum modeling). Moreover, via the uncertainties returned from the Gaussian process regression, the resulting active learning strategy using the Gaussian process regression automatically decides on the required the microscopic calculations given by a specified tolerance. As a result, complex suspension problems can be efficiently solved by the macroscopic continuum model to the desired accuracy with the minimum computational effort (i.e., minimum number of microscopic runs). Four examples are demonstrated, among which two examples involve particle migration as a function of the particle stress, and the resulting nonuniform distribution of particles. This demands a multivariate Gaussian process regression, which is found to be particularly effective in achieving high computational efficiency compared with the off-line cases, i.e., a large number of off-line microscopic calculations to obtain the prescribed constitutive relationships.

Spray behaviour of hydro-treated ester fatty acids fuel made from used cooking oil at low injection pressures

AbstractThe spray characteristics significantly affected the combustion performance. The injection pressure and fuel properties are factors that affect the spray cone angle, penetration, and droplet distribution. Although substantial research has been conducted on spray attributes, understanding the complex biofuel spray dynamics in real nozzles and injectors is crucial. This study examines hydro-processed esters and fatty acid (HEFA) spray characteristics of used cooking oil and palm oil biodiesel in a constant-volume chamber. The study was performed by varying the injection pressures ranging from 30 to 120psi for pure fuels and Jet A-1 blends. Experiments were conducted at standard sea-level atmospheric pressure and an ambient temperature of ±297K using an airblast fuel injector. The initial fuel temperature was set at ±302K. Jet A-1 was used as the baseline fuel for the comparative analysis. Particle image velocimetry (PIV) was employed to visualise the microscopic and macroscopic characteristics of the fuel spray. The results revealed a progressive increase in the penetration length corresponding to an increase in the injection pressure. Lower pressures yielded non-uniform particle distributions across the spray area, whereas cone angle augmentation stabilised at elevated pressures. A blend of 60% Jet A-1 and 40% HEFA closely matched Jet A-1 characteristics, indicating alternative aviation fuel potential. These real-time insights into spray behaviour are critical for enhancing the fuel efficiency and mitigating the generation of particulate emissions resulting from spray combustion.

An adaptive multi-resolution SPH approach for three-dimensional free-surface flow with fluid impacting

This paper presents a new adaptive multi-resolution smoothed particle hydrodynamics (SPH) approach for simulating three-dimensional free surface flows with fluid impacting. The governing equations for multi-resolution SPH are derived based on the finite particle method (FPM) particle interpolation, which is an alternative SPH formulation for improving interpolation accuracy in the non-uniform particle distribution. A variable kernel length related to the volume of neighbor particles is proposed to control the numerical instability due to the multi-resolution scheme. In adaptive particle refinement (APR), a concept named as continuously-variable characteristic length (CCL) for each particle is defined to determine the splitting or merging scheme and then build a smoothing non-uniform particle distribution. Three-dimensional examples are conducted to validate the presented new APR approach, and the capability and numerical efficiency of the proposed numerical model are verified through solving free-surface flow problems with middle and low Weber numbers. After validation, a refined simulation of droplet impacting is conducted using the proposed adaptive multi-resolution SPH approach. The results show that the proposed SPH could well adapt to the large deformation of three-dimensional free-surface flow, and only provide a refinement near the free-surface domain.

A dynamic load balancing algorithm for CFD–DEM simulation with CPU–GPU heterogeneous computing

Computational fluid dynamics combined with discrete element method (CFD–DEM) is widely adopted to study particle–fluid multiphase systems. However, due to the global dynamical evolution and the naturally non-uniform distribution of particles, the computational load in parallel computing varies dynamically in each process, leading to significant load imbalance. Thus, a dynamic load balancing (DLB) algorithm is proposed for parallel CFD–DEM simulations in the CPU–GPU heterogeneous architecture. A type of computing load analyzing grid (AG) is developed which covers the entire simulation region. Then Kalman filtering and linear weighted method is used to solve the weight of fluid grids and particles. The domain decomposition scheme is updated according to the weights of AGs dynamically. This DLB method is tested for gas–solid fluidized beds and periodic gas–solid systems of different scales. Results show that the DLB algorithm can significantly speed up CFD–DEM simulation and the highest performance is 2.28 times faster with DLB.

A Laboratory Study of Use of Crumb Rubber as Partial Replacement of Fine Aggregate in Concrete

Abstract: The aim of this study is to analyse the effect of addition of crumb rubber on the various properties of concrete. The crumb rubber is used as replacement over aggregates in the concrete mix. Crumb rubber is prepared from the scraps of tyres. In this study aggregates are replaced by 0%, 30%, 40%, 60%, 100% crumb rubber. The rubberized concrete is tested for slump and compression strength. It is found that the slump of rubberized concrete increases first but as the amount of rubber is increased the slump starts decreasing. It is observed that initial compression strength of rubberized concrete reduced significantly but the final strength is found to be more than that of ordinary concrete. Waste tyres cause tremendous pressure and ecological issues for the entire tyre industry when accumulated in cultivated land or combusted. Crumb rubber (<5 mm) grinded from waste tyres was introduced to substitute for natural fine aggregate in concrete and effectively solves the consumption challenge. This paper reviews the performance of concrete with crumb rubber as fine aggregate, providing evidence for crumb rubber concrete (CRC) materials design and application. Crumb rubber is characterized with light specific gravity, hydrophobicity and air entrapment in comparison with natural fine aggregate, inducing significant reduction in workability of fresh CRC and exhibiting poor bonding performance with cementitious matrix. By summarizing the compressive/tensile strength, elastic modulus and fracture behaviour of CRC at various rubber content, prediction models of strength reduction factors are proposed. The main reasons for strength deterioration are weak interface transition zone (ITZ) performance and non-uniform distribution of rubber particles. Water/chloride permeability, electrical/carbonation resistance and drying shrinkage of CRC are discussed for durability performance of CRC. Physical/chemical pre-treatments of rubber could alleviate the hydrophobicity and improve the mechanical and durability performance of transition zones between rubber and cement paste. Accordingly, scopes concerning the recycling of crumb rubber and its performance optimization are expected in future studies.

The counterbalance of endothelial glycocalyx and high wall shear stress to low-density lipoprotein concentration polarization in mouse ear skin arterioles

Background and aimsAtherosclerosis preferentially occurs at regions in arterial branching, curvature, and stenosis, which may be explained by the geometric predilection of low-density lipoprotein (LDL) concentration polarization that has been investigated in major arteries in previous studies. Whether this also happens in arterioles remains unknown. MethodsHerein, a radially non-uniform distribution of LDL particles and a heterogeneous endothelial glycocalyx layer in the mouse ear arterioles, as shown by fluorescein isothiocyanate labeled wheat germ agglutinin (WGA-FITC), were successfully observed by a non-invasive two-photon laser-scanning microscopy (TPLSM) technique. The stagnant film theory was applied as the fitting function to evaluate LDL concentration polarization in arterioles. ResultsThe concentration polarization rate (CPR, the ratio of the number of polarized cases to that of total cases) in the inner walls of curved and branched arterioles was 22% and 31% higher than the outer counterparts, respectively. Results from the binary logistic regression and multiple linear regression analysis showed that endothelial glycocalyx thickness increases CPR and the thickness of the concentration polarization layer (CPL). Flow field computation indicates no obvious disturbances or vortex in modeled arterioles with different geometries and the mean wall shear stress is about 7.7–9.0 Pa. ConclusionsThese findings suggest a geometric predilection of LDL concentration polarization in arterioles for the first time, and the existence of an endothelial glycocalyx, acting together with a relatively high wall shear stress in arterioles, may explain to some extent why atherosclerosis rarely occurs in these regions.