THREE-DIMENSIONAL SIMULATIONS OF TRANSIENT RESPONSE OF PEM FUEL CELLS

Transients have utmost importance in the lifetime and performance degradation of PEM fuel cells. Recent studies show that cyclic transients can induce hygro-thermal fatigue. In particular, the amount of water in the membrane varies significantly during transients, and determines the ionic conductivity and the structural properties of the membrane. In this work, we present three-dimensional time-dependent simulations and analysis of the transport in PEM fuel cells. U-sections of anode and cathode serpentine flow channels, anode and cathode gas diffusion layers, and the membrane sandwiched between them are modeled using incompressible Navier-Stokes equations in the gas flow channels, Maxwell-Stefan equations in the channels and gas diffusion layers, advection-diffusion-type equation for water transport in the membrane and Ohm’s law for ionic currents in the membrane and electric currents in gas diffusion electrodes. Transient responses to step changes in load, pressure and the relative humidity of the cathode are obtained from simulations, which are conducted by means of a third party finite-element package, COMSOL

Modeling and simulations of deformation and transport in PEM fuel cells

Performance degradation and durability of PEM fuel cells depend strongly upon transport and deformation characteristics of their components especially the polymer membrane. Physical properties of the membrane, such as its ionic conductivity and Young’s modulus depend on its water content, which varies significantly with operating conditions and during transients. Recent studies indicate that cyclic transients may induce hygro-thermal fatigue that leads to the ultimate failure of the membrane shortening its lifetime, and thus, hindering the reliable use PEM fuel cells for automotive applications. In this work, we present two-dimensional simulations and analysis of coupled deformation and transport in PEM fuel cells. A two-dimensional cross-section of anode and cathode gas diffusion layers, and the membrane sandwiched between them is modeled using Maxwell-Stefan equations in gas diffusion layers, Biot’s poroelasticity and Darcy’s law for deformation and water transport in the membrane and Ohm’s law for ionic currents in the membrane and electric currents in the gas diffusion layers. Steady-state deformation and transport of water in the membrane, transient responses to step changes in load and relative humidity of the anode and cathode are obtained from simulation experiments, which are conducted by means of a commercial finite-element package, COMSOL Multiphysics

Simulations of microflows induced by rotation of spirals in microchannels

In microflows where Reynolds number is much smaller than unity, screwing motion of spirals is an effective mechanism of actuation as proven by microorganisms which propel themselves with the rotation of their helical tails. The main focus of this study is to analyze the flow enabled by means of a rotating spiral inside a rectangular channel, and to identify effects of parameters that control the flow, namely, the frequency and amplitude of rotations and the axial span between the helical rounds, which is the wavelength. The time-dependent three-dimensional flow is modeled by Stokes equation subject to continuity in a time-dependent deforming domain due to the rotation of the spiral. Parametric results are compared with asymptotic results presented in literature to describe the flagellar motion of microorganisms

Validated reduced order models for simulating trajectories of bio-inspired artificial micro-swimmers

Autonomous micro-swimming robots can be utilized to perform specialized procedures such as in vitro or in vivo medical tasks as well as chemical surveillance or micro manipulation. Maneuverability of the robot is one of the requirements that ensure successful completion of its task. In micro fluidic environments, dynamic trajectories of active micro-swimming robots must be predicted reliably and the response of control inputs must be well-understood. In this work, a reduced-order model, which is based on the resistive force theory, is used to predict the transient, coupled rigid body dynamics and hydrodynamic behavior of bio-inspired artificial micro-swimmers. Conceptual design of the micro-swimmer is biologically inspired: it is composed of a body that carries a payload, control and actuation mechanisms, and a long flagellum either such as an inextensible whip like tail-actuator that deforms and propagates sinusoidal planar waves similar to spermatozoa, or of a rotating rigid helix similar to many bacteria, such as E. Coli. In the reduced-order model of the microswimmer, fluid’s resistance to the motion of the body and the tail are computed from resistive force theory, which breaks up the resistance coefficients to local normal and tangential components. Using rotational transformations between a fixed world frame, body frame and the local Frenet-Serret coordinates on the helical tail we obtain the full 6 degrees-of-freedom relationship between the resistive forces and torques and the linear and rotational motions of the swimmer. In the model, only the tail’s frequency (angular velocity for helical tail) is used as a control input in the dynamic equations of the micro-swimming robot. The reduced-order model is validated by means of direct observations of natural micro swimmers presented earlier in the literature and against; results show very good agreement. Three-dimensional, transient CFD simulations of a single degree of freedom swimmer is used to predict resistive force coefficients of a micro-swimmer with a spherical body and flexible tail actuator that uses traveling plane wave deformations for propulsion. Modified coefficients show a very good agreement between the predicted and actual time-dependent swimming speeds, as well as forces and torques along all axes

Numerical analysis of the 3D flow induced by propagation of plane-wave deformations on thin membranes inside microchannels

Propulsion mechanisms of microorganisms are based on either beating or screw-like motion of thin elastic biopolymers. Arguably, this motion is optimal for propulsion at very low Reynolds numbers. Similar actuation mechanisms can be utilized in the design of an autonomous microswimmer or even a micropump. In principle, propagation of plane-wave deformations on a thin-membrane placed inside a channel can lead to a net flow in the direction of the wave propagation. In this study we present effects of the amplitude, frequency, and the width of the membrane on the time-averaged flow rate and the rate of work done on the fluid by the membrane by means of threedimensional transient simulations of flows induced by plane-wave deformations on membranes. Navier- Stokes and continuity equations are used to model the flow on a time-varying domain, which is prescribed with respect to the motion of the membrane. Third party commercial software, COMSOL, is used in to solve the finite-element representation of the 3D time-dependent flow on moving mesh. Numerical simulations show that the flow inside the microchannel depend on the square of the amplitude and is proportional to the excitation frequency. Lastly, characteristic flow rate vs. pressure head curve and efficiency of a typical pump are obtained from 3D transient simulations, and presented here

Numerical simulations of a traveling plane-wave actuator for microfluidic applications

Continuous forming and propagation of large planar deformations on a thin solid elastic film can create propulsion when the film is immersed in a fluid. Microscopic organisms such as spermatozoa use similar mechanisms to propel themselves. In this work, we present a numerical analysis of the effect of traveling plane-wave deformations on an elastic-film actuator within a fluid medium inside a channel. In particular, we analyzed a micropump that consists of a wave actuator, which is placed in a channel to pump the fluid in the direction of the planedeformation waves. The unsteady flow over the moving boundary between the parallel plates has very low Reynolds number, and, hence, is modeled using the two-dimensional time-dependent Stokes equations. The fluid-structure interaction due to moving boundary is modeled with the arbitrary Lagrangian Eulerian (ALE) method incorporating the Winslow smoothing. COMSOL is used to solve two-dimensional timedependent Stokes equations on a deforming mesh, and to carry out simulations of the flow. Effects of the deformation amplitude, wavelength, frequency and channel height on the flow rate are presented

Experiment-based kinematic validation of numeric modeling and simulated control of an untethered biomimetic microrobot in channel

Modeling and control of swimming untethered microrobots are important for future therapeutic medical applications. Bio-inspired propulsion methods emerge as realistic substitutes for hydrodynamic thrust generation in micro realm. Accurate modeling, power supply, and propulsion-means directly affect microrobot motility and maneuverability. In this work, motility of bacteria-like untethered helical microrobots in channels is modeled with the resistive force theory coupled with motor dynamics. Results are validated with private experiments conducted on cm-scale prototypes fully submerged in Si-oil filled glass channel. Li-Po battery is utilized as the onboard power supply. Helical tail rotation is triggered by an IR remote control. It is observed that time-averaged velocities calculated by the model agree well with experimental results. Finally, time-dependent performance of a hypothetical model-based position control scheme is simulated with upstream flow as disturbance