SEARCH · Search NASA

Results for “Complex suspensions”

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.

SourceResearch connectionFrom yearThrough year

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

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Eulerian simulation of complex suspensions and biolocomotion in three dimensions

Significance Fluid–structure interactions are ubiquitous in many natural and man-made environments. They are difficult to study analytically, and therefore accurate and flexible computational methods are an indispensable tool in the field. Typically, fluids are simulated with a fixed background computational mesh, whereas a solid is simulated with a mesh that moves with it, making it challenging to couple the two. Here we develop a three-dimensional computational method where both fluid and solid can be represented on a fixed computational grid, which simplifies the coupling between the two phases considerably. Our method can simulate scenarios that are difficult to do with other methods, such as complex suspensions containing soft particles that are both heavier and lighter than the fluid.

3D fluid-structure interaction↗

DOE OSTI2019

Multiscale method based on coupled lattice‐Boltzmann and Langevin‐dynamics for direct simulation of nanoscale particle/polymer suspensions in complex flows

Summary A hybrid computational method coupling the lattice‐Boltzmann (LB) method and a Langevin‐dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long‐range many‐body hydrodynamic interactions (HIs). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two‐way coupled to the nonfluctuating LB fluid through a discrete LB forcing source distribution to capture the long‐range HI. To ensure intrinsically linear scalability with respect to the number of particles, a Eulerian‐host algorithm for short‐distance particle neighbor search and interaction is developed and embedded to LB‐LD framework. The validity and accuracy of the LB‐LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB‐LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale multicomponent particulate suspension systems such as complex blood suspensions.

Liu, Zixiang↗

DOE OSTI2020

Connecting particle interactions to agglomerate morphology and rheology of boehmite nanocrystal suspensions

Rheology imposes significant challenges on processing of complex suspensions such as nuclear waste slurries at the Hanford and Savannah River sites. Understanding rheology connecting to microstructures and underlying particle interactions in complex slurries is therefore important for both fundamental knowledge and practical applications. Here, we use suspensions of aluminum oxyhydroxide minerals in the form of boehmite as an analog of the radioactive waste slurry to gain physical insights on the correlation between particle interactions, microstructures, and slurry rheology. Specifically, we use a combination of Couette rheometry and small-angle scattering techniques (independently and simultaneously) to understand how the slurry microstructure changes under flow and how these structural changes manifest themselves in the bulk rheology of the suspensions. Our experiments show that the boehmite slurries are thixotropic, with the rheology and structure of the suspensions changing with increasing exposure to flow. In the slurries, particle aggregates begin as loose, system-spanning clusters, but exposure to moderate shear rates causes the aggregates to irreversibly consolidate into denser clusters of finite size. The microstructural changes directly influence the rheological properties of the slurries such as viscosity and viscoelasticity. More importantly, our study shows that solution pH affects the amount of structural rearrangement and the kinetics of the rearrangement process, with an increase in pH leading to faster and more dramatic changes in the bulk rheology. Such dynamic microstructural changes and resultant rheology are understood via correlations between particle interactions and strength of particle network, coupled with surface chemistry and anisotropic nature of particle interactions. Nearly identical structural changes are also observed in Poiseuille flow geometries, implying that the observed changes are relevant in the pipe flow conditions present during waste processing.

Weston, Javen S.↗

DOE OSTI2025

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Recovery of Oxpure 612C-50 Coconut Carbon Particles Using Bump Arrays

Cyanide is used to leach gold from crushed ore (solid matrix) into a gold cyanide solution. The gold is extracted from the cyanide solution by adsorption onto activated carbon. The gold extraction process occurs when the activated coconut carbon is placed into tanks that contain the gold cyanide solution either in a batch process or into a continuous flow circuit. The coconut carbon is then removed from the solution for gold recovery. The use of the mesofluidic separation system (a deterministic lateral displacement system) can significantly reduce costs and waste generated from the process of removing the coconut carbon from the solution. In the coconut carbon recovery process used at many gold mines in Nevada, the gold is attached to carbon particles that are still in the cyanide liquor. The mesofluidic system can rapidly remove the gold bearing carbon particles from the cyanide liquor quickly and with no operating costs. The benefits include: • Reduce operational costs related to filtration and hydro-cyclones • Reduced acid usage for the elution gold stripping phase as only 25% of the fluids will follow the carbon particles to the final express lane when two mesofluidic systems are used in series This alternative particle removal techniques would be advantageous. A promising technique for removing larger particles from slurries is mesofluidic separation, similar to deterministic lateral displacement arrays or “bump” arrays [1] but operates at much larger flow rates. As described by Pease, et al. [2], “Bump arrays in deterministic lateral displacement devices separate large particles from small particles using arrays of staggered posts. Large particles, defined as those with radii larger than the distance between the edge of a post and the stagnation streamline from the next downstream post, must bump toward one side of the device, whereas particles smaller than this distance slalom from entrance to exit without net lateral displacement.” Unlike filters or sieves, the posts in the arrays that cause separation do not block or occlude particle flow but work because particles go around the posts. Unlike hydrocyclones separation, separation is not driven by particle density, because gravitational forces are unimportant to mesofluidic separation. Burns, et al., and Pease, et al., have shown that particles may be separated from complex suspensions, under turbulent flow conditions, and at industrially important flow rates [3-5]. They have also recently shown that mesofluidic devices may be arranged in series to increase separation performance [6]. In this paper, we evaluate the mesofluidic system for the separation of commercial OxPure GR 612 Charred Coconut shell particles as a proof of principle test for the rapid separation using the 1500 micron cut mesofluidic separator. We first describe the experimental system and conditions. We then present the experimental results. [1] Huang, L.R., E.C. Cox, R.H. Austin, J.C. Sturm. 2004. Continuous particle separation through deterministic lateral displacement. Science, 304, 987–990. https://doi.org/10.1126/science.1094567. [2] Pease L.F., J.A. Bamberger, C.A. Burns, and M.J. Minette. 2021. Large Particle Separation from Non-Newtonian Slurries using Bump Arrays. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-65904, V003T08A023. New York, New York: ASME. doi:10.1115/FEDSM2021-65904 [3] Burns C.A., T.G. Veldman, J. Serkowski, R.C. Daniel, X.-Y. Yu, M.J. Minette, L.F. Pease. 2021. Mesofluidic separation versus dead-end filtration. Separation and Purification Technology, 254, 117256. https://doi.org/10.1016/j.seppur.2020.117256. [4] Pease L.F., J.E. Serkowski, T.G. Veldman, J. Williams, X.-Y. Yu, M.J. Minette, J.A. Bamberger, C.A. Burns. 2021. Can Bump Arrays Separate Particles from Turbulent Flows?. In Proceedings of the ASME 2021 Fluids Engineering Division Summer Meeting FEDSM2021-67696, V003T08A024. New York, New York: ASME. doi:10.11

mesofluidic separation, slurry, bump array↗

DOE OSTI2023

Rheological dynamics and structural characteristics of supramolecular assemblies of β-cyclodextrin and sulfonic surfactants

Cyclodextrins are highly functional compounds with a hydrophobic cavity capable of forming supramolecular inclusion complexes with various classes of molecules including surfactants. The resultant rich nanostructures and their dynamics are an interesting research problem in the area of soft condensed matter and related applications. In this paper, we report novel dynamical supramolecular assemblies based on the complexation of β-cyclodextrin with 3 different sulfonic surfactants, which are sodium hexadecylsulfate, sodium dodecylbenzenesulfonate, and myristyl sulfobetaine. It was observed that a β-cyclodextrin : surfactant/2 : 1 molar ratio was ideal for inducing axial growth and imparting large viscosities in the suspensions. Such complexation processes were accompanied by intriguing nanostructural phase behaviors and rheological properties that were very sensitive to the molecular architecture of sulfonic surfactants. The presence of an amino group in the head group of the surfactant allowed for large viscosities that reached 2.4 × 10 4 Pa s which exhibited gel-like behavior. In contrast, smaller viscosity values with a lower consistency index were observed when a bulky aromatic ring was present instead. DIC microscopy was used to visually probe the microstructure of the systems with respect to sulfonate molecular architecture. Additionally, surface tension measurements, and FTIR and NMR spectroscopies were used to gain insights into the nature of interactions that lead to the complexation and nanostructural characteristics. Finally, mechanics correlating the supramolecular morphologies to the rheological properties were proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Experimental evidence for the impact of soil viruses on carbon cycling during surface plant litter decomposition

To date, the potential impact of viral communities on biogeochemical cycles in soil has largely been inferred from correlational evidence, such as virus-driven changes in microbial abundances, viral auxiliary metabolic genes, and links with soil physiochemical properties. To more directly test the impact of soil viruses on carbon cycling during plant litter decomposition, we added concentrated viral community suspensions to complex litter decomposer communities in 40-day microcosm experiments. Microbial communities from two New Mexico alpine soils, Pajarito (PJ) and Santa Fe (SF), were inoculated onto grass litter on sand, and three treatments were applied in triplicate to each set of microcosms: addition of buffer (no added virus), live virus (+virus), or killed-virus (+killed-virus) fractions extracted from the same soil. Significant differences in respiration were observed between the +virus and +killed-virus treatments in the PJ, but not the SF microcosms. Bacterial and fungal community composition differed significantly by treatment in both PJ and SF microcosms. Combining data across both soils, viral addition altered links between bacterial and fungal diversity, dissolved organic carbon and total nitrogen. Overall, we demonstrate that increasing viral pressure in complex microbial communities can impact terrestrial biogeochemical cycling but is context-dependent.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Bridging microscopic dynamics and rheology in the yielding of charged colloidal suspensions

The yielding of soft materials is critical to many natural and industrial processes, yet experimental insights into microscopic aspects of yielding are limited. This study combines angle X-ray scattering, X-ray photon correlation spectroscopy, and in situ rheology (Rheo-SAXS-XPCS) with fast lubrication dynamics simulations to examine how interparticle interactions influence yielding in charged colloidal suspensions. By tuning attraction through salt addition, we compare repulsive and attractive systems under deformation. Repulsive suspensions yield uniformly with Andrade-like creep and minimal structural change. In contrast, attractive suspensions show complex behaviors, including shear banding, delayed yielding, and resolidification, governed by transient dynamics at shear band interfaces. These results directly link microscopic particle dynamics to macroscopic flow and demonstrate how interaction potentials control rheological behavior. This work offers a framework for designing soft materials with tailored properties for applications in coatings, food processing, drug delivery, and other technologies requiring precise mechanical control.

Molecular Dynamics Simulation↗

DOE OSTI2020

Onset of interfacial waves in the terahertz spectrum of a nanoparticle suspension

We used inelastic x-ray scattering to gain insight into the complex terahertz dynamics of a diluted Au-nanoparticle suspension in glycerol. Here, we observe that, albeit sparse, Au nanoparticles leave clear signatures on the dynamic response of the system, the main one being an additional mode propagating at the nanoparticle-glycerol interface. A Bayesian inferential analysis of the line shape reveals that such a mode, at variance with conventional acoustic modes, keeps a hydrodynamiclike behavior well beyond the continuous limit and down to subnanometer distances.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

The Effects of Temperature on the Assembly of Gold Nanoparticle by Interpolymer Complexation

Using synchrotron-based small-angle X-ray scattering techniques, we demonstrate that poly(ethylene glycol)-functionalized gold nanoparticles (PEG-AuNPs) are assembled into close-packed structures that include short-range order with face-centered cubic structure, where crystalline qualities are varied by controlling the electrolyte concentration, pH, and temperature of the suspensions. Here, we show that interpolymer complexation with poly(acrylic acid) (PAA) is induced by lowering the pH level of the PEG-AuNPs suspensions, and furthermore, increasing the temperature of the suspension strengthens interparticle attraction, leading to improved supercrystal structures. Our results indicate that this strategy creates robust nanoparticle superlattices with high thermal stability. The effects of PAA and PEG chain lengths on the assemblies are also investigated, and their optimal conditions for creating improved superlattices are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Manipulation of self-organized multi-vortical states in active magnetic roller suspensions

Ensembles of active magnetic colloids demonstrate complex collective behavior and self-organization when driven out of equilibrium by external magnetic fields. The interplay between magnetic and hydrodynamic interactions often result in spontaneous self-assembly and the emergence of localized vortices in ensembles of ferromagnetic rollers. We report experimental and computational study of the self-organized multi-vortical state in ensembles of magnetic rollers under the imposed external confining potential realized by patterned substrates. We explore the behavior of the system on a substrate patterned with shallow wells ordered in a lattice with square symmetry and lattice spacing smaller than the wells’ diameter. Here we demonstrate, that ensembles of active magnetic rollers evolve in response to changes in the geometry of the confining lattice from a globally correlated state characterized by anti-ferromagnetic vortex ordering to a state with rapidly changing particle flows and self-organized vortical states with a lattice constant larger than one of the confining pattern. Our experimental observations are accompanied by numerical simulations based on phenomenological coarse grained particle dynamics coupled to Navier–Stokes hydrodynamics in shallow water approximation. The reported results provide insights into the collective behavior of active magnetic rollers under confinement and suggest strategies for the manipulation of collective dynamic states in active magnetic liquids.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantifying electron transport in aggregated colloidal suspensions in the strong flow regime

Electron transport in complex fluids, biology, and soft matter is a valuable characteristic in processes ranging from redox reactions to electrochemical energy storage. These processes often employ conductor–insulator composites in which electron transport properties are fundamentally linked to the microstructure and dynamics of the conductive phase. While microstructure and dynamics are well recognized as key determinants of the electrical properties, a unified description of their effect has yet to be determined, especially under flowing conditions. In this work, the conductivity and shear viscosity are measured for conductive colloidal suspensions to build a unified description by exploiting both recent quantification of the effect of flow-induced dynamics on electron transport and well-established relationships between electrical properties, microstructure, and flow. These model suspensions consist of conductive carbon black (CB) particles dispersed in fluids of varying viscosities and dielectric constants. In a stable, well-characterized shear rate regime where all suspensions undergo self-similar agglomerate breakup, competing relationships between conductivity and shear rate were observed. To account for the role of variable agglomerate size, equivalent microstructural states were identified using a dimensionless fluid Mason number, Mn f , which allowed for isolation of the role of dynamics on the flow-induced electron transport rate. At equivalent microstructural states, shear-enhanced particle–particle collisions are found to dominate the electron transport rate. This work rationalizes seemingly contradictory experimental observations in literature concerning the shear-dependent electrical properties of CB suspensions and can be extended to other flowing composite systems.

Science & Technology - Other Topics↗

DOE OSTI2022

Spontaneous self-propulsion and nonequilibrium shape fluctuations of a droplet enclosing active particles

Active particles, such as swimming bacteria or self-propelled colloids, spontaneously assemble into large-scale dynamic structures. Geometric boundaries often enforce different spatio-temporal patterns compared to unconfined environment and thus provide a platform to control the behavior of active matter. Here, we report collective dynamics of active particles enclosed by soft, deformable boundary, that is responsive to the particles’ activity. We reveal that a quasi two-dimensional fluid droplet enclosing motile colloids powered by the Quincke effect (Quincke rollers) exhibits strong shape fluctuations with a power spectrum consistent with active fluctuations driven by particle-interface collisions. A broken detailed balance confirms the nonequilibrium nature of the shape dynamics. We further find that rollers self-organize into a single drop-spanning vortex, which can undergo a spontaneous symmetry breaking and vortex splitting. The droplet acquires motility while the vortex doublet exists. Our findings provide insights into the complex collective behavior of active colloidal suspensions in soft confinement.

42 ENGINEERING↗

DOE OSTI2019

Dispersion and properties of zirconia suspensions for stereolithography

Abstract Stereolithography is an attractive technique for the fabrication of complex‐shaped ceramic components with high dimensional accuracy. One of the challenges in this technology is the development of high solid loading, low viscosity photosensitive ceramic suspension. In this study, the dispersion of zirconia in photocurable resin and the slurry properties were intensively investigated. Rheological measurements showed that DISPERBYK‐103 proved to be an effective dispersant. 42 vol% ZrO 2 suspension was successfully prepared using 3.5 wt% DISPERBYK‐103 as the dispersant, with a suitable viscosity (4.88 Pa·s) below the maximum allowable viscosity value (5 Pa·s) for stereolithography applications. The adsorption behavior of DISPERBYK‐103 on the surface of zirconia powders was characterized by TG and FT‐IR, confirming the dispersion effect of dispersant. Contact angle measurements were also conducted to show that the adsorption of DISPERBYK‐103 could help to improve the wettability between powder and photocurable resin. Results showed that DISPERBYK‐103 was effective for the preparation of suitable slurries for the development of ZrO 2 ceramics through stereolithography.

Li, Xingbang↗

DOE OSTI2020

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Iron precipitation under controlled oxygen flow: Mineralogical implications for BIF precursors in the Archean ocean

Banded iron formation (BIF) deposition in the Archean ocean is thought to have been initiated by the oxidation of dissolved iron into ferric primary precursor phases that descended through the water column and were deposited on the seafloor. Effective interpretation of the trace element and isotopic composition of BIFs in the geologic record requires an understanding of the identity and longevity of those precursor phases. Temporal and spatial variation in oxidant concentration and ocean chemistry may have driven precursor mineralogy. Furthermore, precursor phases may undergo transitions during descent and burial. This experimental study follows the evolution of iron mineralogy, speciation, and redox state during the precipitation of iron under variable oxygen fluxes and in fluids containing different iron-complexing anions (chloride, sulfate, and phosphate). Additionally, suspensions collected after incomplete oxidation were anoxically aged to simulate changes to precursor mineralogy during descent below the seawater redoxcline. Results from X-ray diffraction, sequential dissolution, colorimetric determination of Fe redox state, and Fe K-edge X-ray absorption fine structure spectroscopy showed that under all experimental conditions, intermediate precipitates collected before complete iron uptake contained mineralized ferrous iron. Purely ferric mineral assemblages were not observed until iron was removed from solution. In low-oxygen experiments, intermediate precipitates contained phases with stoichiometric Fe 2+ , including magnetite, vivianite, and green rust (GR). Chloride green rust (GR1) had more ferric iron than sulfate green rust (GR2) and more readily transitioned to magnetite via disproportionation. GR2 was stable over a broader range of iron uptake; magnetite appeared mainly as an additional oxidation product in sulfate solution. In high oxygen experiments, mineralized ferrous iron occurred in amorphous phases or as a non-stoichiometric component in ferric oxides; final ferric assemblages were less crystalline than in low-oxygen equivalents. Low concentrations of phosphate increased the total iron oxidation rate, but also increased the ferrous content of intermediate precipitates. The incorporation of phosphate stabilized GR2 and facilitated GR2 precipitation over a larger range of iron uptake. In higher concentrations of phosphate, vivianite and poorly crystalline ferric hydroxides formed in place of GR. Here, these findings suggest that over a range of possible ocean chemistries, mixed-valence phases may have been short lived but relevant precursors to BIF. Additionally, anionic chemistry and oxidant concentration are shown to influence the crystallinity and chemical resistivity of final ferric assemblages, affecting their reactivity during later anoxic burial.

58 GEOSCIENCES↗

DOE OSTI2021

Temperature and salt controlled tuning of protein clusters

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (c p ) and the concentrations of the trivalent salts YCl 3 and LaCl 3 (c s ). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on c p , c s and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(c p , c s , T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Motility of acoustically powered micro-swimmers in a liquid crystalline environment

Suspensions of microswimmers in liquid crystals demonstrate remarkably complex dynamics and serve as a model system for studying active nematics. So far, experimental realization of microswimmers suspended in liquid crystalline media has relied on biological microorganisms that impose strict limitations on the compatible media and makes it difficult to regulate activity. Here, we demonstrate that acoustically powered bubble microswimmers can efficiently self-propel in a lyotropic liquid crystal. The velocity of the swimmers is controlled by the amplitude of the acoustic field. Histograms of swimming directions with respect to the local nematic field reveal a bimodal distribution: the swimmers tend to either fully align with or swim perpendicular to the director field of the liquid crystal, occasionally switching between these two states. The bubble-induced streaming from a swimmer locally melts the liquid crystal and produces topological defects at the tail of the swimmer. Further, we show that the defect proliferation rate increases with the angle between the swimmer's velocity and the local orientation of the director field.

36 MATERIALS SCIENCE↗