Laser Sheet Research Articles

Sensitivity of wavelet-based optical flow velocimetry (wOFV) to common experimental error sources

Abstract The influence of several potential error sources and non-ideal experimental effects on the accuracy of a wavelet-based optical flow velocimetry (wOFV) method when applied to tracer particle images is evaluated using data from a series of synthetic flows. Out-of-plane particle displacements, severe image noise, laser sheet thickness reduction, and image intensity non-uniformity are shown to decrease the accuracy of wOFV in a similar manner to correlation-based particle image velocimetry (PIV). For the error sources tested, wOFV displays a similar or slightly increased sensitivity compared to PIV, but the wOFV results are still more accurate than PIV when the magnitude of the non-ideal effects remain within expected experimental bounds. For the majority of test cases, the results are significantly improved by using image pre-processing filters and the magnitude of improvement is consistent between wOFV and PIV. Flow divergence does not appear to have an appreciable effect on the accuracy of wOFV velocity estimation, even though the underlying fluid transport equation on which wOFV is based implicitly assumes that the motion is divergence-free. This is a significant finding for the broader applicability of planar velocimetry measurements using wOFV. Finally, it is noted that the accuracy of wOFV is not reduced notably in regions of the image between tracer particles, as long as the overall seeding density is not too sparse i.e. below 0.02 particles per pixel. This explicitly demonstrates that wOFV (when applied to particle images) yields an accurate whole field measurement, and not only at or adjacent to the discrete particle locations.

Three-dimensional dynamics of unstable lean premixed hydrogen-air flames: Intrinsic instabilities and morphological characteristics

The 3D swinging laser sheet technique was employed to study the development and morphological characteristics of premixed hydrogen-air unstable flames in a spherical explosion vessel. Pressure dependencies for laminar flame propagation were sought to exploit the role of the Darrieus-Landau (DL) and Thermal-diffusive (TD) instabilities in the unstable self-accelerating flame regime. A sufficiently low Markstein number, as a consequence of the increased pressure, leads to more cracking and smaller cells over the flame surface. The degree of wrinkling on the flame surface is proportional to the increase in flame burning velocity, a relationship that holds true for low pressures but is not applicable under high pressures. External turbulence can significantly alter the extent of flame surface wrinkling even at low root mean square velocities, producing a more wrinkled flame surface compared to intrinsic cellularity, and distinctly affecting flame dynamics. The increased wrinkling and flame speed due to external turbulence can be attributed to the synergistic effects between thermo-diffusive instabilities and turbulence, resulting in higher fuel consumption rates per flame surface area and the formation of finger-like structures that enhance flame displacement speed in curved segments. The parameters, ϵ, deviation of the Lewis number from a critical value, and ω2, obtained through classical linear stability analysis, display a clear linear relationship with the ratio of the wrinkled surface area observed in planar flames. This study enhances the understanding of hydrogen flame instabilities, which is crucial for preventing explosions in hydrogen storage and utilization, and provides valuable insights into flame dynamics, supporting the design of safer and more efficient hydrogen-fueled engines and turbines.

An investigation on hydrodynamic behavior of horizontal gas–liquid intermittent flow by S-PLIF and parallel-wire conductance sensors

Horizontal gas–liquid two-phase flows widely exist in a variety of industrial fields such as the environment, chemical industry and energy. Intermittent flow, characterized by the quasi-periodic alternation motion of liquid slugs and large-scale Taylor bubbles, is the most prevalent flow pattern in horizontal gas–liquid two-phase flows. Comprehensive research on the hydrodynamic behavior of gas–liquid intermittent flow is essential for revealing the mass and heat transfer characteristics between phases, thereby establishing physical models of flow parameters. In this study, a structured planar laser-induced fluorescence (S-PLIF) system and paired parallel-wire conductance sensors (PWCSs) are designed to detect the gas–liquid interface structures and the hydrodynamic parameters. The S-PLIF system consists of a laser, a visualization box, a Ronchi ruling-plate and a pair of high-speed cameras. An optical linear prism forms a laser sheet in the transverse direction and the middle of the test section, and the cameras are positioned at 70° and 90° respectively to the plane of the laser sheet, referred to as S-PLIF70 and S-PLIF90. Under the laser excitation, bright and dark stripes generated by the Ronchi ruling-plate are deflected at the gas–liquid interface. The influence of total reflection can be effectively avoided by identifying the deflection position of the stripes, allowing for accurate detection of the gas–liquid interfacial structures. The effects of the interface fluctuations and detached bubbles on the S-PLIF70&90 results are investigated. Fluctuation and evolution characteristics of the gas–liquid interface are analyzed by probability density function (PDF) and power spectral density (PSD). Meanwhile, hydrodynamic parameters, i.e., Taylor bubble nose and tail velocity, Taylor bubble and slug length, and slug frequency, are detected by the PWCSs, and compared with previous correlations. The evolution characteristics of the interfacial structures and hydrodynamic parameters are uncovered as the flow condition changes.

Particle Shadow Velocimetry And Its Potential Applications And Advantages Vis-à-Vis Particle Image Velocimetry

Particle image velocimetry (PIV) is an established standard velocimetry technique in experimental fluid mechanics that involves determining a fluid flow velocity field from the motion of tracer particles illuminated by a laser sheet. The necessity of laser illumination poses challenges in certain applications, and is a potential entry barrier due to its complexity and high cost. A laser-free alternative to PIV is particle shadow velocimetry (PSV), which uses images of the shadows cast by the particles on the camera sensor under back-illumination, instead of the Mie scattering signal produced by laser illumination. This study aims to compare various aspects of PSV such as depth of field, seeding density, type of illumination required, particle size, image filtering, cost effectiveness, and limitations with those of PIV. PSV and PIV measurements are conducted in the wake of a flow past a cylinder and in a boundary layer developing over a flat plate. It is found that PSV is capable of achieving equivalent accuracy to PIV and is a viable alternative to PIV in certain applications where light sheet illumination creates experimental challenges.

Planar Concentration Field Measurement In A Transient Jet Using Mie Scattering

An experimental method for time resolved concentration measurements in the near field of a transient starting turbulent air jet is presented. Planar Laser Mie Scattering (PLMS) using DEHS droplets produced by a Laskin nozzle as tracers of concentration is applied on an air jet with Reynolds number of 62.000 and is recorded with 10.000 frames per second. A jet penetration of about 20 nozzle diameters into the ambient environment is resolved with more than 100 frames. A linear response between grey scale values and concentration is found for mixtures up to 13% of the air-particle flow from the Laskin nozzle mixed with clean air. The calibration is based on images with no particles and images with a homogeneous particle concentration. A correction for light scattered by particles outside the laser sheet is proposed. The experimental results agree well with mean concentration measured with other techniques.

PaIRS-UniNa: A Robust And Accurate Free Tool For Digital Particle Image Velocimetry And Optical Camera Calibration

This contribution presents PaIRS-UniNa (PArticle Image Reconstruction Software – University of Naples “Federico II”), a free software application designed to perform digital particle image velocimetry (PIV). PaIRS-UniNa relies on a powerful C library designed for speed, robustness and accuracy, and a user-friendly graphical interface (Figure 1) running in Python environment on all operating systems, including Windows, Linux and MacOS. At the current stage, it features modules to perform two-component (2C) PIV, stereoscopic PIV and optical camera calibration for multi-camera systems. The processing of PIV and stereoscopic PIV data is based on an iterative image deformation method, which can be customized by the user varying different parameters in such a way to find the best compromise between accuracy and computational speed. These parameters include the size, shape and number of interrogation windows, the number of process iterations and the type of interpolating functions for both image deformation and displacement interpolation (such as bilinear, bicubic, methods based on simplexes or the FFT shift theorem and B-splines of varying orders). Additionally, it is possible to select a range of weighting windows for cross-correlation estimation and velocity field filtering, including rectangular, Blackman, Gaussian, Harris, Hann windows, among others, and several vector validation strategies for outlier detection and replacement. Such a level of customization enables to effectively adjust the modulation transfer function of the PIV process to meet specific application requirements. For the stereoscopic PIV process two different approaches are possible: mapping, which involves applying PIV to de-warped images for 2D-2C displacement fields, and warping, which applies PIV to raw images and combine 3C reconstruction and image dewarping in a single subsequent step. In order to accurately estimate the laser sheet plane location and correct errors due to misalignment and inaccurate displacement vector positions, a disparity correction procedure is also available. The calibration module (CalVi) supports accurate calibration with the camera models mostly used in the PIV community: polynomials, rational functions and the pinhole camera model. Moreover, it supports the integration of the pinhole camera model with a refractive correction model for cylindrical geometries and camera calibration procedures working with unknown positions and orientations of the calibration target. In such a way, the target may be moved by hand without losing accuracy in comparison with a conventional calibration, in which the target is displaced with the aid of translational and rotational stages and knowledge of the target positions is strictly required. Already used in numerous experimental studies by the Experimental Thermo-Fluid Dynamics group of University of Naples “Federico II” and international collaborating groups, PaIRS-UniNa offers a rapid and effective free alternative to commercial PIV software for research and industrial use.

High-Speed Schlieren Imaging And PIV Of Unsteady Free Jet Injection With Plume Formation

This work aims to investigate Schlieren image correlation velocimetry for flow field quantification of gaseous jet injection in non-manipulable apparatus like flow chambers. Such flow cases often struggle using tracer-based methods: Seeding particles may locally cause technical problems, or the plant design prevents the principal opportunity of tracer adding, both of which lead to a reduced method repertory. This is specifically the case under flow conditions, where contamination of apparatus internals by seeding particles should be avoided completely. Therefore, in this paper, we explore a seedless, refraction-based method to quantify the velocity field of a jet plume. As an experimental flow case, we examine a transient jet release into a closed compartment under nearly atmospheric conditions. For simplification and greater clarity, the laboratory model geometry is kept essentially two-dimensional, thus we chose a flat-jet nozzle as injection device. The experimental investigations cover high-speed Schlieren imaging and laser sheet visualization, followed by digital image correlation. Combining the traditionally qualitative Schlieren approach for flow visualization with Fast Fourier Transform and cross-correlation algorithms, we calculated the velocity vector field of an unsteady jet plume formation. Furthermore, we determined the axial jet velocity profile at steady-state conditions. Our research findings highlight the applicability of Schlieren image correlation velocimetry to quantify gaseous jet formation at encased flow sections spatially. Concluding, the results coincide with data from particle image velocimetry measurements and point out the potential of expanding flow diagnostics to hitherto hardly explored (industrial) flow cases.

A novel powder sheet laser additive manufacturing method using irregular morphology feedstock

Irregular (i.e. non-spherical) morphology powder is more cost-efficient to produce than spherical shaped powder. However, its reduced levels of flowability limit the wide application ranges of laser beam powder bed fusion (LPBF). To address this issue, a novel powder sheet additive manufacturing concept (MAPS) is proposed. Herein, a pre-manufactured metal particle-polymer binder composite (i.e. powder sheet) feedstock is employed as a raw material. The printability of irregularly shaped powder particles in sheet (MAPS) format was physically investigated using a range of different process parameters. The manufacturing process was observed by high-speed imaging. Microstructural and chemical element characterisations of the irregularly shaped powder particle print were then compared against those prints which were conducted using sheet-based (MAPS) spherical powder morphologies. The results indicated that the geometric accuracy and density of irregular powder morphology sheet printing improved when using a negative defocus strategy of the laser beam. A negative defocusing strategy allowed for the melting mode of the material to be transformed from keyhole to a more favourable conduction state. High speed imaging revealed that more spatter and vapour plume were observed with the increase in the magnitude of the negative defocus. The multi-morphology 304 L stainless steel (SS304) samples were printed in a single printer using an efficient method for the first time, i.e. printing spherical SS304 material on top of irregular SS304. EDX results indicated an insignificant change of chemical elements between the spherical and irregular prints. EBSD results revealed that columnar grains could grow through the irregular-spherical transition zone and similar grain size can form between spherical and irregular prints. The results of this study provide insights into the optimum printing configurations for powder sheet additive manufacturing using a cost-effective solution of irregular morphology material.

In-cylinder flow evolution in the horizontal plane of a motoring compression ignition engine

This study used the two-dimensional particle image velocimetry technique to quantify the in-cylinder horizontal plane velocity field evolution in a swirl-supported light-duty single-cylinder diesel engine. The data were acquired at a constant engine speed of 1600 revolutions per minute. For each case, the distance of the laser sheet from the fire deck was varied (z = 5, 10, and 20 mm) to investigate the axial variations in the flow field during the flow evolution in the compression stroke. A vortex identification algorithm was used to detect the swirl center and its deviation from the rigid body rotation. A Bessel fit was obtained using the experimental data. The result revealed that the in-cylinder flow was not axisymmetric. The swirl center approached the geometrical center as the piston approached the top dead center. The flow evolved at the farthest plane from the fire deck. The axial diffusion of angular momentum resulted in the formation of the swirl flow structure in the plane closer to the fire deck. Angular momentum analysis of a simplified geometry has been presented to explain the swirl amplification. The estimated results were compared with the experimental results to show the momentum stratification in the engine cylinder later in the compression stroke.

Implementation of inverse method for identification of process parameters in a laser heating process

Changes in process parameters can alter the temperature data of a solid body during the laser heating process. This paper proposes an inverse method to determine the laser heating parameters, namely, sheet width, sheet thickness, laser power, and scan speed of the work specimen for the given temperature response which is assigned by the user. The method uses an analytical tool on moving heat source as a forward model capable of real-time temperature prediction of the sheet under the laser heating process. The effectiveness of the present forward model that incorporated the temperature-dependent material properties is tested by experimental findings performed on Al 6061-T6 sheets. Next, an iterative search process based on heuristic method is carried out for satisfying a prescribed temperature response by the optimization of unknown parameters. To illustrate the implementation and assess the practicality of the inverse method, two examples based on laser heating applications are used. The method proposed herein recover the unknowns for the prescribed temperature response after a few iterations, provides an efficient way to optimize the laser heating process. Additionally, the effect of measurement error on the findings of the inverse problem is addressed.

Forced air flow through a rectangular channel with 3D turbulence enhancers: visualization of flow structures by laser sheet scattering

In recent times, the design of heat transfer devices is strongly evolving due to both the diffusion of additive manufacturing techniques, that use a wide variety of materials, and Computational Fluid Dynamics (CFD) -based shape optimization techniques. Considering high-aspect ratio rectangular channels, that represent the geometry of heat transfer devices used, among other applications, for batteries cooling systems and compact automotive heat exchangers, the first step of a project towards the design of optimized turbulence enhancers is the validation of the CFD methodology. In particular, since it is important to have a fast-running 3D simulation with physical meaning, several techniques of increasing complexity and computational requirements and times are currently under scrutiny, and their validation is mainly carried out by experiments that measure both local and global heat transfer enhancement. In this context, flow visualization is a powerful tool to get insights of the physical mechanisms that generate convective heat transfer enhancement. In this paper, the setup of a laser sheet scattering method to visualize turbulent structures in a Plexiglas, 1:10 aspect-ratio, rectangular channel ribbed on one of the main surfaces is presented, and results on two different cutting planes are presented for 90° and V-down squared ribs for Re = 500.

On the errors and their mitigation in planar laser induced fluorescence

In planar laser induced fluorescence (PLIF), scalar concentration fields are obtained from fluorescence images of the flow where the scalar is ‘tagged’ with a suitable fluorescent dye and illuminated with a laser sheet. While the fluorescence occurs, the intensity of the laser sheet is attenuated along its passage through the dye field. This paper discusses the effect of correcting for the attenuation of the laser by the dye while obtaining the scalar concentration field from the fluorescent images of the flow. If the attenuation correction cannot be done, the paper discusses three possible calibration scenarios to minimise the error. It is shown that the error in ignoring the attenuation due to turbulent concentration fluctuations cannot be corrected by any means. However, the remaining error can be exactly corrected if it were possible to calibrate with the mean concentration field. It is shown that a uniform concentration field, centred around the mean field, greatly reduces the error when compared to having no attenuation correction but still is significantly more when compared to the use of the average field. Relying on these results, we propose an iterative technique which ideally achieves the accuracy obtainable when calibrating with the mean concentration field, but which requires no reference cell. In addition, the paper explores the removal striations, of a common artefact during PLIF image acquisition with origins in the image optics. A technique to remove the same is proposed and its advantages discussed when compared to available techniques.

Vortex shedding from a square cylinder interacting with an undular bore wave train

In this work, we studied the wave-induced vortex generation and shedding from a semi-submerged vertical square cylinder interacting with an upcoming undular bore wave train under a shallow water configuration. This unsteady and rapid process was investigated by means of numerical and experimental approaches. A numerical simulation, solving the full turbulent viscous Navier–Stokes equations, was carried out in order to study and characterize both the undular bore wave properties and the vortex dynamics triggered during this unsteady interaction. Starting with fluid at rest, the undular bore was generated by the impulsive translational motion of a piston wavemaker at laboratory scale in both a numerical and an experimental wave tank. When the undular bore impinges on the cylinder, filamentary vortex structures were formed at the four cylinder's edges synchronized with the propagating wave motion, leading to the vortex shedding phenomena at a frequency that matched the wave instantaneous frequency. These vortices extended along the entire cylinder span under the water column, from the free surface to the seabed. At the trailing edge of the cylinder, a pairing process of two shed vortices was observed, similar to a Lamb–Oseen vortex pair. These vortices were present during the whole undular bore wave train dynamic forcing. An overall agreement was found with the experimental version of the bore–cylinder interaction, carried out in a physical wave tank. Laser sheet bore profiling and particle image velocimetry measurements of the velocity field confirmed undular bore properties, the onset of vortex formation, subsequent shedding, and pairing in the experiments performed in similar conditions with the numerical approach.

Particulate dispersion in turbulent wake of Ahmed body and experimental investigation of impact of rear slant angle

This study investigates the particle dispersion characteristics in the turbulent wake of a simplified vehicle model (Ahmed body) for two values of the rear slant angle ϕ (25° and 40°) to study the effect of flow separation. In the experiments (Reynolds number Rel=1.90 ×105), smoke particles were released from a source and visualized with a laser sheet. Concentration fields were analyzed to calculate the vertical (Sy) and lateral (Sz) smoke spread. The findings indicate that the flow topology and concentration fields in the wake are highly dependent on ϕ. In the ϕ=40° case, separation on the rear slant disrupts the trailing vortices originating on the rear slant edges and significantly alters the wake topology. The growth of vertical smoke spread saturates after the recirculation region. The uniform mixing and absence of trailing vortices concentrate the smoke particles in the model midplane. In the ϕ=25° case, the signatures of trailing vortices were observed in the concentration fields behind the model. The vertical smoke spread is less, and the lateral smoke spread is more compared to the ϕ=40° case. The growth of the smoke spread (Sy, Sz) driven by the trailing vortices persists for a long distance, even after the recirculation region. The disruption of trailing vortices brought about by the flow separation appears to be an important effect driving the vertical smoke spread in the wake. The connection between the turbulent velocity structures and concentration structures will need to be explored with combined velocity and concentration measurements in the wake.

Investigation of renal calculi fragmented tracer particles in lithotripsy model by laser speckle technique

We have developed a laser sheet method to evaluate micro-sized fragmented calcium oxalate granular particles created kidney stones suspended in a Newtonian fluid in an in vitro model. These tracer particles are examined by laser optical techniques. An optical setup is configured with a laser sheet to conduct image velocimetry on these calcium oxalate seed particles in the urine environment. The experimental setup involves the application of ultrasonic waves to fragment the calcium oxalate tiny stones of varying sizes and disperse them in random directions in the fluid. The data acquisition process employs double frame-single exposure imaging, which captures images at specified time intervals using a high-resolution CCD camera. This provides information regarding particle displacement and track the flow path within the Newtonian fluid. In addition, the bigger fragmented particles are identified, and their sizes are also measured.

Comparison between Different Optical Configurations of Active-FRAME Setup in Multispectral Imaging of Flames

Snapshot multispectral imaging of chemical species in the flame is essential for improved understanding of the combustion process. In this article, we investigate the different configurations of a structured laser sheet-based multispectral imaging approach called the Frequency Recognition Algorithm for Multiple Exposures (FRAME). Using FRAME, a snapshot of Laser-Induced Fluorescence (LIF) of Polycyclic Aromatic Hydrocarbons (PAH) excited by 283.5 nm laser and Laser-Induced Incandescence (LII) of soot particles excited by 532 nm laser are acquired simultaneously on a single FRAME image. A laminar diffusion flame of acetylene produced by a Gülder burner is used for the experiments. The standard FRAME approach is based on creating two spatially modulated laser sheets and arranging them in a cross-patterned configuration (X). However, the effect of using different configurations (angles) of the two laser sheets on the multispectral planar imaging of the flame has not yet been studied. Therefore, we have compared the FRAME approach in four different configurations while keeping the same flame conditions. First, we have compared the relation between laser fluence and LII signals with and without spatial modulation of the 532 nm laser sheet and found that both detections follow the same curve. When comparing the maps of flame species reconstructed from the standard FRAME configuration and other configurations, there are some dissimilarities. These differences are attributed to minor changes in the imaging plane, optical alignment, laser path length, different modulation frequencies of the laser sheet, laser extinction, laser fluence, etc.

Flow Visualisation by Laser Sheet in a Smoke-Tunnel

The application of wind tunnel for visualizing flow over any bluff or streamline body plays an indispensable role in concepting various aerodynamical phenomenon such as formation of vortices, generating lift and drag, calculation of velocity vector and pressure profiles etc. Flow visualization through smoke and laser optics in a subsonic wind tunnel to capture vortex formation over airfoils, cylindrical and conical bodies is the experimental approach which has been studied and manifested elaborately in the article. Using smoke as a seeding material, the flow of the same is allowed to pass over the model which is illuminated via a double pulse laser beam to visualize the formation of vortex and capture the same for future references. The experimentation of the above models is carried out in a Honeycomb subsonic wind tunnel in the Aerodynamics Lab of Aerospace Engineering in the university. The article will cover Computational Fluid Dynamics (CFD) performed on ANSYS for the particular NACA series, Cross correlation PIV for plotting velocity vector for the captured images of the seed particles and the respective images of the entire experimental set-up. The images captured during this experiment can be later used by the aerodynamicists for optimizing aerodynamic efficiency of the experimented models depending on the vortex formation and flow pattern. The images can also be used for advance calculations required in Particle Image Velocimetry (PIV) for calculating velocity of every speed particle.

FABRICATION OF REALISTIC TRANSPARENT POROUS MEDIA FOR 3D OBSERVATION OF INTERNAL MASS TRANSPORT

A new experimental system was developed to observe the migration of colloidal particles in saturated porous media in 3D, which is important in engineering processes such as contaminant migration and filtration. The system enables the observation of internal migration that is invisible from the outside by fabricating a porous body using a transparent material and a 3D printer. Non-artificial media with arbitrarily controlled internal pore sizes were fabricated by designing the porous media based on the results of numerical analysis on phase-separation phenomena. The refractive indices of the porous media and the interstitial liquid were matched, and a laser sheet was used to extract cross sections of the media. By reconstructing multiple cross-sectional images, it was possible to visualize the 3D migration of suspended particles. The developed system was then applied to the collective gravitational settling of fine particles in saturated porous media. The results showed that the particles selectively settled in the pore networks while repeatedly branching in the lateral direction.

EFFECTS OF REFRACTIVE INDEX AND ASPECT RATIO ON THE PARTICLE IMAGES IN A LEVITATED DROPLET

A technique using droplets suspended by ultrasound has attracted attention as one of the containerless processing methods. While this can avoid contamination from the container, it is known that ultrasonic levitation creates flow fields inside and outside the droplet. For more precise droplet control, it is desirable to elucidate the internal flow of the droplet, and measurements of the internal flow have been performed using the particle image velocimetry (PIV). The aim of this study is to elucidate the internal flow field behavior by solving optical problems and improving the accuracy of velocity field measurements in levitated droplets. The fluid that is to be investigated is scattered with small tracer particles and illuminated by a laser to capture the flow on a created laser sheet. The curvature distortion is successfully visualized using the PIV approach toward distortion correction using calibration methods. The curvature distortion illustration was performed based on the refractive index and aspect ratio of the simulated droplet in acrylic materials. The fluid flow, affected by droplet curvature and refractive index, has been visualized for both levitated and simulated droplets. The experimental results showed that the droplet curvature can be distorted in two types such as radial and tangential distortions and increase as the refractive index and aspect ratio increases.

VARIATION IN THE SETTLING BEHAVIOR OF PARTICULATE SUSPENSION IN SATURATED GRAVEL LAYERS

The gravitational settling of solid particles in saturated porous media, which is common in various engineering processes and natural phenomena, was investigated experimentally. This study examines the settling pattern of fine particles suspended in saturated porous media composed of gravel layers. The experiments were conducted by varying the physical properties of particles and fluid, particle concentration, and the pore scale of porous media. The three-dimensional settling behaviors of suspended particles in an invisible porous media were observed by using refractive index matching together with sheet laser scanning techniques. The experimental results showed that various settling patterns were shown depending on the particle size and concentration. Under conditions of large particle size with a small concentration, suspended particles settle individually through gaps between gravels in the direction of gravity with the velocity similar to that of an isolated particle. However, small particles with large concentrations settle collectively through the pore region as if they are one continuous fluid. In this case, the particles migrate in both vertical and lateral directions with the velocity much larger than the Stokes settling velocity. The settling velocity was quantified by image analysis of reconstructed three-dimensional images in porous media. The results indicated that the settling velocity of particles greatly depends on the settling pattern. Finally, we found that the dimensionless parameter, which describes the collectivity of suspended particles in porous media, quantitatively predicts both the transition of the settling pattern and the settling velocity.