An official website of the United States government
In the study of fluid dynamics and transport phenomena, key quantities of interest are often the force and torque on objects and total rate of heat/mass transfer from them. Conventionally, these integrated quantities are determined by first solving the governing equations for the detailed distribution of the field variables (i.e. velocity, pressure, temperature, concentration, etc.) and then integrating the variables or their derivatives on the surface of the objects. On the other hand, the divergence form of the conservation equations opens the door for establishing integral identities that can be used for directly calculating the integrated quantities without requiring the detailed knowledge of the distribution of the primary variables. This shortcut approach constitutes the idea of the reciprocal theorem, whose closest relative is Green’s second identity, which readers may recall from studies of partial differential equations. Despite its importance and practicality, the theorem may not be so familiar to many in the research community. Ironically, some believe that the extreme simplicity and generality of the theorem are responsible for suppressing its application! In this Perspectives piece, we provide a pedagogical introduction to the concept and application of the reciprocal theorem, with the hope of facilitating its use. Specifically, a brief history on the development of the theorem is given as a background, followed by the discussion of the main ideas in the context of elementary boundary-value problems. After that, we demonstrate how the reciprocal theorem can be utilized to solve fundamental problems in low-Reynolds-number hydrodynamics, aerodynamics, acoustics and heat/mass transfer, including convection. Throughout the article, we strive to make the materials accessible to early career researchers while keeping it interesting for more experienced scientists and engineers.
more »
« less
A reciprocal theorem for biphasic poro-viscoelastic materials
In studying the transport of inclusions in multiphase systems we are often interested in integrated quantities such as the net force and the net velocity of the inclusions. In the reciprocal theorem the known solution to the first and typically easier boundary value problem is used to compute the integrated quantities, such as the net force, in the second problem without the need to solve that problem. Here, we derive a reciprocal theorem for poro-viscoelastic (or biphasic) materials that are composed of a linear compressible solid phase, permeated by a viscous fluid. As an example, we analytically calculate the time-dependent net force on a rigid sphere in response to point forces applied to the elastic network and the Newtonian fluid phases of the biphasic material. We show that when the point force is applied to the fluid phase, the net force on the sphere evolves over time scales that are independent of the distance between the point force and the sphere; in comparison, when the point force is applied to the elastic phase, the time scale for force development increases quadratically with the distance, in line with the scaling of poroelastic relaxation time. Finally, we formulate and discuss how the reciprocal theorem can be applied to other areas, including (i) calculating the network slip on the sphere's surface, (ii) computing the leading-order effects of nonlinearities in the fluid and network forces and stresses, and (iii) calculating self-propulsion in biphasic systems.
more »
« less
- Award ID(s):
- 1944156
- PAR ID:
- 10593109
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 997
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The effects of viscoelasticity have been shown to manifest themselves via symmetry breaking. In this investigation, we show a novel phenomenon that arises from this idea. We observe that when a dense sphere is rotated near a wall (the rotation being aligned with the wall-normal direction and gravity), it levitates to a fixed distance away from the wall. Since the shear is larger in the gap (between the sphere and the wall) than in the open side of the sphere, the shear-induced elastic stresses are thus asymmetric, resulting in a net elastic vertical force that balances the weight of the sphere. We conduct experiments, theoretical models and numerical simulations for rotating spheres of various sizes and densities in a Boger-type fluid. In the small-Deborah-number range, the results are collapsed into a universal trend by considering a dimensionless group of the ratio of elastic to gravitational forces.more » « less
-
The cell cytoskeleton is a dynamic assembly of semi-flexible filaments and motor proteins. The cytoskeleton mechanics is a determining factor in many cellular processes, including cell division, cell motility and migration, mechanotransduction and intracellular transport. Mechanical properties of the cell, which are determined partly by its cytoskeleton, are also used as biomarkers for disease diagnosis and cell sorting. Experimental studies suggest that in whole cell scale, the cell cytoskeleton and its permeating cytosol may be modelled as a two-phase poro-viscoelastic (PVE) material composed of a viscoelastic (VE) network permeated by a viscous cytosol. We present the first general solution to this two-phase system in spherical coordinates, where we assume that both the fluid and network phases are in their linear response regime. Specifically, we use generalized linear incompressible and compressible VE constitutive equations to describe the stress in the fluid and network phases, respectively. We assume a constant permeability that couples the fluid and network displacements. We use these general solutions to study the motion of a rigid sphere moving under a constant force inside a two-phase system, composed of a linear elastic network and a Newtonian fluid. It is shown that the network compressibility introduces a slow relaxation of the sphere and non-monotonic network displacements with time along the direction of the applied force. Our results can be applied to particle-tracking microrheology to differentiate between PVE and VE materials, and to measure the fluid permeability as well as VE properties of the fluid and the network phases.more » « less
-
Curvature mediated elastic interactions between inclusions in lipid membranes have been analyzed using both theoretical and computational methods. Entropic corrections to these interactions have also been studied. Here we show that elastic and entropic forces between inclusions in membranes can compete under certain conditions to a yield a maximum in the free energy at a critical separation. If the distance between the inclusions is less than this critical separation then entropic interactions dominate and there is an attractive force between them, while if the distance is more than the critical separation then elastic interactions dominate and there is a repulsive force between them. We assume the inclusions to be rigid and use a previously developed semi-analytic method based on Gaussian integrals to compute the free energy of a membrane with inclusions.Weshow that the critical separation between inclusions decreases with increasing bending modulus and with increasing tension. We also compute the projected area of a membrane with rigid inclusions under tension and find that the trend of the effective bending modulus as a function of area fraction occupied by inclusions is in agreement with earlier results. Our technique can be extended to account for entropic effects in other methods which rely on quadratic energies to study the interactions of inclusions in membranes.more » « less
-
Curvature mediated elastic interactions between inclusions in lipid membranes have been analyzed using both theoretical and computational methods. Entropic corrections to these interactions have also been studied. Here we show that elastic and entropic forces between inclusions in membranes can compete under certain conditions to a yield a maximum in the free energy at a critical separation. If the distance between the inclusions is less than this critical separation then entropic interactions dominate and there is an attractive force between them, while if the distance is more than the critical separation then elastic interactions dominate and there is a repulsive force between them. We assume the inclusions to be rigid and use a previously developed semi-analytic method based on Gaussian integrals to compute the free energy of a membrane with inclusions.Weshow that the critical separation between inclusions decreases with increasing bending modulus and with increasing tension. We also compute the projected area of a membrane with rigid inclusions under tension and find that the trend of the effective bending modulus as a function of area fraction occupied by inclusions is in agreement with earlier results. Our technique can be extended to account for entropic effects in other methods which rely on quadratic energies to study the interactions of inclusions in membranes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2024.719
