Monday, February 28, 2011 - 14:00 , Location: Skiles 005 , Fangxu Jing , USC Mechanical Engineering , Organizer:
Vortex dynamics and solid-fluid interaction are two of the most important and most studied topics in fluid dynamics for their relevance to a wide range of applications from geophysical flows to locomotion in moving fluids. In this talk, we investigate two problems in these two areas: Part I studies the viscous evolution of point vortex equilibria; Part II studies the effects of body elasticity on the passive stability of submerged bodies.In Part I, we describe the viscous evolution of point vortex configurations that, in the absence of viscosity, are in a state of fixed or relative equilibrium. In particular, we examine four cases, three of them correspond to relative equilibria in the inviscid point vortex model and one corresponds to a fixed equilibrium. Our goal is to elucidate some of the main transient dynamical features of the flow. Using a multi-Gaussian ``core growing" type of model, we show that all four configurations immediately begin to rotate unsteadily, while the shapes of vortex configurations remain unchanged. We then examine in detail the qualitative and quantitative evolution of the structures as they evolve, and for each case show the sequence of topological bifurcations that occur both in a fixed reference frame, and in an appropriately chosen rotating reference frame. Comparisons between the cases help to reveal different features of the viscous evolution for short and intermediate time ! scales of vortex structures. The dynamical evolution of passive particles in the viscously evolving flow associated with the initial fixed equilibrium is shown and interpreted in relation to the evolving streamline patterns. In Part II, we examine the effects of body geometry and elasticity on the passive stability of motion in a perfect fluid. Our main motivation is to understand the role of body elasticity on the stability of fish swimming. The fish is modeled as an articulated body made of multiple links (assumed to be identical ellipses in 2D or identical ellipsoids in 3D) interconnected by hinge joints. It can undergo shape changes by varying the relative angles between the links. Body elasticity is accounted for via the torsional springs at the joints. The unsteadiness of the flow is modeled using the added mass effect. Equations of motion for the body-fluid system are derived using Newtonian and Lagrangian approaches for both hydrodynamically decoupled and coupled models in 2D and 3D. We specifically examine the stability associated with a relative equilibrium of the equations, traditionally referred to as the ``coast motion" (proved to be unstable for a rigid elongated body model), and f! ound that body elasticity does stabilize the system. Stable regions are identified based on linear stability analysis in the parameter space spanned by aspect ratio (body geometry) and spring constants (muscle stiffness), and the findings based on the linear analysis are verified by direct numerical simulations of the nonlinear system.
Monday, February 21, 2011 - 14:00 , Location: Skiles 005 , Enkeleida Lushi , NYU Math Dept. , Organizer:
Micro-organisms are known to respond to certain dissolved chemicalsubstances in their environment by moving preferentially away or towardtheir source in a process called chemotaxis. We study such chemotacticresponses at the population level when the micro-swimmers arehydrodynamically coupled to each-other as well as the chemicalconcentration. We include a chemotactic bias based on the known bacteriarun-and-tumble phenomenon in a kinetic model of motile suspension dynamicsdeveloped recently to study hydrodynamic interactions. The chemicalsubstance can be produced or consumed by the swimmers themselves, as wellas be advected by the fluid flows created by their movement. The linearstability analysis of the system will be discussed, as well as the entropyanalysis. Nonlinear dynamics are investigated using numerical simulationin two dimensions of the full system of equations. We show examples ofaggregation in suspensions of pullers (front-actuated swimmers) anddiscuss how chemotaxis affects the mixing flows in suspensions of pushers(rear-actuated swimmers). Last, I will discuss recent work on numericalsimulations of discrete particle/swimmer suspensions that have achemotactic bias.
Friday, January 21, 2011 - 14:00 , Location: Skiles 006 , Yang Wang , Michigan State University, Department of Mathematics , Organizer: Haomin Zhou
The blind source separation (BSS) problem, also better known as the "cocktail party problem", is a well-known and challenging problem in mathematics and engineering. In this talk we discuss a novel time-frequency technique for the BSS problem. We also discuss a related problem in which foreground audio signal is mixed with strong background noise, and present techniques for suppress the background noise.
[Special Time] A mathematical model for bunching and meandering instabilities during epitaxial growth of a thin filmThursday, January 13, 2011 - 11:00 , Location: Skiles 006 , Michel Jabbour , University of Kentucky , Organizer: Sung Ha Kang
Recent experiments indicate that one- and two-dimensionalinstabilities, bunching and meandering, respectively, coexist duringepitaxial growth of a thin film in the step-flow regime. This is in contrastto the predictions of existing Burton–Cabrera–Frank (BCF) models. Indeed, inthe BCF framework, meandering is predicated on an Ehrlich–Schwoebel (ES)barrier whereas bunching requires an inverse ES effect. Hence, the twoinstabilities appear to be a priori mutually exclusive. In this talk, analternative theory is presented that resolves this apparent paradox. Itsmain ingredient is a generalized Gibbs–Thomson relation for the stepchemical potential resulting in jump conditions along the steps that coupleadatom diffusions on adjacent terraces. Specialization to periodic steptrains reveals a competition between the stabilizing ES kinetics and adestabilizing energetic correction that can lead to step collisions. Theaforementioned instabilities can therefore be understood in terms of thetendency of the crystal to lower, away from equilibrium and in the presenceof dissipation, its total free energy. The presentation will be self-contained and no a priori knowledge of theunderlying physics is needed.
Monday, December 6, 2010 - 13:00 , Location: Skiles 255 , Shawn Walker , LSU Mathematics Dept. , Organizer:
Locomotion at the micro-scale is important in biology and in industrialapplications such as targeted drug delivery and micro-fluidics. Wepresent results on the optimal shape of a rigid body locomoting in 3-DStokes flow. The actuation consists of applying a fixed moment andconstraining the body to only move along the moment axis; this models theeffect of an external magnetic torque on an object made of magneticallysusceptible material. The shape of the object is parametrized by a 3-Dcenterline with a given cross-sectional shape. No a priori assumption ismade on the centerline. We show there exists a minimizer to the infinitedimensional optimization problem in a suitable infinite class ofadmissible shapes. We develop a variational (constrained) descent methodwhich is well-posed for the continuous and discrete versions of theproblem. Sensitivities of the cost and constraints are computedvariationally via shape differential calculus. Computations areaccomplished by a boundary integral method to solve the Stokes equations,and a finite element method to obtain descent directions for theoptimization algorithm. We show examples of locomotor shapes with andwithout different fixed payload/cargo shapes.
Monday, November 29, 2010 - 13:00 , Location: Skiles 255 , Yunho Kim , University of California, Irvine , Organizer: Sung Ha Kang
It has been suggested by Y. Meyer and numerically confirmed by many othersthat dual spaces are good for texture recovery. Among the dual spaces, ourwork focuses on Sobolev spaces of negative differentiability to recovertexture from noisy blurred images. Such Sobolev spaces are good to modeloscillatory component, on the other hand, the spaces themselves hardlydistinguishes texture component from noise component because noise is alsoconsidered to be a highly oscillatory component. In this talk, in additionto oscillatory component recovery, we will further investigate aone-parameter family of Sobolev norms to achieve such a separation task.
A quantitative rigidity result for the cubic to tetragonal phase transition in the geometrically linear theory with interfacial energyMonday, November 22, 2010 - 13:00 , Location: Skiles 255 , Antonio Capella-Kort , Universidad Nacional Autónoma de México (UNAM) , Organizer:
We are interested in the cubic to tetragonal phase transition in a shape memory alloy. We consider geometrically linear elasticity. In this framework, Dolzmann and Mueller have shown the following rigidity result:The only stress-free configurations are (locally) twins (i.e. laminates of just two of the three Martensitic variants).However, configurations with arbitrarily small elastic energy are not necessarily close to these twins: The formation of microstructure allows to mix all three Martensitic variants at arbitrary volume fractions. We take an interfacial energy into account and establish a (local) lower bound on elastic + interfacial energy in terms of the Martensitic volume fractions. The model depends on a non-dimensional parameter that measures the strength of the interfacial energy. Our lower, ansatz-free bound has optimal scaling in this parameter. It is the scaling predicted by a reduced model introduced and analyzed by Kohn and Mueller with the purpose to describe the microstructure near an interface between Austenite and twinned Martensite. The optimal construction features branching of the Martensitic twins when approaching this interface.
A General Framework for a Class of First Order Primal Dual Algorithms for Convex Optimization in Imaging ScienceMonday, November 8, 2010 - 13:00 , Location: Skiles 002 , Ernie Esser , University of California, Irvine , Organizer: Sung Ha Kang
In this talk, based on joint work with Xiaoqun Zhang and Tony Chan, we showhow to generalize the primal dual hybrid gradient (PDHG) algorithm proposedby Zhu and Chan to a broader class of convex optimization problems. A mainfocus will also be to survey several closely related methods and explain theconnections to PDHG. We point out convergence results for some modifiedversions of PDHG that have similarly good empirical convergence rates fortotal variation (TV) minimization problems. We also show how to interpretPDHG applied to TV denoising as a projected averaged gradient method appliedto the dual functional. We present some numerical comparisons of thesealgorithms applied to TV denoising and discuss some novel applications suchas convexified multiphase segmentation.
Monday, October 25, 2010 - 13:00 , Location: Skiles 002 (Ground floor, entrance from Skiles courtyard) , Christopher Larsen , WPI , Organizer:
I will describe a sequence of models for predicting crack paths in brittlematerials, with each model based on some type of variational principleconcerning the energy. These models will cover the natural range ofstatics, quasi-statics, and dynamics. Some existence results will bedescribed, but the emphasis will be on deficiencies of the models and openquestions.
Monday, October 11, 2010 - 13:00 , Location: Skiles 002 , Thorsten Stoesser , Georgia Tech Civil Engineering , Organizer:
In this talk, results of high-resolution numerical simulations of some complex flows that are occurring in the area of hydraulic engineering will be presented. The method of large-eddy simulation is employed to study details of the flow over rough and porous channel beds, flow in an ozone contactor and the flow through idealized emergent vegetation. The main objective of the simulations is to gain insight into physical mechanisms at play. In particular, flow unsteadiness and coherent turbulence structures are important contributors to mass and momentum transfer in open channels. The performed large-eddy simulations allow revealing and quantifying these coherent structures.