Search NASA⌕ Search

SEARCH · Search NASA

Results for “fluid flow systems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quantifying Trapped Powder in Electron Beam Powder Bed Fusion

Abstract Electron beam powder bed fusion (PBF-EB) shows great potential for manufacturing complex parts including those with internal cavities for heat exchanger, manifold systems, or energy absorption purposes. PBF-EB allows for the manufacture of channel geometries without the need for support structures. Due to the nature of the powder spreading process, powder feedstock is often trapped in intentionally manufactured cavities. This trapped powder can often be difficult to remove and can disturb the intended flow of fluid through the cavity or damage downstream components in its use case. These trapped powder particles present a risk of contamination and component failure if not completely evacuated. Ti6Al4V is a choice material for aerospace applications due to its high strength to weight ratio and its composition as a nonferrous metal; however, in weight sensitive applications excess entrapped powders or powders loosely attached to the surface could cause undesirable weight increases. The inherent spreading process of PBF-EB is different than laser powder bed fusion (PBF-LB) in its operational temperature, sintering. In addition, PBF-EB is less commonly studied in literature compared to its PBF-LB counterpart, and as a result the complexity of the semi-sintered powder and its spreading behavior are not well understood. Prior work has investigated the difficulty in removing trapped powder from PBF-EB, but these studies do not address how to quantify the amount of trapped powder in the cavity. Thus, an accurate method to measure the amount of trapped powder in the cavity must be investigated. In this work, Ti6Al4V coupons were manufactured with horizontal and vertical cavities of three different sizes. Archimedes testing allows for the determination of density differences caused by porosity and trapped powders by measuring mass and volumetric dispersion. Computed tomography (CT) is well suited for segmenting the internal structure and features of a part and has been studied for applications including voids, porosity, and dross. Thus, CT was explored as a method for evaluating trapped powder content in this work. The volumetric representation of the segmentation of the reconstructed CT volume can vary greatly depending on the input filter and thresholding methods. In this study, four different types of segmentation approaches were evaluated to determine the best approach for segmenting the volume as compared to an operator labeled ground truth. The percentage density results from the Archimedes testing were compared to the volumetric percent density from the computed tomography approach. Differences in packing density between two different internal channel features were investigated. Overall, this work sought to validate the use of computed tomography for the detection of trapped powders and present a framework for volumetric segmentation.

Johnstone, Brian↗
DOE OSTI2022

Multiscale Interaction of Turbulence, Temperature, and Soot Formation: Measurements for Critical Assessments of Chemical Kinetics and Mechanisms

This program targeted an improved understanding of how turbulence affects soot formation chemistry (directly and indirectly) with a particular focus on detailing the multiscale coupling between flow turbulence, mixing, thermal transport, and soot formation kinetics in gas-phase reacting systems. In this project, simultaneous, quantitative multi-dimensional velocity, soot volume fraction, and gas-phase temperature measurements in turbulent non-premixed sooting flames were developed. These new diagnostics tools can be used to elucidate the relative effects of fluid kinematics, mixing, thermal processes, and turbulence/scalar time/length scales on soot formation chemistry and topology. Simultaneous multi-parameter measurements were acquired within a series of flames with different Reynolds numbers and spatial locations to investigate the inter-dependence of soot, temperature, and turbulence within the soot inception, growth, and oxidation region. Preliminary assessment of the results has been completed at the project end and detailed analysis is planned following the completed program.

02 PETROLEUM↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Thermal hydraulic analysis of the AFIP-7 irradiation test in the Advanced test Reactor – Model correlation and performance evaluation

We report as a system level code, RELAP5 is widely used in the nuclear field for the reactor hydraulic analysis. In the component/experimental level, it is frequently employed as well. This paper demonstrates that the RELAP5 flow simulation of the AFIP-7 irradiation experiment which has a complicated 3D geometry/flow path deviates by ~28% from the flow evaluation through a computational fluid dynamics (CFD) simulation, which has been validated against a flow test performed at Oregon State University. Thermal safety compliance analysis of an experiment planned to be irradiated in the Advanced Test Reactor is usually performed using the finite element analysis code, Abaqus. This research reveals the disadvantages of the Abaqus simulation in the flow instability and departure from nucleate boiling evaluations via a conjugate heat transfer analysis of the AFIP-7 irradiation experiment. As a result, a detailed CFD simulation is suggested for irradiation experiments with complicated flow paths, rather than a simplified RELAP5 simulation or hand calculations. The validated simulation approach will be integrated with the Boehmite correlations to investigate the oxide growth prediction on the fuel cladding in Part II of this research.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Numerical assessment of triply periodic minimal surfaces for direct air capture of carbon dioxide

Direct air capture (DAC) systems often consist of packing material wetted by a capture fluid that reacts with CO 2 in the airstream. The efficiency of the contactor is determined by a complex relationship of fluid dynamics, heat and mass transfer, contactor geometry, and chemical properties. The efficiency of the contactor must be balanced with other factors, primarily pressure drop through the system. Triply periodic minimal surfaces (TPMS) are a class of differential surfaces that have been explored in multiple engineering applications and have been shown to exhibit excellent performance when used in heat exchangers. Their tortuous path provides a high surface-to-volume ratio and favorable trade-off between contact area and pressure drop. In this work, a gyroid-type TPMS contactor was evaluated using computational fluid dynamics for a variety of geometric parameters to explore the potential benefit of TPMS shapes for DAC applications. A thin-film model was employed to model the flow and distribution of the capture solvent, allowing efficient simulations of TPMS structures at scale by eliminating the need for a computationally intensive interface capturing method. A liquid-gas mass transfer model was implemented in the commercial software STAR-CCM+ and used to predict the CO 2 capture efficiency and study the trade-off between capture performance and pressure drop through analysis of capture rates, mass transfer coefficients, and other relevant variables. TPMS contactors with a variety of geometric parameters and two capture solvent options were investigated to determine the effect of design choices on the operational performance of DAC systems. In conclusion, results showed that while contactor geometry is the dominant factor in efficiency and pressure drop, the physiochemical properties of the solvent are an important secondary influence on the contactor performance.

CFD↗
DOE OSTI2025

A comprehensive modeling of falling film evaporators subject to vapor flow, pass arrangements, and refrigerants

Improving the heat transfer performance of falling film evaporators is a crucial step for improving the energy efficiency of the heat pump or refrigeration systems. This study conducts a numerical investigation for the practical-scale falling film evaporator based on the epsilon–NTU method with an updated heat transfer correlation. The algorithm was validated with lab-scale and real-scale falling film evaporator experimental results, and the prediction reaches a mean absolute deviation of 12.5%. Here, the parametric study encompasses eight refrigerants: R-134a, R-410A, R-600a, R-717, R-1270, R-152a, R-1234yf, and R-1234ze(E). The results indicate that vaporization enthalpy of refrigerant is a key property in selecting an appropriate working fluid because it helps minimize severe heat transfer degradation caused by dry-out. Additionally, the vapor-flow–induced heat transfer degradation can be predicted using the critical Weber number. Furthermore, the trade-off between extending the tube length and increasing the number of tubes for heat transfer improvement is discussed. Finally, different two-pass arrangements show deviations of less than 4 %.

Dry-out↗
DOE OSTI2022

Modeling gas breakthrough and flow phenomena through engineered barrier systems using a discrete fracture approach

Compacted clays are being considered to build engineered barrier systems (EBS) intended for the safe isolation of high-level nuclear waste (HLW) and spent nuclear fuel (SNF). The corrosion of the metallic canister containing the HLW/SNF will lead to the generation and buildup of the gas pressure in the more internal part of the clay buffer. This phenomenon would eventually trigger the formation and propagation of fractures in the clay barrier, jeopardizing its safety functions. Here, in this work, we propose to use the fragmentation technique (MFT) to model evolving fractures in clays triggered by gas pressurization. The MFT has been successfully used to model the formation of fractures in concrete, drying cracks in soil, hydraulic and thermo-fractures in rocks. In this work, we extend the MFT to deal with multiphase fluid flow in deformable porous media, and we upgraded a fully coupled computer finite element code using the extended technique. The proposed approach is first verified against analytical solutions and is then applied to model gas breakthrough experiments in clays. A very satisfactory performance of the method is observed in all the analyses, showing the potential of the MFT to tackle multiphase flow problems in deformable porous media with evolving discontinuities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

ECAR-6330 Rev 0 Computational Fluid Dynamics Evaluation of Core Pressure Drop for the Primary Coolant Apparatus Test

PCAT is a thermal-hydraulic test loop for MARVEL reactor. The design gap thickness between the heater rods of PCAT core is 1.0 mm. The heater elements were manufactured within a straightness tolerance for the designed gap thickness, but several elements have been distorted out of tolerance during the welding process. The distortion of heater rods could lead to flow maldistribution and excessive hotspot temperature when the gap between two heater rods is too narrow. Wire-wraps are the spacer devices widely used in the sodium cooled fast reactor system to maintain the design fuel pitch. Generally, a wire spacer can enhance the turbulent heat transfer, but it increases the pressure drop, which is not desirable for a natural-circulation system like PCAT. Therefore, three short wire-wrap spacers per heater have been proposed to contain the core pressure drops and keep the heater-to-heater gap in an acceptable range. The objectives of this ECAR is to evaluate the PCAR core pressure drop with short wire spacers to confirm that the additional pressure drop due to the wire spacer will not disturb the natural circulation of the PCAT.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2023

Hydrodynamic spin-orbit coupling in asynchronous optically driven micro-rotors

Abstract Vortical flows of rotating particles describe interactions ranging from molecular machines to atmospheric dynamics. Yet to date, direct observation of the hydrodynamic coupling between artificial micro-rotors has been restricted by the details of the chosen drive, either through synchronization (using external magnetic fields) or confinement (using optical tweezers). Here we present a new active system that illuminates the interplay of rotation and translation in free rotors. We develop a non-tweezing circularly polarized beam that simultaneously rotates hundreds of silica-coated birefringent colloids. The particles rotate asynchronously in the optical torque field while freely diffusing in the plane. We observe that neighboring particles orbit each other with an angular velocity that depends on their spins. We derive an analytical model in the Stokes limit for pairs of spheres that quantitatively explains the observed dynamics. We then find that the geometrical nature of the low Reynolds fluid flow results in a universal hydrodynamic spin-orbit coupling. Our findings are of significance for the understanding and development of far-from-equilibrium materials.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Extending TOUGH + HYDRATE with a parallel particle transport simulator: numerical investigation of sand production during gas production from hydrate deposits

A new parallel code for simulating particle transport in porous media is integrated with the TOUGH + HYDRATE simulator to investigate sand production associated with gas production from unconsolidated gas hydrate-bearing sediments (HBS). Here, the parallel coupled simulator is named THMPT and uses the integral finite difference method to describe the Darcian and non-Darcian flow of fluids and heat transport, the finite element method to describe the associated geomechanical changes, and the discrete element method to track the trajectory of individual sand particles within the HBS. The THMPT simulator is written in Fortran, incorporates multiple optimized algorithms, and can comprehensively address the coupled flow, thermal, chemical, geomechanical, and particle transport processes that characterize the system behaviors during gas production from HBS. The simulator can capture all processes involved in sand particle transport in porous media, including sand detachment, collision, clogging (i.e., bridging), and migration. A benchmark case study of sand production in the course of depressurization-induced gas production from a representative HBS reveals various distinct microscopic particle migration mechanisms and the adverse impact of sand particle detachment, transport, and clogging. The numerical investigation also examines the effect of bottomhole pressure on mitigating sand production. The simulation results indicate that sand clogging near the wellbore significantly reduces permeability, decreasing gas production by at least 50%. Lastly, the efficiency of gravel packing in mitigating sand production is numerically evaluated, revealing that the structure of the porous media appears to profoundly influence the macroscopic motion behavior of sand particles and sand clogging characteristics.

discrete element method↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2023

Probing turbulence with lasers and coherent light sources

Fluid and plasma turbulence consists of multi-scale interactions where the energy of the large scales is dissipated through a cascade to the smallest scales set by the microphysics. Turbulence and mixing can be detrimental to system performance, such as in inertial confinement fusion (ICF). Our current understanding of turbulence is mostly derived from studies of incompressible flows, where the intermediate-scaled eddies in the “inertial range” bridge the transfer of energy from the large eddies to the dissipative range, and that there is no communication between eddies of disparate scales. Whether a similar cascade process exists for compressible turbulence, however, remains to be seen; experiments have been limited by diagnostics and supersonic flow facilities. This project has focused on the development of techniques to measure the evolution of flows and microphysical transport properties using particle tracking velocimetry. This project has investigated three main areas: (1) particle-tracking experiments with x-ray radiography, (2) shock-particle interactions, (3) quantifying viscous effects on particles in shocked flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Open-Source Tidal Energy Converter (OSTEC) Testbed: Design Basis Report

This report describes the design basis and design details for an instrumented marine turbine system intended to serve as the DOE’s marine tidal turbine test bed for foundational open-source R&D and data generation to advance our understanding and to identify knowledge gaps on the techno-economic performance of tidal energy converters (TEC) under real tidal flow conditions and at sufficiently large scale to enable upscaling of fluid-structure-interactions and component and material load responses

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Selectable inlet for aircraft propulsion system heat exchanger

A propulsion system for an aircraft includes a propulsor rotor, an engine core, a housing, a flowpath and a heat exchange system. The engine core is configured to power rotation of the propulsor rotor. The housing includes an inner structure and an outer structure. The outer structure is radially outboard of and axially overlaps the inner structure. The flowpath is downstream of the propulsor rotor and radially between the inner structure and the outer structure. The heat exchange system includes a duct, an exterior orifice, an interior orifice, a flow regulator and a heat exchanger disposed within the duct. The exterior orifice is formed by the outer structure and disposed along an exterior of the housing. The interior orifice is formed by the outer structure and disposed along the flowpath. The flow regulator is configured to selectively fluidly couple the exterior orifice and the interior orifice to the duct.

Sobanski, Jon E.↗
DOE OSTI2023

Immersive Particle Advection: Through the Scales of Renewable Energy: Preprint

We describe the benefits of immersive flow analysis for three large-scale computational science studies in the field of renewable energy. The studies encompass a range of scales, spanning from the large atmospheric scale of a wind farm to the human scale of an electric vehicle cabin down to the microscopic scale of battery material science. In these studies, users explored the flow patterns and dynamics through immersive particle advection. The integration of high-performance computing with immersive analysis provided a deeper understanding of these systems, helping develop more effective solutions for a sustainable energy future.

computational fluid dynamics↗
DOE OSTI2023

System and method for integrated absorber intercooling via rich solvent stream using 3D-printed active packings

The present disclosure relates to an absorber column apparatus for removing a selected component of a gas. The apparatus may have a first zone, a second zone and a third zone, wherein the first and third zones form a first domain through which a first fluid laden with a select gaseous component to be removed therefrom flows along concurrently with a second fluid. The second fluid at least substantially removes the select gaseous component from the first fluid to create a third fluid. The first fluid leaves the absorber column as a fourth fluid with the select gaseous component at least substantially removed therefrom. The second zone forms an active packing zone including a structure which forms an independent second domain in thermal communication with the first domain. The second receives a quantity of the third fluid and channels it through the second zone to help cool at least one of the first and second fluids.

Moore, Thomas Charles Reading↗
DOE OSTI2021

Clogging, dynamics, and reentrant fluid for active matter on periodic substrates

Here, we examine the collective states of run-and-tumble active matter disks driven over a periodic obstacle array. When the drive is applied along a symmetry direction of the array, we find a clog-free uniform liquid state for low activity, while at higher activity, the density becomes increasingly heterogeneous and an active clogged state emerges in which the mobility is strongly reduced. For driving along nonsymmetry or incommensurate directions, there are two different clogging behaviors consisting of a drive-dependent clogged state in the low activity thermal limit and a drive-independent clogged state at high activity. These regimes are separated by a uniform flowing liquid at intermediate activity. There is a critical activity level above which the thermal clogged state does not occur, as well as an optimal activity level that maximizes the disk mobility. Thermal clogged states are dependent on the driving direction while active clogged states are not. In the low activity regime, diluting the obstacles produces a monotonic increase in the mobility; however, for large activities, the mobility is more robust against obstacle dilution. We also examine the velocity-force curves for driving along nonsymmetry directions and find that they are linear when the activity is low or intermediate but become nonlinear at high activity and show behavior similar to that found for the plastic depinning of solids. At higher drives, the active clustering is lost. For low activity, we also find a reentrant fluid phase, where the system transitions from a high mobility fluid at low drives to a clogged state at higher drives and then back into another fluid phase at very high drives. We map the regions in which the thermally clogged, partially clogged, active uniform fluid, clustered fluid, active clogged, and directionally locked states occur as a function of disk density, drift force, and activity.

36 MATERIALS SCIENCE↗