Hydrodynamic interaction between particles in suspension
Hydrodynamic interaction between particles in suspension
SEARCH · Search NASA
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.
Hydrodynamic interaction between particles in suspension
The structure and dynamics of confined suspensions of particles of arbitrary shape are of interest in multiple disciplines from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an immersed-boundary general geometry Ewald-like method to capture lubrication and long-range hydrodynamics and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction, and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than 10%, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion toward the walls, and a sub-diffusive regime-caused by crowding-in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.
To improve the efficiency of wave farms and achieve maximum power generation, the layout of wave energy converters (WECs) in an array needs to be carefully designed so that the hydrodynamic interactions can be positively exploited. For this, the hydrodynamic characteristics of the WEC array in different layouts need to be calculated. However, such calculations using numerical models usually entail significant computational cost, especially for large arrays of WECs. To address the computational challenge, a physics-constrained Gaussian process (GP) model is proposed to replace the original expensive numerical model and predict the hydrodynamic characteristics of the WECs for any array layout. By exploring the relationship between the WEC array (i.e., the input) and different hydrodynamic characteristics (i.e., the output), here we summarize a set of physical constraints/features, including invariance, symmetry, and additivity. This prior knowledge about the input-output relationship is then directly embedded in the constructed GP model through the design of physics-constrained kernels. In particular, a double-sum invariant kernel is first developed to incorporate the invariance and symmetry features, and then an additive kernel is developed to incorporate the additive feature of the problem. The invariant kernel and the additive kernel are then integrated to construct the physics-constrained GP model. Compared to the standard GP model, the proposed physics-constrained GP models require less training data to achieve the desired accuracy in predicting the hydrodynamic characteristics and are also less vulnerable to the curse of dimensionality (i.e., good scalability for large arrays) due to the use of an additive kernel. The efficiency, accuracy, and scalability of the proposed approach are demonstrated through an application to predict the hydrodynamic characteristics for WEC arrays of different sizes and layouts.
Despite decades of interdisciplinary research on topologically linked ring polymers, their dynamics remain largely unstudied. These systems represent a major scientific challenge as they are often subject to both topological and hydrodynamic interactions (HI), which render dynamical solutions either mathematically intractable or computationally prohibitive. Here we circumvent these limitations by preaveraging the HI of linked rings. We show that the symmetry of ring polymers leads to a hydrodynamic decoupling of ring dynamics. This decoupling is valid even for nonideal polymers and nonequilibrium conditions. Physically, our findings suggest that the effects of topology and HI are nearly independent and do not act cooperatively to influence polymer dynamics. We use this result to develop highly efficient Brownian dynamics algorithms that offer enormous performance improvements over conventional methods and apply these algorithms to simulate catenated ring polymers at equilibrium, confirming the independence of topological effects and HI. The methods developed here can be used to study and simulate large systems of linked rings with HI, including kinetoplast DNA, Olympic gels, and poly[n]catenanes.
The ROSAT HRI Data set on IC 4296 has been reduced and analysed. A draft paper on the small-scale structure of x-ray emission and its relation to the radio emission has been written. Mackie left the Smithsonian Astrophysical Observatory in September 1997 and in January 1998 he joined the staff at Carter Observatory, New Zealand. Mackie is currently (May 1998) upgrading computer software at Carter to run IRAF-PROS and XSPEC x-ray software packages in order to reduce and analyze archival ROSAT PSPC data on IC 4296. The PSPC results are needed to investigate the hot gas temperature and abundance properties that will be used in conjunction with the radio jet properties to determine the hydrodynamical interaction status of IC 4296, and finalize the results of a research paper.
The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.
Not Available
A binary liquid undergoes unidirectional solidification. The one-dimensional steady state is susceptible to morphological instability that causes the solid/liquid interface to change from a planar state to a cellular pattern. This paper examines the effects on this transition of volume-change convection, buoyancy-driven convection or forced flows. It emphasizes how flows alter stability limits, create scale and pattern changes in morphology, and create, through coupling, new instabilities. Emphasis is placed on the physical mechanisms of the interactions.
A typical dense non-Brownian particulate suspension exhibits shear thinning (decreasing viscosity) at low shear rate/stress followed by a Newtonian plateau (constant viscosity) at intermediate shear rate/stress values which transitions to shear thickening (increasing viscosity) beyond a critical shear rate/stress value and finally, undergoes a second shear thinning transition at extremely high shear rate/stress values. In this study, we unify and quantitatively reproduce all the disparate rate dependent regimes and the corresponding transitions for a dense non-Brownian suspension with increasing shear rate/stress. Here, we employ discrete particle dynamics simulations based on the proposed mechanism to elucidate its accuracy. We find that a competition between inter-particle interactions of hydrodynamic and non-hydrodynamic origins and the switching in the dominant stress scale with increasing the shear rate/stress lead to each of the above transitions. Inclusion of traditional hydrodynamic interactions, attractive/repulsive DLVO (Derjaguin and Landau, Verwey and Overbeek) interactions, the inter-particle contact interactions and a constant friction (or other constraint mechanism) reproduce the initial thinning as well as the shear thickening transition. However, to quantitatively capture the intermediate Newtonian plateau and the second shear thinning, an additional non-hydrodynamic interaction of non-DLVO origin and a decreasing coefficient of friction, respectively, are essential; thus, providing the first explanation for the presence of the intermediate Newtonian plateau along with reproducing the second shear thinning in a single framework. Expressions utilized for various interactions and friction are determined from experimental measurements and hence, result in an excellent quantitative agreement between the simulations and previous experiments.
In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the fluctuating immersed-boundary approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both algorithms the immersed-boundary method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long-range electrostatic interactions are computed by solving the Poisson equation, with short-range corrections included using an immersed-boundary variant of the classical particle-particle particle-mesh technique. Also included is a short-range repulsive force based on the Weeks-Chandler-Andersen potential. This methodology is validated by comparison to Debye-Hückel theory for ion-ion pair correlation functions, and Debye-Hückel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In each case, good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid.
Quantitative electron beam-target interaction studies on the Dual-Axis Radiographic Hydrodynamic Test (DARHT) electron linear induction accelerators have only been performed recently. This includes characterization of the temperature, density, pressure, extent, and expansion velocity of the plasma plume. The results presented in this dissertation present a detailed and unified overview of the accelerator systems, target heating physics, beam transport, and diagnostic tools. Additional information includes calibration sources, radiation hydrodynamics, spectroscopic-quality radiation transport modeling, experimental measurements, and analyses of electron beam driven aluminum experiments. The first set of spatially and temporally resolved spectroscopic measurements of electron beam driven aluminum are presented. Contamination quantification analyses are used to understand the origin of the strong Na-I 3p-3s lines that are observed in absorption within the aluminum plasma continuum. These results inform the creation of the first spectroscopic-quality radiation transport model that links several atomic physics codes to interpret the conditions from which the Na-I lines originate. A good agreement is found between the surface analysis results and the model which confirms the concentration of the sodium present within the aluminum alloy foil material. It also demonstrates, for the first time on electron beam-driven target experiments, the ability to interpret plasma conditions from measured absorption lines. In a second experimental campaign which focuses on pure aluminum, the Al-I 3p-4s and 3p3d doublets are both measured in emission. A detailed analysis of the Al-I 3p-4s doublet reveals that the lines undergo moderate self-absorption. A simple model of the self-absorption effect is successfully used to match measured spectra at various temperature/density/plasma scale length combinations. These measurements led to the realization of the minimum density that can be resolved by the spectrometers for the Al-I 3p-4s lines due to the large slit width required to observe a signal on aluminum. These measurements also demonstrate the sensitivity of visible and long wave UV spectroscopy to minor changes in both temperature and density. The simple self-absorption model will be useful for analysis of other beam-target interaction experiments with spectra exhibiting either self-absorption or full self-reversal. Substantial headway has been made on the modeling front by linking together several codes needed to model both the energy deposition, hydrodynamic motion, and atomic kinetics to produce synthetic spectral calculations that are compared with experimental measurements. There exists ample space for improvement, especially with benchmarking the hydrodynamics codes and equation-of-state tables with experimental measurements. The X-ray diagnostics required to make these new measurements along with the simulation capabilities required to interpret the results are under development and will be the subject of future studies.
Highly detailed HST observations of globular-cluster cores and galactic nuclei motivate new theoretical studies of the violent dynamical processes which govern the evolution of these very dense stellar systems. These processes include close stellar encounters and direct physical collisions between stars. Such hydrodynamic stellar interactions are thought to explain the large populations of blue stragglers, millisecond pulsars, X-ray binaries, and other peculiar sources observed in globular clusters. Three-dimensional hydrodynamics techniques now make it possible to perform realistic numerical simulations of these interactions. The results, when combined with those of N-body simulations of stellar dynamics, should provide for the first time a realistic description of dense star clusters. Here I review briefly current theoretical work on hydrodynamic stellar interactions, emphasizing its relevance to recent observations.
Experiments were performed in order to elucidate the effects of hydrodynamic interactions between two drops on their gravity-induced relative motion. The relative trajectories of two drops, their relative velocities, and the travel time for them to flow around each other were measured for different initial horizontal separations. Two size ratios and two viscosity ratios were investigated. Hydrodynamic interactions significantly reduce the relative velocity of two nearby drops and cause them to flow around each other with curved trajectories, resulting in a longer duration of the close encounter, compared with that for two non-interacting drops. These effects increase with decreasing drop separation, decreasing size ratio, and increasing viscosity ratio. Experimental results are in good agreement with theoretical predictions, except when the drops become sufficiently close that interface deformation occurs.
The main objectives of this project are to understand the differing particulate and multiphase flow behaviors that will occur in space and in Earth's gravity. More specifically, the project is concerned with understanding the effect of shear and gravity on two relatively ideal suspensions with significant inertial effects. The first is a gas-solid suspension at small Reynolds numbers and finite Stokes numbers. In this type of suspensions the inertia of the particle phase is significant while the hydrodynamic interactions are dominated by viscous forces in the suspending fluid. The other is a bubble suspension at small Weber and large Reynolds numbers. The hydrodynamic interactions in such suspensions are dominated by the inertial effects in the suspending fluid, but these inertial interactions can be described using potential flow theory. Our main objective is to examine the effects of shear and gravity on the average properties and stability of these two suspensions.
The presentation is an oral defense for the author's PhD.
Explore the source record for details and available documents.
For semidilute polyelectrolyte solutions, it is generally assumed that topological, electrostatic, and hydrodynamic interactions are screened (called triple screening). Despite a large body of research focused on polyelectrolyte solutions, the concept of triple screening has never been rigorously verified. In this work, we test the concept by probing concentration fluctuations in aqueous solutions containing a well-studied polyelectrolyte, sodium poly(styrenesulfonate) (NaPSS) with neutron scattering, theory, and molecular dynamics simulations. Neutron spin–echo (NSE) and small-angle neutron scattering (SANS) data from semidilute solutions of NaPSS are presented at different polymer and salt (NaCl) concentrations. A combined theory for structure (J. Chem. Phys. 105, 5183 (1996)) and dynamics (J. Chem. Phys. 107, 2619 (1997)), which captures effects of hydrodynamic, topological, and electrostatic screening, is used to interpret the experimental results. Here, the theory quantitatively predicts the decay rate obtained from the NSE measurements while capturing the shape and concentration dependencies of the polyelectrolyte peak observed in the SANS spectra. Detailed comparisons of the theory and the experiments reveal that the wavevector-dependent decay rate of concentration fluctuations in semidilute solutions of polyelectrolytes is dictated by the screening of hydrodynamic, topological, and electrostatic interactions. This conclusion is corroborated by coarse-grained molecular dynamics simulations, executed without any hydrodynamic interactions, which fail to capture the correct wavevector dependence of the decay rate. These results highlight that the theories based on the concept of triple screening provide a quantitative framework for predicting a relation between the structure and dynamics of polyelectrolyte solutions.