Non-linear Dynamics Commons^™

Lecture 09: Hierarchically Low Rank And Kronecker Methods, Rio Yokota

Mathematical Sciences Spring Lecture Series

Exploiting structures of matrices goes beyond identifying their non-zero patterns. In many cases, dense full-rank matrices have low-rank submatrices that can be exploited to construct fast approximate algorithms. In other cases, dense matrices can be decomposed into Kronecker factors that are much smaller than the original matrix. Sparsity is a consequence of the connectivity of the underlying geometry (mesh, graph, interaction list, etc.), whereas the rank-deficiency of submatrices is closely related to the distance within this underlying geometry. For high dimensional geometry encountered in data science applications, the curse of dimensionality poses a challenge for rank-structured approaches. On the other …

Lecture 07: Nonlinear Preconditioning Methods And Applications, Xiao-Chuan Cai

Mathematical Sciences Spring Lecture Series

We consider solving system of nonlinear algebraic equations arising from the discretization of partial differential equations. Inexact Newton is a popular technique for such problems. When the nonlinearities in the system are well-balanced, Newton's method works well, but when a small number of nonlinear functions in the system are much more nonlinear than the others, Newton may converge slowly or even stagnate. In such a situation, we introduce some nonlinear preconditioners to balance the nonlinearities in the system. The preconditioners are often constructed using a combination of some domain decomposition methods and nonlinear elimination methods. For the nonlinearly preconditioned problem, …

Lecture 02: Tile Low-Rank Methods And Applications (W/Review), David Keyes

Mathematical Sciences Spring Lecture Series

As simulation and analytics enter the exascale era, numerical algorithms, particularly implicit solvers that couple vast numbers of degrees of freedom, must span a widening gap between ambitious applications and austere architectures to support them. We present fifteen universals for researchers in scalable solvers: imperatives from computer architecture that scalable solvers must respect, strategies towards achieving them that are currently well established, and additional strategies currently being developed for an effective and efficient exascale software ecosystem. We consider recent generalizations of what it means to “solve” a computational problem, which suggest that we have often been “oversolving” them at the …

Five-Wave Resonances In Deep Water Gravity Waves: Integrability, Numerical Simulations And Experiments, Dan Lucas, Marc Perlin, Dian-Yong Liu, Shane Walsh, Rossen Ivanov, Miguel D. Bustamante

Articles

In this work we consider the problem of finding the simplest arrangement of resonant deep water gravity waves in one-dimensional propagation, from three perspectives: Theoretical, numerical and experimental. Theoretically this requires using a normal-form Hamiltonian that focuses on 5-wave resonances. The simplest arrangement is based on a triad of wave vectors K₁ + K₂ = K₃ (satisfying specific ratios) along with their negatives, corresponding to a scenario of encountering wave packets, amenable to experiments and numerical simulations. The normal-form equations for these encountering waves in resonance are shown to be non-integrable, but they admit an integrable reduction …

Numerical Computations Of Vortex Formation Length In Flow Past An Elliptical Cylinder, Matthew Karlson, Bogdan Nita, Ashwin Vaidya

Department of Mathematics Facuty Scholarship and Creative Works

We examine two dimensional properties of vortex shedding past elliptical cylinders through numerical simulations. Specifically, we investigate the vortex formation length in the Reynolds number regime 10 to 100 for elliptical bodies of aspect ratio in the range 0.4 to 1.4. Our computations reveal that in the steady flow regime, the change in the vortex length follows a linear profile with respect to the Reynolds number, while in the unsteady regime, the time averaged vortex length decreases in an exponential manner with increasing Reynolds number. The transition in profile is used to identify the critical Reynolds number which marks the …

Teaching And Learning Of Fluid Mechanics, Ashwin Vaidya

Department of Mathematics Facuty Scholarship and Creative Works

Fluid mechanics occupies a privileged position in the sciences; it is taught in various science departments including physics, mathematics, environmental sciences and mechanical, chemical and civil engineering, with each highlighting a different aspect or interpretation of the foundation and applications of fluids. Doll’s fluid analogy [5] for this idea is especially relevant to this issue: “Emergence of creativity from complex flow of knowledge—example of Benard convection pattern as an analogy—dissipation or dispersal of knowledge (complex knowledge) results in emergent structures, i.e., creativity which in the context of education should be thought of as a unique way to arrange information so …

One-Note-Samba Approach To Cosmology, Florentin Smarandache, Victor Christianto

Branch Mathematics and Statistics Faculty and Staff Publications

Inspired by One Note Samba, a standard jazz repertoire, we present an outline of Bose-Einstein Condensate Cosmology. Although this approach seems awkward and a bit off the wall at first glance, it is not impossible to connect altogether BEC, Scalar Field Cosmology and Feshbach Resonance with Ermakov-Pinney equation. We also briefly discuss possible link with our previous paper which describes Newtonian Universe with Vortex in terms of Ermakov equation.

School Policy Evaluated With Time-Reversible Markov Chain, Trajan Murphy, Iddo Ben-Ari

Honors Scholar Theses

In this work we propose a reversible Markov chain scheme to model for the mobility of students affected by a grade school leveling policy. This model provides unified and mathematically tractable framework in which transition functions are sampled uniformly from the set of {\bf reversible} transition functions. The results from the study appear to confirm the disadvantageous effects of this school policy, on par with the of a previous model on the same policy.

Parameter Estimation And Optimal Design Techniques To Analyze A Mathematical Model In Wound Healing, Nigar Karimli

Masters Theses & Specialist Projects

For this project, we use a modified version of a previously developed mathematical model, which describes the relationships among matrix metalloproteinases (MMPs), their tissue inhibitors (TIMPs), and extracellular matrix (ECM). Our ultimate goal is to quantify and understand differences in parameter estimates between patients in order to predict future responses and individualize treatment for each patient. By analyzing parameter confidence intervals and confidence and prediction intervals for the state variables, we develop a parameter space reduction algorithm that results in better future response predictions for each individual patient. Moreover, use of another subset selection method, namely Structured Covariance Analysis, that …

Surface Waves Over Currents And Uneven Bottom, Alan Compelli, Rossen Ivanov, Calin I. Martin, Michail D. Todorov

Articles

The propagation of surface water waves interacting with a current and an uneven bottom is studied. Such a situation is typical for ocean waves where the winds generate currents in the top layer of the ocean. The role of the bottom topography is taken into account since it also influences the local wave and current patterns. Specific scaling of the variables is selected which leads to approximations of Boussinesq and KdV types. The arising KdV equation with variable coefficients, dependent on the bottom topography, is studied numerically when the initial condition is in the form of the one soliton solution …

Simplicity And Sustainability: Pointers From Ethics And Science, Mehrdad Massoudi, Ashwin Vaidya

Department of Mathematics Facuty Scholarship and Creative Works

In this paper, we explore the notion of simplicity. We use definitions of simplicity proposed by philosophers, scientists, and economists. In an age when the rapidly growing human population faces an equally rapidly declining energy/material resources, there is an urgent need to consider various notions of simplicity, collective and individual, which we believe to be a sensible path to restore our planet to a reasonable state of health. Following the logic of mathematicians and physicists, we suggest that simplicity can be related to sustainability. Our efforts must therefore not be spent so much in pursuit of growth but in achieving …

Iterative Methods To Solve Systems Of Nonlinear Algebraic Equations, Md Shafiful Alam

Masters Theses & Specialist Projects

Iterative methods have been a very important area of study in numerical analysis since the inception of computational science. Their use ranges from solving algebraic equations to systems of differential equations and many more. In this thesis, we discuss several iterative methods, however our main focus is Newton's method. We present a detailed study of Newton's method, its order of convergence and the asymptotic error constant when solving problems of various types as well as analyze several pitfalls, which can affect convergence. We also pose some necessary and sufficient conditions on the function f for higher order of convergence. Different …

Flow Anisotropy Due To Thread-Like Nanoparticle Agglomerations In Dilute Ferrofluids, Alexander Cali, Wah-Keat Lee, A. David Trubatch, Philip Yecko

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Improved knowledge of the magnetic field dependent flow properties of nanoparticle-based magnetic fluids is critical to the design of biomedical applications, including drug delivery and cell sorting. To probe the rheology of ferrofluid on a sub-millimeter scale, we examine the paths of 550 μm diameter glass spheres falling due to gravity in dilute ferrofluid, imposing a uniform magnetic field at an angle with respect to the vertical. Visualization of the spheres’ trajectories is achieved using high resolution X-ray phase-contrast imaging, allowing measurement of a terminal velocity while simultaneously revealing the formation of an array of long thread-like accumulations of magnetic …

Region Of Smooth Functions For Positive Solutions To An Elliptic Biological Model, Joon Hyuk Kang, Timothy Robertson

Faculty Publications

The non-existence and existence of the positive solution to the generalized elliptic model ∆u+g(u v) = 0 in Ω, ∆v+h(u, v) = 0 in Ω, u=v= 0 on∂Ω, were investigated.

Integrability, Recursion Operators And Soliton Interactions, Boyka Aneva, Georgi Grahovski, Rossen Ivanov, Dimitar Mladenov

Book chapter/book

This volume contains selected papers based on the talks,presentedat the Conference Integrability, Recursion Operators and Soliton Interactions, held in Sofia, Bulgaria (29-31 August 2012) at the Institute for Nuclear Research and Nuclear Energy of the Bulgarian Academy of Sciences. Included are also invited papers presenting new research developments in the thematic area. The Conference was dedicated to the 65-th birthday of our esteemed colleague and friend Vladimir Gerdjikov. The event brought together more than 30 scientists, from 6 European countries to celebrate Vladimir's scientific achievements. All participants enjoyed a variety of excellent talks in a friendly and stimulating atmosphere. …

Modeling And Mathematical Analysis Of Plant Models In Ecology, Eric A. Eager

Department of Mathematics: Dissertations, Theses, and Student Research

Population dynamics tries to explain in a simple mechanistic way the variations of the size and structure of biological populations. In this dissertation we use mathematical modeling and analysis to study the various aspects of the dynamics of plant populations and their seed banks.

In Chapter 2 we investigate the impact of structural model uncertainty by considering different nonlinear recruitment functions in an integral projection model for Cirsium canescens. We show that, while having identical equilibrium populations, these two models can elicit drastically different transient dynamics. We then derive a formula for the sensitivity of the equilibrium population to …

Singular Solutions Of Coss-Coupled Epdiff Equations: Waltzing Peakons And Compacton Pairs, Colin Cotter, Darryl Holm, Rossen Ivanov, James Percival

Conference papers

We introduce EPDiff equations as Euler-Poincare´ equations related to Lagrangian provided by a metric, invariant under the Lie Group Diff(Rⁿ). Then we proceed with a particular form of EPDiff equations, a cross coupled two-component system of Camassa-Holm type. The system has a new type of peakon solutions, 'waltzing' peakons and compacton pairs.

Controlling Nanoparticles Formation In Molten Metallic Bilayers By Pulsed-Laser Interference Heating, Mikhail Khenner, Sagar Yadavali, Ramki Kalyanaraman

Mathematics Faculty Publications

The impacts of the two-beam interference heating on the number of core-shell and embedded nanoparticles and on nanostructure coarsening are studied numerically based on the non-linear dynamical model for dewetting of the pulsed-laser irradiated, thin (< 20 nm) metallic bilayers. The model incorporates thermocapillary forces and disjoining pressures, and assumes dewetting from the optically transparent substrate atop of the reflective support layer, which results in the complicated dependence of light reflectivity and absorption on the thicknesses of the layers. Stabilizing thermocapillary effect is due to the local thickness-dependent, steady- state temperature profile in the liquid, which is derived based on the mean substrate temperature estimated from the elaborate thermal model of transient heating and melting/freezing. Linear stability analysis of the model equations set for Ag/Co bilayer predicts the dewetting length scales in the qualitative agreement with experiment.

Formation Of Organized Nanostructures From Unstable Bilayers Of Thin Metallic Liquids, Mikhail Khenner, Sagar Yadavali, Ramki Kalyanaraman

Mathematics Faculty Publications

Dewetting of pulsed-laser irradiated, thin (< 20 nm), optically reflective metallic bilayers on an optically transparent substrate with a reflective support layer is studied within the lubrication equations model. A steady-state bilayer film thickness (h) dependent temperature profile is derived based on the mean substrate temperature estimated from the elaborate thermal model of transient heating and melting/freezing. Large thermocapillary forces are observed along the plane of the liquid-liquid and liquid-gas interfaces due to this h-dependent temperature, which, in turn, is strongly influenced by the h-dependent laser light reflection and absorption. Consequently the dewetting is a result of the competition between thermocapillary and intermolecular forces. A linear analysis of the dewetting length scales established that the non-isothermal calculations better predict the experimental results as compared to the isothermal case within the bounding Hamaker coefficients. Subsequently, a computational non-linear dynamics study of the dewetting pathway was performed for Ag/Co and Co/Ag bilayer systems to predict the morphology evolution. We found that the systems evolve towards formation of different morphologies, including core-shell, embedded, or stacked nanostructure morphologies.

Stability Of A Strongly Anisotropic Thin Epitaxial Film In A Wetting Interaction With Elastic Substrate, Mikhail Khenner, Wondimu T. Tekalign, Margo S. Levine

Mathematics Faculty Publications

The linear dispersion relation for longwave surface perturbations, as derived by Levine et al. Phys. Rev. B 75, 205312 (2007) is extended to include a smooth surface energy anisotropy function with a variable anisotropy strength (from weak to strong, such that sharp corners and slightly curved facets occur on the corresponding Wulff shape). Through detailed parametric studies it is shown that a combination of a wetting interaction and strong anisotropy, and even a wetting interaction alone results in complicated linear stability characteristics of strained and unstrained films.

Stability Of A Strongly Anisotropic Thin Epitaxial Film In A Wetting Interaction With Elastic Substrate, Mikhail Khenner, Wondimu Tekalign, Margo Levine

Mathematics Faculty Publications

The linear dispersion relation for longwave surface perturbations, as derived by Levine et al. Phys. Rev. B 75, 205312 (2007) is extended to include a smooth surface energy anisotropy function with a variable anisotropy strength (from weak to strong, such that sharp corners and slightly curved facets occur on the corresponding Wulff shape). Through detailed parametric studies it is shown that a combination of a wetting interaction and strong anisotropy, and even a wetting interaction alone results in complicated linear stability characteristics of strained and unstrained films.

Modeling Diverse Physics Of Nanoparticle Self-Assembly In Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

Presents physics behind dewetting of thin liquid films and mathematical/computational modeling tools (Educational/Research presentation for senior physics majors).

Morphological Evolution Of Single-Crystal Ultrathin Solid Films, Mikhail Khenner

Mathematics Faculty Publications

An introduction to mathematical modeling of ultrathin solid films and the role of such modeling in nanotechnologies: Educational presentation for senior physics majors

Morphological Evolution Of Single-Crystal Ultrathin Solid Films, Mikhail Khenner

Mathematics Faculty Publications

An introduction to mathematical modeling of ultrathin solid films and the role of such modeling in nanotechnologies: Educational/Research presentation for senior physics majors

Oscillatory And Monotonic Modes Of Long-Wave Marangoni Convection In A Thin Film, Sergey Shklyaev, Mikhail Khenner, Alexei Alabuzhev

Mathematics Faculty Publications

We study long-wave Marangoni convection in a layer heated from below. Using the scaling k=O Bi, where k is the wave number and Bi is the Biot number, we derive a set of amplitude equations. Analysis of this set shows presence of monotonic and oscillatory modes of instability. Oscillatory mode has not been previously found for such direction of heating. Studies of weakly nonlinear dynamics demonstrate that stable steady and oscillatory patterns can be found near the stability threshold.

Oscillatory And Monotonic Modes Of Long-Wave Marangoni Convection In A Thin Film, Sergey Shklyaev, Mikhail Khenner, Alexei Alabuzhev

Mathematics Faculty Publications

We study long-wave Marangoni convection in a layer heated from below. Using the scaling k=O Bi, where k is the wave number and Bi is the Biot number, we derive a set of amplitude equations. Analysis of this set shows presence of monotonic and oscillatory modes of instability. Oscillatory mode has not been previously found for such direction of heating. Studies of weakly nonlinear dynamics demonstrate that stable steady and oscillatory patterns can be found near the stability threshold.

A Generalized Nonlinear Model For The Evolution Of Low Frequency Freak Waves, Jonathan Blackledge

Articles

This paper presents a generalized model for simulating wavefields associated with the sea surface. This includes the case when `freak waves' may occur through an effect compounded in the nonlinear (cubic) Schrodinger equation. After providing brief introductions to linear sea wave models, `freak waves' and the linear and nonlinear Schrodinger equations, we present a unified model that provides for a piecewise continuous transition from a linear to a nonlinear state. This is based on introducing a fractional time derivative to develop a fractional nonlinear partial differential equation with a stochastic source function. In order to explore the characteristics of this …

From Energy Gradient And Natural Selection To Biodiversity And Stability Of Ecosystems, Bo Deng

Department of Mathematics: Faculty Publications

The purpose of this paper is to incorporate well-established ecological principles into a foodweb model consisting of four trophic levels --- abiotic resources, plants, herbivores, and carnivores. The underlining principles include Kimura's neutral theory of genetic evolution, Liebig's Law of the Minimum for plant growth, Holling's functionals for herbivore foraging and carnivore predation, the One-Life Rule for all organisms, and Lotka-Volterra's model for intraand interspecific competitions. Numerical simulations of the model led to the following statistical findings: (a) particular foodwebs can give contradicting observations on biodiversity and productivity, in particular, all known functional forms -- - positive, negative, sigmoidal, and …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

A mathematical model for the evolution of pulsed laser-irradiated, molten metallic films has been developed using the lubrication theory. The heat transfer problem that incorporates the absorbed heat from a single laser beam or the interfering laser beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the reflectivity, the peak laser beam …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

A mathematical model for the evolution of pulsed laser-irradiated, molten metallic films has been developed using the lubrication theory. The heat transfer problem that incorporates the absorbed heat from a single laser beam or the interfering laser beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the reflectivity, the peak laser beam …