Researcher profile

Silas Alben

Silas Alben contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
8topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2025arXiv

Enhancing wall-to-wall heat transport with unsteady flow perturbations

We determine unsteady flow perturbations that are optimal for enhancing the rate of heat transfer between hot and cold walls (i.e. the Nusselt number Nu), under the constraint of fixed flow power (Pe$^2$, where Pe is the Péclet number). The unsteady flows are perturbations of previously computed optimal steady flows and are given by eigenmodes of the Hessian matrix of Nu, the matrix of second derivatives with respect to amplitudes of flow mode coefficients. Positive eigenvalues of the Hessian correspond to increases in Nu by unsteady flows, and occur at Pe $\geq 10^{3.5}$ and within a band of flow periods $τ\sim$ Pe$^{-1}$. For $τ$Pe $\leq 10^{0.5}$, the optimal flows are chains of vortices that move along the walls or along eddies enclosed by flow branches near the walls. At larger $τ$Pe the vorticity distributions are often more complex and extend farther from the walls. The heat flux is enhanced at locations on the walls near the unsteady vorticity. We construct an iterative time-spectral solver for the unsteady temperature field and find increases in Nu of up to 7% at moderate-to-large perturbation amplitudes.

preprint2022arXiv

Efficient bending and lifting patterns in snake locomotion

We optimize three-dimensional snake kinematics for locomotor efficiency. We assume a general space-curve representation of the snake backbone with small-to-moderate lifting off the ground and negligible body inertia. The cost of locomotion includes work against friction and internal viscous dissipation. When restricted to planar kinematics, our population-based optimization method finds the same types of optima as a previous Newton-based method. A few types of optimal motions prevail. We find an s-shaped body with alternating lifting of the middle and ends for small-to-moderate transverse friction. For large transverse friction, curling and sliding motions are typical with small viscous dissipation, replaced by large-amplitude bending with large viscous dissipation. With small viscous dissipation we find local optima that resemble sidewinding motions across friction coefficient space. They are always suboptimal to alternating lifting motions, with average input power 10--100\% higher.

preprint2021arXiv

Packing of elastic rings with friction

We study the deformations of elastic filaments confined within slowly-shrinking circular boundaries, under contact forces with friction. We perform computations with a spring-lattice model that deforms like a thin inextensible filament of uniform bending stiffness. Early in the deformation, two lobes of the filament make contact. If the friction coefficient is small enough, one lobe slides inside the other; otherwise, the lobes move together or one lobe bifurcates the other. There follows a sequence of deformations that is a mixture of spiraling and bifurcations, primarily the former with small friction and the latter with large friction. With zero friction, a simple model predicts that the maximum curvature and the total elastic energy scale as the wall radius to the -3/2 and -2 powers respectively. With nonzero friction, the elastic energy follows a similar scaling but with a prefactor up to 8 times larger, due to delayering and bending with a range of small curvatures. For friction coefficients as large as 1, the deformations are qualitatively similar with and without friction at the outer wall. Above 1, the wall friction case becomes dominated by buckling near the wall.

preprint2020arXiv

Collective locomotion of two-dimensional lattices of flapping plates

We study the propulsive properties of rectangular and rhombic lattices of flapping plates at O(10--100) Reynolds numbers in incompressible flow. We vary five parameters: flapping amplitude, frequency (or Reynolds number), horizontal and vertical spacings between plates, and oncoming fluid stream velocity. Lattices that are closely spaced in the streamwise direction produce intense vortex dipoles between adjacent plates. The lattices transition sharply from drag- to thrust-producing as these dipoles switch from upstream to downstream orientations at critical flow speeds. Near these transitions the flows assume a variety of periodic and nonperiodic states, with and without up-down symmetry, and multiple stable self-propelled speeds can occur. As the streamwise spacing increases, the plates may shed typical vortex wakes that impinge on downstream neighbors. With small lateral spacing, rectangular lattices yield net drag, while rhombic lattices may generate net thrust, sometimes with high efficiency. As lateral spacing increases, rectangular lattices begin to generate thrust, eventually with slightly higher efficiencies than rhombic lattices, as the two types of lattice flows converge. At Re = 70, the maximum Froude efficiencies of time-periodic lattice flows are about twice those of an isolated plate. At lower Re, the lattices' efficiency advantage increases until the isolated flapping plate no longer generates thrust. The mean input power needed to generate the lattice flows can be estimated in the limits of small and large streamwise spacings, with small-gap and Poiseuille-like flows between the plates respectively in the two cases. For both lattices, the mean input power saturates as the lateral spacing becomes large (and thrust occurs). At small lateral spacings, the rhombic lattices' input power may be much larger when the plates overlap, leading to a decrease in Froude efficiency.

preprint2020arXiv

Intermittent sliding locomotion of a two-link body

We study the possibility of efficient intermittent locomotion for two-link bodies that slide by changing their interlink angle periodically in time. We find that the anisotropy ratio of the sliding friction coefficients is a key parameter, while solutions have a simple scaling dependence on the friction coefficients' magnitudes. With very anisotropic friction, efficient motions involve coasting in low-drag states, with rapid and asymmetric power and recovery strokes. As the anisotropy decreases, burst-and-coast motions change to motions with long power strokes and short recovery strokes, and roughly constant interlink angle velocity on each. These motions are seen in the spaces of sinusoidal and power-law motions described by two and five parameters, respectively. Allowing the duty cycle to vary greatly increases the motions' efficiency compared to the case of symmetric power and recovery strokes. Allowing further variations in the concavity of the power and recovery strokes improves the efficiency further only when friction is very anisotropic. Near isotropic friction, a variety of optimally efficient motions are found with more complex waveforms. Many of the optimal sinusoidal and power-law motions are similar to those that we find with an optimization search in the space of more general periodic functions (truncated Fourier series). When we increase the resistive force's power-law dependence on velocity, the optimal motions become smoother, slower, and less efficient, particularly near isotropic friction.

preprint2019arXiv

Large-amplitude membrane flutter in inviscid flow

We study the large-amplitude flutter of membranes (of zero bending rigidity) with vortex-sheet wakes in 2D inviscid fluid flows. We apply small initial deflections and track their exponential decay or growth and subsequent large-amplitude dynamics in the space of three dimensionless parameters: membrane pretension, mass density, and stretching modulus. With both ends fixed, all the membranes converge to steady deflected shapes with single humps that are nearly fore-aft symmetric, except when the deformations are unrealistically large. With leading edges fixed and trailing edges free, the membranes flutter with very small amplitudes and high spatial and temporal frequencies at small mass density. As mass density increases, the membranes transition to periodic and then increasingly aperiodic motions, and the amplitudes increase and spatial and temporal frequencies decrease. With both edges free, the membranes flutter similarly to the fixed-free case but also translate vertically with steady, periodic, or aperiodic trajectories, and with nonzero slopes that lead to small angles of attack with respect to the oncoming flow.

preprint2013arXiv

Optimizing snake locomotion in the plane. I. Computations

We develop a numerical scheme to determine which planar snake motions are optimal for locomotory efficiency, across a wide range of frictional parameter space. For a large coefficient of transverse friction, we find that retrograde traveling waves are optimal. The optimal snake deflection scales as the -1/4 power of the coefficient of transverse friction, in agreement with an asymptotic analysis. At the other extreme, zero coefficient of transverse friction, we propose a triangular direct wave which is optimal. Between these two extremes, a variety of complex, locally optimal motions are found. Some of these can be classified as standing waves (or ratcheting motions).

preprint2013arXiv

Optimizing snake locomotion in the plane. II. Large transverse friction

We determine analytically the form of optimal snake locomotion when the coefficient of transverse friction is large, the typical regime for biological and robotic snakes. We find that the optimal snake motion is a retrograde traveling wave, with a wave amplitude that decays as the -1/4 power of the coefficient of transverse friction. This result agrees well with our numerical computations.

preprint2012arXiv

Color scales that are effective in both color and grayscale

We consider the problem of finding a color scale which performs well when converted to a grayscale. We assume that each color is converted to a shade of gray with the same intensity as the color. We also assume that the color scales have a linear variation of intensity and hue, and find scales which maximize the average chroma (or "colorfulness") of the colors. We find two classes of solutions, which traverse the color wheel in opposite directions. The two classes of scales start with hues near cyan and red. The average chroma of the scales are 65-77% those of the pure colors.

preprint2012arXiv

Optimization of two- and three-link snake-like locomotion

We analyze two- and three-link planar snake-like locomotion and optimize the motion for efficiency. The locomoting system consists of two or three identical inextensible links connected via hinge joints, and the angles between the links are actuated as prescribed periodic functions of time. An essential feature of snake locomotion is frictional anisotropy: the forward, backward and transverse coefficients of friction are different. The dynamics are studied analytically and numerically for small and large amplitudes of the internal angles. Efficiency is defined as the ratio between distance traveled and the energy expended within one period, i.e. the inverse of the cost of locomotion. The optimal set of coefficients of friction to maximize efficiency consists of a large backward coefficient of friction and a small transverse coefficient of friction, compared to the forward coefficient of friction. For the two-link case with a symmetrical motion, efficiency is maximized when the internal angle amplitude is approximately $π/2$, for transverse coefficient sufficiently large. For the three-link case, the efficiency-maximizing paths are triangles in the parameter space of internal angles.

preprint2011arXiv

Bending of bilayers with general initial shapes

We present a simple discrete formula for the elastic energy of a bilayer. The formula is convenient for rapidly computing equilibrium configurations of actuated bilayers of general initial shapes. We use maps of principal curvatures and minimum-curvature direction fields to analyze configurations. We find good agreement between the computations and an approximate analytical solution for the case of a rectangular bilayer. For more general shapes (simple polyiamonds), we find a range of typical bending behaviors: overall bending directions along longest and shortest dimensions, inward bending at corners, curvature intensification near boundaries, and conical bending and partitioned bending zones in some cases.

preprint2010arXiv

Falling Flexible Sheets

We present a fluid dynamics video showing simulations of flexible bodies falling in an inviscid fluid. Vortex sheets are shed from the trailing edges of the bodies according to the Kutta condition. The basic behavior is a repeated series of accelerations to a critical speed at which the sheet buckles, and rapidly decelerates, shedding large vortices. Examples of persistent circling, quasi-periodic flapping, and more complex trajectories are shown.