BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ivan C. Christov

Ivan C. Christov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynmath.APcond-mat.softmath.DSnlin.PSmath-phmath.MPphysics.comp-phComputational Engineering, Finance, and Sciencecond-mat.mtrl-scihep-thmath.NAnlin.CDnlin.SIphysics.class-ph

Catalog footprint

What is connected

28works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

28 published item(s)

preprint2025arXiv

Oscillatory flows in three-dimensional deformable microchannels

Deformable microchannels emulate a key characteristic of soft biological systems and flexible engineering devices: the flow-induced deformation of the conduit due to slow viscous flow within. Elucidating the two-way coupling between oscillatory flow and deformation of a three-dimensional (3D) rectangular channel is crucial for designing lab- and organ-on-a-chip microsystems and eventually understanding flow-structure instabilities that can enhance mixing and transport. To this end, we determine the axial variations of the primary flow, pressure, and deformation for Newtonian fluids in the canonical geometry of a slender (long) and shallow (wide) 3D rectangular channel with a deformable top wall under the assumption of weak compliance and without restriction on the oscillation frequency (\textit{i.e.}, on the Womersley number). Unlike rigid conduits, the pressure distribution is not linear with the axial coordinate. To validate this prediction, we design a PDMS-based experimental platform with a speaker-based flow-generation apparatus and a pressure acquisition system with multiple ports along the axial length of the channel. The experimental measurements show good agreement with the predicted pressure profiles across a wide range of the key dimensionless quantities: the Womersley number, the compliance number, and the elastoviscous number. Finally, we explore how the nonlinear flow-deformation coupling leads to self-induced streaming (rectification of the oscillatory flow). Following Zhang and Rallabandi (\textit{J.\ Fluid Mech.}, vol.~996, 2024, A16), we develop a theory for the cycle-averaged pressure based on the primary problem's solution, and we validate the predictions for the axial distribution of the streaming pressure against the experimental measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hydrodynamic Bulge Testing: Materials Characterization without Measuring Deformation

Characterizing the elastic properties of soft materials through bulge testing relies on accurate measurement of deformation, which is experimentally challenging. To avoid measuring deformation, we propose a hydrodynamic bulge test for characterizing the material properties of thick, pre-stressed elastic sheets via their fluid--structure interaction with a steady viscous fluid flow. Specifically, the hydrodynamic bulge test relies on a pressure drop measurement across a rectangular microchannel with a deformable top wall. We develop a mathematical model using first-order shear-deformation theory of plates with stretching, and the lubrication approximation for Newtonian fluid flow. Specifically, a relationship is derived between the imposed flow rate and the total pressure drop. Then, this relationship is inverted numerically to yield estimates of the Young's modulus (given the Poisson ratio), if the pressure drop is measured (given the steady flow rate). Direct numerical simulations of two-way-coupled fluid--structure interaction are carried out in ANSYS to determine the cross-sectional membrane deformation and the hydrodynamic pressure distribution. Taking the simulations as ``ground truth,'' a hydrodynamic bulge test is performed using the simulation data to ascertain the accuracy and validity of the proposed methodology for estimating material properties. An error propagation analysis is performed via Monte Carlo simulation to characterize the susceptibility of the hydrodynamic bulge test estimates to noise. We find that, while a hydrodynamic bulge test is less accurate in characterizing material properties, it is less susceptible to noise, in the input (measured) variable, than a hydrostatic bulge test.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Long-time asymptotics of non-degenerate non-linear diffusion equations

We study the long-time asymptotics of prototypical non-linear diffusion equations. Specifically, we consider the case of a non-degenerate diffusivity function that is a (non-negative) polynomial of the dependent variable of the problem. We motivate these types of equations using Einstein's random walk paradigm, leading to a partial differential equation in non-divergence form. On the other hand, using conservation principles leads to a partial differential equation in divergence form. A transformation is derived to handle both cases. Then, a maximum principle (on both an unbounded and a bounded domain) is proved, in order to obtain bounds above and below for the time-evolution of the solutions to the non-linear diffusion problem. Specifically, these bounds are based on the fundamental solution of the linear problem (the so-called Aronson's Green function). Having thus sandwiched the long-time asymptotics of solutions to the non-linear problems between two fundamental solutions of the linear problem, we prove that, unlike the case of degenerate diffusion, a non-degenerate diffusion equation's solution converges onto the linear diffusion solution at long times. Select numerical examples support the mathematical theorems and illustrate the convergence process. Our results have implications on how to interpret asymptotic scalings of potentially anomalous diffusion processes (such as in the flow of particulate materials) that have been discussed in the applied physics literature.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revisiting steady viscous flow of a generalized Newtonian fluid through a slender elastic tube using shell theory

A flow vessel with an elastic wall can deform significantly due to viscous fluid flow within it, even at vanishing Reynolds number (no fluid inertia). Deformation leads to an enhancement of throughput due to the change in cross-sectional area. The latter gives rise to a non-constant pressure gradient in the flow-wise direction and, hence, to a nonlinear flow rate--pressure drop relation (unlike the Hagen--Poiseuille law for a rigid tube). Many biofluids are non-Newtonian, and are well approximated by generalized Newtonian (say, power-law) rheological models. Consequently, we analyze the problem of steady low Reynolds number flow of a generalized Newtonian fluid through a slender elastic tube by coupling fluid lubrication theory to a structural problem posed in terms of Donnell shell theory. A perturbative approach (in the slenderness parameter) yields analytical solutions for both the flow and the deformation. Using matched asymptotics, we obtain a uniformly valid solution for the tube's radial displacement, which features both a boundary layer and a corner layer caused by localized bending near the clamped ends. In doing so, we obtain a ``generalized Hagen--Poiseuille law'' for soft microtubes. We benchmark the mathematical predictions against three-dimensional two-way coupled direct numerical simulations (DNS) of flow and deformation performed using the commercial computational engineering platform by ANSYS. The simulations show good agreement and establish the range of validity of the theory. Finally, we discuss the implications of the theory on the problem of the flow-induced deformation of a blood vessel, which is featured in some textbooks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Static response of deformable microchannels: A comparative modelling study

We present a comparative modelling study of fluid-structure interactions in microchannels. Through a mathematical analysis based on plate theory and the lubrication approximation for low-Reynolds-number flow, we derive models for the flow rate-pressure drop relation for long shallow microchannels with both thin and thick deformable top walls. These relations are tested against full three-dimensional two-way-coupled fluid-structure interaction simulations. Three types of microchannels, representing different elasticity regimes and having been experimentally characterized previously, are chosen as benchmarks for our theory and simulations. Good agreement is found in most cases for the predicted, simulated and measured flow rate-pressure drop relationships. The numerical simulations performed allow us to also carefully examine the deformation profile of the top wall of the microchannel in any cross section, showing good agreement with the theory. Specifically, the prediction that span-wise displacement in a long shallow microchannel decouples from the flow-wise deformation is confirmed, and the predicted scaling of the maximum displacement with the hydrodynamic pressure and the various material and geometric parameters is validated.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-averaged transport in oscillatory squeeze flow of a viscoelastic fluid

Periodically-driven flows are known to generate non-zero, time-averaged fluxes of heat or solute species, due to the interactions of out-of-phase velocity and temperature/concentration fields, respectively. Herein, we investigate such transport (a form of the well-known Taylor--Aris dispersion) in the gap between two parallel plates, one of which oscillates vertically, generating a time-periodic squeeze flow of either a newtonian or Maxwellian fluid. Using the method of multiple time-scale homogenization, the mass/heat balance equation describing transport in this flow is reduced to a one-dimensional advection--diffusion--reaction equation. This result indicates three effective mechanisms in the mass/heat transfer in the system: an effective diffusion that spreads mass/heat along the concentration/temperature gradient, an effective advective flux, and an effective reaction that releases or absorbs mass/heat - in the time-averaged frame. Our results demonstrate that there exist resonant modes under which the velocity peaks when the dimensionless plate oscillation frequency (embodied by the Womersley number, the ratio of the transient inertia to viscous forces) approaches specific values. As a result, transport in this flow is significantly influenced by the dimensionless frequency. On the one hand, the effective, time-averaged dispersion coefficient is always larger than the molecular diffusivity, and is sharply enhanced near resonance. The interaction between fluid elasticity and the oscillatory forcing enhances the efficiency of transport in the system. On the other hand, the identified effective advection and reaction mechanisms may transport mass/heat from regions of high concentration/temperature to those of low concentration/temperature, or vice versa, depending on the value of dimensionless frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Tuning a magnetic field to generate spinning ferrofluid droplets with controllable speed via nonlinear periodic interfacial waves

Two dimensional free surface flows in Hele-Shaw configurations are a fertile ground for exploring nonlinear physics. Since Saffman and Taylor's work on linear instability of fluid--fluid interfaces, significant effort has been expended to determining the physics and forcing that set the linear growth rate. However, linear stability does not always imply nonlinear stability. We demonstrate how the combination of a radial and an azimuthal external magnetic field can manipulate the interfacial shape of a linearly unstable ferrofluid droplet in a Hele-Shaw configuration. We show that weakly nonlinear theory can be used to tune the initial unstable growth. Then, nonlinearity arrests the instability, and leads to a permanent deformed droplet shape. Specifically, we show that the deformed droplet can be set into motion with a predictable rotation speed, demonstrating nonlinear traveling waves on the fluid-fluid interface. The most linearly unstable wavenumber and the combined strength of the applied external magnetic fields determine the traveling wave shape, which can be asymmetric.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Unsteady fluid--structure interactions in a soft-walled microchannel: A one-dimensional lubrication model for finite Reynolds number

We develop a one-dimensional model for the unsteady fluid--structure interaction (FSI) between a soft-walled microchannel and viscous fluid flow within it. A beam equation, which accounts for both transverse bending rigidity and nonlinear axial tension, is coupled to a one-dimensional fluid model obtained from depth-averaging the two-dimensional incompressible Navier--Stokes equations across the channel height. Specifically, the Navier--Stokes equations are scaled in the viscous lubrication limit relevant to microfluidics. The resulting set of coupled nonlinear partial differential equations is solved numerically through a segregated approach employing fully-implicit time stepping. We explore both the static and dynamic FSI behavior of this example microchannel system by varying a reduced Reynolds number $Re$, which necessarily changes the Strouhal number $St$, while we keep the geometry and a modified dimensionless Young's modulus $Σ$ fixed. At steady state, an order-of-magnitude analysis (balancing argument) shows that the axially-averaged pressure in the flow, $\langle P\rangle$, exhibits two different scaling regimes, while the maximum deformation of the top wall of the channel, $H_{\mathrm{max}}$, can fall into four different regimes, depending on the magnitudes of $Re$ and $Σ$. These regimes are physically explained as resulting from the competition between the inertial and viscous forces in the fluid flow as well as the bending resistance and tension in the elastic wall. Finally, the linear stability of the steady inflated microchannel shape is assessed via a modal analysis, showing the existence of many highly oscillatory but stable modes, which further highlights the computational challenge of simulating unsteady FSIs.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Computational Analysis of Interfacial Dynamics in Angled Hele-Shaw Cells: Instability Regimes

We present a theoretical and numerical study on the (in)stability of the interface between two immiscible liquids, i.e., viscous fingering, in angled Hele-Shaw cells across a range of capillary numbers ($Ca$). We consider two types of angled Hele-Shaw cells: diverging cells with a positive depth gradient and converging cells with a negative depth gradient, and compare those against parallel cells without a depth gradient. A modified linear stability analysis is employed to derive an expression for the growth rate of perturbations on the interface and for the critical capillary number ($Ca_c$) for such tapered Hele-Shaw cells with small gap gradients. Based on this new expression for $Ca_c$, a three-regime theory is formulated to describe the interface (in)stability: (i) in Regime I, the growth rate is always negative, thus the interface is stable; (ii) in Regime II, the growth rate remains zero (parallel cells), changes from negative to positive (converging cells), or from positive to negative (diverging cells), thus the interface (in)stability possibly changes type at some location in the cell; (iii) in Regime III, the growth rate is always positive, thus the interface is unstable. We conduct three-dimensional direct numerical simulations of the full Navier--Stokes equations, using a phase field method to enforce surface tension at the interface, to verify the theory and explore the effect of depth gradient on the interface (in)stability. We demonstrate that the depth gradient has only a slight influence in Regime I, and its effect is most pronounced in Regime III. Finally, we provide a critical discussion of the stability diagram derived from theoretical considerations versus the one obtained from direct numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Solving Nonlinear Parabolic Equations by a Strongly Implicit Finite-Difference Scheme

We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error $\mathcal{O}\left[(Δx)^2 + (Δt)^2 \right]$. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape $h(x,t)$. To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear Acoustics and Shock Formation in Lossless Barotropic Green--Naghdi Fluids

The equations of motion of lossless compressible nonclassical fluids under the so-called Green--Naghdi theory are considered for two classes of barotropic fluids: (\textit{i}) perfect gases and (\textit{ii}) liquids obeying a quadratic equation of state. An exact reduction in terms of a scalar acoustic potential and the (scalar) thermal displacement is achieved. Properties and simplifications of these model nonlinear acoustic equations for unidirectional flows are noted. Specifically, the requirement that the governing system of equations for such flows remain hyperbolic is shown to lead to restrictions on the physical parameters and/or applicability of the model. A weakly nonlinear model is proposed on the basis of neglecting only terms proportional to the square of the Mach number in the governing equations, without any further approximation or modification of the nonlinear terms. Shock formation via acceleration wave blowup is studied numerically in a one-dimensional context using a high-resolution Godunov-type finite-volume scheme, thereby verifying prior analytical results on the blowup time and contrasting these results with the corresponding ones for classical (Euler) fluids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear waves in electromigration dispersion in a capillary

We construct exact solutions to an unusual nonlinear advection--diffusion equation arising in the study of Taylor--Aris (also known as shear) dispersion due to electroosmotic flow during electromigration in a capillary. An exact reduction to a Darboux equation is found under a traveling-wave anzats. The equilibria of this ordinary differential equation are analyzed, showing that their stability is determined solely by the (dimensionless) wave speed without regard to any (dimensionless) physical parameters. Integral curves, connecting the appropriate equilibria of the Darboux equation that governs traveling waves, are constructed, which in turn are shown to be asymmetric kink solutions ({\it i.e.}, non-Taylor shocks). Furthermore, it is shown that the governing Darboux equation exhibits bistability, which leads to two coexisting non-negative kink solutions for (dimensionless) wave speeds greater than unity. Finally, we give some remarks on other types of traveling-wave solutions and a discussion of some approximations of the governing partial differential equation of electromigration dispersion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stress retardation versus stress relaxation in linear viscoelasticity

We present a preliminary examination of a new approach to a long-standing problem in non-Newtonian fluid mechanics. First, we summarize how a general implicit functional relation between stress and rate of strain of a continuum with memory is reduced to the well-known linear differential constitutive relations that account for "relaxation" and "retardation." Then, we show that relaxation and retardation are asymptotically equivalent for small Deborah numbers, whence causal pure relaxation models necessarily correspond to ill-posed pure retardation models. We suggest that this dichotomy could be a possible way to reconcile the discrepancy between the theory of and certain experiments on viscoelastic liquids that are conjectured to exhibit only stress retardation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Parametric Study of Mixing in a Granular Flow a Bi-Axial Spherical Tumbler

We report on a computational parameter space study of mixing protocols for a half-full bi-axial spherical granular tumbler. The quality of mixing is quantified via the intensity of segregation (concentration variance) and computed as a function of three system parameters: angles of rotation about each tumbler axis and the flowing layer depth. Only the symmetric case is considered in which the flowing layer depth is the same for each rotation. We also consider the dependence on $\bar{R}$, which parametrizes the concentric spheroids ("shells") that comprise the volume of the tumbler. The intensity of segregation is computed over 100 periods of the mixing protocol for each choice of parameters. Each curve is classified via a time constant, $τ$, and an asymptotic mixing value, $bias$. We find that most choices of angles and most shells throughout the tumbler volume mix well, with mixing near the center of the tumbler being consistently faster (small $τ$) and more complete (small $bias$). We conclude with examples and discussion of the pathological mixing behaviors of the outliers in the so-called $τ$-$bias$ scatterplots.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comment on "The velocity field due to an oscillating plate in an Oldroyd-B fluid" by C.C. Hopkins and J.R. de Bruyn [Can. J. Phys. 92, 533 (2014)]

We correct certain errors and ambiguities in the recent pedagogical article by Hopkins and de Bruyn. The early-time asymptotics of the solution to the transient version of Stokes' second problem for an Oldroyd-B fluid in a half-space is presented, as an Appendix, to complement the late-time asymptotics given by Hopkins and de Bruyn.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On a hierarchy of nonlinearly dispersive generalized KdV equations

We propose a hierarchy of nonlinearly dispersive generalized Korteweg--de Vries (KdV) evolution equations based on a modification of the Lagrangian density whose induced action functional the KdV equation extremizes. It is shown that two recent nonlinear evolution equations describing wave propagation in certain generalized continua with an inherent material length scale are members of the proposed hierarchy. Like KdV, the equations from the proposed hierarchy possess Hamiltonian structure. Unlike KdV, however, the solutions to these equations can be compact (i.e., they vanish outside of some open interval) and, in addition, peaked. Implicit solutions for these peaked, compact traveling waves ("peakompactons") are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Study in Three-Dimensional Chaotic Dynamics: Granular Flow and Transport in a Bi-Axial Spherical Tumbler

We study 3D chaotic dynamics through an analysis of transport in a granular flow in a half-full spherical tumbler rotated sequentially about two orthogonal axes (a bi-axial "blinking" tumbler). The flow is essentially quasi-2D in any vertical slice of the sphere during rotation about a single axis, and we provide an explicit exact solution to the model in this case. Hence, the cross-sectional flow can be represented by a twist map, allowing us to express the 3D flow as a linked twist map (LTM). We prove that if the rates of rotation about each axis are equal, then (in the absence of stochasticity) particle trajectories are restricted to 2D surfaces consisting of a portion of a hemispherical shell closed by a "cap"; if the rotation rates are unequal, then particles can leave the surface they start on and traverse a volume of the tumbler. The period-one structures of the governing LTM are examined in detail: analytical expressions are provided for the location of period-one curves, their extent into the bulk of the granular material, and their dependence on the protocol parameters (rates and durations of rotations). Exploiting the restriction of trajectories to 2D surfaces in the case of equal rotation rates about the axes, a method is proposed for identifying and constructing 3D Kolmogorov-Arnold-Moser (KAM) tubes around the normally elliptic period-one curves. The invariant manifold structure arising from the normally hyperbolic period-one curves is also examined. When the motion is restricted to 2D surfaces, the structure of manifolds of the hyperbolic points in the bulk differs from that corresponding to hyperbolic points in the flowing layer. Each is reminiscent of a template provided by a non-integrable perturbation to a Hamiltonian system, though the governing LTM is not. This highlights the novel 3D chaotic behaviors observed in this model dynamical system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Shear dispersion in dense granular flows

We formulate and solve a model problem of dispersion of dense granular materials in rapid shear flow down an incline. The effective dispersivity of the depth-averaged concentration of the dispersing powder is shown to vary as the Péclet number squared, as in classical Taylor--Aris dispersion of molecular solutes. An extensions to generic shear profiles is presented, and possible applications to industrial and geological granular flows are noted.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Successive phase transitions and kink solutions in $ϕ^{8}$, $ϕ^{10}$, and $ϕ^{12}$ field theories

We obtain exact solutions for kinks in $ϕ^{8}$, $ϕ^{10}$ and $ϕ^{12}$ field theories with degenerate minima, which can describe a second-order phase transition followed by a first-order one, a succession of two first-order phase transitions and a second-order phase transition followed by two first-order phase transitions, respectively. Such phase transitions are known to occur in ferroelastic and ferroelectric crystals and in meson physics. In particular, we find that the higher-order field theories have kink solutions with algebraically-decaying tails and also asymmetric cases with mixed exponential-algebraic tail decay, unlike the lower-order $ϕ^4$ and $ϕ^6$ theories. Additionally, we construct distinct kinks with equal energies in all three field theories considered, and we show the co-existence of up to three distinct kinks (for a $ϕ^{12}$ potential with six degenerate minima). We also summarize phonon dispersion relations for these systems, showing that the higher-order field theories have specific cases in which only nonlinear phonons are allowed. For the $ϕ^{10}$ field theory, which is a quasi-exactly solvable (QES) model akin to $ϕ^6$, we are also able to obtain three analytical solutions for the classical free energy as well as the probability distribution function in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On a difficulty in the formulation of initial and boundary conditions for eigenfunction expansion solutions for the start-up of fluid flow

Most mathematics and engineering textbooks describe the process of "subtracting off" the steady state of a linear parabolic partial differential equation as a technique for obtaining a boundary-value problem with homogeneous boundary conditions that can be solved by separation of variables (i.e., eigenfunction expansions). While this method produces the correct solution for the start-up of the flow of, e.g., a Newtonian fluid between parallel plates, it can lead to erroneous solutions to the corresponding problem for a class of non-Newtonian fluids. We show that the reason for this is the non-rigorous enforcement of the start-up condition in the textbook approach, which leads to a violation of the principle of causality. Nevertheless, these boundary-value problems can be solved correctly using eigenfunction expansions, and we present the formulation that makes this possible (in essence, an application of Duhamel's principle). The solutions obtained by this new approach are shown to agree identically with those obtained by using the Laplace transform in time only, a technique that enforces the proper start-up condition implicitly (hence, the same error cannot be committed).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cutting and Shuffling a Line Segment: Mixing by Interval Exchange Transformations

We present a computational study of finite-time mixing of a line segment by cutting and shuffling. A family of one-dimensional interval exchange transformations is constructed as a model system in which to study these types of mixing processes. Illustrative examples of the mixing behaviors, including pathological cases that violate the assumptions of the known governing theorems and lead to poor mixing, are shown. Since the mathematical theory applies as the number of iterations of the map goes to infinity, we introduce practical measures of mixing (the percent unmixed and the number of intermaterial interfaces) that can be computed over given (finite) numbers of iterations. We find that good mixing can be achieved after a finite number of iterations of a one-dimensional cutting and shuffling map, even though such a map cannot be considered chaotic in the usual sense and/or it may not fulfill the conditions of the ergodic theorems for interval exchange transformations. Specifically, good shuffling can occur with only six or seven intervals of roughly the same length, as long as the rearrangement order is an irreducible permutation. This study has implications for a number of mixing processes in which discontinuities arise either by construction or due to the underlying physics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Wave Solutions

In classical continuum physics, a wave is a mechanical disturbance. Whether the disturbance is stationary or traveling and whether it is caused by the motion of atoms and molecules or the vibration of a lattice structure, a wave can be understood as a specific type of solution of an appropriate mathematical equation modeling the underlying physics. Typical models consist of partial differential equations that exhibit certain general properties, e.g., hyperbolicity. This, in turn, leads to the possibility of wave solutions. Various analytical techniques (integral transforms, complex variables, reduction to ordinary differential equations, etc.) are available to find wave solutions of linear partial differential equations. Furthermore, linear hyperbolic equations with higher-order derivatives provide the mathematical underpinning of the phenomenon of dispersion, i.e., the dependence of a wave's phase speed on its wavenumber. For systems of nonlinear first-order hyperbolic equations, there also exists a general theory for finding wave solutions. In addition, nonlinear parabolic partial differential equations are sometimes said to posses wave solutions, though they lack hyperbolicity, because it may be possible to find solutions that translate in space with time. Unfortunately, an all-encompassing methodology for solution of partial differential equations with any possible combination of nonlinearities does not exist. Thus, nonlinear wave solutions must be sought on a case-by-case basis depending on the governing equation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comments on: "Energetic balance for the Rayleigh--Stokes problem of an Oldroyd-B fluid" [Nonlinear Anal. RWA 12 (2011) 1]

We point out that an erroneous derivation in the recent paper [Fetecau et al., Nonlinear Anal. RWA 12 (2011) 1] yields a correct solution by accident. Additionally, a number of misrepresentations and inaccuracies in the latter recent paper are identified, corrected and/or clarified in this Comment. Finally, a listing of recent papers in this journal that make a mistake applying the Fourier sine transform, and thus present erroneous solutions, is given as an Appendix.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Comments on: "Starting solutions for some unsteady unidirectional flows of a second grade fluid," [Int. J. Eng. Sci. 43 (2005) 781]

A significant mathematical error is identified and corrected in a recent highly-cited paper on oscillatory flows of second-grade fluids [Fetecau & Fetecau (2005). Int. J. Eng. Sci., 43, 781--789]. The corrected solutions are shown to agree identically with numerical ones generated by a finite-difference scheme, while the original ones of Fetecau & Fetecau do not. A list of other recent papers in the literature that commit the error corrected in this Comment is compiled. Finally, a summary of related erroneous papers in this journal is presented as an Appendix.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stretching and folding versus cutting and shuffling: An illustrated perspective on mixing and deformations of continua

We compare and contrast two types of deformations inspired by mixing applications -- one from the mixing of fluids (stretching and folding), the other from the mixing of granular matter (cutting and shuffling). The connection between mechanics and dynamical systems is discussed in the context of the kinematics of deformation, emphasizing the equivalence between stretches and Lyapunov exponents. The stretching and folding motion exemplified by the baker's map is shown to give rise to a dynamical system with a positive Lyapunov exponent, the hallmark of chaotic mixing. On the other hand, cutting and shuffling does not stretch. When an interval exchange transformation is used as the basis for cutting and shuffling, we establish that all of the map's Lyapunov exponents are zero. Mixing, as quantified by the interfacial area per unit volume, is shown to be exponentially fast when there is stretching and folding, but linear when there is only cutting and shuffling. We also discuss how a simple computational approach can discern stretching in discrete data.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comment on "On a class of exact solutions of the equations of motion of a second grade fluid" by C. Fetecau and J. Zierep (Acta Mech. 150, 135-138, 2001)

In a 2001 article, Fetecau and Zierep [1] considered Stokes' first problem for a second grade (SG) fluid, unaware that it had already been solved (correctly) by Puri [2] in 1984. These authors used the Fourier sine transform to obtain, what they believed to be, the exact solution of the resulting initial-boundary-value problem (IBVP). Unfortunately, due to their incorrect computation of the distributional derivative of the Heaviside function, an omission of a critical term in the subsidiary equation occurred, an elementary mistake (gracefully explained in [3]) apparently not uncommon in the literature. Therefore, the solution given in [1] is, generally speaking, incorrect.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Hidden solitons in the Zabusky-Kruskal experiment: Analysis using the periodic, inverse scattering transform

Recent numerical work on the Zabusky--Kruskal experiment has revealed, amongst other things, the existence of hidden solitons in the wave profile. Here, using Osborne's nonlinear Fourier analysis, which is based on the periodic, inverse scattering transform, the hidden soliton hypothesis is corroborated, and the \emph{exact} number of solitons, their amplitudes and their reference level is computed. Other "less nonlinear" oscillation modes, which are not solitons, are also found to have nontrivial energy contributions over certain ranges of the dispersion parameter. In addition, the reference level is found to be a non-monotone function of the dispersion parameter. Finally, in the case of large dispersion, we show that the one-term nonlinear Fourier series yields a very accurate approximate solution in terms of Jacobian elliptic functions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stokes' first problem for some non-Newtonian fluids: Results and mistakes

The well-known problem of unidirectional plane flow of a fluid in a half-space due to the impulsive motion of the plate it rests upon is discussed in the context of the second-grade and the Oldroyd-B non-Newtonian fluids. The governing equations are derived from the conservation laws of mass and momentum and three correct known representations of their exact solutions given. Common mistakes made in the literature are identified. Simple numerical schemes that corroborate the analytical solutions are constructed.

0 citations0 reviewsNo code linkedOpen access