4 AN10 Abstracts
IC1 IC4 On Dispersive Equations and Their Importance in Kinematics and Numerical Algebraic Geometry Mathematics Kinematics underlies applications ranging from the design Dispersive equations, like the Schr¨odinger equation for ex- and control of mechanical devices, especially robots, to the ample, have been used to model several wave phenom- biomechanical modelling of human motion. The major- ena with the distinct property that if no boundary con- ity of kinematic problems can be formulated as systems ditions are imposed then in time the wave spreads out of polynomial equations to be solved and so fall within spatially. In the last fifteen years this field has seen an the domain of algebraic geometry. While symbolic meth- incredible amount of new and sophisticated results proved ods from computer algebra have a role to play, numerical with the aid of mathematics coming from different fields: methods such as polynomial continuation that make strong Fourier analysis, differential and symplectic geometry, an- use of algebraic-geometric properties offer advantages in alytic number theory, and now also probability and a bit of efficiency and parallelizability. Although these methods, dynamical systems. In this talk it is my intention to present collectively called Numerical Algebraic Geometry, are ap- few simple, but still representative examples in which one plicable wherever polynomials arise, e.g., chemistry, biol- can see how these different kinds of mathematics are used ogy, statistics, and economics, this talk will concentrate on in this context. applications in mechanical engineering. A brief review of the main algorithms of the field will indicate their broad Gigliola Staffilani applicability. MIT [email protected] Charles Wampler General Motors Research Laboratories, Enterprise Systems Lab IC2 30500 Mound Road, Warren, MI 48090-9055, USA On the Role of Error and Uncertainty in the Nu- [email protected] merical Simulation of Complex Fluid Flows
The failure of numerical simulation to predict physical real- IC5 ity is often a direct consequence of the compounding effects Cloaking and Transformation Optics of numerical error arising from finite-dimensional approx- imation and physical model uncertainty resulting from in- We describe recent theoretical and experimental progress exact knowledge and/or statistical representation. In this on making objects invisible to detection by electromagnetic topical lecture, we briefly review systematic theories for waves, acoustic waves and matter waves. For the case of quantifying numerical errors and restricted forms of model electromagnetic waves, Maxwell’s equations have transfor- uncertainty occurring in simulations of fluid flow. A goal mation laws that allow for design of electromagnetic ma- of this lecture is to elucidate both positive and negative terials that steer light around a hidden region, returning aspects of applying these theories to practical fluid flow it to its original path on the far side. Not only would ob- problems. Finite-element and finite-volume calculations of servers be unaware of the contents of the hidden region, subsonic and hypersonic fluid flow are presented to contrast they would not even be aware that something was being the differing roles of numerical error and model uncertainty hidden. The object, which would have no shadow, is said for these problems. to be cloaked. We recount some of the history of the sub- ject and discuss some of the mathematical issues involved. Timothy J. Barth NASA Ames Research Center [email protected] Gunther Uhlmann University of Washington Department of Mathematics IC3 [email protected] Communication Complexity of Algorithms Algorithms have two kinds of costs: arithmetic and com- IC6 munication, which means either moving data between lev- Algebraic Geometric Algorithms in Discrete Opti- els of a memory hierarchy (in the sequential case) or over a mization network connecting processors (in the parallel case). The costs of communication often dominate the cost of arith- It is common knowledge that the understanding of the ge- metic. Thus, minimizing communication is often highly de- ometry of convex bodies has helped speed up algorithms sirable. The main goal of this talk is to describe a general in discrete optimization. For example, cutting planes and method for deriving upper and lower bounds on the amount facet-description of polyhedra have been crucial in the suc- of data moved (also known as bandwidth) for a very gen- cess of branch-and-bound algorithms for mixed integer lin- eral class of algorithms, including most dense and sparse ear programming. Another example is how the ellipsoid linear algebra algorithms. The method involves combinato- method can be used to prove polynomiality results in com- rial analysis of computation graphs, in particular, of their binatorial optimization. For the future, the importance expansion properties. This is joint work with Grey Ballard, of algebra and geometry in optimization is even greater James Demmel, and Oded Schwartz. since applications now demand non-linearity constraints together with discrete variables. In the past 5 years two Olga Holtz beautiful algebraic geometric algorithms on polyhedra have University of California, Berkeley been used to prove unexpected new results on the compu- Technische Universitat Berlin tation of integer programs with non-linearly objective func- [email protected] tions. The first is Barvinok’s algorithm for polytopes, and the second is Graver’s bases method on polyhedral cones. I will describe these two algorithms and explain why we AN10 Abstracts 5
can now prove theorems that were beyond our reach be- understanding of the options available to reduce energy fore. I will also describe attempts to turn these two al- use in buildings and will highlight the role of mathematics gorithms into practical computation, not just theoretical and particularly computational science in delivering low results. This is a nice story collecting results contained energy buildings to the market in cost-effective ways. A in several papers, joint work with various subsets of the systems approach to designing and operating low energy following people: R. Hemmecke, M. Koeppe, S. Onn, U. buildings will be described and the opportunities to de- Rothblum, and R. Weismantel. velop and deploy methods for uncertainty quantification, dynamical systems and embedded systems will be high- Jesus De Loera lighted. University of California, Davis [email protected] Clas Jacobson United Technologies Research Center IC7 [email protected] Scalable Tensor Factorizations with Incomplete Data IC10 Incomplete data is ubiquitous in biomedical signal pro- Semidefinite Programming: Algorithms and Appli- cessing, network traffic analysis, bibliometrics, social net- cations work analysis, chemometrics, computer vision, communi- cation networks, etc. We explain how factor analysis can The last two decades have seen dramatic advances in the identify the underlying latent structure even when signif- theory and practice of semidefinite programming (SDP), icant portions of data are missing and the remainder is stimulated in part by the realization that SDP is an ex- contaminated with noise. Though much of what we will tremely powerful modeling tool. Applications of SDP in- say is also applicable to matrix factorizations, we focus clude signal processing, relaxations of combinatorial and on the CANDECOMP/PARAFAC (CP) tensor decompo- polynomial optimization problems, covariance matrix es- sition. We demonstrate that it is possible to factorize in- timation, and sensor network localization. The first part complete tensors that have an underlying low-rank struc- of the talk will describe several applications, such as ma- ture. Further, there are approaches that can scale to sparse trix completion, that have attracted recent interest in large-scale data, e.g., 1000 x 1000 x 1000 with 99.5% miss- large scale SDP. The second part will focus on algorithms ing data. Real-world applicability is demonstrated in two for SDP. First, we present interior-point methods (IPMs) examples: EEG (electroencephalogram) applications and for medium scale SDP, follow by inexact IPMs (with lin- network traffic data. ear systems solved by iterative solvers) for large scale SDP and discuss their inherent limitations. Finally, we Tamara G. Kolda present recent algorithmic advances for large scale SDP Sandia National Laboratories and demonstrate that semismooth Newton-CG augmented [email protected] Lagrangian methods can be very efficient. Kim-Chuan Toh IC8 National University of Singapore Factorization-based Sparse Solvers and Precondi- Department of Mathematics tioners [email protected]
Efficient solution of large-scale, ill-conditioned and highly- indefinite algebraic equations often relies on high quality IP0 preconditioners together with iterative solvers. Because W. T. and Idalia Reid Prize in Mathematics Lec- of their robustness, factorization-based algorithms often ture: William T. and Idalia Reid: His Mathematics play a significant role in developing scalable solvers. We and Her Mathematical Family present our recent work of using state-of-the-art sparse fac- torization techniques to build domain-decomposition type The W. T. and Idelia Reid Prize in Mathematics was es- direct/iterative hybrid solvers and efficient incomplete fac- tablished in 1993 by SIAM and is funded by an endowment torization preconditioners. In addition to algorithmic prin- from the late Mrs. Idalia Reid to honor the memory of her ciples, we also address many practical aspects that need to husband Dr. W. T. Reid and his love of mathematics. Over be taken into consideration in order to deliver high speed time we often forget the people whose names are attached and robustness to the users of today’s sophisticated high to the prizes that are periodically awarded for outstanding performance computers. research or other contributions to our profession. Since I am the first of Professor Reid’s students to be honored by Xiaoye S. Li this Prize, it is only fitting that I take this opportunity to Computational Research Division review some of Dr. Reid’s contributions to mathematics Lawrence Berkeley National Laboratory and Mrs. Reid’s support of his mathematical colleagues [email protected] and students.
John A. Burns IC9 Interdisciplinary Center for Applied Mathematics Energy Efficiency in the Built Environment: Sys- Virginia Tech tems Approaches and Mathematical Challenges [email protected] Buildings consume nearly 40% of the world’s energy, sig- nificantly more than either the transportation or industrial IP0 sectors. Any comprehensive plan to reduce atmospheric Julian Cole Prize Lecture: Mathematical Mod- carbon must include actions to reduce energy consumption in the building sector. This talk will focus on the current 6 AN10 Abstracts
elling of Tissue Growth models that allow us to see and hear music simultaneously.
Approaches to the continuum modelling of biological tis- sue growth will be described, with specific emphasis on (i) Dmitri Tymoczko novel aspects of the resulting PDE formulations that arise Princeton University from the biological contexts to which they are intended [email protected] to apply and (ii) the role of singular-perturbation meth- ods in elucidating model properties. The talk will focus on recent developments in the application of multiphase IP1 continuum-mechanics descriptions, particularly in describ- Mathematical Challenges in Climate Change Sci- ing the growth of engineered tissue within a porous scaffold. ence John King Climate models solve the equations for the conservation University of Nottingham of momentum, mass and energy in the atmosphere and School of Mathematical Science oceans, the equations of state of air and for sea water, as [email protected];[email protected] well as equations for energy and water exchange with the land and cryosphere. This talk will present a brief review of the progress of climate modeling and will emphasis new IP0 challenges in projecting future climate change. Examples The AWM-SIAM Sonia Kovalevsky Lecture: Mix- will include prediction of co-evolution of atmospheric CO2 ing It Up: Discrete and Continuous Optimal Con- and climate, carbon data assimilation, and chaotic transi- trol for Biological Models tions. This presentation will illustrate optimal control methods Inez Fung applied to several types of models, including a mixture of University of California, Berkeley discrete and continuous features. The applications range [email protected] from a discrete model for cardiopulmonary resuscitation to partial differential equation models for rabies in raccoons. Detailed results will be given for harvesting in a PDE fish- IP2 ery model that answer the question: Does a marine reserve Compressed Sensing: Survey and Applications occur when maximizing harvest yield? Compressed sensing is a very active area of recent research Suzanne M. Lenhart attracting researchers from approximation theory, informa- University of Tennessee tion theory, harmonic analysis, statistics and signal pro- Department of Mathematics cessing. It has applications in areas ranging from speed- [email protected] ing up medical imaging to enabling more ambitious deep- sky astronomy. In addition the topic combines a bit of paradox – seeming to show that traditional ’barriers’ such IP0 as the Shannon-Nyquist sampling theorem don’t apply – The John von Neumann Lecture - Algebra: From while also involving some pretty advanced ideas from high- Linear to Non-Linear dimensional geometry and even statistical physics. In my talk I will try to engage newcomers with a survey some This lecture will discuss recent applications of methods of the excitement of this area, while also mentioning the from abstract algebra in modeling and solving non-linear sophisticated new ideas appearing almost daily. I will also problems across the mathematical sciences. Techniques try to bridge my talk to the associated minisymposia on that are familiar from linear algebra have natural exten- Compressed Sensing. sions in the non-linear world of algebraic geometry, whose recent advances are now transforming our thinking about David L. Donoho problems arising in domains such as convex optimization, Stanford University statistical inference and computational molecular biology. Department of Statistics [email protected]
Bernd Sturmfels University of California , Berkeley IP3 [email protected] Closing the Loop: Computational Physics and the Industrial Design Process IP0 Over the last two decades, improvements in computer hard- I.E. Block Community Lecture: The Geometry of ware and the development of nearly linear scaling methods Music have brought a variety of large-scale calculations within practical reach. While these methods have had substantial In my talk, I show how to translate the language of elemen- impact in industry, a number of improvements are needed tary music theory into that of contemporary geometry. It in order for these methods to be fully integrated into acous- turns out that concepts such as ”chord” and ”chord type” tic, electromagnetic, thermal and mechanical design envi- are naturally represented by points in geometrical spaces ronments. In this talk, we will describe the current state known as ”orbifolds.” Understanding these spaces can help of the art, as well as some recent ideas that may lead to a us to understand general constraints on musical style, as new generation of fast solvers, enabling robust, automatic well as the inner workings of specific pieces. For example, and geometrically flexible design by simulation. we will see that Mozart, Chopin, and Schubert made very sophisticated use of a necklace of four-dimensional cubes Leslie Greengard representing four-note chords. The talk will be accessible Courant Institute to non-musicians and will exploit interactive 3D computer New York University AN10 Abstracts 7
[email protected] sumes a narrow plume of hot gas rising and entraining the surrounding air. The surrounding air is assumed to have constant density and is irrotational far from the fire line. IP4 Then the flow outside the plume is described by a Biot- Transaction Costs Savart integral with jump conditions across the position of the plume. The plume model describes the unsteady Stochastic calculus models for financial markets have evolution of the mass, momentum and energy inside the played a fundamental role in finance in the past 30 years. plume, with sources derived to model mixing in the style The first application of these models is to determine the of Morton, et al. [1956]. The model is implemented in the price of a derivative security by working out the strategy same manner as structure-fluid interaction models [Alben by which the security payoff can be replicated by trading and Shelley 2008], except that we include a sink term in in primary assets. The second application is to solve the the Biot-Savart integral to couple the entrainment into the stochastic control problem of maximizing the utility of in- plume to the potential flow around the plume. The results vestment. If there are no transaction costs for trading in show that this model is capable of capturing a complicated assets, these problems have elegant solutions. This talk interaction of the plume with the surrounding air. Morton, is discusses what can be done when transaction costs are B. R., Taylor, Sir G. I., and Turner, J. S., Proc. Roy. Soc. present. The principal tool for analysis is a variational in- London, A 234, 1-23 (1956). Alben, S. and Shelley, M. J., equality version of the Hamilton-Jacobi-Bellman equation. Phys. Rev. Letters, 100, 074301, (2008). Although this equation does not have an analytical solu- tion, the initial terms in an asymptotic expansion of the Donald A. Drew solution can be obtained. Rensselaer Polytech Institute Dept of Mathematical Sciences Steven E. Shreve [email protected] Carnegie Mellon University Dept of Mathematical Sciences Kevin DelBene [email protected] Rensselaer Polytech Institute Department of Mathematical Sciences IP5 [email protected] Extension of Functions and Interpolation of Data CP1 We discuss the problem of finding a Cm extension of a given real-valued function defined on an arbitrary given subset E Asymptotic Structure of Diffusion Flames at High of Rn. If E is infinite, then we ask how to decide whether Pressure with Soret Transport an extension exists. If E is finite, then we ask for efficient Mathematical models of nonpremixed flames usually as- algorithms to compute an extension with close-to-minimal m sume Fick’s law applies to the molecular diffusivities. How- C norm. ever, at elevated pressures, it is known that the Soret ef- Charles Fefferman fects, i.e. mass transport via temperature gradients, be- Princeton University, USA come important. Moreover, at very high pressures, includ- [email protected] ing supercritical conditions, Soret effects may become dom- inant. We investigate the asymptotic structure of diffusion flames with Soret mass transport. Our theory provides ex- JP1 plicit expressions for flame temperature, location and ex- The Dynamics of Obesity tinction conditions in terms of physico-thermal parameters affecting combustion. The past two decades have seen a surge in the incidence of obesity in the developed world. Changes in body weight Daniel Fong,JohnK.Bechtold that can lead to obesity are known to result from imbal- New Jersey Institute of Technology ances between the energy derived from food and the en- [email protected], [email protected] ergy expended to maintain life and perform physical work. However, quantifying this relationship has proved difficult. CP1 Here, I will show how simple concepts from dynamical sys- tems can be used to provide a general quantitative descrip- Modes of Buckling of Viscous Fluid Threads tion of how body weight will change over time. The model The ‘coiling’ of a thin thread of viscous fluid falling onto a can then be used to answer open questions (and dispel surface is a familiar example of a buckling instability. Us- some myths) regarding weight loss and gain and provide ing an asymptotic ‘slender thread’ model embedded in a an explanation for the American obesity epidemic. numerical continuation procedure, we show that the onset Carson C. Chow of coiling can occur in three distinct modes involving dif- National Institutes of Health ferent balances of viscous forces, gravity, and inertia. We [email protected] confirm these numerical predictions using laboratory ex- periments, which also reveal a previously unobserved ‘ro- tatory folding’ mode of finite-amplitude buckling. Hydro- CP1 dynamic stability; turbulence (1704) Computational fluid Mathematical Model for the Behavior of Wildfires dynamics (1707) Incompressible fluids (1701)
Wildfires have been a long-standing problem in today’s so- Neil Ribe ciety. In this paper, we derive and solve a fluid dynam- Laboratoire FAST, France ics model to study a specific type of wildfire, namely, a [email protected] two dimensional flow around rising plume above a con- centrated heat source modeling a fire line. This flow as- 8 AN10 Abstracts
CP1 terials Bifurcation and Chaos in Rotating Driven Cavity Flow Ostwald ripening is the self-organization of components of two phase mixtures through a diffusive mechanism. This Numerical simulations of the incompressible Navier Stokes phenomenon has been modeled using stochastic partial equations in a two-dimensional driven cavity are performed differential equations derived from the underlying micro- using high-order numerical methods. The effect of solid physics. In this talk, results from simulations using spec- body rotation rate and Reynolds on flow stability is stud- tral schemes for stochastic partial differential equations are ied. The results show that the flow converged to a station- described and new results for improving the efficiency of ary state up to a certain rotation rate for a given Reynolds such schemes will be described. These simulation results number. Above this critical value periodic oscillations in- are compared with theoretical results such as the Lifshitz- dicative of Hopf bifurcation and eventually chaos are ob- Slyozov growth law. served and analyzed. David J. Horntrop Dinesh A. Shetty Dept of Mathematical Sciences, Center for Applied Math School of Mechanical Engineering New Jersey Institute of Technology Purdue University [email protected] [email protected] CP2 Steven Frankel School of Mechanical Engineering A Highly Effective Finite Difference Approach for Purdue University Oscillatory Optical Beam Equations [email protected] This talk discusses a class of efficient and effective fi- nite difference schemes for solving highly oscillated optical CP2 wave equations. To achieve the computational efficiency, mathematical transformations are used to map the com- Weak-Material Approximations of Natural Bound- plex model equations to coupled real differential equations. ary Conditions Then operator splitting strategies are introduced to de- The technique of Topology Optimization using the material compose the acquired nonlinear differential equations. An distribution approach is increasingly used in the design of Crank-Nicolson type discretization incorporates with par- advanced mechanical components. This approach approx- allel computation architectures is finally implemented un- imates natural boundary conditions at material interfaces der proper adaptations. The finite difference schemes con- by a fictitious-domain approach: a very weak material re- structed are numerically stable. Simulation experiments places regions of non-material (“void’). The energy-norm are given to illustrate our results, conclusions and expec- error associated with this approximation is shown to scale tations. linearly with the weak-material density and with the norm Qin Sheng of the boundary flux of a continuously extended solution Department of Mathematics to the original problem. Baylor University Martin Berggren Qin [email protected] Department of Computing Science, Umea University [email protected] Shekhar Guha U.S. Air Force Research Laboratory Fotios Kasolis [email protected] Department of Computing Science Umea University, Sweden Leonel Gonzalez [email protected] Air Force Research Laboratory Wright-Patterson AFB, OH 45433-7702 [email protected] CP2 The Effective Nonlinear Schrodinger Equation for Dark Spatial Plasmon-Polariton Solitons CP2 Shocks Versus Kinks in a Discrete Model of Dis- We derive an effective Nonlinear Schrodinger Equation for placive Phase Transitions propagation an optical dark spatial soliton at the subwave- length with surface plasmonic interaction. Starting with We consider dynamics of phase boundaries in a bistable Maxwell’s Equations we derive TM polarized type spatial one-dimensional lattice with harmonic long-range inter- solitons on a metal dielectric interface in which the dielec- actions. Using Fourier transform and Wiener-Hopf tech- tric is a Kerr medium that has self-defocusing. We numer- nique, we construct traveling wave solutions that represent ically and theoretically study the beam dynamics of this both subsonic phase boundaries (kinks) and intersonic ones nano-waveguide. (shocks). We derive the kinetic relation for kinks that pro- vides a needed closure for the continuum theory. We show Sihon H. Crutcher that the different structure of the roots of the dispersion U. S. Army, Research, Development,& Engineering Center relation in the case of shocks introduces an additional free Redstone Arsenal, Alabama parameter in these solutions, which thus do not require [email protected] a kinetic relation on the macroscopic level. The case of ferromagnetic second-neighbor interactions is analyzed in detail. We show that the model parameters have a signif- CP2 icant effect on the existence, structure and stability of the Stochastic Simulation of Self-Organization in Ma- traveling waves, as well as their behavior near the sonic AN10 Abstracts 9
limit. case there is measurement error, the Tikhonov regulariza- tion method is still convergent. The study is motivated by Evgeni Trofimov physical considerations. Department of Mathematics University of Pittburgh Xinfu Chen [email protected] Mathematics University of Pittsburgh Anna Vainchtein [email protected] Department of Mathematics University of Pittsburgh Yan-Hsiou Cheng [email protected] nsysu [email protected]
CP2 C.K. Law Structure Preserving Tau Methods for Stochastic NSYSU Chemical Systems [email protected] Tau leaping methods can efficiently simulate models of stiff stochastic chemical systems. However, most existing meth- CP3 ods do not naturally preserve some chemical structures, ASolutiontoAnAmbarzumyanProblemonTrees such as integer-valued and nonnegative molecular popula- tion states. In this talk, we present structure preserving We consider the Neumann Sturm-Liouville problem defined tau methods for simulating stochastic chemical systems. on trees such that the ratios of lengths of edges are not nec- We illustrate the new methods through a number of bio- essarily rational. It is shown that the potential function of chemically motivated examples, and provide comparisons the operator must be zero if the spectrum is equal to that with existing implicit tau methods. for zero potential. This extends a previous result which states that if the edgelengths of a tree are in rational ra- Yushu Yang tio and the spectrum contains {(nπ)2}, then the potential University of Maryland Baltimore County function is zero. Our result gives a complete solution to this [email protected] Ambarzumyan problem. For a proof, we compute approxi- mated eigenvalues for zero potential by using a generalized Muruhan Rathinam pigeon hole argument, and make use of a recursive formula University of Maryland, Baltimore County for characteristic functions. [email protected] Chun-Kong Law Jinglai Shen Dept of Applied Mathematics University of Maryland Baltimore County National Sun Yat-sen University [email protected] [email protected]
Eiji Yanagida CP3 Mathematical Institute Hopf Bifurcation of a Ratio-Dependent Predator Tohoku University Prey System with Time Delay [email protected]
In this study, we consider a predator-prey system with time delay where the predator dynamics is logistic with the car- CP3 rying capacity proportional to prey population. We study Existence of Solutions and Asymptotic Analysis of the impact of the time delay on the stability of the model a Singularly Perturbed BVP and by choosing the delay time t as a bifurcation parame- ter, we show that Hopf bifurcation can occur as the delay I will outline a proof of the existence of two solutions of time t passes some critical values. Using normal form the- a singularly perturbed BVP εy +2y + ey =0,y(0) = ory and central manifold argument, we also establish the 0,y(1) = 0 for 0 <ε<1/10. Both the solutions have direction and the stability of Hopf bifurcation. Finally, we the same outer solution on (0, 1], but different boundary perform numerical simulations to support our theoretical behavior near x = 0. One of them is bounded, the other results. unbounded as ε → 0. If time permits, I will briefly discuss a rigorous proof that the conjectured formal asymptotic Canan Celik expansion for the smaller solution is correct. Bahcesehir University, Turkey [email protected] Susmita Sadhu University Of Pittsburgh [email protected] CP3 Reconstructing Potentials from Zeros of One Eigenfunction CP3 Value Functions and Transversality Conditions for We study an inverse nodal problem, concerning the re- Infinite-Horizon Optimal Control Problems construction of a potential of a Sturm-Liouville operator by using zeros of one eigenfunction as input. We pro- The purpose of this paper is to establish the transversality pose three methods for the reconstruction, one of which condition at infinity as necessary and sufficient conditions is the Tikhonov regularization method. The explicit er- for optimality in optimal control problems with an infi- ror bounds are calculated for all the three methods. In nite horizon under nonsmoothness assumptions, dispens- 10 AN10 Abstracts
ing with any controllability condition on the velocity. We CP4 introduce the coercivity of the integrand, which plays a sig- A Mixed Implicit-Explicit Multirate Numerical nificant role in verifying the sufficiency of the transversality Scheme for Time-Dependent Equations condition at infinity together with convexity assumptions. Consequently, it is possible to derive a new characteriza- We develop a multirate time integration method for sys- tion of optimality in terms of the adjoint inclusions for tems of time-dependent equations that present two signif- the value function and the Hamiltonian, without using any icantly different scales within the model. We use an it- transversality condition at infinity. eration scheme to decouple the two time scales. At each iteration, we use an implicit Galerkin method to solve for Nobusumi Sagara the fast scale variable and an explicit method to solve for Faculty of Economics, Hosei University the slow variable. The error equation consists of a com- [email protected] putable leading order term and a provably higher order expression. CP3 Brandon Chabaud The Inverse Nodal Problem and Ambarzumyan Pennsylvania State University Problem for the p-Laplacian with Various Bound- [email protected] ary Conditions Qiang Du We study the issues of the reconstruction and stability Penn State University of the inverse nodal problem for the one-dimentional p- Department of Mathematics Laplacain eigenvalue problem. A key step is the appli- [email protected] cation of a modified Prufer substitution. Two associated Ambarzumyan problems are also solved. We consider the 1 case of C potentials with the Dirichlet boundary condi- CP4 tion and extend to the case of integrable potentials with Domain Decomposition Methods for Coupled the periodic or anti-periodic boundary conditions. Stokes-Darcy Model C.K. Law The Stokes-Darcy model arises in many applications such NSYSU as surface water flows, groundwater flows in karst aquifers, [email protected] and petroleum extraction. It has higher fidelity than either the Darcy or Stokes systems on their own. However, cou- W.C. Lian pling the two constituent models leads to a very complex NKMU system. This presentation discusses the domain decompo- [email protected] sition methods for solving the coupled Stokes-Darcy sys- tem. Convergence of these algorithms is demonstrated and Wei-Chuan Wang the results of computational experiments are presented to [email protected] illustrate their features. [email protected] Xiaoming He Y.H. Cheng Department of Scientific Computing NSYSU Florida State University [email protected] [email protected]
CP4 CP4 Effect Of Abrupt Expansion On Nanofluid Flow Solving High Speed Flow Problems with An Im- And Heat Transfer Characteristics: A Numerical proved Version of Ice. Study The Implicit Continuous-fluid Eulerian(ICE), a semi- Using the standard k-e turbulence model, incompressible implicit finite-volume solver, is used for simulating prob- nanofluid flow is simulated. Effect of alumina volume lems in multiphase flow which span a wide area of science concentration in water flow and heat transfer character- and engineering. ICE is utilized by the C-SAFE code Uin- istics downstream of an axisymmetric sudden expansion tah at the University of Utah to simulate explosions, fires is investigated. Assuming homogenous nanofluid, different and other fluid and fluid-structure interaction phenomena. alumina volume fractions in water, ranging from 0.01 to The implementation of ICE used in Uintah is given in many 0.05 are considered. The simulation revealed recirculation papers by Kashiwa at Los Alamos and extended to solve downstream of the expansion. Wall shear stress and heat multifield cases by Harman at Utah. In its original form transfer coefficient are shown to increase linearly with in- the ICE algorithm does not perform as well as the best creasing alumina volume fraction, up to 14% for alumina current methods for compressible flow problems. An im- volume fraction of 0.05. It is shown that the boundary proved version of ICE method using limiters is discussed layer reattachment length is insensitive to alumina volume and the obtained numerical results for several test prob- fraction under the assumptions invoked. lems are shown.
Khalid N. Alammar Lethuy T. Tran,MartinBerzins Department of Mechanical Engineering University of Utah King Saud University, Riyadh, SA [email protected], [email protected] [email protected] CP4 A Particle-in-Cell Method with Remapping and AN10 Abstracts 11
Local Mesh Refinement for Plasma Physics equations arising from motion within a sphere with a small gate for escape. Computational methods will be used to The Vlasov-Poisson equation describes the kinetic behavior compare via statistical analysis with the theoretical results of a collisionless plasma. We present an accurate and sta- to determine whether they are consistent, and also to ex- ble particle-in-cell based algorithm for solving the Vlasov- plore the limitations of the theory. Poisson equation. The method overcomes the numerical noise inherent in the usual particle-in-cell methods by pe- Carey Caginalp riodically remapping the distribution function on a hier- University of Pittsburgh archical of regularized grids. Positivity can be preserved [email protected] through the novel use of a limiter. The method has suc- cessfully been applied to a set of one-dimensional problems, e.g., strong Landau damping, two stream instability. Sec- CP5 ond order accuracy is observed over short times. Over long Predictability in Stochastic Dynamical Systems times, the test problems gradually evolve into a weakly collisional domain. There, filamentation in velocity space This talk presents sensitivity and predictability analysis of develops which leads to instability and loss of accuracy. stochastic reaction networks. Polynomial chaos expansions We solve this issue by introducing a collisional term to the are used to efficiently perform non-intrusive parametric un- Vlasov-Poisson equation. The collision equation is solved certainty propagation in the presence of intrinsic noise. on the remapping grids periodically, and coupled to the The effects of strong nonlinearities and bimodalities in the system using operator splitting. Second order accuracy is forward model are relieved with an adaptive, dimension- obtained for the coupled system. specific domain decomposition technique. The curse of di- mensionality is addressed using sparse quadrature methods Bei Wang for orthogonal projections. Department of Applied Science University of California at Davis Bert J. Debusschere [email protected] Energy Transportation Center Sandia National Laboratories, Livermore CA Greg Miller [email protected] University of California Department of Applied Science Khachik Sargsyan [email protected] Sandia National Laboratories [email protected] Phillip Colella Lawrence Berkeley National Laboratory Habib N. Najm [email protected] Sandia National Laboratories Livermore, CA, USA [email protected] CP4 A Lagrangian Vortex Method for Barotropic Vor- ticity Equation on a Rotating Sphere CP5 On the Maximum Entropy Solution to PDE-Based We present a Lagrangian vortex method for the barotropic Inverse Problems vorticity equation on a rotating sphere. The method tracks the flow map and absolute vorticity using Lagrangian par- We present a novel approach, based on the principle of ticles and panels. The velocity is computed from the Biot- maximum entropy, to statistical inverse problems con- Savart integral on the sphere. An adaptive refinement strained by partial differential equations (PDEs). Indi- strategy is implemented to resolve small-scale features and rectly observable parameters are estimated by a probabil- a treecode is used for efficient computation. Results are ity distribution that best represents our current state of presented for Rossby-Haurwitz waves and vortex interac- knowledge and does not impose additional synthetic con- tions. straints. Uncertainty is quantified naturally by information entropy, a by-product of our inversion method. Results are Lei Wang compared to the standard Bayesian approach for a thermal University of Michigan conduction problem. [email protected] Chad E. Lieberman John P. Boyd MIT University of Michigan [email protected] Ann Arbor, MI 48109-2143 USAUSA [email protected] Karen Willcox Massachusetts Institute of Technology Robert Krasny [email protected] University of Michigan Department of Mathematics CP5 [email protected] Multiscale Methods for Statistical Inference in El- liptic Problems CP5 Estimating the coefficients of an elliptic or parabolic PDE Asymptotics and Computation on Brownian Mo- from observations of the solution is typically an ill-posed tion Escape problem. Finite data resolution and the smoothing char- This talk will focus on asymptotics of partial differential acter of the forward operator limit one’s ability to recover 12 AN10 Abstracts
fine-scale information. We use the multiscale finite ele- models and numerical methods. We treat American op- ment method (MsFEM) to formulate a low-dimensional tion models as a special class of obstacle problems. Fi- and computationally efficient inference problem in this con- nite element formulation is introduced together with error text. A fully Bayesian treatment of the problem conditions analysis of numerical solutions. Some interesting arbitrage both small-scale and large-scale variability in the coeffi- properties about sensitivity of the option price to the pay- cients on available data. off function are proved. We also give a criterion for the convergence of numerical free boundaries (optimal exercise Youssef M. Marzouk, Matthew Parno boundaries) under mesh refinement. Some future research Massachusetts Institute of Technology plans will be discussed. [email protected], [email protected] Yongmin Zhang Xi’an Jiaotong - Liverpool University CP5 [email protected] A Multi-Scale Random Basis Method for Stochas- tic PDEs CP6 Uncertainty quantification is important in the field of engi- The Two-Machine Flowshop Scheduling Problem neering and applied mathematics. In this talk I will intro- to Minimize Maximum Lateness duce a multi-scale random basis method for elliptic PDE with random coefficients. We first derive the multi-scale This study addresses a two-machine flowshop scheduling random basis offline based on the Karhunen-Lo`eve decom- problem to minimize maximum lateness where processing position using Monte-Carlo simulations for one particular times are random variables with lower and upper bounds. non-zero force term. Then we use the multi-scale random Given that the problem is NP-hard, we propose several basis to solve the elliptic PDE for a general force term. algorithms and heuristics. The proposed algorithms and One important advantage of this method is that only very heuristics are compared through randomly generated data. small number of random bases is required to accurately The computational analysis has shown that one of the pro- represent the stochastic solution, thus providing significant posed heuristics performs very well with an overall average saving in the on-line computational effort compared with percentage error of less than one. existing methods. Numerical results will be provided to demonstrate the effectiveness of the method. Ali Allahverdi Kuwait University Mike Yan [email protected] Applied and computational mathematics California Institute of Technology Harun Aydilek [email protected] Gulf University for Science and Technology [email protected] Tom Hou California Institute of Technology [email protected] CP6 Matrix Splitting Method Combined with the Gra- Mulin Chen dient Projection Method for Support Vector Ma- Caltech chines [email protected] We propose the matrix splitting method combined with the gradient projection method can be applied for the CP5 quadratic programming problem with a single linear equal- A Stochastic Heterogeneous Multiscale Method for ity constraint and box constraints arising in training sup- Porous Media Flow port vector machines. The support vector machine prob- lem has a dense Hessian matrix. The matrix splitting al- A stochastic heterogeneous multiscale method is intro- gorithm transforms the problem into a sequence of approx- duced for modeling flows in heterogeneous media. To ac- imate subproblems which have a diagonal Hessian matrix count for the high stochastic dimensionality of the perme- and can be solved by the efficient pegging algorithm. The ability field, we employ the recently developed stochastic gradient projection method is combined with the matrix high dimensional model representation (HDMR) technique splitting method to reduce its sensitiveness. Some tech- for the solution of stochastic PDEs. HDMR decomposes niques for calculating the stepsize and non-monotone line the original high-dimensional problem into several lower search are applied for the proposed method to compare dimensional sub-problems which are efficiently solved us- their performances. In addition, we use the incomplete ing adaptive sparse grid collocation. The methodology is Cholesky decomposition method for large scale problems. demonstrated through a number of benchmark examples Our experimental results show the viability of the proposed of flows in porous media. algorithm.
Nicholas Zabaras,XiangMa Gitae Kim, Chih-Hang Wu Cornell University Kansas State University [email protected], [email protected] [email protected], [email protected]
CP5 CP6 American Option Pricing Models and Obstacle On the Extremal Problem of Polya Problems The notion of transfinite diameter of planar sets was in- We first give a brief overview of American option pricing troduced by M. Fekete around 1920s. This concept plays AN10 Abstracts 13
an important role in the classical complex analysis and is constrained problems. Compared to traditional SQP algo- related to other well-known concepts such as the logarith- rithms, this algorithm does not require the exact evaluation mic capacity and Chebyshev polynomials. The transfinite of constraint Jacobian in each optimization step but uses diameter of a compact set in the complex plain is the limit an approximation of the first-order derivative information. of n-diameters of the set. For each n ≥ 3, the n-diameter Hence the proposed algorithm is well suited to solve opti- dn(E)ofE is given by mization problems where the constraint Jacobian is dense i.e. the time required for the computation of the Jaco- 2 bian and its factorization dominates the overall optimiza- dn(E)=max /zi − zj / n(n−1) , (1) tion process. The quality of the approximated constraint 1≤i 1 Carnegie Mellon University dn(E) ≤ n n−1 (2) [email protected] and the equality holds for n-tuples of equally spaced points on the boundary of D. While investigating the transfinite Biegler Lorenz diameter of sets of constant width, Prof. Zair Ibragimov Carnegie Mellon University was led to the following weaker version of Polya’s problem: [email protected] among all n-tuples E = {z1,z2, ···,zn} with |zi − zj |≤2, find one with the largest n-diameter. He conjectured that Walther Andrea the extremal configuration will also be the vertices of a Institute of Mathematics regular n-gon, at least when n is odd. In this paper, I will University of Paderborn present my solution of Ibragimov’s problem for the case [email protected] n = 5 and approach for n = 7. I will also discuss why other cases can not be proved using the same method as for n =5. CP6 Image Space Analysis and Scalarization for Opti- Tuan N. Le mization of Multifunctions Fairmont High School [email protected] Using a new method based on generalized sections of fea- sible sets, we obtain optimality conditions for vector op- timization of objective multifunction with multifunction CP6 constraints. In particular, necessary and sufficient condi- A Binary Tree Based Search Algorithm for Global tions for scalarization of ε-optimization for multifunctions Optimization are deduced. We describe a binary tree based search method for solving Jafar Zafarani global optimization problems in the framework of a recently Dept of Maths. , Sheikhbahaee University developed memoryless interval-based optimization method. University of Isfahan, Isfahan, 81745-163, Iran It attempts to improve the reliability of the memoryless al- [email protected] gorithm while sacrificing little both in memory space usage and in overall speed of convergence. A binary tree data structure is introduced as a means of memory that would CP7 guide the algorithm to locate a new approximate solution, Mortar Mixed Finite Element Method for Curved while the overall algorithm employs the memoryless algo- Domains rithm iteratively. We present a mortaring (or a domain decomposition) Min Sun method for the Darcy problem that can handle curved do- Univ. of Alabama mains. The meshes are allowed to be non-matching over Math Dept the subdomain interfaces and the continuity is enforced [email protected] weakly via mortars, i.e. Lagrange multipliers. We show error estimates and confirm the results with numerical ex- Laura Ingram amples. The University of Alabama Mika Juntunen [email protected] University of Texas at Austin, ICES / CSM [email protected] CP6 An Inexact Trust-Region Based SQP Method Mary F. Wheeler for Nonlinear Programming Problems with Dense Center for Subsurface Modeling Constraint Jacobians University of Texas at Austin [email protected] We present a trust-region based SQP algorithm for the solution of optimization problems with general non-linear 14 AN10 Abstracts CP7 includedintwospacedimensions(d = 1 or 2). By applying Numerical Solution of Nonlinear Two-Dimensional the interpolated postprocessing technique, similar results Parabolic Partial Differential Equations by Branch- are also obtained on the error between the interpolation of ing Stochastic Processes the approximate solution and the exact solution. A new parallel numerical algorithm based on generating Kening Wang suitable random trees by Monte Carlo has been developed University of North Florida for solving nonlinear parabolic partial differential equa- [email protected] tions. While classical techniques based on a deterministic domain decomposition exhibits strong limitations, proba- Shuang Li bilistic methods are capable of exploit massively parallel Ernst & Young LLP architectures since the problem can be fully decoupled. [email protected] New examples of nonlinear equations have been run on a high performance supercomputer, showing a remarkable scalability and performance. CP7 Numerical Methods for Elliptic and Parabolic Angel Rodriguez-Rozas Equations in Perforated Domains Center for Mathematics and its Applications Department of Mathematics - Instituto Superior T´ecnico Let denote the size ratio of the holes of some perforated [email protected] domain to the whole domain. As closes to 0, a direct numerical simulation of the solutions of elliptic equations Juan Acebron in perforated domains can be very expensive. It is known Center for Mathematics and its Applications that the elliptic solutions approach a solution of some ho- Department of Mathematics - Instituto Superior Tecnico mogenized equation when is very small. It is reasonable [email protected] to expect that the numerical approximation of the solution of the homogenized equation is a good approximation for Renato Spigler the elliptic solutions in perforated domains when is small Dipartimento di Matematica, Universita Roma tre. enough. We use standard numerical methods to obtain the approximation for the homogenized elliptic solution and to Rome (Italy) ∞ 2 [email protected] derive the L estimate, L gradient estimate, and Lips- chitz estimate for the difference between the exact elliptic solutions and the numerical approximation of the homoge- CP7 nized solution. Counterpart results for parabolic equations Multiple Scattering from Complexly Shaped are also obtained. Coated Obstacles Li-Ming Yeh A numerical technique consisting of the Helmholtz equa- Department of Applied Mathematics tion in curvilinear coordinates coupled with an almost ex- National Chiao Tung University, Taiwan act Dirichlet-to-Neumann boundary condition is applied [email protected] to multiple scattering from coated obstacles of arbitrary shape. Novel elliptic grids conforming to complex geomet- CP7 rical configurations of several two-dimensional coated ob- stacles are constructed and approximations of the scattered New Immersed Finite Element Methods For Planar field supported by them are obtained. The numerical re- Elasticity Interface Problem sults illustrate the ability of the proposed technique to deal We will discuss immersed finite element (IFE) methods for with heterogeneous media. solving the planar elasticity interface problems with Carte- Vianey R. Villamizar sian meshes. The new IFE functions are formed on Carte- Department of Mathematics, Brigham Young University, sian meshes independent of the interface between different Provo, Utah 84602-6539 elastic materials. Some basic properties of the IFE func- [email protected] tions will be presented including the optimal approxima- tion capability of these IFE functions spaces. Furthermore, numerical examples indicate that these IFE methods con- Sebastian Acosta verge optimally in both L2 and semi-H 1 norms. Department of Mathematics, Brigham Young University [email protected] Xu Zhang Virginia Tech [email protected] CP7 Strong Superconvergence of Finite Element Meth- Tao Lin ods for Linear Parabolic Problems Department of Mathematics We study the strong superconvergence of a semi-discrete Virginia Tech finite element scheme for linear parabolic problems on Q = [email protected] Ω×(0,T], where Ω is a bounded domain in R(≤)with piecewise smooth boundary. We establish the global two CP8 order superconvergence results for the error between the Improved Bounds on Restricted Isometry Con- approximate solution and the Ritz projection of the exact 1,p stants for Gaussian Matrices solution of our model problem in W (Ω) and Lp(Q)with ≤ ∞ 2 p< and the almost two order superconvergence in The Restricted Isometry Constants (RIC) of a matrix A 1,∞ ∞ W (Ω) and L∞(Q). Results of the p = case are also measures how close to an isometry is the action of A on AN10 Abstracts 15 vectors with few nonzero entries, measured in the 2 norm. Networks Using Slime Mold Specifically, the upper and lower RIC of a matrix A of size n × N is the maximum and the minimum deviation from We present a new protocol for wireless sensor networks unity (one) of the largest and smallest, respectively, square (WSNs) based on models of network self-assembly in N Physarum polycephalum, a true slime mold. Slime mold of singular values of all k matrices formed by taking k columns from A. Calculation of the RIC is intractable assembles and modifies networks of tubes to distribute re- for most matrices due to its combinatorial nature; how- sources throughout its cell body. This system offers an ever, many random matrices typically have bounded RIC elegant solution to resource allocation problems similar to in some range of problem sizes (k, n, N). We provide the extracting data from WSNs. We present modeling, anal- best known bound on the RIC for Gaussian matrices, which ysis and simulation of this protocol for WSNs based on is also the smallest known bound on the RIC for any large models of Physarum polycephalum. rectangular matrix. Improvements over prior bounds are Louis F. Rossi achieved by exploiting similarity of singular values for ma- University of Delaware trices which share a substantial number of columns. [email protected] Bubacarr Bah, Jared Tanner University of Edinburgh Ke Li [email protected], [email protected] University of Delaware Computer and Information Sci. [email protected] CP8 A Mathematical Model of Communications Net- Claudio Torres works University of Delaware [email protected] The behavior of a communication network can be mod- eled at the node/link scale as traffic units flowing along Kyle Thomas the links connected by nodes. We present an approach to University of Delaware packet flow that derives the node/link network model and Chemical Engineering connects it to a fluid-like model of traffic flow. The discrete [email protected] node/link model describes packet queuing and the flow of packets from spatial point to spatial point. The model as- Chien-Chung Shen sumes that packets reside in buffers at each node, and are University of Delaware identified by their destination and the length of time they Computer and Information Sci. have resided in the buffer. An algorithm was created for [email protected] packets to exit the buffer at each node according to their age (“first in, first out”) and travel to the next node along a predetermined path to their destination. This algorithm CP8 calculates the rate at which packets distribute themselves An Iterative Thresholding Algorithm Using En- to the next link in the route to their destination, assumes hanced Sparsity a source of packets originating at the node, and subtracts packets arriving at their destination. A continuum model In this paper we propose a variation of the iterative soft- is derived from this discrete flow model by associating each thresholding algorithm that is used to find sparse (approx- node with a spatial area defined by Voronoi polygons. The imate) solutions to the equation Ax = b. In this variation, association with spatial areas leads to a flow continuity the sparsity of the iterate vector xn is used to redefine the equation. The continuity equation describes the density regularizing penalty function. Minimization of the newly of packets as a function of time and space, so that we are defined penalty function, then leads to the next iterate. able to predict changes in global flow patterns on a macro- Numerical data showing effectiveness of this approach will scopic scale due to aberrant flow behavior, such as outages be presented as well. or heavy traffic. Exact (non-classical) solutions are derived for one dimensional flows using the method of character- Hugo J. Woerdeman istics. These solutions show that if a combination of the Drexel University packet sources and network flow saturate the network flow hugo @ math.drexel.edu capacity, the packet density grows at the nearest upstream node. When the source strength is reduced, or when flow Ingrid Daubechies is restored, the buffered packets flow at capacity until the Princeton University density has been reduced. Multiple sources and destina- Department of Mathematics tions and multidimensional effects are discussed [email protected] Donald A. Drew Rensselaer Polytech Institute CP9 Dept of Mathematical Sciences A-Posteriori Estimator for Parabolic Coupled Sys- [email protected] tems of Pdes James R. Gatewood We present an extension of recent results by [Bernardi et al] Rensselear Polytechnic Institute and [Verfurth] on residual a-posteriori estimators for scalar [email protected] parabolic PDEs. Our results concern multilevel discretiza- tion of parabolic systems in which each component can be computed on a different grid. We analyze contribution of CP8 each term and show convergence results. The estimator is Modeling Adaptive Protocols for Wireless Sensor robust i.e. the efficiency index remains essentially constant 16 AN10 Abstracts when the coefficients of the system vary even by several or- refining a mesh locally. We present a residual-based a pos- ders of magnitude. teriori error estimate for the finite element discretization of elliptic equations using PHT-splines basis functions. Then Viviane Klein numerical experiments are presented to verify the theoret- Oregon State University ical results and demonstrate the robustness of the error [email protected] estimate. Malgorzata Peszynska Li Tian Department of Mathematics Penn State University Oregon State University [email protected] [email protected] Falai Chen Department of Mathematics CP9 University of Science and Technology of China Preconditioning for Implicit Ocean Models chenfl@ustc.edu.cn One challenge in ocean modeling is the spin-up problem, Qiang Du which requires integrations over centuries of time. Using Penn State University implicit methods to take large timesteps requires effective Department of Mathematics preconditioners. We examine thin stratified fluid problems [email protected] as modeled by the Parallel Ocean Program(POP), consid- ering multilevel preconditioners for the velocity-salinity- temperature subproblem and block preconditioners fully CP9 coupled system. We present theoretical results to guide Robust Multigrid Preconditioners for the High- the choice of preconditioners as well as numerical results Contrast Diffusion and Biharmonic-Plate Equa- illustrating the effect of those choices in practice. tions Chris Siefert, Andy Salinger We study finite volume discretization of high-contrast dif- Sandia National Laboratories fusion equation and HCT and Morley discretizations of [email protected], [email protected] high-contrast plate equation. We construct precondition- ers that are robust with respect to the magnitude of the CP9 coefficient contrast and mesh size simultaneously. For that, we prove robustness of the preconditioner proposed by Ak- ADI As Preconditioner for Krylov Subspace Meth- soylu et al. (2008) by extending the devised singular per- ods in Serial and Parallel turbation analysis to the above discretizations leading to The alternating directions implicit (ADI) method is a a same family of preconditioners used for different PDEs classical iterative method for numerically solving linear and discretizations. systems arising from discretizations of partial differential Zuhal Yeter equations. We use ADI as a preconditioner for Krylov sub- Department of Mathematics, Louisiana State University space methods for a linear system from a finite difference [email protected] approximation of an elliptic test problem in three dimen- sions. The method is attractive because it allows highly efficient matrix-free implementations both in serial Matlab Burak Aksoylu and in parallel C with MPI. Louisiana State University, Department of Mathematics Center for Computation & Technology Kyle Stern [email protected] University of Maryland, Baltimore County [email protected] CP10 \√ Matthias K. Gobbert New Upper Bounds on the Complexity of - University of Maryland, Baltimore County Complete Algorithms Department of Mathematics and Statistics [email protected] In this talk we derive some new upper bounds on the com- plexity of algorithms which can solve NP-Complete prob- Andrei Draganescu lems. Also upper bounds on the complexity of algorithms University of Maryland, Baltimore County which have a high likelihood of solving NP-Complete prob- [email protected] lems within the given complexity bounds will be discussed. Emphasis will be placed on one or two NP-Complete prob- lems. CP9 Adaptive Finite Element Method for Elliptic Equa- Michael Laidacker tions over Hierarchical T-Meshes Lamar University Beaumont, Texas Isogeometric analysis based on NURBS (Non-Uniform Ra- [email protected] tional B-Splines) as basis functions preserves the exact ge- ometry but suffers from the inconvenience of a purely local refinement. In this paper, we use the so called polyno- CP10 mial splines over hierarchical T-meshes (PHT-splines) to Production Lot Sizing with a Secondary Outsourc- construct basis functions which share the similar proper- ties as B-splines, and allow us to overcome the difficulty of AN10 Abstracts 17 ing Facility lem instances quickly on a PC. An extended Economic Production Quantity (EPQ) model Hossein M. Soroush is investigated in this paper where a fixed lot sizing pol- Kuwait University icy is implemented to reduce the complexity of produc- P.O. Box 5969 tion planning and inventory control, and outsourcing with [email protected] a secondary facility is used to supplement the lot sizing policy and to cope with the random demand. The consid- ered cost includes: set up cost for the batch production, CP10 inventory carrying cost, backorder cost when the demand A Polynomial Arc-Search Path-Following Algo- cannot be met during the production period, and outsourc- rithm for Linear Programming ing cost when the total demand is greater than the lot size in one production cycle. Under some mild conditions, the A simple analytic arc is proposed to approximate the cen- expected cost per unit time can be shown to be convex tral path of the linear programming. A primal-dual path- with respect to the lot size. The average cost reduction of following interior-point algorithm is developed to search the proposed model is 57.5%, when compared with that of the optimal solution along the analytic curve. The algo- rithm is proved to be polynomial with complexity bound the classical lot sizing policy. Outsourcing in the produc- 1 tion lot sizing policy contributes to a significant portion O(n 2 log(1/)). Numerical test is conducted for prob- of this cost savings when the mean demand rate is high. lems in Netlib. The result shows that the new algo- The numerical results also demonstrate that randomness rithm is promising compared to Matlab optimization tool- of demand has significant impact on the lot sizing policy. box linprog which implements the state-of-art Mehrotra’s predictor-corrector algorithm. Shine-Der Lee National Cheng Kung University Yaguang Yang Dept. of Industrial & Information Management NRC [email protected] [email protected] Shu-Chuan Lan, Chin-Ming Yang National Cheng Kung University CP11 Dept. of Industrial & Inform Management Nonlinear Asymptotic Stability in the Semi-Strong [email protected], [email protected] Pulse Regime for the Gierer-Meinhardt Equation with N-Pulse Positions CP10 This paper shows the nonlinear asymptotic stability of an Static Vertex Reordering Schemes for Local Mesh N-pulse solution to the Gierer-Meinhardt equation. Us- Quality Improvement ing semigroup and resolvent estimates, we are able to prove stability with renormalization group methods. In Vertex reordering can be performed within the context of the semi-strong limit the localized activator pulses interact local mesh quality improvement with the goal of decreas- strongly through the slowly varying inhibitor. The interac- ing the amount of time required for the mesh optimization. tion causes pulse amplitudes and speeds to change as the In this talk, we investigate various static vertex reordering pulse separation evolves. In addition the point spectrum schemes, based on mesh quality and the performance of the of the associated linearized operator evolves with the pulse optimization algorithm, and the trade-offs between order- dynamics. ing and overall performance of the optimization algorithm. The study uses the Laplace smoother within the Mesquite Thomas Bellsky package to optimize hexahedral meshes. Michigan State University Department of Mathematics Jeonghyung Park, Suzanne M. Shontz [email protected] The Pennsylvania State University [email protected], [email protected] Keith Promislow Michigan State University Patrick M. Knupp [email protected] Sandia National Laboratories [email protected] CP11 Nonsmooth Systems and Their Application CP10 Minimizing Weighted Quadratic Functions of Job For many years, nonsmooth systems have been used as Lateness in the Single Machine Scheduling System simplifications for certain types of smooth systems. In this with Inserted Idle Time talk we will show how this approach can lead to significant differences between solutions of both systems. We will il- This paper studies the problem of finding optimal sched- lustrate the talk using the example of a genetic regulatory ules that minimize the weighted sum of some quadratic network. functions of job lateness on a single processor where an idle time is inserted before the processing of the first job John Hogan begins. We introduce a novel exact two-stage algorithm Bristol Centre for Applied Nonlinear Mathematics for this NP-hard problem based on a precedence relation Department of Engineering Mathematics, University of structure among adjacent jobs. Our computational results Bristol demonstrate that the algorithm can solve very large prob- [email protected] 18 AN10 Abstracts CP11 Cornell University Traveling Waves Driven by Spatio-Temporally [email protected] Varying Stimuli in Neural Field Equations Bernd Krauskopf We examine the existence of traveling waves for continuum University of Bristol neuronal networks modeled by integro-differential equa- Department of Engineering Mathematics tions. For a scalar field model with a general firing rate [email protected] function and spatio-temporally varying stimulus, we show stimulus-locked fronts exist for a certain interval of stim- Hinke M. Osinga ulus speeds. We also add a slow adaptation equation and University of Bristol obtain a formula, involving an adjoint solution, for stimu- Department of Engineering Mathematics lus speeds that induce locked pulses and perform a singular [email protected] perturbation analysis to approximate the adjoint. Jozsi Z. Jalics Martin Wechselberger Department of Mathematics and Statistics University of Sydney Youngstown State Univeristy [email protected] [email protected] CP11 G. Bard Ermentrout, Jonathan E. Rubin University of Pittsburgh Parallel Time Integration for a Membrane Problem Department of Mathematics A Ratio-based Parallel Time Integration method is applied [email protected], [email protected] to a membrane second order initial value problem, given by y −b|y |q−1y +|y|p−1y =0,wherep ≤ q ≤ 2p/(p+1). The CP11 solution y(t) exhibits an oscillatory behavior and blows- On the Stability of An Operator-Difference Scheme up as t tends to infinity. After transforming the initial for Non-Linear Hyperbolic Differential Equation value problem into a first order system of ordinary dif- ferential equations, we use a method that automatically We study the stability of second-order of accuracy differ- generates time-slices and rescales the time variable over ence scheme for the nonlinear hyperbolic equation u (t)+ every slice. Theoretically, in the case when q =2p/(p +1), A(t)u(t)=f (t, u, u )(0≤ t ≤ T ),u(0) = ϕ, u (0) = ψ in a such method leads to invariance with respect to the time Hilbert space H with the self-adjoint positive-definite oper- slices in the sense that the computation of y(t) is reduced ator A(t) with the domain of D (A (t)) . A new second-order to its finding over one slice. Whereas in the case when of accuracy difference scheme with operator coefficient in p ≤ q<2p/(p + 1), the method yields asymptotic simi- a Hilbert space for the approximate solution of the initial- larity to a limit model, in the sense that the the rescaled value problem is contracted. The stability estimates for systems admit a limit system. In both cases, this leads to the solution of this difference scheme are established. an interesting Ratio based Parallel Time Integration algo- rithm - RaPTI. Its implementation on a cluster of 2 to 8 Mehmet Emir Koksal,RaviP.Agarwal processors proves to be extremely efficient (decreasing ex- Florida Institute of Technology ecution time when increasing the number of processors). mkoksal@fit.edu, agarwal@fit.edu ALLABEREN Ashyralyev Nabil R. Nassif Fatih University Mathematics Department [email protected] American University of Beirut [email protected] CP11 CP12 Mixed-Mode Oscillations in the Koper Model Shift-Invert Arnoldi Approximation to the Toeplitz Mixed-mode oscillations (MMOs) appear in a wide variety Matrix Exponential of applied dynamical systems arising in chemistry, biology and physics. The model investigated by Koper in 1995 pro- The shift-invert Arnoldi method is employed to generate an vides an illustration how fast-slow system mechanisms can orthonormal basis from the Krylov subspace corresponding explain the complicated MMO patterns. We shall illustrate to a real Toeplitz matrix and an initial vector. The vec- the techniques required for the analysis which include ana- tors and recurrence coefficients produced by this method lytical as well as numerical ideas. An outlook on a possible are exploited to approximate the Toeplitz matrix exponen- general classification of MMOs will be given as well. tial. Toeplitz matrix inversion formula and rapid Toeplitz matrix-vector multiplications are utilized to lower the com- Christian Kuehn putational costs. For convergence analysis, a sufficient con- Center for Applied Mathematics dition is established to guarantee that the error bound is Cornell University independent of the norm of the matrix. Numerical results [email protected] are given to demonstrate the efficiency of the method. Spike T. Lee, Hong-Kui Pang, Hai-Wei Sun Mathieu Desroches Department of Mathematics University of Bristol University of Macau Engineering Mathematics [email protected], [email protected] [email protected], [email protected], [email protected] John Guckenheimer AN10 Abstracts 19 CP12 method. Inexact Newton and Krylov Methods for Riccati Equations Hai-Wei Sun Department of Mathematics We explore the use of Inexact Newton methods to obtain a University of Macau low-rank solution of the Continuous-time algebraic Riccati [email protected] T equation. That is, we are looking for X∗ = YY ,withY ∈ n×p R , p n, such that F (X∗)=0whereF (X)=AX + T T T n×n n×m CP12 AX − XBB XA + C C, A ∈ R ,B ∈ R ,C ∈ Rq×n,andm, q n. The standard formulation of the Interpreting IDR As a Petrov-Galerkin Method Newton method for this problem implies the solution of a In this talk, we show that the IDR method of Sonn- (different) Lyapunov equation at each iteration. We use a eveld and van Gijzen [SIAM J. Sci. Comput., 31:1035– Krylov projection method to find an approximate low-rank 1062, 2008] can be interpreted as a Petrov-Galerkin (pro- solution of this Lyapunov equation. As is well known,the jection) method with a particular choice of left Krylov sub- accuracy of these approximations do not need to be too spaces; these left subspaces are rational Krylov spaces. precise at first, but need to improve as the (outer) Newton Consequently, other methods, such as BiCGStab and method proceeds. We also explore the use of some recycling ML(s)BiCG, which are mathematically equivalent to some strategies to reduce the computational cost. versions of IDR, can also be interpreted as Petrov-Galerkin Marlliny Monsalve methods. The connection with rational Krylov spaces in- Department of Mathematics spired a new version of IDR, called Ritz-IDR, where the Temple University poles of the rational function are chosen as certain Ritz [email protected] values. Experiments are presented illustrating the effec- tiveness of this new version. CP12 Valeria Simoncini Univeristy of Bologna The Iterative Solution of General Finite Linear Sys- and IMATI-CNR, Pavia tems Via Flexible New Dual Variational Principles [email protected] These principles provide promising new approaches to numerically solving general systems of linear equations Daniel B. Szyld and/or linear inequalities – including general linear opti- Temple University mization problems. The resulting methodologies require Department of Mathematics only the numerical computation of a critical solution for [email protected] an unconstrained objective function – carefully chosen for each individual problem from an infinite number of candi- date functions, with the intent of employing parallel proces- CP13 sors and/or exploiting any special system structure (such Charge Retrieval Analysis for Lightning Flashes in as sparsity). This infinite flexibility makes these new it- a Mountain Thunderstorm erative solution methodologies competitive with, as well as much more general than, the previously developed iter- On 24 August 2007, a balloon was launched from Langmuir ative solution methodologies for solving large-scale linear Laboratory, NM. An electric field sonde (Esonde) measured systems – including both the conjugate-gradient method the electric field during the flight, while simultaneously for symmetric positive-definite systems of linear equations Lightning Mapping Array recorded the location of VHF and the interior-point methodologies for linear optimiza- pulses generated during lightning. The Esonde and LMA tion. Finally, this presentation provides an introduction data are the basis for an inverse problem whose solution to the basic ideas and fundamental theory of generalized yields the charge transport due to a lightning stroke. The geometric programming, within the familiar context of el- analysis of several interesting flashes provides new insight ementary linear algebra, using only elementary convexity into charge transport processes in a thunderstorm. theory and multi-variable differential calculus. Beyza C. Aslan Elmor L. Peterson University of North Florida Systems Science Consulting [email protected] [email protected] CP13 CP12 An Adaptive Mesh Refinement Strategy for Finite Fast Exponential Time Integration for Pricing Op- Volume Methods tions in Jump-Diffusion Models In this talk we will present a goal-oriented adaptive mesh We consider pricing options in jump-diffusion models which refinement strategy for finite volume methods. The strat- require solving partial integro-differential equations. Dis- egy is based on the a posteriori error analysis of the numer- cretizing by the central spatial finite difference scheme ical solution using the duality between the primal equation leads to linear systems of ordinary differential equations and the adjoint equation. This approach has two major with Toeplitz structure. A fast exponential time integra- advantages: one is that it does not require the numerical tion scheme, where the shift-invert Arnoldi method is em- method to be in a variational form, for which most FVMs ployed for the Toeplitz matrix exponential, is proposed to are not; the other is that the well-posedness of the ad- approximate the solutions of those systems. Numerical re- joint problem can be studied using standard PDE theories. sults are given to demonstrate the efficiency of the proposed Examples in 1D and 2D will be shown. When time permit- ting, the application of such strategy to regional modeling 20 AN10 Abstracts will be discussed. consider how they behave in the ADI iteration and connect this to projection methods for solving Lyapunov equations. Qingshan Chen Florida State University [email protected] Garret M. Flagg, Serkan Gugercin, Chris Beattie Virginia Tech Max Gunzburger [email protected], [email protected], [email protected] Dept. of Scientific Computing Florida State University CP14 [email protected] 2 A Characterizaton of Dn(q) by Order of Normal- izer of Sylow Subgroups CP13 Equilibrium Conditions and Sound Velocities in Let (V,f) be a orthogonal space, where V = V2n(q), f is a nondegenerate orthogonal form and there are maximum Two-Phase Flows − (n − 1) distinct hyperplanes in V. Define SO2n(q)={A ∈ We consider a hierarchy of hyperbolic models describing SL2n(q) | f(Av, Aw)=f(v, w) for all v, w ∈ V } and 2 − − − single-component two-phase flows in pipelines, with ap- Dn(q)=Ω2n(q)/Z,whereΩ2n(q)=(SO2n(q)) and Z − plications to CO2 capture and storage. The hierarchy is is the center of Ω2n(q). In this paper, we characterize the 2 characterized by the number of equilibrium assumptions simple group Dn(q) by the order of normalizer of its Sylow made. We present a formal proof that every additional level subgroups. of enforced equilibrium lowers the propagation velocity of pressure waves. This subcharacteristic condition holds for A. Iranmanesh arbitrary thermodynamic state equations. We present nu- Modares University, Iran merical examples relevant for CO2 transport, and argue [email protected] the importance for pipeline integrity simulations. Halvor Lund,ToreFl˚atten CP14 SINTEF Energy Research Measuring the Connectedness of Geographically [email protected], tore.fl[email protected] Embedded Graphs Many networks surround us share the property that there CP13 exists underlying metric space in which the nodes are em- Nonlinear Homogenization Problem in the Acous- bedded in. For example, neurons in the brain, or flight tics of a Two-Phase Medium connections. To measure the efficiency one can navigate in such networks with only local information, we introduce Effective equations modeling the acoustics of a two-phase the concept of greedy connectivity for such geographical medium with periodic microstucture are derived. The networks. By using a chemical potential like parameter, medium (such as cancellous bone) is composed of a linear the greedy connectivity accounts also for imperfect trans- viscoelastic matrix filled with non-Newtonian fluid. Non- mission across established links. linear governing equations are obtained via two-scale ho- mogenization and other weak convergence techniques. The Jie Sun effective model is a two-velocity system for the effective ve- Clarkson University locity v and a corrector velocity w. The effective stress is [email protected] explicitly dependent on the sum of the strain rates. Ana Vasilic CP14 UAE University On the Largest Eigenvalue of Distance Matrix of [email protected] Some Nanotubes Robert Gilbert Let G be a connected simple graph with set of vertices { } Department of Mathematical Sciences V (G)= v1,...,vn . The distance matrix D(G)ofG is University of Delaware a square matrix of order n,whoseentrydij the between [email protected] the vertices vi and vj in G. The eigenvalues of D(G)are denoted by μ1(G),μ2(G),...,μn(G). Since the distance matrix is symmetric, its eigenvalues are real numbers and Alexander Panchenko can be ordered as μ1(G) ≥ μ2(G) ≥ ...≥ μn(G). Washington State University [email protected] A. Tehranian Science and Research Branch, Islamic Azad University, Tehran CP14 [email protected] H2-Optimal Interpolation: New Properties and Applications CP14 H2-optimal interpolation points have received great atten- You Can’t Beat Gibbs and Runge tion recently due to their success in producing accurate reduced-order models in a numerically efficient way. In Suppose you sample an analytic function f in n equispaced this talk, first, we will show how H2 optimal interpola- points in [−1, 1]. Can you use this data to approximate f in tion points can be applied to optimal H∞ model reduction. a manner that converges exponentially as n →∞?Many Second we consider their potential-theoretic properties and algorithms have been proposed that are effective for prac- connection to the rational Zolotarev problem. Thirdly, we tical values of n, but we prove that they must all fail in the AN10 Abstracts 21 limit n →∞: exponential convergence implies exponential measure theoretic approach to approximate a probability ill-conditioning. measure on the input space. Nick Trefethen Troy Butler Computing Laboratory Institute for Computational Engineering and Sciences University of Oxford University of Texas at Austin [email protected] [email protected] CP14 CP15 Application of Fractional Calculus to the Analysis Large Deviation of Quantum Markov Semigroups. of Laplace Transformed Data It has been found recently in a series of papers by F. Fag- This paper describes a novel method using fractional cal- nola and R. Rebolledo that the time averages of a Markov culus to estimate non-integer moments of a random vari- semigroup of completely positive linear maps acting on a able from the measured Laplace transform of its probability von Neumann algebra converges in weak* topology to a density function. We demonstrate that the ω-th moment stationary quantum state and that necessary and sufficient (ω ∈ R) of the random variable can be directly obtained conditions have been established for the faithfulness of the by a linear transformation of the data. When ω>0, com- stationary state. In this paper we obtain a large deviation putation of moments corresponds to fractional integration principle for this convergent sequence. The rate function of the data. When ω ≤ 0, computation of moments corre- for large deviation is also established. sponds to fractional differentiation. Mou-Hsiung Chang Lalitha Venkataramanan, Tarek Habashy, Denise Freed U. S. Army Research Office Schlumberger Mathematics Division [email protected], [email protected], [email protected] [email protected] CP15 CP15 Analysis of Centralized Sequential Probability Ra- Optimal Control of a Stochastic Volterra Integral tio Test (SPRT)-Based Models for Group Decision- Equation in Hilbert Space with Fractional Brown- Making ian Motion Input Centralized group decision-making models are of interest This paper investigates stochastic Volterra integral equa- because important decisions are frequently made by a com- tions in Hilbert spaces. The existence and uniqueness of mittee. Our goal is to analyze the relative benefits of differ- their adapted solutions is established. The regularity of the ent group decision rules, in terms of speed and accuracy. adapted solutions to these equations is proved by means of We model the individuals in our group with the SPRT, Malliavin calculus. For an application, we study an optimal an optimal and biologically-relevant method, then combine control problem for a stochastic Volterra integral equation the individual decisions to arrive at the group decision. We with Hilbert space-valued fractional Brownian motion in- will present both analytical and computational results. put. A Pontryagin-type maximum principle is formulated for the problem, and an example is presented. Margot Kimura, Jeffrey Moehlis UCSB Vo V. Anh [email protected], [email protected] Discipline of Mathematical Sciences Queensland University of Technology, Brisbane, Qld. 4001, Au CP15 [email protected] Evaluating Expectations of Functionals of Brown- ian Motions: a Multilevel Idea Wilfried Grecksch Pricing a path-dependent financial derivative, such as an Institute of Mathematics, Martin-Luther-University · Halle-Wittenberg 06099 Halle (Saale), Germany Asian option, requires the computation of E[g(B( ))], the expectation of a payoff functional, g, that depends on a [email protected] T Brownian motion, (B(t))t=0. This problem turns out to be an infinite dimensional problem. A multilevel algorithm Jiongmin Yong with low discrepancy designs is introduced to approximate University of Central Florida the infinite dimensional integral. The worst case error as a [email protected] functional of each level’s sample size and truncated dimen- sion is minimized. Numerical examples in computational CP15 finance will be presented. A Computational Measure Theoretic Approach to Ben Niu Inverse Sensitivity Analysis Department of Applied Mathematics Illinois Institute of Technology We consider the inverse sensitivity analysis problem of [email protected] quantifying the uncertainty of inputs to a finite dimensional map given specified uncertainty in a linear functional of the output of the map. Our formulation is based on the Law of Fred J. Hickernell Total Probability that represents a direct inversion of the Illinois Institute of Technology forward stochastic sensitivity problem for a deterministic Department of Applied Mathematics model. We derive and analyze an efficient computational [email protected] 22 AN10 Abstracts Thomas Muller-Gronbach tem in order to enhance detection of faults. Unlike previous U Passau studies where it is usually assumed that there is only one [email protected] fault occurring at a time, we allow for multiple faults to occur simultaneously which is a natural assumption in the Klaus Ritter incipient case. A computational method for the construc- TU Darmstadt tion of an optimal input signal for achieving guaranteed Fachbereich Mathematik detection with specified precision is presented. This work [email protected] is an extension of the multi-model approach used for the construction of auxiliary signals for incipient failure detec- tion. There are both discrete time and continuous time CP15 versions of this problem. The discrete time case with a Benchmark Approach for Defaultable Claims modified noise bound is examined. The continuous time case with a similar model that requires solving a linear We study the benchmark approach for defaultable claims. quadratic regulator problem is also studied. Both models The local risk minimizing strategy is derived in the case involve only additive uncertainty. when the agent information takes into account the possibil- ity of a default event. Here we do the pricing and hedging Martene L. Fair of the defaultable claim under the real world probability North Carolina State University measure only. martene [email protected] Indukuri Raju Stephen L. Campbell Indian Institure of Tehcnology Guwahati North Carolina State Univ [email protected] Department of Mathematics [email protected] CP15 Convergence of Tau Leaping Schemes for Markov CP16 Jump Systems on a Lattice PageRank Optimization Through Ergodic Control Tau leaping methods provide an efficient way to simulate We study a general class of pagerank optimization prob- Markov jump systems on a lattice. Existing convergence lems, which consist in finding an optimal outlink strategy results on tau leap methods apply only to systems that re- for a web site subject to design constraints. We model main in a bounded region and/or are specific to the explicit these problems by constrained Markov decision processes, and the implicit tau leap schemes with Poisson random in which the webmaster determines the transition probabil- variates. We present new convergence results that deal ities of websurfers. We identify assumptions under which with fairly general tau leap schemes applied to unbounded there exists a “master’ page to which all controlled pages systems that possess certain moment growth bounds. should point. We report numerical results exploiting dy- namic programming and convex programming techniques. Muruhan Rathinam University of Maryland, Baltimore County [email protected] Olivier Fercoq INRIA Saclay and CMAP Ecole Polytechnique [email protected] CP16 Convergence of Discontinuous Time-Stepping Marianne Akian Schemes for the Gradient Minimization Boundary INRIA Saclay–Ile-de-France and CMAP, Ecole Control Problem Polytechnique The minimization of the energy functional having states [email protected] constrained to semi-linear parabolic PDE’s is considered. The controls act on the boundary and are of Robin type. Mustapha Bouhtou The discrete schemes under consideration are discontinu- France T´el´ecom R & D ous in time but conforming in space. Stability estimates [email protected] are presented at the energy norm and at arbitrary times for the state, and adjoint variables. The estimates are de- St´ephane Gaubert rived under minimal regularity assumptions and allow us INRIA Saclay – Ile-de-France & to prove onvergence of the discrete solutions. Centre de Math´ematiques appliqu´ees, Ecole Polytechnique Konstantinos Chrysafinos [email protected] Department of Mathematics National Technical University of Athens, Greece chrysafi[email protected] CP16 Adaptive and Higher Order Numerical Solution for Optimal Control of Monodomain Equations in Car- CP16 diac Electrophysiology Active Incipient Fault Detection With Multiple Si- multaneous Faults The focus of this work is on the development and im- plementation of an efficient numerical technique to solve The problem of detecting small parameter variations in an optimal control problem related to a reaction-diffusions linear uncertain systems due to incipient (slowly develop- system arising in cardiac electrophysiology. A Newton-type ing) faults is considered. Using an active fault detection method for the monodomain model, which is a well estab- approach, an input signal is injected into the observed sys- AN10 Abstracts 23 lished model for describing wave propagation of the action Ume{ potential in the heart, is developed. The numerical treat- aa} Univsersity ment is enhanced by using a second order time stepping [email protected] method, adaptive grid refinement techniques and receding horizon methods. A super linear convergence is achieved Martin Berggren for Newton’s optimization algorithm. Department of Computing Science, Umea University [email protected] Chamakuri Nagaiah Institute of Mathematics and Scientific Computing, University of Graz, Austria CP17 [email protected] Modeling Complex Fluid Flows in Porous Media Karl Kunisch Flow of complex fluids in porous media occurs in many Universitat Graz processes, some naturally occurring and some by design. Institut fuer Mathematik Modeling of such flows accurately and efficiently is a non- [email protected] trivial task. We will discuss in this talk many non-trivial issues that come up in the context of setting the correct Gernot Plank mathematical model that includes myriad effects that play Institute of Biophysics a role. In this context we may take some simple model Medical University of Graz, Austria problems and show some exact and computational results. [email protected] This is an ongoing work. Prabir Daripa CP16 Texas A&M University Why Walk and Run: Towards a Predictive Theory Department of Mathematics of Legged Locomotion Based on Energy Optimality [email protected] Healthy human legs are capable of great behavioral com- plexity – some human legs are capable of even ”break- CP17 dancing”. But mostly, people use stereotypical gaits – Multiscale Modeling of Heterogenous Flow with ”walking” at low speeds and ”running” at high speeds. Non-Newtonian Behavior Why these particular gaits from among the millions our legs are capable of? Using optimal control calculations us- A new multiscale approach to modeling non-Newtonian ing a simple mathematical model and classic prior experi- /heterogenous flow will be discussed. The computations ments, I will show evidence that the patterns of movement are carried out by combining an Immersed-Boundary- that we use roughly minimize energy cost of locomotion. I Method with Brownian Dynamics. The exchange of in- will show how to extend this simple theory to obtain better formation between both computational techniques is done quantitative predictions of people walk and run. through the velocity gradient and stress state tensors. Sev- eral examples will be presented including the response of Manoj Srinivasan a Newtonian drop suspended in a viscoelastic fluid and Movement Lab, Mechanical Engineering the dynamics of Newtonian bubbles rising in a viscoelastic Ohio State University fluid. [email protected] Arturo Fernandez Department of Mechanical Engineering CP16 Worcester Polytechnic Institute Efficient Loudspeaker Horn Design: from Mathe- [email protected] matical Model to Experimental Validation During the last few years, we have designed numerous CP17 acoustic devices using gradient-based boundary shape and Three-Dimensional Multispecies Nonlinear Tumor topology optimization. This talk outlines strategies to op- Growth: Tumor Invasion and Angiogenesis timize acoustic devises with respect to efficiency and di- rectivity properties. The acoustical properties are calcu- We study 3D non-linear tumor growth by tracking multiple lated through numerical solutions of the Helmholtz equa- viable cell species populations, discrete angiogenesis vessels tion. Recently, we manufactured one of our designs and an- and individual cell movement using a hybrid continuum- alyzed it in an anechoic chamber. In the design frequency discrete approach. We investigate disease progression as band, the acoustic impedance is close to ideal. a function of cellular-scale parameters such as prolifera- tion and nutrient uptake rates in both avascular and vas- Eddie Wadbro cular regimes. We find that heterogeneity in the physio- Department of Computing Science logically complex tumor microenvironment can be quanti- Ume˚aUniversity tatively linked to the tumor macro-scale as a mechanism [email protected] that promotes morphological instability. Fang Jin Rajitha Udawalpola University of California, Irvine Department of Information Technology [email protected] Uppsala University [email protected] Herman Frieboes, Yaoli Chuang School of Health Information Science Daniel Noreland The University of Texas Health Science Center, Houston, Department of Computing Science 24 AN10 Abstracts TX resolved by introducing a singular perturbation reduction [email protected], [email protected] of the transition layer into full numerical solution of the free boundary problem. Numerical results are presented Steven Wise for two-phase flows in the Stokes flow limit. Department of Mathematics University of Tennessee, Knoxville, TN Kuan Xu [email protected] Department of Mathematical Sciences New Jersey Institute of Technology [email protected] John Lowengrub Department of Mathematics University of California at Irvine Michael R. Booty [email protected] Department of Mathematical Sciences, NJIT [email protected] Vittorio Cristini School of Health Information Science Michael Siegel The University of Texas Health Science Center, Houston, New Jersey Institute of Technology TX [email protected] [email protected] CP18 CP17 On Accurate Reduced-Order Modeling Using Eulerian Formulation of Fluid/structure Interac- Proper Orthogonal Decomposition tion Reduced-order models based on the proper orthogonal de- The interaction of an elastic structure and a fluid occurs in composition (POD) are often used to represent a complex many phenomena in physics. To avoid the difficulty of cou- dynamical system in fluid flows. However, for Navier- pling lagrangian elasticity and eulerian fluid we consider a Stokes equations, these models have been successfully de- whole eulerian formulation. The elasticity of the structure veloped mainly for 2-D flows and are limited to a reference is computed with retrograde caracteristics which satisfy a simulation. We present improvements in the reduced-order transport equation and a level set function is introduced to model to expand its applicability and encounter various capture the interface between the fluid and the structure. stability issues in the model. We then compare results of In application, we will present some numerical simulation reduced-order models with direct numerical simulations. relative to biolocomotion. Imran Akhtar Thomas M. Milcent Virginia Tech University of Bordeaux 1, INRIA [email protected] [email protected] Jeff Borggaard Virginia Tech CP17 Department of Mathematics A Parallel Monolithic Approach for Blood Flow [email protected] Modeling Traian Iliescu, Zhu Wang We study a parallel algorithm for the simulation of blood Department of Mathematics flows in arteries using a coupled system of PDEs consist- Virginia Tech ing of incompressible Navier-Stokes equations and a linear [email protected], [email protected] elastic equation. The coupled nonlinear system is solved by a Newton-Krylov algorithm preconditioned by an over- lapping Schwarz method. We focus on the investigation CP18 of the algorithm for different outflow boundary conditions, Semiclassical Asymptotics in Norm for Non- and also on the scalability of the software on computers Smooth and for Nonlinear Problems with a large number of processors. Quantum problems, formulated appropriately, are well Yuqi Wu known to converge to their classical counterparts as Dept. of Applied Math., University of Colorado Planck’s constant ¯h tends to zero. On a mathematical level, [email protected] this is traditionally seen as a weak-∗ proximity between the solutions of the quantum and classical phase-space kinetic Xiao-Chuan Cai equations (Wigner and Liouville respectively). Indeed a University of Colorado, Boulder lot is known for the weak-∗ convergence to the semiclas- Dept. of Computer Science sical limit, both for linear and nonlinear problems with [email protected] sufficient smoothness. However, much less is known for strong topology approximations; the first systematic ex- ploration of H s convergence, for linear problems, appeared CP17 in 2006. In this work we show semiclassical approximation A Hybrid Numerical Method for Two-Phase Flow results, for the linear case (for pure and mixed states), in with Soluble Surfactant the Hs, s ∈ Z, norms, with weaker assumptions than the state of the art. Moreover, in certain cases we can lower the We present a hybrid numerical method for computation of smoothness assumptions enough for our results to be appli- fluid interfaces with soluble surfactant at large bulk Peclet cable to C1 \ C1,1 potentials. In such problems the weak-∗ number. A narrow transition layer adjacent to the interface approximation is known to hold, but the limit Liouville which is important in controlling the interface dynamics is AN10 Abstracts 25 equation has multiple solutions: the semiclassical limit is known to be a(weak)solution of the Liouville equation, but which one is not known. In a family of such problems, Anita Mareno we can choose the solution of the Liouville equation which Mathematics corresponds to the limit of the Wigner equation. Moreover, Penn State University for the first time to the best of our knowledge, we show an [email protected] L2 semiclassical limit for mixed states evolving under the Hartree equation, i.e. we show that the solutions of the Wigner-Poisson equation and the nonlinear Vlasov equa- CP18 tion are close in L2 as ¯h tends to zero. These results are The Effect of Upscaling on Seismic Inversion joint with T. Paul, M. Pulvirenti and F. Pezzoti, and use heavily the technique of the smoothed Wigner transform. Small-scale geologic features may have a large impact on estimates of material parameters determined by seismic Agisilaos G. Athanasoulis inversion. Unfortunately there is no magic scale known CMLS a priori on which the important features of the subsur- Ecole Polytechnique face reside. Two-scale forward modeling leads to interest- [email protected] ing questions about model parameter estimation. In the case of constant density acoustics, studies indicate the es- timated sound velocity should correspond to an average, CP18 but numerical experimentation is required in the case of On the Right Hand Side Identification Problem for variable-density acoustics. a Parabolic Equation with Nonlocal Conditions Susan E. Minkoff We consider the inverse problem of reconstructing the right University of Maryland, Baltimore County ∂u(t,x) ∂u(t,x) side of a parabolic equation = ∂ k (x) − Dept of Mathematics and Statistics ∂t ∂x ∂x sminkoff@umbc.edu σu (t, x)+η (t) ψ (x) , 0 [email protected] MATLAB for sparse rectangular matrices. Timothy A. Davis CP19 University of Florida Storm Surge Simulations Using Adaptive Mesh Re- Computer and Information Science and Engineering finement [email protected] We present preliminary results of storm surge simulations with idealized bathymetry and hurricane model data with CP20 adaptive mesh refinement. The goal of this study is to Fast Computation of Schur Complements of Large gain an understanding of how adaptive mesh refinement Finite Difference Matrices could provide greater computational efficiency and higher resolution than traditional methods used to model these It is well-known that the linear equations associated with phenomena. Multi-layer methods are also investigated in finite difference operators on structured grids admit very order to include basic three dimensional affects in the sim- fast solution techniques. Recently developed fast direct ulation. solvers require much less memory than FFT for constant coefficient operators, and are more robust than iterative Kyle T. Mandli techniques in the general case. The direct solvers are also University of Washington extremely fast for problems with multiple right-hand sides. Dept. of Applied Mathematics We will present fast methods for computing Schur comple- [email protected] ments of a finite difference matrix. Randall J. LeVeque Adrianna Gillman Applied Mathematics University of Colorado at Boulder University of Washington (Seattle) Department of Applied Mathematics [email protected] [email protected] Per-Gunnar Martinsson CP19 University of Colorado at Boulder Accuracy of Tail Regions in Uncertain Climate [email protected] Model Predictions Conventional spectral methods for uncertainty quantifica- CP20 tion are generally challenged in the “tails” of input distri- Parallel Solution of Large Symmetric Eigenvalue butions. Extensive sampling in the tail regions is especially Problems Via Trace Minimization costly in climate models. We will illustrate how the usage of non-conventional basis functions and surrogate model The Trace Minimization algorithm (TraceMIN) is a sym- analysis improves the convergence of the spectral expan- metric eigenvalue problem solver that minimizes the trace sions in the tail regions. Preliminary results will be shown over a subspace. TraceMIN requires the solution of a set for the Atlantic meridional overturning circulation shut- of linear systems with modest accuracy. In this talk we off case as a result of abnormally high climate sensitivity present the convergence properties and the parallel scal- values. ability of a new version of TraceMIN. In addition, we compare the parallel performance of our algorithm with Khachik Sargsyan, Cosmin Safta PARPACK for solving eigenvalue problems that arise in Sandia National Laboratories MEMS device simulation and other applications. [email protected], [email protected] Murat Manguoglu Habib N. Najm Purdue University Department of Computer Science Sandia National Laboratories [email protected] Livermore, CA, USA [email protected] Faisal Saied Purdue University Bert J. Debusschere [email protected] Energy Transportation Center Sandia National Laboratories, Livermore CA Ahmed Sameh [email protected] Department of Computer Science Purdue University [email protected] CP20 Multifrontal Multithreaded Rank-Revealing Sparse QR Factorization CP20 New Results in Numerically Solving Large-Scale SuiteSparseQR is a sparse QR factorization package based Algebraic Sylvester Equations on the multifrontal method. LAPACK and the multi- threaded BLAS are used within each frontal matrix. Par- Two approaches are investigated for constructing numer- allelism across different frontal matrices is handled with ical solutions of large-scale Sylvester equations. The first Intel’s Threading Building Blocks library. Rank-detection is to transform the Sylvester equation into an invariant is performed within each frontal matrix using Heath’s subspace problem, which can be solved efficiently in real method. The resulting sparse QR factorization obtains a arithmetic. Optimal shift selection is presented to enhance substantial fraction of the theoretical peak performance of convergence of iterative eigensolvers. This optimal shift se- a multicore computer. The method is used in backslash in lection analysis is then generalized for the computation of AN10 Abstracts 27 multiple optimal real shifts for the ADI approach, resulting CP21 in a more compact search region for the optimal shifts. System Identification in Tumor Growth Modeling Hung Nong Several theoretical efforts were devoted to setting up Sandia National Laboratories phenomenological models describing cancer propagation. [email protected] Among them models based on partial differential equa- tions describe macroscopic and mechanical features of tu- Danny C. Sorensen mor propagation on the basis of mixture theory. Such mod- Rice University els are parametric but parameters cannot be fixed apri- [email protected] ori because they do not derive from first principles. Also, they cannot be measured directly because they represent mesoscale properties. In order to approach realistic appli- CP20 cations, i.e., provide a prognosis or a therapy evaluation Piro sky: Pipelined Rotations for Sparse Skyline tool, realistic values for these parameters must be found. Svd The aim of this work is to present an efficient procedure to solve data assimilation problems, the data coming from PIRO SKY is a library for computing all the singular val- medical imagery. Several realistic cases are analyzed, show- ues of a sparse matrix efficiently. PIRO SKY uses sparse ing that the prognosis can be accurate on a significant time data structures and blocked algorithms to compute the scale. singular values. The impact of various profile reduction ordering strategies on computing the sparse SVD will be Angelo Iollo, Thierry Colin, Damiano Lombardi, Olivier discussed. We show PIRO SKY is competitive against it- Saut erative methods even when only a reasonable number of Institut de Math´ematiques de Bordeaux singular values are required. We also present the software [email protected], [email protected] engineering aspects of PIRO SKY’s design. bordeaux1.fr, [email protected], [email protected] Sivasankaran Rajamanickam Sandia National Laboratories [email protected] CP21 A Lagrangian Scheme to Solve the Monge- Timothy Davis Kantorovich Problem University of Florida Dept of CISE The Monge-Kantorovich problem was recently re- [email protected]fl.edu discovered and studied, having applications in statistics, meteorology, numerical optimization, physics, medical im- agery. In the literature three family of solution methods CP21 were proposed: the first class is based on the Kantorovich Numerical Tests of An Algorithm for Seismic Imag- duality, the second one on the continuous formulation and ing and Inversion the third one on the critical points of a dynamical flow (AHT). In this work we propose a different approach: we We present numerical tests for an acoustic inverse scat- implicitly satisfy the transport equation by propagating tering algorithm for simultaneous geophysical imaging and mass particles on rays. Particle discretization makes the amplitude correction from measured data. No knowledge integration simpler and the parallelization efficient. Sev- about the medium under investigation is assumed. We eral 2D and 3D cases are presented. model the data from several one-dimensional earth config- urations and show how the algorithm can find the precise Damiano Lombardi, Angelo Iollo location and a good estimate of the layers’ parameters from Institut de Math´ematiques de Bordeaux data only. Our tests will include different number of layers, [email protected], high/low contrasts, velocity inversions and noisy data. [email protected] Ashley Ciesla, Bogdan G. Nita Montclair State University CP21 [email protected], [email protected] Simultaneous Imaging and Inversion Using An In- verse Scattering Algorithm CP21 We present a method, derived from inverse scattering the- Reduced Order Modeling to Enhance Thermal ory, for geophysical imaging and amplitude correction from Imaging for Crack Detection measured data. No knowledge about the medium under in- vestigation is assumed. Although derived as a series, the Proper orthogonal decomposition is employed as part of an algorithm is shown to converge to a closed form indepen- algorithmic means to identify micro-cracks in metal plates. dent of the parameters involved in the problem. Analytic As model-updating approaches are used to solve the re- and numerical one dimensional examples show excellent sulting inverse problem, POD offers relief from the ensuing results in finding both the location of interfaces and the computational challenges. A new approach to incorporat- amplitude of acoustic reflections. ing POD modes, as an enrichment to the finite element ba- sis functions used in the forward models, will be treated. Bogdan G. Nita, Ashley Ciesla Montclair State University [email protected], [email protected] Christopher J. Earls Cornell University [email protected] 28 AN10 Abstracts CP21 CGS2 algorithm. Semi-Blind Deconvolution in 4Pi-Microscopy Jesse L. Barlow 4Pi-microscopy has been invented by Stefan Hell to im- Penn State University prove the resolution of confocal fluorescence microscopy Dept of Computer Science & Eng along the optical axis. The convolution kernel (or point [email protected] spread function, short psf) of a 4Pi-microscope differs from the psf of a corresponding confocal microscope approxi- Alicja Smoktunowicz mately by a squared cosine factor resulting in a main peak Warsaw University of Technology of much smaller band-width and several (typically at least Warsaw, 00-0661 Poland two) side lobes. The size and position of the main peak [email protected] and the side lobes of the 4Pi-psf depend on a phase pa- rameter, which has been assumed to be space-invariant so far, yielding standard deconvolution problems with posi- CP22 tivity constraint. However, this assumption is violated e.g. Polynomial Transform Decomposition Using Mod- for inhomogeneous refractive indices, and the space depen- ule Induction dence of the psf is considered as one of the main problems in 4Pi-microscopy. In this project the joint recovery of the Polynomial transforms, or generalized Vandermonde ma- three-dimensional density of fluorescent markers and the trices, are an important class of matrix-vector products slowly varying phase function is considered as a nonlinear used in mathematics and engineering. Important exam- inverse problem, which is tackled by a variant of the it- ples include the discrete Fourier transform, discrete cosine eratively regularized Gauss-Newton method (IRGNM) re- and sine transforms, and transforms related to orthogonal specting the positivity constraint. We present a conver- polynomials. We present a novel method to derive fast gence proof for this constraint IRGNM and reconstruction transform algorithms by decomposing the associated poly- results for artificial and real data. nomial algebra using a technique from representation the- ory called induction. The method algebraically character- Robert St¨uck, Thorsten Hohage izes and generalizes transform algorithms reported in the Georg-August-Universit¨at G¨ottingen literature. [email protected], [email protected] Aliaksei Sandryhaila, Markus Pueschel Carnegie Mellon University [email protected], [email protected] CP21 Attenuation Compensation and Boundary Segmen- tation in Ultrasound Images CP22 A Schur Complement Approach for Solving a Spatial variations of attenuation across tissue can result Nearly Hermitian System in shadowing and enhancement in ultrasound B-scan im- ages. These artifacts affect the underlying signal backscat- We discuss an approach for solving an n × n linear system ter which is the main component of ultrasound images and with rank s skew-Hermitian part (s n). This can be has clinical significance in detecting diseases and tumors. interpreted as the Schur complement of a larger (n + s) × We present a joint estimation method based on the varia- (n+s) system. We can solve the original system by solving tional principle to compensate for attenuation artifacts via s + 1 Hermitian systems, e.g., using MINRES, and directly functional minimization and regularization. Pair of patho- solving one s × s system. We present a bound for the logical useful backscatter and attenuation fields is recon- outer residual norm in terms of residuals of the Hermitian structed along with segmented anatomic structures. systems and develop a suitable stopping criteria. Jue Wang Kirk M. Soodhalter,DanielB.Szyld Union College Temple University [email protected] Department of Mathematics [email protected], [email protected] Yongjian Yu InfiMed Inc. Mark Embree yuy@infimed.com Department of Computational and Applied Mathematics Rice University [email protected] CP22 The Numerical Stability of Reorthogonalized Block Josef Sifuentes Classical Gram-Schmidt Rice University [email protected] A reorthogonalized block classical Gram–Schmidt algo- rithm is given that factorizes a full column rank matrix A into A = QR where Q is left orthogonal (has orthonor- CP22 mal columns) and R is upper triangular and nonsingular. Numerical Stability of a Complete CS Decomposi- With appropriate assumptions on the diagonal blocks of R, tion Algorithm the algorithm, when implemented in floating point arith- metic with machine unit εM , produces Q and R such that A proof of backward stability for the simultaneous bidiago- T I − Q Q F = O(εM )and A − QR F = O(εM A F ). nalization of all four blocks of a partitioned unitary matrix The resulting bounds also improve a previous bound by will be outlined. The reduction serves as the first phase in Giraud et al. [Numer. Math. 101:87–100,2005] on their AN10 Abstracts 29 an algorithm for the complete CS decomposition. Northwestern University [email protected] Brian D. Sutton Randolph-Macon College Michael Miksis [email protected] Department of Engineering Science and Applied Mathematics CP22 Northwestern University [email protected] Fast Inexact Implicitly Restarted Arnoldi Method for Generalized Eigenvalue Problems with Spectral Transformation Stephen H. Davis Dept. of Engineering Sciences and Applied Math We study an inexact implicitly restarted Arnoldi (IRA) Northwestern University method for computing a few eigenpairs of generalized non- [email protected] Hermitian eigenvalue problems with spectral transforma- tion, where in each Arnoldi step (outer iteration) the matrix-vector product involving the transformed operator CP23 is performed by iterative solution (inner iteration) of the A Singularly Perturbed Formulation of the corresponding linear system of equations. We provide new Quantum-Corrected Energy Transport Model perspectives and analysis of two major strategies that help reduce the inner iteration cost: a special type of precon- A scaling reformulation of the quantum-corrected energy ditioner with “tuning’, and gradually relaxed tolerances transport model is presented. We find that the scaled De- for the solution of the linear systems. The effectiveness of bye length as the standard singular perturbation parameter these strategies is demonstrated by numerical experiments. for the classical models is no longer valid for the nano-scale semiconductor devices and that the scaled Planck constant Fei Xue is moderately small. We introduce here a new parame- Temple University ter called the scaled thermal conductivity and modify the [email protected] standard scaled intrinsic density to include a dependence on potentials. It is numerically shown that the modified intrinsic density and the scaled thermal conductivity are CP23 strongly singular for the continuity and energy transport Uzawa Algorithms for Non-Symmetric Saddle equations, respectively. Point Problems Jinn-Liang Liu This is joint work with Brendan McCracken (REU -2009), National Hsinchu University of Education and Lu Shu (GEMS-UDEL 2009). We introduce and ana- Dept. of Applied Mathematics lyze new Uzawa algorithms for non-symmetric saddle point [email protected] systems. Typical problems where such systems appear are, for example, in finite element or finite difference discretiza- tion of steady state Navier-Stokes systems. Convergence CP23 for the algorithms are established based on new spectral A Gradient-Augmented Level Set Approach with results about Schur complements. A new Uzawa type al- Subgrid Resolution gorithm with variable relaxation parameters is introduced and analyzed in a general framework. Numerical results Level set methods are versatile interface tracking ap- supporting the efficiency of the algorithms are presented for proaches. Common difficulties are the loss of small struc- finite element discretization of the the steady state Navier- tures, an accurate approximation of curvature, and the Stokes equations. Thanks to NSF. combination with adaptive mesh refinement. We present a gradient-augmented approach, which advects interfaces Constantin Bacuta with third order accuracy, and admits a second order ac- University of Delaware curate approximation of curvature, both with optimally Department of Mathematics local stencils. In addition, the presence of gradient infor- [email protected] mation yields a certain level of subgrid resolution, which is beneficial in the tracking of thin films or droplets. We demonstrate this advantage in fluid dynamics simulations. CP23 Self-Similar Solutions and Their Instability for Benjamin Seibold Liquid-to-Glass Phase Transitions in Vitrification Temple University [email protected] Some theoretical aspects of vitrification phase transitions are addressed by modeling the process as an isothermal Jean-Christophe Nave moving-boundary problem in which the two phases have Massachusetts Institute of Technology differing diffusivities. A similarity solution is obtained, [email protected] its linear stability is examined, and an analytic dispersion relation is investigated, yielding the stability boundaries. Rodolfo Ruben Rosales Extensions of this work in which the vitrified phase demon- Department of Mathematics strates subdiffusive behavior also will be discussed and the Massachusetts Institute of Technology dynamics and instabilities of this adjusted problem exam- [email protected] ined. Christopher A. Gruber CP23 Department of Engineering Sciences and Applied The Immersed Interface Method for Two-Fluid Mathematics 30 AN10 Abstracts Problems The osmotic effect is modeled by the relative sliding be- tween fluid and membrane normal to the membrane, which Many problems of fluid mechanics involve the interaction of is mediated by the interaction between solute and mem- two immiscible fluids. It is generally difficult or inefficient brane through immersed chemical potentials. The influ- to simulate each fluid separately using an interface-fitted ence of internal versus external ion concentration in deter- grid method. In the immersed interface method, a two- mining the volume of an elastic semi-permeable vesicle is fluid problem is formulated as one set of governing equa- studied with and without electrical effects and/or exchange tions and simulated on a fixed Cartesian grid. The effect of pumping. the two-fluid interface enters the formulation as a singular force and a numerical scheme as jump conditions. In this Pilhwa Lee talk, we will present the principal jump conditions, discuss University of Connecticut Health Center the difficulties in implementing them, and provide a few [email protected] options to overcome the difficulties. We will use the two- fluid circular Couette flow to demonstrate the accuracy and efficiency of these options. CP24 Mathematical Model of Self-sustained Firing in Miguel A. Uh Motoneurons After Spinal Cord Injury Southern Methodist University [email protected] The self-sustained firing in motoneurons caused by plateau potentials provide a mechanism to translate short lasting synaptic inputs into long lasting motor output. During the CP24 acute stage of spinal cord injury the endogenous ability to Transport Effects on Surface Reaction Arrays: Ap- generate plateaus is lost; however, during the chronic stage plications to Biosensors the plateau potentials reappear. Using two-compartment conductance based differential equation model we investi- The ubiquity of surface-volume reactions makes knowledge gate mechanisms that might contribute to spasticity due of their kinetics critical. Several biosensors designed to to the self-sustained firing. measure rate constants, such as the Flexchip and the dot- Lab, have multiple reacting zones in a single flow chan- Mini P. Kurian nel. A basic unidirectional flow model is developed and School of Mathematical and Statistical Sciences solved for typical experimental parameters using perturba- Arizona State University tion methods. The effect of zone placement along the chan- [email protected] nel can be quantified in terms of an effective Damk¨ohler number based upon position. Moreover, it is established Sharon M. Crook that zone placement across the channel does not affect the Dept of Mathematics and Statistics & School of Life measurements. Sciences Arizona State University David A. Edwards [email protected] University of Delaware Newark, DE [email protected] CP24 Mathematical Modeling of Cell Biopreservation by CP24 Desiccation Learning-Rate-Dependent Clustering and Self- We consider the model of a cell with an incompressible Development in a Network of Coupled Oscillators water-sugar mixture inside and outside the cell. The cell membrane has resistance to bending and local incompress- We investigate the role of the learning rate in a Kuramoto ibility. As water leaves through the surface, a glassy state model of coupled phase oscillators in which the coupling co- can be formed inside the cell, whose volume decreases while efficients dynamically vary according to the Hebbian learn- local area of the membrane is fixed. The hydrodynamics ing rule; here synapses between neurons are strengthened is governed by the Stokes equation. We investigate linear when they fire in synchrony. We show that either one or stability and dynamics of an initial spherical membrane. two synchronized clusters form, depending on the learning rate relative to a critical value. We also comment on the Alla Podolny possibility of self-development of neural networks within Northwestern University this model. [email protected] Lars English Dickinson College, Pennsylvania Michael Miksis USA Department of Engineering Science and Applied [email protected] Mathematics Northwestern University [email protected] Ritwik Niyogi Dickinson College [email protected] Stephen H. Davis Dept. of Engineering Sciences and Applied Math Northwestern University CP24 [email protected] On the Cellular Volume Regulation in 2-D The electrical and osmotic effects with ion exchange trans- CP24 port mechanism for the volume regulation are investigated. Frequency Modulation of the Excitable Fitzhugh- AN10 Abstracts 31 Nagumo System phase. The continuous phase is determined by large-eddy simulation in the Eulerian frame of reference whereas the Using a generalized Fitzhugh-Nagumo model we studied dispersed phase is simulated in a Lagrangian frame of ref- how temporal characteristics of an external signal affect the erence. dynamics of excitable systems. In particular, we choose a subset of parameters for the Fitzhugh-Nagumo system so Marcel Ilie that the solution contains (1) no limit cycle or (2) a unique University of Central Florida limit cycle. For this set of parameters, we investigate the Department of Mechanical, Material and Aerospace dynamic consequences of injecting the system a variable Engineering current I(t,w) with different frequency w. [email protected] Xujing Wang University of Alabama at Birmingham CP25 [email protected] Cluster Computing with Intel Nehalem and Infini- Band Interconnect CP25 High performance parallel computing depends on the in- Finite Difference Method for Smoothed Gradually teraction of a number of factors including the proces- Varied Surfaces sors, the architecture of the compute nodes, their inter- connect network, and the numerical code. In this talk,we A digital method for obtaining continuous functions with present performance and scalability studies on a new 86- (n) smoothness to a certain order (C ) from sample data node distributed-memory cluster with two quad-core Intel is designed. This method is based on gradually varied Nehalem processors per node and a quad-data rate Infini- functions discovered by Chen in 1989 and finite difference Band interconnect. These components are currently the method. This new method has been applied to real ground- most popular for general-purpose clusters and they exhibit water data and the results have validated the method. This excellent performance. method is independent from existing popular methods such as the cubic spline method and the finite element method. David Trott The new digital-discrete method has considerable advan- University of Maryland, Baltimore County tages for a large amount of real data applications. This [email protected] digital method differs from the classical discrete method that usually uses triangulations. Matthias K. Gobbert University of Maryland, Baltimore County Li Chen Department of Mathematics and Statistics University of the District of Columbia [email protected] [email protected] CP25 CP25 Nonlinear Evolution of Electrified Liquid Threads Lime: Lightweight Integrated Multiphysics Envi- and Annular Layers ronemnt The nonlinear dynamics of an electrified liquid thread in a LIME is a multiphysics coupling environment that enables cylindrical tube (electrode) which is concentrically placed coupling of both legacy application codes as well as novel is studied. The radial electric field is generated by the po- codes. It is lightweight by requiring application codes to tential difference between the fluid interface and the outer only provide a response given requisite inputs. It provides electrode. We solve the system in the zero Reynolds num- value-added coupling by incorporating any additional in- ber limit by the boundary-integral method. Pinch-off so- formation to effect more robust coupling algorithms based lutions and satellite formation are discussed. In addition, on variations of Jacobian-Free Newton-Krylov. LIME is different breakup behaviors are found at other parameter currently implemented as extensions to the Trilinos Solver values. Framework thereby leveraging state-of-the-art solution al- gorithms to drive multiphysics coupling. Qiming Wang Department of Mathematical Sciences Russell W. Hooper, Rodney Schmidt, Kenneth Belcourt New Jersey Institute of Technology Sandia National Laboratories [email protected] [email protected], [email protected], [email protected] Demetrios Papageorgiou Department of Mathematics CP25 Imperial College London, UK [email protected] Numerical Computations of Particle-Laden Turbu- lent Flows CP25 An efficient Eulerian-Lagrangian particle dispersion algo- rithm for the prediction of particle-laden flows is proposed. Mathematical and Numerical Analysis of Peridy- The volume fraction of the dispersed phase is assumed to namic Models with Nonlocal Boundary Conditions be small enough such that particle-particle collisions are We study the linear bond-based peridynamic models with a negligible and properties of the carrier flow are not modi- particular focus on problems associated with nonstandard fied. With the examination of dilute systems only the effect nonlocal displacement loading conditions. Both stationary of turbulence on particle motion has to be taken into ac- and time dependent problems are considered for a one di- count (one-way coupling). With this assumption the con- mensional scalar equation defined on a finite bar and for a tinuous phase can be treated separate from the particulate 32 AN10 Abstracts two dimensional system defined on a square. The related Dept of Chemistry and Mathematics peridynamic operators and associated functional spaces are [email protected] analyzed. Applications to the analysis of the finite dimen- sional approximations to peridynamic models are also il- lustrated. CP26 A Spatiotemporal Model of Barley Yellow Dwarf Qiang Du Virus in Competing Plant Species with Seasonality Penn State University and Age Structure Department of Mathematics [email protected] Barley yellow dwarf virus infects over 100 grass species in agricultural and natural systems, mitigating success of in- Kun Zhou vasive species. We model BYDV transmission within a Penn State University metapopulation framework incorporating the movement of [email protected] aphid vectors between discrete host patches. We exam- ine BYDV persistence and alteration of pathogen-mediated interactions between perennial and annual competitors CP26 across a range of host community configurations. Sensi- A Differential Equation/agent-Based Hy- tivity analysis of the basic reproduction number for two brid Model of Antibiotic-Resistant Infection in a simplified models identifies how key parameters influence Hospital Ward pathogen invasion. When modeling the spread of bacterial infections in a popu- Carrie A. Manore, Sean Moore, Vrushali A. Bokil lation, dynamics occur on two levels: Interaction between Oregon State University the people in the ward, and the bacterial growth taking [email protected], place inside infected ward members. We propose a method [email protected], for combining both into a single model, where each person [email protected] is treated as an agent in an agent-based system, with the agent-specific dynamic parameters driven by a differential Elizabeth Borer equations system that models the within-host bacterial dy- University of Minnesota namics for that agent. Ecology, Evolution, and Behavior [email protected] Lester Caudill University of Richmond [email protected] CP26 Control of Mosquito Populations Through the In- troduction of Sterile Mosquitoes CP26 Optimal Efficacy of Ribavirin in the Treatment of We construct an ODE model for the release of sterile Hepatitis C mosquitoes into a population of wild mosquitoes. Assum- ing Holling-II dynamics, we examine the optimal control A mathematical model to represent hepatitis C is pre- of the wild mosquito population through several objec- sented. An objective functional keeping biomedical goals in tive function formulations. Numerical simulations are dis- mind is formulated. The optimal control problem is solved cussed. numerically to obtain an optimal efficacy of ribavirin in a combination treatment of ribavirin with interferon, where Maeve L. McCarthy, K. Renee Fister the efficacy of the latter has a clinically validated functional Murray State University form [email protected], renee.fi[email protected] Siddhartha P. Chakrabarty Dept. of Mathematics Seth F. Oppenheimer Indian Institute of Technology Guwahati Mississippi State University [email protected] [email protected] CP26 CP26 Drug Resistance Effects on Weighted Optimal Con- PhaseLockingBetweentheExpressionProfilesof trol for Chemotherapy Transcription Factors and Their Target Genes Dur- A differential equations model for cancer and its treatment ing the Yeast Cell Cycle is examined. The governing equations are developed from Increasing evidence now suggest the role of frequency mod- a four-compartmental model that shows the development ulation in gene expression. We examined this question of resistance over time. Optimal control techniques are using the time course gene expression data of yeast cell then used to investigate the effects of drug resistance (both cycle. We found that between transcription factors and natural and induced) on cell-cycle specific chemotherapy their binding targets, the expressions are significantly more regimes. Optimal strategies are suggested that take into likely to be phase locked than random gene pairs. In ad- account factors such as toxicity level. Once the existence of dition, ∼ 5–10% show higher order, more than 1 : 1, phase the optimal control is established, the control is then char- locking. acterized through the Hamiltonian and the adjoint system to governing equations. Xujing Wang University of Alabama at Birmingham Daniel L. Kern [email protected] Florida Gulf Coast University AN10 Abstracts 33 CP27 pose a generic model of a reaction-diffusion system in the Quantitative Dynamics of Spiral Waves in Active presence of a morphogen gradient. This morphogen gra- Media. dient is established independently of the reaction-diffusion system and acts by increasing the production of the acti- If perturbed, spiral waves in active media slowly move their vator proportional to the morphogen concentration. The core. In linear approximation, the drift velocity is pro- model is motivated by several existing models in develop- portional to convolution of perturbation with translational mental biology in which a Turing patterning mechanism response function (RF), the eigenfunction of adjoint lin- is proposed and various chemical gradients are known to earized operator corresponding to translational symmetry. exist. We investigate how the Turing pattern is affected, For variety of perturbations, we compute drift velocities if it exists. We also apply our results to a model of bone using RFs in Barkley and FitzHugh-Nagumo models and pattern formation in vertebrate limbs and show how they show good quantitative agreement with direct numerical may shed light on some experimental findings. simulations. We also demonstrate orbital movement of spi- ral due to interaction with localized inhomogeneities. Tilmann Glimm Western Washington University Irina V. Biktasheva [email protected] University of Liverpool [email protected] CP27 Dwight Barkley Spatial Modeling of Bacterial Antibiotic Tolerance Mathematics Institute, The University of Warwick, UK Persisters are phenotypic bacterial variants which exhibit [email protected] ”antibiotic tolerance”: persisters survive when exposed to bactericidal antibiotics, but the bacterial colonies re- grown from persister cells remain susceptible to killing by Vadim N. Biktashev further antibiotic treatments.Persister formation has been Dept of Applied Maths, University of Liverpool, linked to quorum sensing, the ability of bacteria to deter- UK mine their local population density, and with the formation [email protected] of biofilms, surface-growing bacterial communities particu- larly difficult to eradicate with antibiotic treatments. Ex- Andrew J. Foulkes tending a previous model of bacterial colony pattern for- Dept of Computer Science, mation, we present a reaction-diffusion model of bacterial University of Liverpool, UK growth which includes density-dependent persister forma- [email protected] tion. The model also exhibits separation of temporal and spatial scales between persister and regular bacteria. We investigate the pattern formation properties of the model CP27 as a number of experimentally accesible parameters are Complex Media, Complex Fluids and Complex varied. Flows William E. Sherwood Flows involving simple fluids are relatively well understood Center for Biodynamics in comparison to complex fluids such as blood, petroleum, Boston University jam, honey, just to cite a few. In fact, many of the biologic, [email protected] edible, and chemical fluids are complex in nature. Many such fluids are product of nature and many are man made John Burke by design in order to achieve some specific objectives. For Boston University example, complex fluid such as aqueous solution containing [email protected] polymer and surfactants are used in improved oil recovery processes.It is necessary to understand at the fundamen- tal level the origin of many physical phenomena that can CP27 be triggered by the complex rheological properties of these Mathematical Analysis of a System of Pde Arising complex fluids, in particular in complex media. Such un- From Epidermal Wound Healing derstanding can lead to better design of many processes and better numerical methods to simulate these processes. Systems of reaction-diffusion equations have been proposed to study cell migration across the surface of an epidermal Prabir Daripa wound. Approximation of traveling waves which describe Texas A&M University the cell migration into the wound has been studied by a Department of Mathematics number of researchers. We will investigate asymptotic be- [email protected] havior and speed of the travel wave solutions. Analysis of the model system reveals biological interactions in the CP27 wound healing process and the role of the model parame- ters in determining the speed. Turing Patterns with External Morphogen Gradi- ents Haiyan Wang Arizona State University In embryonic development, cells self-organize into complex [email protected] distributions giving rise to tissues and organs. The basic mechanisms of how this happens are still not very well un- derstood. One experimentally verified mechanism is the action of gradients of signalling molecules. Another one is the Turing mechanism in reaction-diffusion systems. To in- vestigate how these two mechanisms may interact, we pro- 34 AN10 Abstracts CP27 tors Analysis of Equations with Nonlinear Dispersion Anetworkofn coupled spherical-cap droplet oscillators is At this time, there is little existence theory for equations considered. The coupling by inviscid flow between droplets in which the mechanism which generates dispersion is itself is modeled as a system of n − 1 second-order differential nonlinear. In this talk, I will present a new equation which equations under a constant-volume constraint. For Sn sym- possesses compactly supported traveling wave solutions (a metric coupling, analytic stability and bifurcation results hallmark of degenerate dispersive equations) and for which are obtained for an arbitrary number of droplets. Nonlin- the Cauchy problem has a weak solution. This work is joint ear dynamics are explored numerically in both undamped with D. Ambrose. and weakly-damped cases. In addition, the influence of a variety of coupling laws is analyzed. J. Douglas Wright Drexel University David M. Slater Department of Mathematics Cornell University [email protected] [email protected] Paul Steen CP28 School of Chemical Engineering Mathematical Analysis of Swarming: Integrat- Cornell University ing Agent-Based Simulation with Population-Level [email protected] Modeling and Analysis An agent-based model is used to inform the design of a CP28 biologically relevant continuum model for swarming. Lin- Robustness and Efficiency of Ant-Based Routing ear stability analysis reveals the existence of low-frequency and Forwarding Protocols. instabilities as key behavioral parameters are varied. The results of the stability analysis are verified via numerical Ant-based protocols provide elegant solutions to routing simulation of the continuum model. Simulations of the and forwarding problems. These protocols use control agent-based model demonstrate a correspondence for finite packets or ”ants” to discover new routes and reinforce ef- swarms and are used to evaluate the effects of noise and ficient routes by adjusting pheromone values. We revisit changes in communication on stability. the new nonlinear stochastic models and compare them to network simulations with high-fidelity protocol, and phys- Allison Kolpas ical layer models. Through analysis of the models, we can Department of Mathematical Sciences understand how network behavior depends upon routing University of Delaware exponents and pheromone deposition/evaporation param- [email protected] eters. Jennifer Miller, Joao Juchem Neto, Louis F. Rossi Claudio Torres, Louis F. Rossi, Chien-Chung Shen, Ke Li, University of Delaware Jeremy Keffer, Rui Fang [email protected], [email protected], University of Delaware [email protected] [email protected], [email protected], [email protected], [email protected], jeremy@yekeffer.net, [email protected] CP28 Mathematical Modeling and Analysis of a Contin- uum Model for Three-Zone Swarming Behavior CP28 Dynamics of Light Beam Switching at the Interface Swarms in nature have been modeled with individuals but of Two Nonlinear Optical Media a continuum model, in which individuals are replaced with fields, may be better suited to theoretical analysis and scal- The dynamics of two light beams interacting at an interface ing up to larger swarms. Models including zones of repul- separating tow nonlinear optical media is studied numeri- sion, orientation, and attraction are popular in ecological cally and analytically. The beam dynamics is simulated by modeling. Our integro-differential equations describe the the beam propagation method. The numerical simulations reactions to the varying density in these three zones. We agree with the results of analytical model. The trapped use linear stability analysis to explore first and second or- beam at the interface acts as a power controllable switch der models of constant density swarms. to reflect or transmit the incident beam at the interface. Some new results will be presented in this talk. Jennifer Miller, Louis F. Rossi University of Delaware Rajah P. Varatharajah [email protected], [email protected] North Carolina A &T State University Allison Kolpas [email protected] Department of Mathematical Sciences University of Delaware CP28 [email protected] High-Order Bisection Method for Computing In- variant Manifolds of 2-D Maps CP28 We describe an efficient and accurate numerical method for Nonlinear Dynamics of Coupled Droplet Oscilla- computing smooth approximations to invariant manifolds of planar maps, based on geometric modeling ideas from Computer Aided Geometric Design (CAGD). The unstable AN10 Abstracts 35 manifold of a hyperbolic fixed point is modeled by B´ezier we explore both novel and traditional supporting compo- interpolant (a Catmull-Rom spline) and properties of such nents of high-order linear and nonlinear wave solvers for curves are used to define a rule for adaptively adding points their suitability for modern massively parallel computer ar- to ensure that the approximation resolves the manifold to chitectures. Components examined range from time step- within a specified tolerance. Numerical tests on a variety pers to preconditioners, linear solvers and shock capturing of example mappings demonstrate that the new method schemes, for which we examine design considerations, al- produces a manifold of a given accuracy with far fewer gorithms, and performance data. calls to the map, compared with previous methods. A brief introduction to the relevant ideas from CAGD is provided. Andreas Kloeckner Additionally, we will discuss the difficulties of extending Brown University these ideas to higher dimensions. [email protected] Jacek Wrobel, Roy H. Goodman New Jersey Institute of Technology MS1 Department of Mathematical Sciences An Accelerated Discontinuous Galerkin Method [email protected], [email protected] for Time-Domain Wave Propagation on Curvilin- ear Domains MS1 We will introduce the low storage curvilinear discontinu- Absorbing Boundary Condition Taking into Ac- ous Galerkin (DG) method for solving time-dependent par- count the Grazing Modes tial differential equations on curvilinear domains. We will present a proof of stability and show results that indicate We propose a new Absorbing Boundary Condition for the the method is typically as accurate as standard DG dis- wave equation, which is adapted to arbitrarily-shaped con- cretizations on curvilinear domains. Finally we will discuss vex boundary. After having decomposed the exact solution how this approach has enabled three dimensional electro- into propagating, evanescent and grazing fields, this con- magnetic scattering and also fluid flows to be computed on dition is obtained by combining the approximations of the a workstation using commodity graphics processing units. transparent condition in the three corresponding regions. The absorption of the grazing field involves a fractional derivative in space which requires a specific discretization Tim Warburton in order to implement the condition in a finite element Computational and Applied Mathematics scheme. Rice University [email protected] H´el`ene Barucq, Julien Diaz Team-Project Magique-3D INRIA Bordeaux Sud-Ouest MS2 [email protected], [email protected] Model Reduction of Large Scale System in Oil In- dustries V´eronique Duprat LMA, CNRS UMR 5142 Large-scale numerical models for suburface hydrocarbon Universit flow are frequently used during the design phase of oil ’e de Pau et des Pays de l’Adou fields. Moreover, an emerging technique is the use of such [email protected] models during the operational phase. To speed-up simula- tion of these models there appears to be a large scope for simplifying them using model-order reduction techniques. MS1 In particular we present an example of the use of Proper A Orthogonal Decomposition applied to a simple reservoir High Order Runge-Kutta Discontinuous Galerkin model for two-phase (oil-water) flow. Scheme with Time Accurate Local Time Stepping Patricia Astrid A new explicit RK based method for the time integration Shell Projects & Technology of the DG scheme is presented. The considered class of RK The Netherlands schemes have the property that a time polynomial of the [email protected] solution is available, which can be used to introduce local time stepping. Despite the local time steps, the scheme is Jan-Dirk Jansen high order accurate in time, conservative and ideally suited Delft University of Technology for massively parallel computations, allowing an efficient The Netherlands discretization of large scale time dependent problems. [email protected] Gregor Gassner Institut fuer Aerodynamik und Gasdynamik MS2 [email protected] Structure-preserving Model Reduction for MEMS Modeling MS1 Abstract unavailable at time of publication. Machine-adapted Methods: High-order Wave Propagation on GPUs David Bindel Cornell University Having recently shown that high-order unstructured dis- [email protected] continuous Galerkin methods are a spatial discretization that is well-matched to execution on GPUs, in this talk 36 AN10 Abstracts MS2 methods. Balanced Truncation for Switched Dynamical Sys- tems Amit Hochman, Jacob White Research Lab. of Electronics A hybrid dynamical system is a system described by MIT both differential equations (continuous flows) and differ- [email protected], [email protected] ence equations (discrete transitions). It has the benefit of allowing more flexible modeling of dynamic phenomena, including physical systems with impact such as a bouncing MS3 ball, switched systems such as a thermostat, and even in- Stochastic and Kinetic Approaches in Mesoscale ternet congestion. Hybrid dynamical systems pose a chal- Modeling of Grain Growth in Polycrystalline Ma- lenge since almost all reduction methods cannot be di- terials rectly applied. Here we introduce a balanced truncation- like method for switched linear dynamical systems. For Microstructure evolution is a problem of central impor- this we define the controllability and observability Grami- tance in materials science. Grain boundary engineering ans for either time-transmission switching or time-invariant techniques allow to design and manufacture materials pre- switching. For the first case, the Gramians are only the disposed to certain behavior under external conditions. sum of the Gramians of the subsystems, and so the bal- This type of work relies heavily on statistical analysis anced truncated switched system is the sum of the balanced of grain boundary evolution and of an interplay between truncated subsystems. For the second case the method is micro- and macroscopic properties. We present an account tricky as a notion of uncertainty should be added in order of recent developments in mesoscale theory of texture evo- to deal with the special structure of the switched system. lution in polycrystalline materials. The talk will demon- We also introduce two new definitions of the Gramians, strate how a carefully chosen combination of simulation, solutions of either special LMIs or Lyapunov equations to theory and modeling techniques can bring a deeper under- come up with the balanced truncated-like reduced model. standing to coarsening processes governed by the motion With this approach we will preserve the stability for the of triple junctions. By means of modulated stochastic pro- reduced model. Here one should notice that we supposed cess characterization and kinetic evolution equations, we implicitly that each subsystem is stable. aim at quantifying the rates of grain boundary disappear- ance events and predicting the influence of these events on Younes Chahlaoui macro-level materials properties. Comparison with large- CICADA scale 2D and 3D simulations will provide the necessary vali- The University of Manchester, UK dation and reveal similarities and differences between these [email protected] approaches. Maria Emelianenko MS2 George Mason University Interpolatory Projection Methods for Optimal [email protected] Model Reduction of Large-scale Dynamical Sys- tems MS3 We investigate interpolatory optimal model reduction of Predictive Theory for the Grain Boundary Char- large-scale dynamical systems. First, we present a trust- acter Distribution region approach for optimal H∈ approximation of MIMO systems. The approach generates a sequence of reduced or- Most technologically useful materials are polycrystalline microstructures composed of a myriad of small grains sepa- der models producing monotone improving H∈ error norms and is globally convergent to a reduced order model guar- rated by grain boundaries. The energetics and connectivity anteed to satisfy first-order optimality conditions. Then, of the grain boundary network plays a crucial role in deter- we extend this approach to optimal model reduction of pa- mining the properties of a material across a wide range of scales. A central problem is to develop technologies capable rameterized systems using an H∈ ⊗L∈ joint error measure. of producing an arrangement of grains—a texture—that provides for a desired set of material properties. Since the Serkan Gugercin dawn of man, coarsening has been the principal feature of Virginia Tech. these technologies. Cellular structures coarsen according to Department of Mathematics a local evolution law, a gradient flow or curvature driven [email protected] growth, for example, limited by space filling constraints, which give rise to changes in the configuration. There are two aspects of coarsening, geometric growth and texture MS2 development. The main objective of this presentation is Stability Issues and DEIM-based Nonlinear Model to show that they may be decoupled and characterized by Reduction different types of evolution processes, leading to statisti- cally different types of coarsening rates. We consider the Reduced order model stability is examined for the discrete grain boundary character distribution, the GBCD, a ba- empirical interpolation based model reduction method sic texture measure, and establish an entropy based theory (DEIM). A micro-electromechanical and a circuit exam- for it which suggests that GBCD satisfies Fokker-Planck ple are used to demonstrate the stability issue, and type kinetics. For this, a simplified critical event model is stability-enhancing modifications of the algorithm are de- introduced and studied. scribed. We show that extracting an equilibrium point Jacobian guarantees local DIEM stability, and that multi- Yekaterina Epshteyn, David Kinderlehrer, Katayun point methods for extracting stabilizing projection matri- Barmak ces yields superior global stability. The examples are also Carnegie Mellon University used to demonstrate the failure characteristics of DEIM [email protected], [email protected], AN10 Abstracts 37 [email protected] [email protected] Maria Emelianenko George Mason University MS4 [email protected] Modularity and Graph Algorithms A number of graph partitioning algorithms are based on Eva Eggeling the concept of modularity. In particular Clauset, Newman Fraunhofer InstitutComputer Graphics Research IGD and Moore (CNM) have developed a greedy agglomerative [email protected] graph partitioning algorithm that scales well but is known to have several flaws. Fortunato and Barthelemy have per- Shlomo Ta’asan formed a rigorous analysis of the CNM algorithm that elu- Department of Mathematical Sciences cidates it problems. More recently Berry, Hendrickson, Carnegie Mellon University Laviolette, and Phillips have derived a weighted variant of [email protected] CNM that performs much better in practice. This talk will focus on a different version of the parent CNM algorithm based on a statistical re-interpretation of CNM that also MS3 addresses some of the issues with the original algorithm. The Janossy Effect and Hybrid Variational Princi- ples Joe McCloskey National Security Agency Light changes the orientation of a liquid crystal. This is the [email protected] optical Freeder- icksz transition. In the Janossy effect, the threshold intensity for the optical Freedericksz transition David A. Bader is dramatically reduced by the additon of a small amount College of Computing of dye to the sample. This has been interpreted as an Georgia Institute of Technology optically pumped orientational rachet mechanism, similar [email protected] to the rachet mechanism in bi- ological molecular motors and requires an innovative hybrid gradient flow. Joint work with Michal Kowalczyk. MS4 David Kinderlehrer Scalable Methods for Representing, Characteriz- Carnegie Mellon University ing, and Generating Large Graphs [email protected] Abstract unavailable at time of publication. Ali Pinar MS3 Sandia National Labs Kinetic Theory and Lax Equations for Shock Clus- [email protected] tering and Burgers Turbulence We study shock statistics in scalar conservation laws with MS4 random initial data by deriving kinetic equations, which Exploiting Sparsity in the Statistical Analysis of have the form of a Lax pair. This, in addition to some Gene Expression Data exact solutions to the equation, suggest that the model is a completely integrable system in 2+1 dimensions. Abstract unavailable at time of publication. Ravi Srinivasan Padma Raghavan University of Texas at Austin The Pennsylvania State Univ. [email protected] Dept of Computer Science Engr. [email protected] Govind Menon Brown University Anirban Chatterjee, Francesca Chiaromonte [email protected] Pennsylvania State University [email protected], [email protected] MS4 People You May Know MS5 Multiscale Simulation Techniques for Fluid Struc- Facebook’s friend recommendation system helps people ture Interaction connect with their friends. Our system, called People You May Know, uses a combination of results from sociology In this talk, I will discuss multiscale simulation techniques and machine learning to make the best suggestions pos- for fluid-structure interaction. Our goal is to develop simu- sible. We will look at some of the challenges involved in lation techniques for solving fluid-structure problems with building a system that can handle the scale of Facebook complex microgeometries on a coarse grid. A particular ap- and provide high quality recommendations. In this talk plication is fluid flows in deformable fracture porous media. I will discuss both the algorithmic and machine learning Due to changes in pore pressure, the deformation can sig- challenges that we have faced and overcome in building nificantly alter pore geometry and change permeabilities. this system. Coarse-scale models will be discussed and numerical results will be presented. This is a joint work with Peter Popov, Lars Backstrom Yuliya Gorb, and Donald Brown. Facebook Palo Alto, CA Yalchin Efendiev 38 AN10 Abstracts Dept of Mathematics tracts key nonlinearities, thereby providing a basis for un- Texas A&M University derstanding competing motivations that drive investor de- [email protected] cisions. One of the motivations involved is trend, or mo- mentum, a key component of the system of ordinary differ- ential equations modeling price dynamics developed since MS5 1990 (see www.pitt.edu/∼caginalp). We show equilibrium On a Nonlinear Problem Arising in Poroelasticity and stability results for this system. In this work we study a model for steady fluid flow in a Mark DeSantis saturated deformable (elastic) porous medium. The model Mathematics which was developed in the context of modeling one of the Univ of Pittsburgh processes in paper production is distinguished in its use of [email protected] a Kozeny-Carman type relation for the fluid’s permeabil- ity, resulting in a quasilinear elliptic system. We describe and analyze a finite element method for approximating so- MS6 lutions of this nonlinear problem. TheDynamicsofaDurableConsumptionGood: Modeling for Experimental Design and Testing Yanzhao Cao Department of Mathematics & Statistics An unprecedented housing bubble, 1997 to 2007 subse- Auburn University quently engulfed the economy from 2007 to the present. [email protected] Economic theory and experiments have modeled the two polar extremes: (1) Services and other non-durables that A. J. Meir perish on consumption, and (2) Durable goods with infinite Auburn University life. We build on the previous literature to model a simple [email protected] consumer durable asset, that has a finite expected life, with agent producers who can vary their output of the durable asset depending on its market price and production cost. MS5 We characterize the long run stationary equilibrium of the Coupling of Elastic and Electromagnetic Waves and system, and show how it can be combined with existing Related Mathematical Problems models of asset pricing to enable the study of short-term dynamic behavior. Part 2. There are considered theoretical and numerical aspects of the interaction of elastic and electromagnetic waves in com- Steven Gjerstad plex media. The models considered are different combina- Chapman University tions of the Lam´e and Maxwell equations. The coupling [email protected] mechanism is based on the seismomagnetic and electroki- netic effects. Theoretical results of the analytical solution are discussed for various direct and inverse problems. Fi- MS6 nally, we give some results of the numerical solution of some On the Specialization and Division of Exchange, direct and inverse problems describing linear processes of Durability and Market Performance interaction of elastic and electromagnetic waves. Robust laboratory studies have long established contrast- Viatcheslav I. Priimenko, Mikhail Vishnevsky ingly different behavior depending whether the items ex- North Fluminense State University Darcy Ribeiro changed are final use goods, with specialized buyers, and [email protected], [email protected], sellers, or assets that can be re-traded over multiple peri- [email protected] ods. The former perform with high efficiency, the latter poorly, tending to bubble and crash. New experiments combine these features in simple markets where durability, MS5 and specialization, is varied by imposing or relaxing a re- Fluid Flows in Karst Aquifer striction on re-trade within each period. Since about 60 percent of national output is consumer non-durable goods We discuss two types of models for flows in karst aquifer. and services, a component that resists change in good times The first is the coupled Navier-Stokes-Darcy model and its or recessions, while economic troubles originate in asset variants. The second is the so-called coupled continuum markets, the results illuminate the much greater stability pipe flow model (CCPF) or conduit flow process model of expenditure patterns in perishable final goods markets (CFP). Mathematical well-posedness of the models, their relative to durable goods, credit and financial markets. numerical approximation, calibration of key physical pa- rameters will be discussed. David Porter Economics and Mathematics Xiaoming Wang Chapman University Department of Mathematics [email protected] Florida State University [email protected] MS6 Housing Bubbles and the Economy: Challenge to MS6 Modeling Nonlinearity and Instability in the Dynamics of Fi- nancial Markets An unprecedented housing bubble, 1997 to 2007 subse- quently engulfed the economy from 2007 to the present. Two recent studies in collaboration with Prof. Caginalp are Economic theory and experiments have modeled the two presented. A large-scale statistical methodology demon- polar extremes: (1) Services and other non-durables that strates a significant role for numerous variables, and ex- perish on consumption, and (2) Durable goods with infinite AN10 Abstracts 39 life. We build on the previous literature to model a simple MS7 consumer durable asset, that has a finite expected life, with Topics in Kinetic Simulation: A Brief Look at the agent producers who can vary their output of the durable Landscape asset depending on its market price and production cost. We characterize the long run stationary equilibrium of the This talk will be a brief survey of numerical methods used system, and show how it can be combined with existing to solve kinetic equations and the current challenges faced models of asset pricing to enable the study of short-term by computational scientists in various applications. This dynamic behavior. Part 1. includes an overview of topics to be addressed in more de- tail during the symposium. Vernon Smith Economics and Law Cory Hauck Chapman University Oak Ridge National Laboratory [email protected] [email protected] MS6 MS7 A Profit Keeping and Loss Protection Strategy in A Nonlinear Spatial Closure for the 2-D Sn Equa- Portfolio Management tions The widely studied power-utility-maximization strategy is We present a strictly positive non-linear Petrov-Galerkin optimal only in the sense of ensemble averaging. However, method that reduces to the standard linear-discontinuous in reality, only one random path will be realized and the Galerkin spatial discretization when that method produces value of the portfolio at the end of the investment horizon a positive solution. In addition to strictly positive angu- could be dramatically lower than its historical high. This lar flux solutions, our method offers preservation of the P1 is evident in the recent financial crisis. We will present a spatial moments of the transport equation. Computational new strategy to overcome this problem. results for several test problemsinbothslabandrectan- gular geometries are presented. Qiang Zhang Mathematics Jim Morel, Peter Maginot City Univ of Hong Kong Texas A&M [email protected] [email protected], [email protected] Jean Ragusa MS7 Department of Nuclear Engineering Positivity-preserving Discontinu- Texas A&M University ous Galerkin Schemes for Linear Vlasov-Boltzmann [email protected] Transport Equations We develop a high-order positivity-preserving discontinu- MS7 ous Galerkin (DG) schemes for linear Vlasov-Boltzmann Conservative High Order Semi-Lagrangian Finite transport equations (BTE) under the action of quadrati- Difference Weno Methods for Advection in Incom- cally confined electrostatic potentials. The solutions of the pressible Flow and Kinetic Equations BTEs are positive probability distribution functions. It is very challenging to have a mass-conservative, high-order We propose a novel semi-Lagrangian finite difference for- accurate scheme that preserves positivity of the numeri- mulation for approximating conservative form of advec- cal solutions in high dimensions. Our work extends the tion equations with general variable coefficients. Com- maximum-principle-satisfying scheme for scalar conserva- pared with the traditional semi-Lagrangian finite differ- tion laws by Zhang and Shu (2010) to include the linear ence schemes, which are designed by approximating the Boltzmann collision term. The DG schemes we developed advective form of equation via direct characteristics trac- conserve mass and preserve the positivity of the solution ing, the scheme we proposed approximates the conservative without sacrificing accuracy. A discussion of the standard form of equation. This essential difference makes the pro- semi-discrete DG schemes for the BTE are included as a posed scheme conservative by nature, and extendable to foundation for the stability and error estimates for this equations with variable coefficients. The proposed semi- new scheme. Potential applications of this solver include Lagrangian finite difference framework is coupled with high the modeling of charge transport in nano-scale semiconduc- order essentially non-oscillatory (ENO) or weighted ENO tor and solar devices and Vlasov-Poisson, Vlasov-Maxwell (WENO) reconstructions to achieve high order accuracy in systems in plasma physics. smooth parts of the solution and capture sharp interfaces without introducing oscillations. The scheme is extended Yingda Cheng to high dimensional problem by Strang splitting. The per- Dept. of Math and ICES formance of the proposed schemes is demonstrated by lin- University of Texas at Austin ear advection, several challenging examples of rigid body [email protected] rotation and swirling deformation in multi-dimensions. As the information is propagating along characteristics, the Irene M. Gamba semi-Lagrangian scheme does not have CFL time step re- Department of Mathematics and ICES striction, allowing for a cheaper and more flexible numeri- University of Texas cal realization than the regular finite difference scheme for [email protected] some problems. Jennifer Proft Jing-Mei Qiu University of Texas at Austin Mathematical and Computer Sciences [email protected] Colorado School of Mines 40 AN10 Abstracts [email protected] ties of the method, as well as numerical examples illustrat- ing its high-order convergence properties. Chi-Wang Shu Division of Applied Mathematics Jia Li, Dongbin Xiu Brown University Purdue University [email protected] [email protected], [email protected] MS8 MS8 Anchor Point in ANOVA Decomposition with Ap- A Stochastic Heterogeneous Multiscale Method for plications to Stochastic Problems Porous Media Flow Analysis of variance (ANOVA) is presented to lift the curse Abstract not available at time of publication. of dimensionality for high dimensional problems. We focus Nicholas Zabaras,XiangMa on anchored-ANOVA representation employing the Dirac Cornell University measure that converges exponentially for certain classes [email protected], [email protected] of functions. However, the error depends on the anchor points. We investigate which anchor points give the best approximation. We then present examples of a function ap- MS9 proximation as well as numerical solutions of the stochas- Krylov Subspace Enhanced Parareal Algorithms tic advection equation using a combination of anchored- ANOVA and polynomial chaos expansions. We are interested in a variant of the parareal algorithm in- troduced by Farhat et al., which uses previous iterations in Minseok Choi, Zhongqiang Zhang, George Karniadakis order to create a subspace in which a more accurate coarse Brown University approximation of the solution is sought at each iteration. minseok [email protected], zhongqiang [email protected], We show that this approach can be interpreted as a Krylov george [email protected] subspace method with an incomplete subspace: there are missing powers of the iteration matrix. This leads to a new MS8 and interesting polynomial approximation problem. A Sparse Approximation of Partial Differential Stefan Guettel Equations with Random Inputs Mathematics University of Geneva In this work, using techniques based on concentration of [email protected] measure theory we derive a method for approximation of sparse solution of stochastic PDEs. The proposed method is: Martin Gander Department of Mathematics • Non-intrusive: It is based on the direct solution sam- Universite’ de Geneve, Switzerland pling, [email protected] • Provably convergent: We obtain probabilistic bounds on the approximation error proving the convergence MS9 of the method, Parallel Deferred Corrections and Parareal • Well suited for high dimensional problems. A method for temporal parallelization which combines fea- Several numerical experiments will be presented to explore tures of the parareal algorithm and spectral deferred cor- different aspects of the proposed method. rection methods will be presented. This approach can be viewed either as a way to parallelize deferred corrections, or Alireza Doostan as a technique for improving the efficiency of the parareal Department of Aerospace Engineering Sciences algorithm. University of Colorado, Boulder [email protected] Michael Minion University of North Carolina-Chapel Hill Houman Owhadi Mathematics Department Applied Mathematics [email protected] Caltech [email protected] Jan Prins University of North Carolina Computer Science MS8 [email protected] Stochastic Collocation Methods on Arbitrary Nodes Alex Keng Computer Science We present a stochastic collocation strategy based on mul- University of North Carolina tivariate interpolation theory. The key feature of the [email protected] method is that it allows simulations on collocation nodes (sampling points) of arbitrary number and at arbitrary locations in arbitrary dimensional spaces. Therefore the method results in great flexibility for practical stochastic simulations, particularly in high dimensions. We present mathematical analysis on the well-posedness and proper- AN10 Abstracts 41 MS9 in studying colloids is that there are very few rigorous re- Parallel Implicit Deferred Corrections sults and it even appears difficult to efficiently simulate these systems to test hypotheses experimentally. We study We present an extension of revisionist integral deferred cor- several discrete models of colloids with large squares sus- rection, a parallel high order time integrator, to solve stiff pended in a sea of small diamonds and rigorously show that systems implicitly. Based on the corrector-predictor model, the squares will cluster together beyond a certain density, the main idea is to decouple the prediction and correction extending a result of Frenkel and Louis. We also study steps so that they can be run in parallel. Time permitting, heuristics due to Buhat, Dress and Krauth for sampling we apply our algorithms to stiff rate equations. configurations of colloids and show that their algorithms are provably efficient in certain settings and inefficient in Andrew J. Christlieb others. (This is joint work with Sarah Miracle and Amanda Michigan State Univerity Pascoe.) Department of Mathematics [email protected] Dana Randall Georgia Institute of Technology Ben Ong [email protected] michigan state university [email protected] MS10 Phase Transitions for Random Graph Processes MS9 A Parallel in Time Solver for Optimal Control The classic random graph evolution of Paul Erd˝os and Al- fred R´enyi undergoes a phase transition. Parametrizing the Problems n number of edges as t 2 there is a critical point tcr =1.For In this talk, we present a general methodology that cou- t MS10 This work attempts to apply radial basis collocation Discrete Models of Colloids method to fracture mechanics by introducing a domain decomposition technique with proper interface conditions. Colloids are mixtures of substances that tend to separate The proposed method allows (1) natural representation of so that like particles are clustered together. Unlike re- discontinuity across the crack surfaces and (2) enrichment lated models with enthalpic forces causing like particles to of crack-tip solution in a local subdomain. Exponential attract or disparate particles to repel, colloids are driven convergence rate can be achieved while keeping the dis- solely by entropic forces. It seems that in these models crete system well-conditioned. The analytical prediction clustering occurs above a certain density if configurations and numerical results demonstrate that an optimal dimen- are chosen from the uniform distribution. The challenge 42 AN10 Abstracts sion of the near-tip subdomain exists. time dependent expansion coefficients. Using rotational and translational transformations, the Euler equations are Jiun-Shyan Chen solved as exact time dependent differential equations in a University of California, Los Angeles (UCLA) moving frame. Riemann solvers are used to calculate vari- Dept. of Civil & Environmental Eng. ous wave interactions. The eigenvalues of the time amplifi- [email protected] cation matrix are all unity. No artificial viscosity, upwind- ing or entropy violating methods were needed. Lihua Wang Tongji University, China Edward Kansa [email protected] University of California, Davis (UC Davis) [email protected] or [email protected] Hsin-Yun Hu Tunghai University, Taiwan MS12 [email protected] On the Stochastic Navier-Stokes Equations Abstract unavailable at time of publication. MS11 The Gibbs Phenomenon for Radial Basis Functions Alexander Labovsky Department of Scientific Computing The Gibbs phenomenon was first observed in the con- Florida State University text of truncated Fourier expansions, but it arises also in [email protected] cases such as truncated integral transforms and different interpolation methods. Radial basis functions (RBFs) in- clude both splines and trigonometric interpolation as spe- MS12 cial cases in 1-D, and they generalize these methodologies Filter Based Stabilization for Evolution Equations to scattered node layouts in any number of dimensions. We investigate here the Gibbs phenomenon in the case of RBF This project, which describes joint work with Vince Ervin interpolation on an equispaced grid in 1-D. (Clemson) and Monika Neda (UNLV), considers the follow- ing idea: Step 1: Take one step of any numerical method for Bengt Fornberg advancing an evolution equation in time, Step 2: take the University of Colorado result of Step 1, Filter it and use that for the following step. Boulder We show how this process Evolve then Filter can overdif- [email protected] fuse an approximate solution. We then show how to fix it to produce a highly accurate, stable and non-overdiffused Natasha Flyer approximation. One advantage of this approach to stabi- National Center for Atmospheric Research lization is that Steps 1 and 2 can be performed by calls to Division of Scientific Computing independent, black box programs. fl[email protected] William Layton University of Pittsburgh MS11 [email protected] An RBF-Gegenbauer Method for Time Dependent Discontinuous Problems MS12 Radial basis function(RBF) methods have been actively Analysis of Generalized Alpha-models of Turbu- developed in the last decades. The advantages of RBF lence methods are that these methods are mesh-free and yield high order accuracy if the function is smooth enough. The In this talk I will discuss the global existence and unique- RBF approximation for discontinuous problems, however, ness results for the general family of regularized Navier- deteriorates its high order accuracy due to the Gibbs phe- Stokes and Magnetohydrodynamics (MHD) models . This nomenon. In this talk we will show that high order accu- family captures most of the specific regularized models racy can be recovered from the RBF approximation con- that have been proposed and analyzed in the literature, in- taminated with the Gibbs phenomenon if a proper recon- cluding the Navier-Stokes equations, the Navier-Stokes-α struction method is applied. In this work, the Gegenbauer model, the Leray-α model, the Modified Leray-α model, reconstruction method is used to reconstruct the RBF ap- the Simplified Bardina model, the Navier-Stokes-Voight proximation for the recovery of high order accuracy. Nu- model, the Navier-Stokes-α-like models, and certain MHD merical examples including the linear and nonlinear hyper- models, in addition to representing a larger 3-parameter bolic partial differential equations are presented. family of models not previously analyzed. We give a unified analysis of the entire three-parameter family of models us- Sigal Gottlieb ing only abstract mapping properties of the principal dissi- Department of Mathematics pation and smoothing operators, and then use assumptions University of Massachusetts at Dartmouth about the specific form of the parameterizations, leading to [email protected] specific models, only when necessary to obtain the sharpest results. MS11 Evelyn Lunasin Meshfree Solutions of One-dimensional Inviscid University of California San Diego Euler Equations [email protected] The inviscid Euler equations are solved by expanding any dependent variable in terms of radial basis functions and AN10 Abstracts 43 MS12 Oscar Bruno Drag and Lift Computations of High Order Accu- Division of Engineering and Applied Science racy Fluid Flow Models California Institute of Technology [email protected] Numerical studies based on finite element discretizations of a higher order time relaxation models of fluid motion will be presented. The family of models is based on filtering MS13 and deconvolution techniques for driving the unresolved Complete Radiation Boundary Conditions: Theory fluctuations to zero. Optimal error estimates for velocity and Application convergence analysis are derived. The family of models is tested on a benchmark problem where drag and lift esti- Complete radiation boundary conditions (CRBCs) are ge- mates and pressure drop on a body immersed in a fluid ometrically flexible, applicable to all standard wave equa- were investigated. tions in homogeneous, isotropic media as well as to a lim- ited class of anisotropic models, and satisfy optimal com- Monika Neda plexity estimates. They provide provably accurate long- University of Nevada Las Vegas time results on small domains using parameters which are [email protected] determinable a priori. We review the CRBC theory and describe our ongoing efforts to develop a generally usable implementation and to extend their applicability to more MS12 complex models. Solving the Navier-Stokes Equations for Velocity, Vorticity and Helicity Thomas M. Hagstrom Southern Methodist University For the three-dimensional incompressible Navier-Stokes Department of Mathematics equations, we present a formulation featuring velocity, vor- [email protected] ticity and helical density as independent variables. We show that mathematically the helical density can be ob- served as a Lagrange multiplier corresponding to the MS13 divergence-free constraint on the vorticity variable. As one Stable Fourier Based Embedding Methods for possible practical application of this formulation, we con- Wave Propagation on Complex Domains sider a time-splitting numerical scheme based on a simple alternating procedure between vorticity - helical density The FC-AD methodology provides efficient new PDE and velocity - Bernoulli pressure systems of equations. We solvers based on the use of a certain “Fourier continua- discuss the relation of the numerical method to some well tion’ (FC) method in conjunction with alternating direc- known turbulence models. tion (AD) strategies and can be applied directly to com- plex computational domains with unconditional stability Maxim Olshanskii and high-order accuracy. FC-AD techniques are uniquely Department of Mechanics and Mathematics suited for wave propagation problems as the Fourier basis Moscow State M.V.Lomonosov University, Moscow, essentially eliminates pollution error. New extensions in- Russia cluding those to handle Neumann- and mixed-type bound- [email protected] ary conditions will be presented. Mark Lyon MS13 University of New Hampshire High Order Accurate Discretization of the Wave [email protected] Equation on Structured Grids Enforcement of boundary conditions in high-order MS13 embedded-boundary methods that removes small-cell stiff- Boundary and Interface Conditions for Non-Linear ness for both Dirichlet and Neumann b.c. and guarantee Wave Propagation Problems single valued solutions for slender bodies is presented. The b.c. enforcement is paired with Taylor-series-time-stepping Most wave propagation problems are linear and the re- and a high-order spatial discretization to accurately solve lated boundary and interface conditions that lead to well- the wave equation. Three spatial discretizations are consid- posedness and stability are relatively well known. For ered: Pade (implicit), summation by parts finite difference non-linear problems that is not the case. We will discuss and Fourier continuation. Numerical experiments illustrat- boundary treatments where the non-linearity adds addi- ing the properties of the methods are presented. tional complications. Typical applications that will be dis- ussed are related to shocks in fluid flow, propagation of Daniel Appelo discontinuities in stochastic PDE’s and earthquake prob- Division of Engineering and Applied Science lems. California Institute of Technology [email protected] Jan Nordstrom University of Witwatersrand Anders Petersson South Africa Center for Applied Scientific Computing [email protected] Lawrence Livermore National Lab [email protected] MS14 Randomization in the Default Boundary Problem Nathan Albin California Institute of Technology We consider an inverse problem for the first hitting time [email protected] distribution of a Wiener process and its application for 44 AN10 Abstracts credit risk pricing. We demonstrate that randomization of transparent and robust methods for valuing and hedging the initial state of the process lead to analytically tractable structured credit portfolios. First-generation models, such solution. as the Gaussian copula-based methods have well docu- mented practical and theoretical limitations. In this paper, Alexander Kreinin we demonstrate the practical application of the weighted University of Toronto and Monte Carlo methodology to value and compute sensitivi- Algorithmics Inc. ties and risk statistics for ABS and CLO structures. The [email protected] model extends the full bottom-up approach of Rosen and Saunders (2009) for pricing bespoke CDOs to include can- Sebastian Jaimungal, Angel Valov cellability and stochastic LGDs in a natural way. The University of Toronto performance of the model is analyzed throughout a three- [email protected], [email protected] month period during the credit crisis in 2008. The model calibrates very well to observed prices for the CDXHY and LCDX indices, and provides stable implied distributions MS14 for the systematic factor. Furthermore, it gives robust, Modeling Default Correlation and Clustering: A consistent prices and sensitivities , even during this very Time Change Approach volatile period. We present a novel framework for modeling correlated de- Dan Rosen faults where it is possible to capture the so-called “con- R2 Financial Technologies tagion effect’. Time changing multi-parameter Markov dan.rosen@r2-financial.com processes with multivariate subordinators leads to jump- diffusion processes that are correlated through their jump measures. When unpredictable shocks arrive, the default MS15 intensity of multiple firms will shift simultaneously, which Linear Feedback Control Strategies for the Navier- can trigger multiple defaults. Stokes Equations Rafael Mendoza-Arriaga One systematic control strategy often proposed for fluid University of Texas at Austin flows is optimal feedback control of the linearized system. [email protected] If the linearized portion of the Navier-Stokes equations are stabilized, then the full nonlinear system is stabilized. The Vadim Linetsky natural questions are practical-how does this strategy per- Northwestern University form in practice. This talk will focus on a numerical proce- [email protected] dure for solving the LQR control problem associated with the Oseen equations and numerical experiments to deter- mine the effectiveness of the feedback control strategy. MS14 Pricing Counterpary Risk at the Trade Level and Imran Akhtar CVA Allocation Virginia Tech [email protected] We address the problem of allocating the counterparty- level credit valuation adjustment (CVA) to the individual Jeff Borggaard trades composing the portfolio. We show that this problem Virginia Tech can be reduced to calculating contributions of the trades to Department of Mathematics the counterparty-level expected exposure (EE) conditional [email protected] on the counterparty’s default. We propose a methodology for calculating conditional EE contributions for both col- Miroslav Stoyanov lateralized and non-collateralized counterparties. Calcu- Florida State University lation of EE contributions can be easily incorporated into Department of Scientific Computing exposure simulation processes that already exist in a finan- [email protected] cial institution. We also derive closed-form expressions for EE contributions under the assumption that trade values are normally distributed. Analytical results are obtained MS15 for the case when the trade values and the counterparty’s Bridging the Boussinesq and Primitive Equations credit quality are independent as well as when there is a Through Spatio-Temporal Filtering dependence between them (wrong-way risk). We propose a novel approach for bridging the Boussinesq Michael Pykhtin equations and the primitive equations. This approach uses Federal Reserve Board spatio-temporal filtering as an alternative to traditional [email protected] scaling arguments used in the derivation of simplified mod- els for the ocean and atmosphere. Dan Rosen R2 Financial Technologies Traian Iliescu dan.rosen@r2-financial.com Department of Mathematics Virginia Tech [email protected] MS14 Pricing Structured Credit Products with Implied Jinqiao Duan Factor Models Illinois Institute of Technology [email protected] The current financial crisis has underscored the need for AN10 Abstracts 45 Paul F. Fischer ing Models Argonne National Laboratory fi[email protected] In this talk I’ll focus on fundamental properties of two spe- cific car-following (CF) models of highway traffic; namely Tamay Ozgokmen he Intelligent Driver Model, IDM, and the Aw-Rascle- University of Miami/RSMAS Greenberg Model, ARG. I shall show that for both of these [email protected] models there are no car crashes and no velocity rever- sals. These results concern well-posedness of the under- lying ODE system and are addressed by a-priori estimates MS15 which are independent of the number of cars in the system. Improving Mass Conservation in FE Computations I shall also establish that both the IDM and ARG models of NSE with C0 Velocities support large amplitude Stop and Go waves; these waves are traveling wave solutions of the underlying systems. This talk will begin by introducing the Scott-Vogelius el- ement pair, then show a connection between it and the James M. Greenberg grad-div stabilized Taylor-Hood element. Implications of Carnegie Mellon University this connection will be discussed, and numerical experi- Dept of Mathematical Sciences ments will also be given. [email protected] Michael Case Dept. of Mathematical Sciences MS16 Clemson University Impact of Protein and Lipid Binding on Drug Dis- [email protected] position Abstract unavailable at time of publication. Leo Rebholz Clemson University L. Peletier Department of Mathematical Sciences Mathematical Institute [email protected] Leiden University [email protected] Vince Ervin Department of Mathematical Sciences Clemson University MS16 [email protected] Solutions of Semilinear Equations Alexander Linke We study solutions of a second order semilinear pde, with WIAS Berlin particular focus on singular solutions. For both the fo- [email protected] cussing and non-focussing versions of the equation we find new singular solutions. MS15 William Troy University of Pittsburgh Two-Level Strategies for Large Eddy Simulation Department of Mathematics Proper Orthogonal Decomposition Models [email protected] One of the main hurdles in the development of modern clo- sure models in Proper Orthogonal Decomposition of tur- MS16 bulent flows has been the lack of efficient computational strategies for the discretization of the associated nonlinear- A Weak KAM Theorem for the Nonlinear Vlasov ities. We introduce a two-level discretization methodology System that reduces dramatically the computational cost without The space L2(0, 1) has a natural Riemannian structure on compromising the accuracy of the closure model. This is the basis of which we introduce an L2(0, 1)–infinite dimen- demonstrated numerically in the 3D simulation of a flow sional torus T. For a class of Hamiltonians defined on its past a cylinder at Re=1,000. cotangent bundle we establish existence of a viscosity so- Zhu Wang, Traian Iliescu lution for the cell problem on T or, equivalently, we prove Department of Mathematics a Weak KAM theorem. As an application, we obtain exis- Virginia Tech tence of absolute action-minimizing solutions of prescribed [email protected], [email protected] rotation number for the one-dimensional nonlinear Vlasov system with periodic potential. Imran Akhtar Wilfrid Gangbo Virginia Tech Georgia Institute of Technology [email protected] [email protected] Jeff Borggaard Adrian Tudorascu Virginia Tech Department of Mathematics Department of Mathematics University of Wisconsin, Madison [email protected] [email protected] MS16 MS17 Collective Behavior of Two Classes of Car Follow- Practical Heuristics for Inexact Subgraph Isomor- 46 AN10 Abstracts phism MS18 Taking the Road Not [Usually] Taken Abstract unavailable at time of publication. There are commonalities among research career paths; Jon Berry however, internally motivated and externally mandated Sandia Natinal Laboratories transitions help mold each experience into our own unique [email protected] journey. I will discuss major transition points in my profes- sional career ranging from graduate school to early career MS17 and personal as well as externally motivated events that brought me to career crossroads. My experience and pos- Spectral Methods for Subgraph Detection sible lessons learned in taking less traveled paths will be We describe a statistical test for subgraph detection and examined. localization using spectral properties of the so-called mod- Elebeoba May ularity matrix, a type of residual under the Chung-Lu ran- Sandia National Laboratories dom graph model. We show that the resultant algorith- [email protected] mic procedure can be applied to very large graphs (< 106 vertices), with complexity dominated by that of standard sparse eigensolver methods, and can successfully isolate MS18 anomalous vertices in real data examples. Panel Discussion: Success through Transitions Nadya Bliss, Benjamin A. Miller This panel session provides an opportunity for AWM work- MIT Lincoln Laboratory shop members to ask questions and discuss issues with pro- [email protected], [email protected] fessional development speakers. Patrick J. Wolfe Elebeoba May Harvard University Sandia National Laboratories [email protected] [email protected] Gigliola Staffilani MS17 MIT Tools and Primitives for High-performance Graph [email protected] Computation Abstract unavailable at time of publication. Mary Ann Horn National Science Foundaton John Gilbert [email protected] U. California at Santa Barbara [email protected] MS18 Title Not Available at Time of Publication MS17 Hybrid Parallel Programming for Massive Graph Abstract unavailable at time of publication. Analysis Gigliola Staffilani Abstract unavailable at time of publication. MIT [email protected] Kamesh Madduri Lawrence Berkeley National Lab [email protected] MS19 FETI Methods for Mortar Coupling of Stokes- Darcy Systems MS18 Deciding to Give Up Tenure: Surprising Decisions Abstract unavailable at time of publication. Along the Path Juan Galvis In a world where a successful research program and tenure Texas A&M University often seem to be the main goals of a Ph.D. scientist, one [email protected] may find it hard to see the other options available. After receiving my degree, my choices were not unusual and I MS19 spent three years as a postdoc before accepting a tenure- track position and, later, receiving tenure. But life some- An Efficient Method for Solving Multiscale Elliptic times throws unexpected opportunities in your path and Problems with Randomly Perturbed Data I will describe a number of my own experiences, includ- We propose a method for efficient solution of elliptic prob- ing how I eventually reached a decision to resign a tenured lems with multiscale features and randomly perturbed co- position. efficients. We use the multiscale finite element method Mary Ann Horn (MsFEM) as a starting point and derive an algorithm for National Science Foundaton solving a large number of problems with randomly per- [email protected] turbed coefficients. We show that the method converges to the multiscale finite element solution in the limit. A set of numerical examples is presented to illustrate theoretical AN10 Abstracts 47 findings. tested on both single-machine and multi-machine system with various random disturbances. Our numerical results Victor E. Ginting demonstrate the validity and performance of the proposed Department of Mathematics dynamic state estimation approach with different random University of Wyoming disturbance. The numerical results also indicate the pro- [email protected] posed approach can include the full dynamics of power sys- tems and ensure an accurate representation of the changing Axel Malqvist states in the power system. Department of Information Technology Uppsala University Guang Lin [email protected] Pacific Northwest National Laboratory [email protected] Michael Presho Department of Mathematics MS20 University of Wyoming [email protected] Capturing All Scales of Random Modes in Stochas- tic Pdes Based on Polynomial Chaos and a Dynam- ically Adaptive Wavelet Method MS19 We consider the numerical solution of stochastic partial dif- Some Computational Issues in the Modeling and ferential equations based on polynomial chaos expansion Simulation of Coupled Surface and Subsurface exemplified by the convection-diffusion equation model. Fluid Flow This typically produces a coupled system of determinis- Abstract unavailable at time of publication. tic equations for the coefficients of the expansion whose size can become very large with an increasing number of Junping Wang independent random variables and/or order of polynomial National Science Foundation chaos expansion (dubbed the dimensionality curse). We [email protected] tackle this challenge using a dynamically adaptive wavelet- based methodology utilizing space-refinement that allows capturing all scales of each random/unceretainty mode on MS20 indepenent grids. We present numerical results illustrating Adaptive Strategies for Uncertain Hyperbolic Sys- the salient features of the proposed method. tems Xiaoan Ren Abstract not available at time of publication. University of Shanghai for Science and Technology [email protected] Olivier Le Maitre LIMSI-CNRS Wenquan Wu [email protected] Univ of Shanghai Sci & Technology, Shanghai 200093, PR China Julie Tryoen xiao ann [email protected] Universit´eParis-Est CERMICS - ENPC Leonidas S. Xanthis [email protected] Department of Aeronautics, Imperial College [email protected] Alexandre Ern Universite Paris-Est CERMICS, Ecole des Ponts MS20 [email protected] Multi-element Uncertainty Quantification Based on Simplex Elements Stochastic Collocation MS20 A Simplex Elements Stochastic Collocation (SESC) A Multi-element Probabilistic Collocation Based method is presented for robust and efficient multi-element Kalman Filter for Power System Dynamic State uncertainty quantification. The method is equally robust Estimation as Monte Carlo simulation in terms of the Extremum Di- minishing (ED) robustness concept extended to probabil- In this talk, we propose an efficient multi-element proba- ity space. Its efficient extension to multiple uncertainties bilistic collocation method (MEPCM) based Kalman filter is considered based on higher degree polynomial interpola- (KF) to include dynamic state variables in the power sys- tion, randomized refinement sampling, and Essentially Ex- tem dynamic state estimation process. Comparing with tremum Diminishing (EED) extrapolation. These strate- Monte Carlo (MC) based KF approach, the proposed gies employ the flexibility of simplex elements for han- MEPCM based KF implementation can solve the system dling randomized samples and non-hypercube probability of stochastic state equations much more efficient. Ad- spaces. ditionally, the proposed approach can sample the gener- alized polynomial chaos approximation of the stochastic Jeroen Witteveen solution with an arbitrarily large number of samples, at Stanford University virtually no additional computational cost. The proposed Center for Turbulence Research algorithm thus can drastically reduce the sampling errors [email protected] and achieve a high accuracy at reduced computational cost, compared to the classical MC implementations of KF. The Gianluca Iaccarino proposed MEPCM-KF based dynamic state estimation is Stanford University 48 AN10 Abstracts [email protected] [email protected] MS21 MS21 Anatomy of a Young Giant Component in the Ran- Structure of Random r-SAT below the Pure Literal dom Graph Threshold In this talk, we will completely describe the structures of gi- It is well known that there is a sharp density threshold for ant components in random graphs G(n, p)withn−1/3 << a random r-SAT formula to be resolved by the pure literal pn − 1 << n−1/4. The description can be made by us- rule, and a higher threshold for satisfiability. We analyze ing random 3-regular graphs and Galton-Watson Poisson the structure of the (rare) unsatisfiable instances below the branching processes. The proof uses the Poisson cloning pure literal threshold, contrasting this with the situation model and interesting computation arguments. We will above the satisfiability threshold. It remains open what also present some results regarding the diameter and the happens between the two thresholds. mixing time of the giant component. Gregory Sorkin Jeong Han Kim IBM Research Yonsei University Mathematical Sciences National Institute for Mathematical Sciences [email protected] [email protected] MS21 Jian Ding Phase Transition for the Mixing Time of the University of California, Berkeley Glauber Dynamics for Coloring Regular Trees [email protected] We prove that the mixing time of the Glauber dynamics Eyal Lubetzky, Yuval Peres for random k-colorings of the complete tree with branch- Microsoft Research ing factor b undergoes a phase transition at k = b(1 + [email protected], [email protected] ob(1))/ ln b. Our main result shows nearly sharp bounds on the mixing time of the dynamics on the complete tree with n vertices for k = Cb/ln b colors with constant C.ForC ≥ MS21 1+o (1) 2 1 we prove the mixing time is O(n b ln n). On the How Frequently is a System of 2-linear Boolean other side, for C<1 the mixing time experiences a slowing Equations Solvable? down, in particular, we prove it is O(n1/C+ob(1) ln2 n)and 1/C−ob(1) We consider a random system of equations xi + xj = Ω(n ). The critical point C=1 is interesting since it b(i,j)(mod2), (xu ∈{0, 1},b(u,v) = b(v,u) ∈{0, 1}), with coincides (at least up to first order) to the so-called recon- the pairs (i, j)fromE, a symmetric subset of [n] × [n]. E struction threshold which was recently established by Sly. is chosen uniformly at random among all such subsets of a The reconstruction threshold has been of considerable in- given cardinality m; alternatively (i, j) ∈ E with a given terest recently since it appears to have close connections to probability p, independently of all other pairs. Also, given the efficiency of certain local algorithms, and this work was E, Pr{be =0} = Pr{be =1} for each e ∈ E, indepen- inspired by our attempt to understand these connections dently of all other be .Itiswellknownthat,asm passes in this particular setting. This is joint work with Prasad through n/2(p passes through 1/n, resp.), the underlying Tetali, Juan C. Vera, and Linji Yang which appeared at random graph G(n, #edges = m), (G(n, P r(edge) = p), SODA’10. resp.) undergoes a rapid transition, from essentially a for- est of many small trees to a graph with one large, mul- Eric Vigoda ticyclic, component in a sea of small tree components. Georgia Institute of Technology We should expect then that the solvability probability de- [email protected] creases precipitously in the vicinity of m ∼ n/2(p ∼ 1/n), 1/4 and indeed this probability is of order (1 − 2m/n) , MS22 1/4 for m Jae-Hun Jung Department of Mathematics Department of Mathematics University of Massachusetts at Dartmouth SUNY at Buffalo [email protected] jaehun@buffalo.edu MS23 MS22 Regulatory Network Analysis and Approximation Anti-Gibbs and Anti-Runge Strategies in Spectral Using Tensor Products Methods We consider the problem of recovering a system of func- We discuss the relationship between the Gibbs and Runge tions governing a regulatory network using time course Phenomena and compare and contrast single-interval and data, such as that obtained from microarrays. The num- three-interval strategies for combating them. RBF, Fourier ber of components in the regulatory networks and the cost and polynomial bases are analyzed. The limitations of of data collection yield an ugly showing of the curse of these tactics in both accuracy and ill-conditioning and the dimensionality. In [Gregory Beylkin, Jochen Garcke, and trade-offs are discussed. Martin J. Mohlenkamp. ’Multivariate regression and Ma- chine Learning with Sums of Separable Functions,” SIAM John P. Boyd Journal on Scientific Computing (2009)] it was shown that University of Michigan regression with tensor products provides a way to effec- Ann Arbor, MI 48109-2143 USAUSA tively overcome this and obtain good approximations. We [email protected] will incorporate time derivatives into a novel ALS algo- rithm for approximating time course data. We will then consider the benefits of using tensor product in the analy- MS22 sis of the approximated network. Kernel-based Approximation for Ill-posed and Dis- continuous Problems Ryan Botts Ohio University In this talk we present the idea of meshless computational Department of Mathematics methods based on the use of kernel-based functions for solv- [email protected] ing various ill-posed problems and problems with disconti- nuity. In particular, the recent development of the method Martin J. Mohlenkamp of fundamental solutions combined with various regular- Ohio University ization techniques to solve Cauchy problems and inverse [email protected] heat conduction/source identification problems will be dis- cussed. The main idea is to approximate the unknown solution by a linear combination of fundamental solutions MS23 whose singularities are located outside the solution domain. Solution of High-Dimensional Uncertain Systems More recent works on using the particular solutions and Using Separated Tensor Decomposition Green’s functions will be highlighted. Separated representations have proven efficient to relax the Benny Hon curse-of-dimensionality associated with approximation of City University of Hong Kong high-dimensional functions and tensors. In this talk, we [email protected] will demonstrate how such techniques can be employed to efficiently compute the solution to linear elliptic PDEs with MS22 high-dimensional random inputs. In particular, we illus- trate the existence of an approximate low-rank separated Detecting Edges in Two Dimensional Functions solution to this class of stochastic PDEs and further explore with Radial Basis Functions the efficiency and accuracy of alternating least-squares and We extend the iterative adaptive multiquadric radial Rayleigh-Ritz methods in computing the solution. basis function (RBF) method for the detection of lo- Alireza Doostan cal jump discontinuities in one-dimensional problems to Department of Aerospace Engineering Sciences two-dimensional problems. The iterative edge detection University of Colorado, Boulder method is based on the observations that the absolute val- [email protected] ues of the expansion coefficients of RBF approximation have different growth rates. We consider two approaches: the dimension-by-dimension technique and the global ex- Gianluca Iaccarinno tension approach. Numerical examples demonstrate that Stanford University the two-dimensional iterative adaptive RBF method yields Mechanical Engineering accurate results. [email protected] Saeja Kim University of Massachusetts Dartmouth MS23 [email protected] Source Apportionment of Particulate Matter Based on Segregated Composition Size Data Jae-Hun Jung Source apportionment on bulk airborne particle composi- Department of Mathematics tion data resolves the data into a product of two matrices, SUNY at Buffalo source profiles and source contributions. However, size- jaehun@buffalo.edu segregated data requires a higher dimension model since source emissions are size dependent. Source profiles are Sigal Gottlieb a composition matrix as a function of size. The model 50 AN10 Abstracts sums outer products of matrices times vectors. The ap- have implemented this scheme to model hot electron-hole plication of this model to data collected using a multiple transport on Kane and full energy bands in the simulation stage DRUM impactor analyzed by synchrotron XRF will of one and two-physical and three-phase space dimensional be presented. nano meter devices with a base channels of the order of 50 nm. The energy bands (conduction and valence) are Philip Hopke, Jihong Wang, Punith Nallathamby computed by eingenvalue solvers for electronic structures Clarkson University and coupled to the Boltzmann Poisson solver. Center for Air Resources Engineering and Science [email protected], [email protected], James R. Chelikowsky [email protected] University of Texas at Austin Institute for Computational Engineering and Sciences [email protected] MS23 Convergence Results of a Regularized Tensor Al- Yingda Cheng ternating Least-Squares Dept. of Math and ICES University of Texas at Austin The Alternating Least-Squares (ALS) for tensor decom- [email protected] position has been widely applied to many applications, for example, in signal processing, chemometrics, and data min- ing. However, ALS can be prohibitively slow for some data. Irene M. Gamba We improved ALS with a regularization term. Our regu- Department of Mathematics and ICES larized method dramatically accelerates ALS. In addition University of Texas to providing numerical results for several data sets, we also [email protected] analyze the convergence of the regularized ALS. Armando Majorana Na Li Dipartimento di Matematica e Informatica Clarkson University University of Catania - Italy Dept. of Mathematics [email protected] [email protected] Chi-Wang Shu Stefan Kindermann Division of Applied Mathematics RICAM Institute, Linz Brown University Austria [email protected] [email protected] Carmeliza Navasca MS24 Clarkson University New Models of Chemotaxis: Analysis and Numer- [email protected] ics Patlak-Keller-Segel (PKS) system is a classical PDE model MS24 of the chemotaxis. The system admits solutions that de- Diffeomorphic Deformation of Textured Shapes velop delta-type singularities within a finite time. Even though such blowing up solutions model a concentration This work proposes a new method to match shapes with in- phenomenon, they are not realistic since biological cells terior textures through the large deformation diffeomorphic do not converge to one point (while the cell density grows metric mapping, resulting in optimizing for geodesics on sharply, it must remain bounded at all times). I will present the space of diffeomorphisms connecting textured shapes. a new chemotaxis model, which can be viewed as a regular- This study is to address the problems caused by boundary ized PKS system. The proposed regularization is based on conditions and response to the need of matching regions of a basic physical principle: boundedness of the chemotac- interest directly in medical applications. The new method tic convective flux, which should depend on the gradient is illustrated by numerical experiments on MRI images. of the chemoattractant concentration in a nonlinear way. Solutions of the new system may develop spiky structures Yan Cao that model the concentration phenomenon. However, both The University of Texas at Dallas cell density and chemotattractant concentration remain Department of Mathematical Sciences bounded. The proposed model is studied both analytically [email protected] and numerically. I will first prove a global existence result and then use the bifurcation theory to investigate existence of nontrivial steady states and their stability properties. MS24 Finally, I will persent one- and two-dimensional numerical Discontinuous Galerkin Solvers for Boltzmann- examples that support the analytical findings and demon- Poisson Systems in Semiconductor Device Simu- strate the formation and stability of the spiky solutions. lations The Boltzmann transport equation (BTE) is an integro- Alina Chertock differential model describing the evolution of a single point North Carolina State University probability distribution function. The mathematical and Department of Mathematics computational difficulties associated to this equation is due [email protected] to the balance of a first order transport term to an integral operator accounting for their interactions and the coupling Alexander Kurganov to the Poisson equation. We will present a discontinuous Tulane University Galerkin solver for the Boltzmann Poisson systems. We Department of Mathematics AN10 Abstracts 51 [email protected] Enhanced heavy oil recovery methods are an intensive re- search area in the oil industry, and have recently generated Xuefeng Wang a battery of recovery methodsinwhatisthelargestgrow- Tulane University ing sector of this industry, such as steam assisted grav- [email protected] ity drainage (SAGD) and cyclic steam stimulation (CSS). However, the environmental impacts of these processes and Yaping Wu the use of a high volume of water and natural gas suggest Department of Mathematics that extensive research is required for economic and en- Capital Normal University, Beijing vironmentally friendly development of heavy oil reserves. [email protected] This presentation will give an overview on current research in heavy oil modeling, and the presenter will also describe his current research program. MS24 Zhangxin Chen Well-Balanced University of Calgary Positivity Preserving Central-Upwind Scheme on [email protected] Triangular Grids for the Saint-Venant System We introduce a new second-order central-upwind scheme MS25 for the Saint-Venant system of shallow water equations on triangular grids. We prove that the scheme both preserves Rate and Interfacial Tension Dependent CO2 Rel- stationary steady states (lake at rest) and guarantees the ative Permeability positivity of the computed fluid depth. Moreover, it can We must develop models that better quantify and predict be applied to models with discontinuous bottom topogra- the influence of important processes and small-scale fea- phy and irregular channel widths. We demonstrate these tures on fluid behavior and CO2 containment potential. features of the new scheme, as well as its high resolution Recent relative permeability experiments indicate the de- and robustness in a number of numerical examples. pendency of supercritical CO2 on interfacial tension as re- Yekaterina Epshteyn sults of changes in insitu pressure and temperature. We Carnegie Mellon University investigate the roles of capillary and gravitational forces on [email protected] trapped CO2 taking into the account the density, viscosity, and interfacial tension computed from the compositional phase behavior model. The combined capillary number Alexander Kurganov and bond number referred to as trapping number is suc- Tulane University cessfully implemented in IPARS simulator. Carbon storage Department of Mathematics simulation results with this new modeling capability will be [email protected] presented. Steve Bryson Mojdeh Delshad NASA Ames Research Center Department of Petroleum and Geosystems Engineering [email protected] The University of Texas at Austin [email protected] Guergana Petrova Texas A&M University X. Kong, M.F. Wheeler [email protected] The University of Texas at Austin [email protected], [email protected] MS24 Stencil Choosing Sensitivity in WENO MS25 Towards Accurate and Efficient Computational The weighted essentially non-oscillatory method is used to Modeling of Compositional Multi-Phase Flow approximate the numerical flux in hyperbolic PDEs. It automatically centers the stencil in the case of a smooth Compositional multiphase flow modeling is demanded in solution, but places larger weights on the smoother stencils a number of important reservoir applications. One im- when the solution has a discontinuity pr sharp gradient. A portant application is the injection of CO2 in hydrocarbon significant parameter in determining the WENO weights, reservoirs, which is becoming increasingly attractive due to epsilon, is examined in this talk. issues related to global warming, and urgent need for im- proved hydrocarbon recovery. CO2 injection may be one of Sigal Gottlieb the most attractive options for water-flooded reservoirs ei- Department of Mathematics ther from natural water drive or from external water injec- University of Massachusetts at Dartmouth tion. Coinjection of water with CO2 for improved mobility [email protected] and increased sweep is also another process of high interest. To evaluate the efficiency of CO2 injection in water-flooded reservoirs, three-phase compositional numerical models are MS25 required. Such simulators have to be equipped with unique Mathematical and Computational Techniques for capabilities due to special features of CO2 mixing in hydro- Modeling Heavy Oil Reservoirs carbon liquids and gases for CO2 applications. Specifically, when dissolved in liquids, CO2 generally increases the den- While mathematical modeling has been successful in the sity. The density increase may have a significant effect on recovery of conventional oil, it is still in the early stage flow path and may result in convective mixing. Moreover, of heavy oil modeling. As conventional oil reserves dwin- CO2 mixing may also lower the viscosity; it may also in- dle and oil prices rise, heavy oil is now the center stage. crease the viscosity when there is substantial vaporization 52 AN10 Abstracts of liquid components in the gas phase. In this talk, we MS26 present three-phase compositional flow using higher-order On the Geometry of Discrete Exponential Families finite element methods. Our numerical model is based on with Application to Exponential Random Graph an iterative IMPEC coupling of the pressure equation and Models species transport equations, which are solved by mixed finite element (MFE) and discontinuous Galerkin (DG) There has been an explosion of interest in statistical models methods, respectively. Compared with lower order finite for analyzing network data, and considerable interest in the volume methods, the proposed MFE-DG method has low class of exponential random graph (ERG) models. In this numerical diffusion and can capture the solution discon- talk I will relate the properties of ERG models to the prop- tinuities and yield accurate prediction of shock locations erties of the broader class of discrete exponential families. arising in computational three-phase flow. I will describe a general geometric result about discrete ex- ponential families with polyhedral support. Specifically, I Shuyu Sun will show how the statistical properties of these families can Clemson University and KAUST be well captured by the normal fan of the convex support. [email protected] I will discuss the relevance of such results to maximum like- lihood estimation and apply them to the analysis of ERG models. By means of a detailed example, I will provide MS25 some characterization and a partial explanation of certain Compositional Flow Modeling in Porous Media pathological features of ERG models known as degeneracy. Joint work with S.E. Fienberg and Y. Zhou. Abstract unavailable at time of publication. Alessandro Rinaldo Mary F. Wheeler Carnegie Mellon University Center for Subsurface Modeling [email protected] University of Texas at Austin [email protected] MS26 MS26 Introduction to Algebraic Statistics Methods and Extensions for Parametric Inference I will provide an introduction to one of the main principles of algebraic statistics– that statistical models are semial- In computational biology, many dynamic programming al- gebraic sets– and illustrate this idea with examples. gorithms can be interpreted as sum-product algorithms. Examples include Needleman–Wunsch and Viterbi algo- Seth Sullivant rithms. Equivalently, these algorithms find the leading North Carolina State University term of a given multivariate polynomial F presented by a [email protected] short and structured factorization. The term-order is part of the input to the sum-product algorithm – and usually depends on biological parameters which are estimated or MS27 guessed. Given this unavoidable uncertainty in the term- Improved Accuracy and Efficiency in Leray-alpha order, it is desirable to compute all leading terms of F Computations under all term-orders – or equivalently, compute the New- ton polytope NP(F ). This is called parametric inference. In this talk I will discuss how some simple, yet fundamen- In this talk, we discuss a general framework for parametric tally important changes can be made to numerical algo- inference, and review existing results. Then, motivated by rithms for the Leray-alpha model, which lead to significant biological applications, we propose two key extensions: 1) improvements in accuracy and efficiency. parametric k-best inference, and 2) constrained parametric inference. We finish by presenting recent results for these Abigail Bowers extensions. Department of Mathematical Sciences Clemson University Peter Huggins [email protected] Carnegie Mellon University [email protected] Leo Rebholz Clemson University Department of Mathematical Sciences MS26 [email protected] Hermite Polynomial Aliasing A representation of Hermite polynomials of degree 2n +1, MS27 as sum of an element in the polynomial ideal generated by Uncertainty Quantification for a Two Domain Nat- the roots of the Hermite polynomial of degree n and of a re- ural Convection Problem minder, suggests a folding of multivariate polynomials over a finite set of points. From this, the expectation of some Numerical algorithms are studied for a Boussinesq model of polynomial combinations of random variables normally dis- natural heat convection in two domains, motivated by the tributed is computed. This is related to quadrature formu- dynamic core of climate models. One (monolithic) algo- las and has strong links with designs of experiments. This rithm is coupled across the fluid-fluid interface. Another is is joint work with G. Pistone. decoupled for parallel implementation using a partitioned time stepping approach which retains unconditional sta- Eva Riccomagno bility, a property not enjoyed by climate codes. Reliability Politecnico di Torino of computations is critical in climate applications. Here, [email protected] stochastic noise is introduced into two nonlinear coupling AN10 Abstracts 53 terms that play an important role in stability, by modeling x-ray beam entering the object, resulting in reconstruction them with random variables. Uncertainty quantification inaccuracies. In this talk, we discuss a novel mathemat- is performed for each numerical model via the stochastic ical model based on a polyenergetic x-ray spectrum and collocation method. Comparison is drawn in predicting develop statistically based iterative optimization methods average surface temperature with a “one-way” coupled al- for image reconstruction. gorithm, based on a common climate modeling assumption. Uncertainty in average surface temperature is smallest us- Julianne Chung ing the monolithic algorithm and largest using the one-way University of Maryland coupled algorithm. [email protected] Jeff Conners Department of Mathematics MS28 University of Pittsburgh Modeling Combustion Reactions with Step- [email protected] function Kinetics Here we develop reaction-diffusion models for self- MS27 propagating reactions where the reaction rates are reduced Unconditional Convergence of Extrap- to have a step-function dependence on temperature. This olated Crank-Nicolson, Finite Element Method for approach is first studied in the case of a single reaction in the Navier-Stokes Equations both adiabatic and non-adiabatic environments. Of partic- ular interest, however, are systems of two reactions which Error estimates for the Crank-Nicolson in time, finite ele- are thermodynamically coupled. This includes parallel, ment in space (CN-FE) discretization of the Navier-Stokes competing and sequential reactions with both analytical equations require a discrete version of the Gronwall in- and numerical studies. equality, which leads to a time-step restriction. Previous convergence analyses of CN-FE with linear extrapolation Erin Lennon rely on a similar time-step restriction as the full CN-FE. Northwestern University We show that CN-FE with linear extrapolation is uncon- [email protected] ditionally convergent in the energy norm. We also show optimal convergence of CN-FE with extrapolation in a dis- crete L∞(H1)-norm and convergence of the corresponding MS28 discrete time derivative in a discrete L2(L2)-norm. Non-linear Wave Interactions in Rotating Strati- fied Fluid Flow Ross Ingram Univeristy of Pittsburgh By utilizing the asymptotic renormalization theory of Department of Mathematics Wirosoetisno et. al. (2002), we propose a dynamical ex- [email protected] planation of the origin, nature, and energetics of the spon- taneously generated inertia gravity waves which appear in the trough of the baroclinic mode in the rotating annular MS27 experiment of Williams et. al. (1995). We compute the Continuum-microscopic Interaction Computation O() Wirosoetisno et. al. correction term to the under- of Viscoelastic Turbulent Flow lying quasi-geostrophic dynamics and compare our results with the location, amplitude, and morphology of the IGWs A continuum-microscopic approach is presented for com- observed in the lab. We also present energy results which puting dilute polymer flow at high Reynolds number. o?er a strong indication that in a real ?uid, spontaneously The finitely extensible nonlinear elastic (FENE) model is emitted IGWs are an extremely efficient way to transfer adopted to describe the polymer. The continuum conser- geostrophic energy to ageostrophic energy, ultimately to vation equations require a closure based upon the micro- be dissipated by viscosity scopic configuration state of the polymer. A fast compu- tational method to obtain the closure is presented based Dawn Ring upon a multiphysics PARAREAL approach that simulta- Wentworth Institute of Technology neously solves stochastic differential equations for the poly- [email protected] mer state and a Fokker-Planck equation for the probability distribution function of the states. The PARAREAL algo- rithm is implemented on graphical processing units. The MS28 microscopic configuration of the polymer is investigated at The Mechanical Stability of Growing Arteries maximum drag reduction conditions. In many cylindrical structures in biology, residual stress Sorin Mitran fields are created through differential growth. The possible University of North Carolina Chapel Hill role of axial residual stress in regulating stress in arteries [email protected] and preventing buckling instabilities is investigated. It is shown that axial residual stress lowers the critical internal pressure leading to buckling and that a reduction of ax- MS28 ial loading may lead to a buckling instability which may Numerical Methods for a Problem arising in 3D eventually lead to arterial tortuosity. Breast Image Reconstruction Rebecca Vandiver Digital tomosynthesis imaging, the process of reconstruct- Bryn Mawr College ing a 3D object from a few 2D projection images, is a vi- [email protected] able alternative to standard mammography in breast can- cer imaging. However, current algorithms for image recon- struction do not incorporate the polyenergetic nature of the 54 AN10 Abstracts MS29 bility. Our numerical results suggest the use of high-order An Implicit Asymptotic Preserving Maxwell Solver schemes in both time and space is advantageous when con- sidering the Vlasov-Maxwell system. We present an extension of the Boundary Integral Treecode (BIT), a grid free electrostatic O(N log N) field solver, to Guang Lin an implicit electromagnetic field solver. The new method Pacific Northwest National Laboratory computes a magnetic field that is by construction diver- [email protected] gence free in the computational domain. The implicit Maxwell solver is Asymptotic Preserving, recovering the Jingmei Qiu Darwin limit of Maxwell’s equations in the long time limit. Colorado School of Mines The method is used to simulate the wave equation for a Department of Mathematical and Computer Sciences range of geometries. [email protected] Andrew J. Christlieb Fengyan Li Michigan State Univerity Rensselaer Polytechnic Institute Department of Mathematics [email protected] [email protected] Yingda Cheng Benjamin Ong Dept. of Math and ICES Department of Mathematics University of Texas at Austin Michigan State University [email protected] [email protected] MS29 MS29 A Hybrid Method for Particle Transport Based on Quadrature-based Moment Methods for Gas- a Collided/ Uncollided Split particle Flows We present methods for solving particle transport prob- Gas-particle flows can be modeled by a kinetic equation lems by splitting the phase space density into particles that (KE) for the particle velocity distribution function. The have collided during a time step and particles that have not KE contains terms for acceleration due to fluid drag and collided in the step. Using this paradigm we can exploit gravity, transport due to the particle velocity, and inelastic benefits of particular methods: moment based methods be- particle-particle collisions. The KE is coupled to the conti- have well for particles that have undergone many collisions nuity and momentum equations for the gas phase through whereas other methods such as Monte Carlo and integral the particle volume fraction and fluid drag, respectively. methods are inefficient in the absence of collisions. We an- In most practical gas-particle flows, the particle Knudsen alyze our splitting technique and apply it to several prob- and Mach numbers are far from the equilibrium limit so lems of interest in linear transport and radiative transfer. that hydrodynamic descriptions of the KE are invalid. In this talk we describe how quadrature-based moment meth- Ryan G. McClarren ods can be used to capture non-equilibrium solutions to the Texas A & M KE. Examples for gas-particle riser flows ranging from very [email protected] dilute to moderately dense particle hold up will be used to illustrate the quadrature-based numerical methods. Cory Hauck Alberto Passalacqua, Cansheng Yuan Oak Ridge National Laboratory Iowa State University [email protected] [email protected], [email protected] MS30 Rodney O. Fox A Least Squares Interpolation Method for Param- Iowa State University eterized Systems of Equations Department of Chemical and Biological Engr. [email protected] Exhaustive parameter studies for engineering models are typically infeasible, particularly for high dimensional pa- rameter spaces; cheap and accurate surrogate models are MS29 necessary. We propose a reduced basis method for approx- A High-order WENO Method for the Vlasov- imating the model output given a set of snapshots. The Maxwell System method requires only the residual of the model and poses a least squares problem to compute the coefficients of the In this talk, we propose a high-order Vlasov-Maxwell solver approximation. We derive the method in the linear case based on high-order Runge-Kutta scheme in time and high- and extend its application to nonlinear problems order WENO (weighted essentially nonoscillatory) recon- struction in space. The spatial WENO reconstruction de- Paul Constantine veloped for this method is conservative and Lax- Friedrichs Sandia National Labs flux is employed. While the third, fifth, seventh and ninth [email protected] order reconstructions are presented in this talk, the scheme can be extended to arbitrarily high order. WENO recon- Qiqi Wang struction is able to achieve high-order accuracy in smooth Massachusetts Institute of Technology parts of the solution while being able to capture sharp in- [email protected] terfaces without introducing oscillations. The quality of the proposed method is demonstrated by applying the ap- proach to classical plasma problems, such as Weibel’s insta- AN10 Abstracts 55 MS30 MS31 A Geometric Interpretation of Stochastic Processes Self-correcting Estimates from the Differential Equations Method The time-dependent stochastic harmonic oscillator can ge- ometrically be described as the Lie-derivative of a 1-form In this talk we discuss applications of the differential equa- on a cylinder, [a, b] × S1 tions method for random graph processes in which the (1) bounds on variation around the expected trajectory de- LX α =0, crease as the process evolves. These methods are illus- (1) trated in the context of the random greedy algorithm for where α is the 1-form which describes the probability constructing a large partial Steiner-Triple-System. This distribution and X is the vector field on the cylinder. This (1) stochastic graph process is defined as follows: We begin equations means that α does not change along the flow with a complete graph on n vertices and proceed to remove generated by X. In this talk a discrete polynomial repre- the edges of triangles one at a time, where each triangle re- sentation of differential forms will be given which possesses moved is chosen uniformly at random from the collection the invariance under general transformations. As a conse- of all remaining triangles. The process terminates when it quence, integrals are preserved. arrives at a triangle-free graph. We show that with high probability the number of edges in the final graph is at Mark Gerritsma 7/4 5/4 Delft University of Technology most O(n log n). [email protected] Tom Bohman Carnegie Mellon University MS30 Department of Mathematics [email protected] Proper Generalized Decomposition and Separated Representations for Uncertainty Propagation Alan Frieze A family of model reduction methods, called Proper Gen- Carnegie Mellon University eralized Decompositions methods, has been recently pro- [email protected] posed for the propagation of parametric uncertainties in very high dimensional stochastic models. It is based on the Eyal Lubetzky a priori construction of separated representations of the so- Microsoft Research lution of stochastic partial differential equations. They can [email protected] be seen as generalizations of Karhunen-Lo`eve decomposi- tions. Here, we review basic definitions of decompositions and propose new definitions in order to improve conver- MS31 gence properties of separated representations. A Probabilistic Technique for Finding Almost Pe- riods in Additive Combinatorics Anthony Nouy GeM, Ecole Centrale Nantes We (E. Croot and Olof Sisask) introduce a new probabilis- University of Nantes, CNRS tic technique for finding ‘almost periods’ of convolutions [email protected] of subsets of finite groups. This allows us to give: a new probabilistic proof of Roth’s theorem; a new way to ap- Mathilde Chevreuil proach the 2D corners problem; a new result on the exis- GeM - CNRS tence of long arithmetic progressions in sumsets A+B; a Ecole Centrale Nantes, University of Nantes translation-invariance result for “discontinuous sets’ (sets [email protected] A whose convolution function A*A is somewhat discontin- uous); and several non-abelian analogues of classical theo- rems. In many cases, the proofs we give are the shortest MS30 known, and in some cases, like the case of long APs in Practical Rare Event Simulation in High Dimen- A+B, we obtain the strongest bounds to date. sions Ernie Croot I will discuss an importance sampling method for certain Georgia Institute of Technology rare event problems involving small noise diffusions. Stan- [email protected] dard Monte Carlo schemes for these problems behave expo- nentially poorly in the small noise limit. Previous work in rare event simulation has focused on developing, in specific MS31 situations, estimators with optimal exponential variance Maximal Independent Sets in Graphs with No decay rates. I will introduce an estimator related to a de- Quadrilaterals terministic control problem that not only has an optimal variance decay rate under certain conditions, but that can There has recently been substantial attention on studying even have vanishingly small statistical relative error in the the problem of studying the size of a maximal independent small noise limit. The method can be seen as the limit of set in a d-regular graph of girth g. The most precise results a well known zero variance importance sampling scheme for large girth and fixed d are recent, due to Gamarnik and for diffusions which requires the solution of a second order Goldberg, following from an analysis of a natural random- partial differential equation. I will also report on progress ized greedy algorithm, and generalizing earlier results of toward applying the algorithm within the design of mag- Lauer and Wormald. In this talk I will present a new re- netic memory devices. sult stating that every d-regular n-vertex graph containing no cycles of length four contains a maximal independent Jonathan Weare set I such that NYU |I| log d + O(1) = [email protected] n d 56 AN10 Abstracts and this is best possible in light of random regular graphs nonlinear biharmonic equation describing the stream func- containing no cycles of length four, and in view of unions of tion. For the Stokes problem (Re = 0), the accuracy dete- complete bipartite graphs with d vertices in each part. In riorates in the neighborhood of the singularities in the two particular, this shows that it is the lack of cycles of length upper corners. Enlarging the space spanned by the RBF four in a random d-regular graph which gives maximal in- basis functions with additional functions that capture the dependent sets of this size. The degree ratio of a graph singular behavior of the solution restores the spectral con- G is defined to be the maximum degree divided by the vergence of the RBF method. minimum degree. We show that there are n-vertex graphs with degree ratio c>1 and maximum degree d such that Manuel Kindelan every maximal independent set has size at least roughly Universidad Carlos III de Madrid n/2d2/(c−1) , which shows that the condition of regularity [email protected] is necessary in the above-mentioned result. We conjecture that such graphs should have maximal independent sets Francisco Bernal of size at most n/2d1/c up to logarithmic factors in d for UT Dresden each c>1. Some possible applications of the result will be [email protected] mentioned. Jacques Verstraete MS32 University of California-San Diego Approximation of Linear PDEs by Gaussian Pro- [email protected] cesses via Mat´ern Functions N M Given discrete data X1 = {xi}i=1 ⊂ Ω, X2 = {xN+j }j=1 ⊂ Paul Horn N,M,n,m Emory University ∂ΩandY = {fk(xi),gl(xN+j )}i,j,k,l=1 generated from a [email protected] linear partial differential equations model ··· ∈ MS31 L1u(x)=f1(x), ,Lnu(x)=fn(x),xΩ, B1u(x)=g1(x), ···,Bmu(x)=gm(x),x∈ ∂Ω, Optimal Inverse Littlewood-Offord Theorems Let ηi,i=1,...,n be iid Bernoulli random variables, tak- where Lk are linear differential operators, Bl are linear ± 1 ing values 1 with probability 2 .GivenamultisetV of n boundary operators and Ω is a bounded domain in R . integers v1,...,vn, we define the concentration probability A Gaussian process (Gaussian field) is firstly constructed as on a reproducing-kernel Hilbert space corresponding to ρ(V ):=supPr(v1η1 + ...vnηn = x). aMat´ern function. Combining this with a Bayesian ap- x proach, we are able to obtain confidence interval for u(x) A classical result of Littlewood-Offord and Erd˝os from (instead of an explicit approximating solution) through a the 1940s asserts that if the vi are non-zero, then ρ(V ) covariance matrix for the above differential operators and is O(n−1/2). Since then, many researchers obtained im- boundary operators on the data. proved bounds by assuming various extra restrictions on V . About 5 years ago, motivated by problems concern- Qi Ye, Greg Fasshauer ing random matrices, Tao and Vu introduced the Inverse Illinois Institute of Technology Littlewood-Offord problem. In the inverse problem, one [email protected], [email protected] would like to give a characterization of the set V ,given that ρ(V ) is relatively large. In this talk, I describe a new method to attack the inverse problem. As an application, MS32 we strengthen a previous resultof Tao and Vu, obtaining an RBF Interpolation using Gaussians with Domain optimal characterization for V . This characterization im- Decomposition on GPUs mediately implies several classical theorems, such as those of S´ark¨ozy-Szemer´edi and Hal´asz. As another application, We have developed a parallel algorithm for radial basis we obtain an asymptotic, stable version of a famous theo- function (RBF) interpolation with Gaussians that uses a rem of Stanley that shows that under the assumption that GMRES iterative solver with a restricted additive Schwarz method (RASM) as a preconditioner and a fast matrix- the vi are different, ρ(V ) attains its maximum value when V is a symmetric arithmetic progression. All results ex- vector algorithm. The fast decay of the Gaussian basis tend to the general case when V is a subset of an abelian function allows rapid convergence of the iterative solver, while retaining the interpolation accuracy. The precon- torsion-free group and ηi are independent variables satis- fying some weak conditions. The method also works well ditioning and matrix-vector product have been ported to in the continuous setting, and gives a simple proof for the CUDA to accelerate the most computationally intensive β-net theorem of Tao and Vu, parts of the algorithm. Van Vu Rio Yokota Rutgers University University of Bristol [email protected] [email protected] Lorena Barba MS32 Department of Mechanical Engineering RBF Solution of Lid-Driven Cavity Problem using Boston University the Stream Function Formulation. [email protected] An RBF meshless formulation of the lid-driven cavity prob- lem is presented. It uses the streamfunction-vorticity for- MS33 mulation so that the problem reduces to the solution of a Nonlinear Model Reduction Using Petrov-Galerkin AN10 Abstracts 57 Projection and Data Reconstruction MS33 Nonlinear Model Reduction for Porous Media Flow A computational method for constructing a stable non- linear reduced-order model and analyzing it as fast as A Discrete Empirical Interpolation Method (DEIM) is ap- possible is presented. This method operates on a semi- plied in conjunction with Proper Orthogonal Decomposi- discretized partial differential equation, relies on a time- tion (POD) to construct a nonlinear reduced-order model variant Petrov-Galerkin projection for performing model of finite difference discretized system used in the simula- reduction, and on a reconstruction of gappy data for ap- tion of nonlinear miscible viscous fingering in a 2-D porous proximating the reduced system in order to accelerate its medium. POD is first applied to extract a low-dimensional analysis. Its performance is highlighted for sample prob- basis that optimally captures the dominant characteris- lems in nonlinear structural dynamics and computational tics of the system trajectory. This basis is then used in fluid dynamics a Galerkin projection scheme to construct a reduced-order system. DEIM is then applied to greatly improve the effi- Kevin T. Carlberg, Charbel Farhat ciency in computing the projected nonlinear terms in the Stanford University POD reduced system. DEIM achieves a complexity reduc- [email protected], [email protected] tion of the nonlinearities which is proportional to the num- ber of reduced variables while POD retains a complexity MS33 proportional to the original number of variables. Numer- ical results demonstrate that the dynamics of the viscous Interpolation for Adaptation of Parameterized fingering in the full-order system of dimension 15000 can Reduced-Order Models be captured accurately by the POD-DEIM reduced system Differential geometry-based interpolation methods are pre- of dimension 40 with the computational time reduced by sented for adapting parameterized reduced-order models factor of 1000. to new configurations. Two different approaches are dis- Danny C. Sorensen, Saifon Chaturantabut cussed. The first one is based on the interpolation of the Rice University underlying reduced-order bases, and the second one on the [email protected], [email protected] interpolation of the reduced-order models themselves. It is shown that the first approach privileges robustness and stability at the expense of CPU time, whereas the second MS34 approach delivers real-time performance at the expense of Correctors and Field Fluctations for the p(x)- generality and robustness. Laplacian With Rough Exponents Charbel Farhat, David Amsallem We provide a corrector theory for the strong approximation Stanford University of fields inside composites made from two materials with [email protected], amsallem @stanford.edu different power law behavior. The correctors are used to develop bounds on the local singularity strength for gradi- MS33 ent fields inside micro-structured media. The bounds are multi-scale in nature and can be used to measure the am- A Projection-based Moment-matching Interpola- plification of applied macroscopic fields by the microstruc- tion for Large-scale Frequency Response Problems ture. Frequency response problems appear in many computa- Silvia Jim´enez Bola˜nos tional engineering applications. Our attention is focused − 2 −1 Louisiana State University on evaluating functions H(σ)=(K + iσD σ M) f on [email protected] an interval [σl,σr]. The matrices K, D,andM are very large and they may have complex entries. We will present a method that exploits values of H and its derivatives at mul- Robert P. Lipton tiple frequencies to construct a reduced-order model for H, Department of Mathematics with a rational Krylov algorithm. The computations of val- Louisiana State University ues and derivatives of H at multiple frequencies will use a [email protected] domain-decomposition based iterative solver for large-scale indefinite systems. Numerical experiments from structural MS34 dynamics and fluid-structure interaction problems will il- lustrate the method. Simultaneous Seismic Imaging and Inversion Using an Inverse Scattering Algorithm for One Dimen- Charbel Farhat sional Media Stanford University [email protected] We present the theoretical development and numerical tests for an acoustic inverse scattering algorithm for si- multaneous geophysical imaging and amplitude correction Ulrich Hetmaniuk directly from measured data. No knowledge about the University of Washington medium under investigation is assumed. We model several Department of Applied Mathematics one-dimensional earth configurations and show how the al- [email protected] gorithm can find the precise location and a good estimate of the layers’ parameters. Our tests will include different Radek Tezaur number of layers, high/low contrasts, velocity inversions Stanford University and noisy data. [email protected] Ashley Ciesla, Bogdan G. Nita Montclair State University [email protected], [email protected] 58 AN10 Abstracts MS34 MS34 A Split Sparse/Dense Matrix-vector Multiplication A Novel Pseudospectral Method for Solving a Non- Routine linear Volterra Partial Integro-differential Equa- tion on a Polar Geometry A new approach for performing sparse matvecs consists of two stages: (i) reordering the original general sparse In this talk we study a nonlinear Volterra partial integrod- matrix A, and (ii) tearing it into two parts – a dense banded ifferential equation used to model swelling porous mate- matrix B and a sparse matrix E of much lower rank. The rials where the problem domain is a unit disk. We show multiplication of B by a block of vectors is accomplished well-posedness is established under a given set of assump- via dense BLAS level-3. Comparisons are made against tions and introduce a novel approach to constructing pseu- Intel’s MKL matvec routine, mkl csrmm, and an extension dospectral differentiation matrices in a polar geometry for using multi-level splitting will be discussed. computing the spatial derivatives. An exponential time- differencing scheme proposed by Cox and Matthews and Eric Cox, Murat Manguoglu, Faisal Saied stabilized by Kassam and Trefethen is employed for the Purdue University time-stepping. [email protected], [email protected], [email protected] Keith Wojciechowski Department of Mathematical Sciences Ahmed Sameh University of Colorado Denver Department of Computer Science [email protected] Purdue University [email protected] Lynn S. Bennethum University of Colorado at Denver [email protected] MS34 Understanding the Impact of Low Permeability Jinhai Chen Media on Groundwater Remediation University of Colorado Denver [email protected] In groundwater remediation, contaminants trapped in low permeability soil after clean-up can cause re-contamination in the same area after a long period of time. This study fo- Kristian Sandberg cuses on analyzing and improving existing transport mod- GeoEnergy, Corp. els to understand the contaminant behavior in terms of [email protected] the rate and extent of release from its storage at the HPM/LPM interface. Current approaches include a sensi- MS35 tivity analysis of MT3D model, constructing a numerical simulator using dual-resistance model, and lab cell experi- Tensor Decomposition and Application in Signal ments. and Image Processing Xiaojing (Ruby) Fu We describe a module of n-by-n matrices upon which every linear transformation can be represented by tensor-matrix Clarkson University multiplication by a third-order tensor. An inner product is [email protected] defined upon this space, which leads to natural extensions of many familiar linear algebra concepts including length Kathleen Fowler and orthogonality of matrices and symmetry of tensors. Clarkson University Using these tools, we investigate applications in the areas Department of Mathematics of image reconstruction and pose estimation. [email protected] Karen S. Braman, Randy Hoover Michelle Crimi, Guannan Jiang South Dakota School of Mines Clarkson University [email protected], [email protected] [email protected], [email protected] Misha E. Kilmer Tufts University MS34 [email protected] Option Pricing Using Monte Carlo Methods One main problem in mathematical finance is calculating MS35 the fair price of derivatives. The fair price of these deriva- Tensor-SVD with Applications in Image Processing tives can be complicated to obtain analytically. Monte Carlo simulation offers a way to compute the fair price, We define a new type of multiplication between two ten- relying on the approximation of an expected value by the sors that allows for the representation of tensors as prod- average of the simulated values. We created a random num- ucts of tensors that is reminiscent of matrix factorizations. ber generator and used it in the Monte Carlo simulations In particular, a tensor-SVD is introduced which gives a to compute the fair price of a European Call Option. representation of tensors as the sum of outer-products of matrices. This alternate representation of tensors shows Steven Melanson, Yered Pita-Juarez promise with respect to the tensor approximation problem California State University in the context of image processing. We test our approach [email protected], [email protected] on an image deblurring application. Ning Hao Tufts University AN10 Abstracts 59 Department of Mathematics Carmeliza Navasca [email protected] Clarkson University [email protected] Misha E. Kilmer Tufts University [email protected] MS36 On Multirate SSP Methods for Hyperbolic PDEs Carla Martin We discuss a set of explicit multirate time discretization James Madison University methods for hyperbolic conservation laws based on multi- Department of Mathematics step and Runge-Kutta methods. Multirate methods allow [email protected] different time steps to be used in different parts of the spa- tial domain while preserving the consistency and conser- MS35 vation properties of the ”classical” methods. Linear and nonlinear stability are guaranteed only under local CFL Iterative Methods for Multilinear Systems conditions. The necessity to take small global time steps We have shown a tensor group endowed with the Einstein restricted by the largest CFL number is thus avoided. The (contracted) product is isomorphic to the general linear theoretical results are illustrated on advection and Burgers group of degree n. In consequence, higher order tensor equations. inversion is possible for even order. Although odd order Emil M. Constantinescu tensors are not invertible under the Einstein product, their Argonne National Laboratory inversions are still possible through pseudo-inversion tech- Mathematics and Computer Science Division niques. With the notion of tensor inversion, multilinear [email protected] systems are solvable. Numerically we solve multilinear sys- tems arising from partial differential equations and quan- tum mechanics. Iterative methods for tensors are the main Adrian Sandu computational techniques. We demonstrate the efficacy of Virginia Tech these iterative methods through some motivating exam- [email protected] ples. Carmeliza Navasca MS36 Clarkson University The History of SSP Methods [email protected] Abstract not available at time of publication. Michael Brazell Sigal Gottlieb Clarkson University Department of Mathematics Dept. of Mechanical Engineering University of Massachusetts, Dartmouth [email protected] [email protected] Na Li Clarkson University MS36 Dept. of Mathematics On the Sharp SSP Timestep for Runge-Kutta Dis- [email protected] cretizations of IVPs Christino Tamon The product of the absolute monotonicity radius of the Clarkson University Runge-Kutta (RK) method and the strong stability pre- Department of Computer Science serving (SSP) timestep of the explicit Euler method is of- [email protected] ten used as a timestep limit that guarantees strong stability preservation of the RK method. In this talk we relax this limit in two ways. First, for certain subsets of the initial MS35 vector we prove that in the above limit a much larger factor Tensors as Module Homomorphisms over Group than absolute monotonicity rules the SSP property. Sec- Rings ond, for some classical RK methods we prove SSP property also in the case when the explicit Euler method is not SSP Abstract not available at time of publication. for any positive timestep. Finally, we demonstrate our re- sults by computational experiments. Michael Opperman Clarkson University Zoltan Horvath Department of Mathematics Sz´echenyi Istv´an University [email protected] [email protected] Timothy Penderghest MS36 Clarkson University [email protected] Multi-step Multi-stage SSP Runge–Kutta Methods High-order strong stability preserving (SSP) time integra- Christino Tamon tors of are often desirable in the numerical solution of hy- Clarkson University perbolic PDEs. Multistep Runge–Kutta methods general- Department of Computer Science ize both linear multistep methods and multistage Runge– [email protected] Kutta methods. We present some recent results on SSP multistep Runge–Kutta methods, including numerically 60 AN10 Abstracts optimal methods, both implicit and explicit. The meth- location. The spatial domain is decomposed into a series of ods are tested with numerical experiments. small subdomains (coarse grid) of either Stokes or Darcy type. The flow solution is resolved locally (on each coarse David I. Ketcheson element) on a fine grid, allowing for non-matching grids Mathematical and Computer Sciences & Engineering across subdomain interfaces. Coarse scale mortar finite King Abdullah University of Science & Technology elements are introduced on the interfaces to approximate [email protected] the normal stress and impose weakly continuity of flux. The transport equation is discretized via a local discon- Sigal Gottlieb tinuous Galerkin method. Computational experiments are Department of Mathematics presented. University of Massachusetts at Dartmouth [email protected] Benjamin Ganis University of Pittsburgh Colin Macdonald Department of Mathematics Oxford University [email protected] [email protected] Danail Vassilev University of Pittsburgh MS37 [email protected] Hydrodynamics Coupled with Bed Morphology for Modeling Sediment Scour and Deposition Ivan Yotov Univeristy of Pittsburgh Abstract unavailable at time of publication. Department of Mathematics [email protected] Clint Dawson Institute for Computational Engineering and Sciences University of Texas at Austin Ming Zhong [email protected] University of Pittsburgh [email protected] MS37 MS37 Adaptive Strategies in the Multilevel Multiscale Mimetic (M3) Method for Two-phase Flows in Stochastic Collocation Data Assimilation Porous Media In this work, a stochastic data assimilation approach is pre- The Multilevel Multiscale Mimetic (M3) method was de- sented for estimating model parameters such as permeabil- veloped to simulate accurately two-phase flows in highly ity and porosity from direct and indirect measurements. heterogeneous media with long correlation lengths. A mul- This approach combines the advantages of the ensemble tilevel approximation strategy is used to build a hierarchy Kalman filter (EnKF) for dynamic data assimilation and of models that are locally conservative at each level. The the polynomial chaos expansion (PCE) for efficient uncer- mimetic finite difference method is used to handle arbi- tainty quantification. In the latter, the model parameters trary meshes and tensor permeability fields. We describe are represented by the Karhunen-Loeve expansions and the adaptive strategies that reduce further the cost of the M3 model responses such as pressure and saturation are ex- method and illustrate them with numerical simulations of pressed by the PCE. The coefficients of PCE are solved well-driven flows. with a collocation technique. The approach is non-intrusive in that such realizations are solved forward in time via ex- Daniil Svyatskiy, Konstantin Lipnikov isting deterministic solver independently. The needed en- Los Alamos National Laboratory tries of the state covariance matrix are approximated with [email protected], [email protected] the PCE coefficients. It is shown that the approach is com- putationally efficient and provides satisfactory estimations David Moulton of the model parameters with dynamic measurements. Los Alamos National Laboratory Dongxiao Zhang Applied Mathematics and Plasma Physics USC [email protected] [email protected] MS37 MS38 Stochastic Multiscale Modeling of Coupled Adaptive Rational Krylov Subspaces for Time- Surface-subsurface Flow and Contaminant Trans- invariant Dynamical Systems port The rational Krylov space is recognized as a powerful tool We discuss a multiscale stochastic framework for uncer- within Model Order Reduction techniques for linear dy- tainty quantification in modeling flow and transport in namical systems. However, its success has been hindered surface-subsurface hydrological systems. The governing by the lack of procedures, which would generate the se- flow equations are the Stokes-Darcy system with Beavers- quence of shifts used to build the space with good approx- Joseph-Saffman interface conditions. The permeability in imation properties. We start with a-priori and adaptive the Darcy region is stochastic and it is represented with a approaches generating such shifts for the first order prob- Karhunen-Lo`eve (KL) expansion. The porous media can lems. Then a-priori approach is generalized to structure be statistically non-stationary, which is modeled by differ- preserving algorithms for high order systems. ent KL expansions in different regions. Statistical moments of the solution are computed via sparse grid stochastic col- Vladimir Druskin,L.Knizhnerman,C.Lieberman,V. AN10 Abstracts 61 Simoncini, M. Zaslavsky TU Braunschweig Schlumberger Doll Research [email protected] [email protected], [email protected], [email protected], [email protected], [email protected] Alexander Vendl TU Braunschweig, Germany Institut Computational Mathematics, AG Numerik MS38 [email protected] Robust Low-order Models for Flow Control The main issue discussed will be how to build a low-order Stefan G¨ortz, Ralf Zimmermann model that is robust to parameter variation. To do that an DLR Braunschweig, Germany optimal sampling strategy of the parameter space based on [email protected], [email protected] the full model residuals is proposed. In addition to that, the reduced model is calibrated using Tikhonov regulariza- MS39 tion. This system identification approach can accurately describe unsteady phenomena encountered in flow control. A Geometric Approach to Partition Functions Computing partition functions to get normalizing con- Angelo Iollo stants or marginals is an ubiquitous, yet generally in- Institut de Math´ematiques de Bordeaux tractable, task in statistical modelling. Many tricks and [email protected] workarounds have been developed to make such computa- tions feasible. I discuss how some of these strategies can be understood, and new ones manufactured, using tensor M Bergmann geometry. Universite Bordeaux [email protected] Jason Morton Pennsylvania State University E. Lombard [email protected] Universit´e Bordeaux 1 [email protected] MS39 Combinatorial Commutative Algebra in Action: MS38 Markov Moves from Classical Algebraic Construc- Model Reduction Algorithm for Robust Control of tions Linear PDE Systems Markov bases, as a fundamental object of study in algebraic We consider the computation of a robust control law for statistics, can often be difficult to construct and compute. large-scale finite dimensional linear systems and a class of It is often the case that a sub-basis may be sufficient for linear distributed parameter systems. The controller is ro- the problem at hand; however, theoretically, one would like bust with respect to left coprime factor perturbations of the to really understand the variety of the model, its algebraic nominal system. We present a convergent algorithm based and geometric properties. In this talk I will describe chal- on balanced proper orthogonal decomposition to compute lenges in Markov bases constructions, posed by the poorly the nonstandard features of this robust control law. Nu- understood geometry of the models, illustrated on an ex- merical results are presented for a convection diffusion par- ample. tial differential equation. Sonja Petrovic John Singler University of Illinois, Chicago Missouri S & T [email protected] Mathematics and Statistics [email protected] MS39 Belinda Batten Algebraic Statistics Framework for Causal Infer- Mechanical, Industrial, and Manufacturing Engineering ence and Data Privacy with Discrete Data Oregon State University [email protected] We present an algebraic computational framework that handles special cases of latent class analyses. Specifically, we consider discrete data problems with unobserved vari- MS38 ables such that arbitrary linear constraints are imposed on Model Order Reduction for Steady Aerodynamic the possible realizations of the complete data, and thus on Applications the possible states of the joint distribution of all the vari- ables (observed and unobserved) in the analysis. The con- An approach combining proper orthogonal decomposition straints are imposed either by the modeling assumptions, (POD) with linear regression, which is called gappy POD, the structure of the latent variables or for consistency rea- is used to obtain complete flow solutions in steady aero- sons. We illustrate our methods by applying them to two dynamic applications from knowledge of a (suitable) POD important related problems. The first problem pertains basis and the solution at very few points. In practice the to the assessment of disclosure risk of releasing potentially partial or gappy data can easily be obtained by experimen- sensitive information from a latent class analysis in the tally measuring the flow variables. Therefore this method form of class membership probabilities and probability dis- is effective in combining experimental with computational tribution of covariates conditional on the classes. The sec- data. ond problem pertains to estimation of average causal effect in presence of unobserved confounders , under the Neyman- Heike Fassbender Rubin framework of potential outcomes. The code is im- 62 AN10 Abstracts plemented in R, but interfaces with 4ti2. the approach to the Boussinesq system is also discussed. The resulting schemes are well suited for the simulation Aleksandra Slavkovic of moderate to high Reynolds and Rayleigh number flows. Penn State University Accuracy checks, simulations of the lid-driven cavity flow, [email protected] of a differentially heated cavity flow, and of a Rayleigh- Bernard convection problem for Rayleigh number up to 1010 are presented to demonstrate that the schemes are MS39 capable of producing accurate results at a reasonable com- Betti Numbers of Experimental Designs putational cost. For polynomial regression models in the context of factorial Hans Johnston experimental design there have been both formal and infor- Department of Mathematics & Statistics mal way of representing estimable effects and interactions. University of Massachusetts Amherst We generalize this to a larger study of the combinatorial [email protected] structure of models derived by the algebraic method, that is using ideals-of-points theory. It turns out, by exploit- Cheng Wang ing known results on monomial ideals, that the complexity University of Massachusetts of models can be measured by the size of the Betti num- Dartmouth, MA bers of the Stanley-Reisner complex. It transpires that the [email protected] most complex models are available for designs which are generic in the sense of ’corner cut’ theory. An complete enumeration for Plackett-Burman designs gives 19 equiva- Jian-guo Liu lence classes of complex. Dept. of Mathematics and Physics Duke University Henry Wynn [email protected] London School of Economics [email protected] MS40 Numerical Approxima- MS40 tion for Generalized-Newtonian Flows with Flow Dual-mixed Finite Element Methods for the Rate Boundary Conditions Navier-Stokes Equations We consider a generalized-Newtonian fluid with defective In many applications within mechanical, materials, and boundary conditions where only flow rate are prescribed on biomedical engineering, it is important to accurately pre- a part of boundary. The defect boundary condition prob- dict fluid stresses. However, most existing numerical lem is formulated as an minimization problem in which schemes for fluids are formulated with velocity as the pri- a Neumann or Dirichlet boundary control is used for the mary unknown of interest. In this talk, a dual-mixed finite flow rate matching. We will discuss the derivation of an element method for the Navier-Stokes equations, in which optimality system based on the first necessary condition the stress is a primary unknown of interest, is derived and and examine the corresponding adjoint problem. Com- analyzed. The method employs symmetric tensor finite putational algorithms and numerical results will be also elements for stress and is well-suited for non-Newtonian presented. fluids. Hyesuk Lee Jason Howell Clemson University Carnegie Mellon University Dep. of Mathematical Sciences Department of Mathematical Sciences [email protected] [email protected] MS40 Noel Walkington Carnegie Mellon University Optimal Control of Flow Past a Rectangle [email protected] We address optimal control of the frequency and amplitude of two-dimensional vortex-shedding behind a rectangular MS40 obstacle. Control is effected by small periodic oscillation of the rigid obstacle. Numerical results are presented. Local Boundary Condition based Spectral Colloca- tion Methods for 2D and 3D Naiver-Stokes Equa- Mike Sussman tions Department of Mathematics, University of Pittsburgh [email protected] We present a simple approach to accurately computing lo- cal boundary conditions in spectral collocation schemes for the Navier-Stokes equations in 2D and 3D. Access to these MS41 local boundary values makes possible the decoupling of Using Sequence Coverage Statistics to Determine the computation of the primary flow variables, resulting in Protein Binding Sites in a Genome highly efficient schemes. In 2D the local vorticity boundary values are employed in the vorticity-stream function for- Inspired by the notion of persistence in topological data mulation, while in 3D local values of a Neumann boundary analysis, we introduce a tree depicting sequence coverage condition for the pressure in the velocity-pressure formu- via fragment placement on a genome. We then describe lation are used. In both cases these boundary conditions statistically the trees that correspond to random fragment are approximated by differentiating a local Lagrange inter- placement and use this theory to determine the binding polant at the boundary. The straightforward extension of sites for a given protein in a genome. Our method for AN10 Abstracts 63 calling statistically signi?cant protein binding sites reduces introduction to some of these methods and applications, in- to the study of certain tree-based statistics derived from cluding some open problems. This talk (in fact, the entire the data. minisymposium) is intended to support the plenary lecture of Charles Wampler. Valerie Hower University of California, Berkeley Daniel J. Bates [email protected] Colorado State University Department of Mathematics [email protected] MS41 Lower and Upper Bounds on the Probability Dis- tributions of the Wasted Spaces of a Processor- MS42 Sharing Storage Allocation Model Regeneration and Numerical Algebraic Geometry We consider a storage allocation model with m primary Large-scale polynomial systems that often arise in kine- holding spaces, infinitely many secondary ones, and one matics and other engineering applications typically have processor servicing customers. We define the traffic inten- fewer solutions than their total degree root count. Regen- sity ρ to be λ/μ where λ is the customer arrival rate and eration can solve such large-scale polynomial systems by μ is the service rate of the processor. We study the lower algorithmically revealing the underlying structure of the and upper bounds on the probability distributions of the polynomial system during the equation-by-equation solv- wasted spaces and the largest index of the occupied spaces ing process. After describing regeneration, we will apply for fixed rho,0<ρ<1, and ρ =1− (heavy traffic case). it to polynomial systems arising in kinematics and other engineering applications. Eunju Sohn Jonathan Hauenstein University of Georgia Texas A&M University [email protected] [email protected] Andrew Sommese MS41 Department of Mathematics, University of Notre Dame Mechanisms of Simple Perceptual Decision-making Notre Dame, IN 46556-4618, USA Processes [email protected] Perceptual decision-making, an omnipresent component of everyday life, plays a pivotal role in cognitive tasks. In Charles Wampler this presentation, I will talk about mechanisms underlying General Motors Research Laboratories, Enterprise simple two-option perceptual decision-making processes by Systems Lab studying a biological-realistic reduced two-variable model 30500 Mound Road, Warren, MI 48090-9055, USA and phenomenological drift-diffusion models. [email protected] Xueying Wang Statistical and Applied Mathematical Science Institute MS42 [email protected] Basic Algebraic Geometry of Acoustic Arrays This talk is on the basic algebraic geometry of acoustic MS41 arrays. Given an array of microphones and a point P there Oscillations in NFkB Signaling Pathway is a set of time delays that, when applied to the signals collected by the microphones, lead to a focusing of the Upon stimulation, oscillations of NF-kB localization are microphones towards P. The locus of such coherent time observed at both single cell and population levels. A re- delays determine an algebraic variety. We will discuss this cent work reported that different frequencies of the oscil- variety, an associated moduli space and will discuss themes lations leads to different gene expression. Many authors in common with kinematic varieties. point out that NF-kB may interact with other pathways. However, the existence and mechanism of those potential Chris Peterson interactions are not clear. In this talk, we study this is- Colorado State University sue by considering the pathway subjected to two types of [email protected] putative signals: sinusoid and pulsatile signals. Michael Kirby Yunjiao Wang Colorado State University Ohio State University Department of Mathematics [email protected] [email protected] Josh Thompson MS42 Colorado State University Numerical Methods for Solving Polynomial Sys- [email protected] tems The term numerical algebraic geometry describes a variety MS42 of numerical methods for finding and manipulating the so- Numerical Matrix Computation in Algebraic Ge- lution sets of polynomial systems. There are numerous ap- ometry plications throughout the sciences and engineering as well as many within mathematics. This talk will serve as a brief Numerical polynomial algebra and numerical algebraic ge- 64 AN10 Abstracts ometry become fast growing fields of studies where tremen- MS43 dous progress has been achieved and new challenges emerge Numerical Simulations of Radiative Shock Tube (e.g. handling ill-posedness). On the other hand, numeri- Experiments: Challenges and Approaches cal polynomial algebra is a natural extension and an appli- cation area of numerical linear algebra, as matrix compu- The Center for Radiative Shock Hydrodynamics (CRASH) tation plays an indispensible role. This talk introduces the investigates ways to improve the predictive capability of matrix computation strategies in numerical polynomial al- models for shock waves produced when a laser is used gebra and numerical algebraic geometry using several case to shock, ionize, and accelerate a beryllium plate into a studies along with computing results. xenon-filled shock tube. These shocks, when driven above a threshold velocity of about 100 km/s, become strongly Zhonggang Zeng radiative and convert most of the incoming energy flux to Northeastern Illinois University radiation. What results is a complex, evolving, radiation- Department of Mathematics hydrodynamic, multiphysics environment. I will discuss [email protected] the challenges and approaches we use to model this system as a tractable problem in large-scale parallel computing. MS43 Eric S. Myra Hybrid Monte Carlo Methods for Kinetic Equa- University of Michigan tions [email protected] Probabilistic techniques such as DSMC (Direct Simulation Monte Carlo) are extensively used in real simulations of MS43 the Boltzmann equation for their great flexibility, ability High-order Discontinuous Galerkin Schemes for in treating different collision terms and low computational 2+2 Vlasov Models on Unstructured Grids cost compared to any type of deterministic scheme. On the other hand, solutions are affected by large fluctuations The purpose of this work is to explore mesh-based alterna- and, in non stationary situations, the impossibility of aver- tives to the Particle-In-Cell (PIC) methods that are cur- aging quantities leads to, or low accurate solutions, or high rently favored in many plasma physics applications. In costly simulations. Moreover, close to thermodynamical particular, we present here a technique for solving the equilibrium, the cost of Monte Carlo methods increases. In Vlasov-Poisson system that is based on operator splitting this talk, emphasis will be addressed to a recent developed in time and high-order discontinuous Galerkin discretiza- new family of Hybrid Monte Carlo methods which permits tions in space. The efficiency of the method is increased by to treat both accurately and fast transitional regimes and considering semi-Lagrangian approaches for the advection multiscale problems. The key aspects, on which methods in phase space. The flexibility of the approach is enhanced rely, are the choice of a suitable hybrid representation of by allowing the mesh in physical space to be unstructured. the solution and a merging of Monte Carlo methods in non- equilibrium regimes with deterministic methods in equilib- rium ones. James A. Rossmanith,DavidSeal University of Wisconsin Giacomo Dimarco Department of Mathematics Institut de Math´ematiques de Toulouse [email protected], [email protected] [email protected] Lorenzo Pareschi MS44 Department of Mathematics Asynchronous Time Integration for Periodic Dy- University of Ferrara namical Systems with Uncertain Parameters [email protected] Application of PC expansions to stochastic dynamical sys- tems exhibiting cyclic oscillations frequently leads to a MS43 broadening solution spectrum, and consequently to exces- An Asymptotic Preserving Scheme for Non- sive computational requirements. This talk outlines an classical Particle Transport asynchronous integration framework that offers to substan- tially mitigate this problem. The approach is based on a We present a numerical scheme to solve the nonclassical stochastic time rescaling which is implemented so that the particle transport model suggested by Larsen [Larsen, A spectrum of the rescaled solution remains narrow-banded. generalized Boltzmann equation for non-classical particle Different variants of this approach are illustrated for simple transport, ANS M&C 2007]. This model has applications systems having almost surely a stable limit cycle, including in radiative transfer through clouds and novel nuclear reac- a stochastic linear oscillator, and a stiff nonlinear chemical tor types. While basically any scheme for linear transport system with uncertain rate parameters. has a diffusion limit (albeit with the wrong diffusion coeffi- cient), we show here that the discretization of the transport Omar M. Knio part has to be chosen carefully so that the scheme has a Dept. of Mech. Eng. (correct) diffusion limit. Johns Hopkins University [email protected] Martin Frank RWTH Aachen University Olivier LeMaitre Center for Computational Engineering Science LIMSI-CNRS [email protected] [email protected] Lionel Mathelin LIMSI - CNRS AN10 Abstracts 65 [email protected] Sandia National Laboratories, Livermore CA [email protected] M. Yousuff Hussaini Florida State University Sonday/Najm Sonday/Najm [email protected] Princeton [email protected] MS44 Phase Conditions for Autonomous Oscillators with MS44 Random Parameters Uncertainty Quantification via Codimension One Domain Decomposition and a New Concentration We consider autonomous systems of ordinary differen- Inequality tial equations (ODEs) or differential algebraic equations (DAEs), which exhibit periodic solutions. Since a contin- We propose a localized variant of McDiarmid’s inequality uum of periodic solutions exists, we require a phase condi- in order to obtain upper bounds on the probability that a tion to isolate a particular solution. We assume uncertain- system of interest assumes certain values (“fails’). By par- ties in some parameters of the autonomous systems. Hence titioning the input parameter space appropriately, much the relevant parameters are replaced by random variables. sharper bounds than the usual McDiarmid bound are ob- We apply the technique of the generalized polynomial chaos tained. We prove an error estimate for the method, define to resolve the stochastic model. A Galerkin approach yields a codimension one recursive partitioning scheme and prove a larger coupled system of ODEs or DAEs. The phase con- its convergence properties, and use a new concentration in- ditions of the original systems imply additional boundary equality to give confidence levels when empirical means are conditions in the periodic problems of the larger coupled used in place of exact ones. systems. However, the expected value and the variance of the solution depend on the choice of the phase condition. Timothy J. Sullivan We construct an alternative constraint, which minimizes California Institute of Technology the variance of the solution. Numerical simulations of a [email protected] test example are presented. Houman Owhadi Roland Pulch Applied Mathematics University of Wuppertal Caltech [email protected] [email protected] MS44 MS45 Eigenvalue Analysis of Uncertain ODE Systems Operator Splitting Methods for Maxwell’s Equa- tions in Dispersive Media The polynomial chaos expansion provides a means of rep- 2 resenting any L random variable as a sum of polynomi- We consider Maxwell’s equations in dispersive media of De- als that are orthogonal with respect to a chosen measure. bye type. In such relaxing dielectric media, the presence Examples include the Hermite polynomials with Gaussian of different wave speeds leads to stiffness in the temporal measure on the real line and the Legendre polynomials with domain. We present an operator splitting scheme that de- uniform measure on an interval. Polynomial chaos can couples fast and slow moving processes in the problem to be used to reformulate an uncertain ODE system, using develop separate sub-problems. This alleviates the strin- Galerkin projection, as a new, higher-dimensional, deter- gent requirements on the time-step which along with sta- ministic ODE system which describes the evolution of each bility conditions requires small spatial steps and hence ex- mode of the polynomial chaos expansion. It is of interest to cessive computations for long-time integration of Maxwell’s explore the eigenstructure of the original and reformulated equations. ODE systems, by studying the eigenvalues and eigenvectors of their Jacobians. In this talk, we study the distribution Vrushali A. Bokil of the eigenvalues of the two Jacobians. We outline in gen- Oregon State University eral the location of the eigenvalues of the new system with [email protected] respect to those of the original system, and examine the effect of expansion order on this distributio MS45 Benjamin Sonday Overlapping Domain Decomposition Method for Princeton University the Helmholtz Equation [email protected] The model problem addressed in this talk concerns the Robert Berry analysis and computation of a radiated or scattered time- Sandia National Laboratories harmonic acoustic solution, where the obstacle is composed [email protected] of dielectric and metal. A new adaptive radiation condition method, that localizes the artificial interface only around Habib N. Najm the dielectric object, is described. An appropriate algo- Sandia National Laboratories rithm coupling finite and boundary element methods is Livermore, CA, USA solved iteratively using an overlapping domain decompo- [email protected] sition method. Numerical results are presented validating the theoretical results. Bert J. Debusschere Yassine Boubendir Energy Transportation Center Department of Mathematical Sciences 66 AN10 Abstracts [email protected] very ill-conditioned. Based on our analysis, we propose a regularization of the local model problem that leads to much more effective steps. We demonstrate our approach MS45 on nonlinear least-squares problems arising in shape-based A Numerical Method for Simulating EM Wave regularization for diffuse optical tomography. Propagation in Dielectrics that Exhibit Fractional Relaxation Eric De Sturler Virginia Tech The frequency-dependent Havriliak-Negami dielectric per- [email protected] s−∞ mittivity model, (ω)=∞+ (1+(iωτ)α)β ,where0<α,β< 1, generates the entire class of dielectric models used in nu- Misha E. Kilmer merical simulations of time-domain propagation and scat- Tufts University tering of electromagnetic waves. We present a numerical [email protected] method to incorporate this model in a solver for Maxwell’s equations which now contain fractional differential opera- tors. We give a complete stability and error analysis of the MS46 method and a validation against the exact solution avail- A Numerical Scheme for Monge-Kantorovich Mass able in some simple cases. Transfer Problem, in the Context of Image Regis- tration Matthew F. Causley New Jersey Institute of Technology (NJIT) The Monge-Kantorovich mass transfer problem addresses [email protected] the question of how to move a pile of soil from one lo- cation to another with minimal cost. In this work we Peter Petropoulos introduce computationally efficient numerical method for Department of Mathematical Sciences non-rigid image registration based on the optimal mass New Jersey Institute of Technology transport theory. We will consider the formulation of [email protected] Banemou and Brenier (2000), which requires the solution of a time-dependent PDE constrained optimization, with mass-preserving PDE. MS45 Raya Horesh Polynomial Chaos Approach for Approximating Emory University Cole-Cole Dispersive Media Atlanta, GA Time-domain simulations involving the Cole-Cole model [email protected] are not straight-forward as the model corresponds to a frac- tional order ODE. We introduce an alternative approach Eldad Haber based on using the first order linear ODE Debye model, Department of Mathematics but with distributions of relaxation times. We apply gen- The University of British Columbia eralized Polynomial Chaos, and then discretize the system [email protected] coupled with Maxwell’s equations. We present a stability and dispersion analysis of the overall method and discuss the impact of the variance of relaxation times. MS46 Iteratively Regularized BFGS-type Algorithm for Nathan L. Gibson an Inverse Problem in Groundwater Hydrology Oregon State University Department of Mathematics A novel iteratively regularized BFGS-type algorithm with [email protected] simultaneous updates of the operator