Georgia Institute of Technology College of Sciences

Fluid membranes are area-preserving interfaces that resist bending. They are models of cell membranes, intracellular organelles, and viral particles. We are interested in developing simulation tools for dilute suspensions of deformable vesicles. These tools should be computationally efficient, that is, they should scale well as the number of vesicles increases. For very low Reynolds numbers, as it is often the case in mesoscopic length scales, the Stokes approximation can be used for the background fluid. We use a boundary integral formulation for the fluid that results in a set of nonlinear integro-differential equations for the vesicle dynamics. The motion of the vesicles is determined by balancing the nonlocal hydrodynamic forces with the elastic forces due to bending and tension. Numerical simulations of such vesicle motions are quite challenging. On one hand, explicit time-stepping schemes suffer from a severe stability constraint due to the stiffness related to high-order spatial derivatives and a milder constraint due to a transport-like stability condition. On the other hand, an implicit scheme can be expensive because it requires the solution of a set of nonlinear equations at each time step. We present two semi-implicit schemes that circumvent the severe stability constraints on the time step and whose computational cost per time step is comparable to that of an explicit scheme. We discretize the equations by using a spectral method in space, and a multistep third-order accurate scheme in time. We use the fast multipole method to efficiently compute vesicle-vesicle interaction forces in a suspension with a large number of vesicles. We report results from numerical experiments that demonstrate the convergence and algorithmic complexity properties of our scheme. Joint work with: Shravan K. Veerapaneni, Denis Gueyffier, and Denis Zorin.

We consider the three-dimensional gravity-capillary waves on water of finite-depth which are uniformly translating in a horizontal propagating direction and periodic in a transverse direction. The exact Euler equations are formulated as a spatial dynamical system in stead of using Hamiltonian formulation method. A center-manifold reduction technique and a normal form analysis are applied to show that the dynamical system can be reduced to a system of ordinary differential equations. Using the existence of a homoclinic orbit connecting to a two-dimensional periodic solution for the reduced system, it is shown that such a generalized solitary-wave solution persists for the original system by applying a perturbation method and adjusting some appropriate constants.

There has been growing interest in developing numerical methods for stochastic computations. This is motivated by the need to conduct uncertainty quantification in simulations, where uncertainty is ubiquitous and exists in parameter values, initial and boundary conditions, geometry, etc. In order to obtain simulation results with high fidelity, it is imperative to conduct stochastic computations to incorporate uncertainty from the beginning of the simulations. In this talk we review and discuss a class of fast numerical algorithms based on generalized polynomial chaos (gPC) expansion.The methods are highly efficient, compared to other traditional In addition to the forward stochastic problem solvers, we also discuss gPC-based methods for addressing "modeling uncertainty", i.e., deficiency in mathematical models, and solving inverse problems such as parameter estimation. ones, and suitable for stochastic simulations of complex systems.