Physics Commons^™

Year-2 Progress Report On Numerical Methods For Bgk-Type Kinetic Equations, Steven M. Wise, Evan Habbershaw

Faculty Publications and Other Works -- Mathematics

In this second progress report we expand upon our previous report and preliminary work. Specifically, we review some work on the numerical solution of single- and multi-species BGK-type kinetic equations of particle transport. Such equations model the motion of fluid particles via a density field when the kinetic theory of rarefied gases must be used in place of the continuum limit Navier-Stokes and Euler equations. The BGK-type equations describe the fluid in terms of phase space variables, and, in three space dimensions, require 6 independent phase-space variables (3 for space and 3 for velocity) for each species for accurate simulation. …

The Lagrangian Formulation For Wave Motion With A Shear Current And Surface Tension, Conor Curtin, Rossen Ivanov

Articles

The Lagrangian formulation for the irrotational wave motion is straightforward and follows from a Lagrangian functional which is the difference between the kinetic and the potential energy of the system. In the case of fluid with constant vorticity, which arises for example when a shear current is present, the separation of the energy into kinetic and potential is not at all obvious and neither is the Lagrangian formulation of the problem. Nevertheless, we use the known Hamiltonian formulation of the problem in this case to obtain the Lagrangian density function, and utilising the Euler-Lagrange equations we proceed to derive some …

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 …

Swirling Fluid Flow In Flexible, Expandable Elastic Tubes: Variational Approach, Reductions And Integrability, Rossen Ivanov, Vakhtang Putkaradze

Articles

Many engineering and physiological applications deal with situations when a fluid is moving in flexible tubes with elastic walls. In real-life applications like blood flow, a swirl in the fluid often plays an important role, presenting an additional complexity not described by previous theoretical models. We present a theory for the dynamics of the interaction between elastic tubes and swirling fluid flow. The equations are derived using a variational principle, with the incompressibility constraint of the fluid giving rise to a pressure-like term. In order to connect this work with the previous literature, we consider the case of inextensible and …

On The Intermediate Long Wave Propagation For Internal Waves In The Presence Of Currents, Joseph Cullen, Rossen Ivanov

Articles

A model for the wave motion of an internal wave in the presence of current in the case of intermediate long wave approximation is studied. The lower layer is considerably deeper, with a higher density than the upper layer. The flat surface approximation is assumed. The fluids are incompressible and inviscid. The model equations are obtained from the Hamiltonian formulation of the dynamics in the presence of a depth-varying current. It is shown that an appropriate scaling leads to the integrable Intermediate Long Wave Equation (ILWE). Two limits of the ILWE leading to the integrable Benjamin-Ono and KdV equations are …

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 …

Equatorial Wave–Current Interactions, Adrian Constantin, Rossen Ivanov

Articles

We study the nonlinear equations of motion for equatorial wave–current interactions in the physically realistic setting of azimuthal two-dimensional inviscid flows with piecewise constant vorticity in a two-layer fluid with a flat bed and a free surface. We derive a Hamiltonian formulation for the nonlinear governing equations that is adequate for structure-preserving perturbations, at the linear and at the nonlinear level. Linear theory reveals some important features of the dynamics, highlighting differences between the short- and long-wave regimes. The fact that ocean energy is concentrated in the long-wave propagation modes motivates the pursuit of in-depth nonlinear analysis in the long-wave …

Godunov-Type Upwind Flux Schemes Of The Two-Dimensional Finite Volume Discrete Boltzmann Method, Leitao Chen, Laura Schaefer

Publications

A simple unified Godunov-type upwind approach that does not need Riemann solvers for the flux calculation is developed for the finite volume discrete Boltzmann method (FVDBM) on an unstructured cell-centered triangular mesh. With piecewise-constant (PC), piecewise-linear (PL) and piecewise-parabolic (PP) reconstructions, three Godunov-type upwind flux schemes with different orders of accuracy are subsequently derived. After developing both a semi-implicit time marching scheme tailored for the developed flux schemes, and a versatile boundary treatment that is compatible with all of the flux schemes presented in this paper, numerical tests are conducted on spatial accuracy for several single-phase flow problems. Four major …

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 …

Introduction To The Usu Library Of Solutions To The Einstein Field Equations, Ian M. Anderson, Charles G. Torre

Tutorials on... in 1 hour or less

This is a Maple worksheet providing an introduction to the USU Library of Solutions to the Einstein Field Equations. The library is part of the DifferentialGeometry software project and is a collection of symbolic data and metadata describing solutions to the Einstein equations.

Evolution Of Superoscillations For Schrödinger Equation In A Uniform Magnetic Field, Fabrizio Colombo, Jonathan Gantner, Daniele C. Struppa

Mathematics, Physics, and Computer Science Faculty Articles and Research

Aharonov-Berry superoscillations are band-limited functions that oscillate faster than their fastest Fourier component. Superoscillations appear in several fields of science and technology, such as Aharonov’s weak measurement in quantum mechanics, in optics, and in signal processing. An important issue is the study of the evolution of superoscillations using the Schrödinger equation when the initial datum is a weak value. Some superoscillatory functions are not square integrable, but they are real analytic functions that can be extended to entire holomorphic functions. This fact leads to the study of the continuity of a class of convolution operators acting on suitable spaces of …

Interplay Of Quantum Size Effect, Anisotropy And Surface Stress Shapes The Instability Of Thin Metal Films, Mikhail Khenner

Mathematics Faculty Publications

Morphological instability of a planar surface ([111], [011], or [001]) of an ultra-thin metal film is studied in a parameter space formed by three major effects (the quantum size effect, the surface energy anisotropy and the surface stress) that influence a film dewetting. The analysis is based on the extended Mullins equation, where the effects are cast as functions of the film thickness. The formulation of the quantum size effect (Z. Zhang et al., PRL 80, 5381 (1998)) includes the oscillation of the surface energy with thickness caused by electrons confinement. By systematically comparing the effects, their contributions into the …

Numerically Solving A System Of Pdes Modeling Chronic Wounds Treated With Oxygen Therapy, Stefan Stryker

Mahurin Honors College Capstone Experience/Thesis Projects

Chronic wounds such as diabetic foot ulcers are the leading cause of non-traumatic amputations in developed countries. For researchers to better understand the physiology of these wounds, a mathematical model describing oxygen levels at the wound site can be used to help predict healing responses. The model utilizes equations that are modified from work by Guffey (2015) that consists of four variables – oxygen, bacteria, neutrophils, and chemoattractant within a system of partial differential equations. Our research focuses on numerically solving these partial differential equations using a finite volume approach. This numerical solver will be important for future research in …

An Rbf Interpolation Blending Scheme For Effective Shock-Capturing, M. Harris, Eduardo Divo, Alain J. Kassab

Publications

In recent years significant focus has been given to the study of Radial basis functions (RBF), especially in their use on solving partial differential equations (PDE). RBF have an impressive capability of inter- polating scattered data, even when this data presents localized discontinuities. However, for infinitely smooth RBF such as the Multiquadrics, inverse Multiquadrics, and Gaussian, the shape parameter must be chosen properly to obtain accurate approximations while avoiding ill-conditioning of the interpolating matrices. The optimum shape parameter can vary significantly depending on the field, particularly in locations of steep gradients, shocks, or discontinuities. Typically, the shape parameter is chosen …

A Coupled Localized Rbf Meshless/Drbem Formulation For Accurate Modeling Of Incompressible Fluid Flows, Leonardo Bueno, Eduardo Divo, Alain J. Kassab

Publications

Velocity-pressure coupling schemes for the solution of incompressible fluid flow problems in Computational Fluid Dynamics (CFD) rely on the formulation of Poisson-like equations through projection methods. The solution of these Poisson-like equations represent the pressure correction and the velocity correction to ensure proper satisfaction of the conservation of mass equation at each step of a time-marching scheme or at each level of an iteration process. Inaccurate solutions of these Poisson-like equations result in meaningless instantaneous or intermediate approximations that do not represent the proper time-accurate behavior of the flow. The fact that these equations must be solved to convergence at …

Model For Computing Kinetics Of The Graphene Edge Epitaxial Growth On Copper, Mikhail Khenner

Mathematics Faculty Publications

A basic kinetic model that incorporates a coupled dynamics of the carbon atoms and dimers ona copper surface is used to compute growth of a single-layer graphene island. The speed of theisland's edge advancement on Cu[111] and Cu[100] surfaces is computed as a function of the growthtemperature and pressure. Spatially resolved concentration pro les of the atoms and dimers aredetermined, and the contributions provided by these species to the growth speed are discussed.Island growth in the conditions of a thermal cycling is studied.

Theoretical Modeling Of The Effect Of Polymer Chain Immobilization Rates On Holographic Recording In Photopolymers, Dana Mackey, Paul O'Reilly, Izabela Naydenova

Articles

This paper introduces an improved mathematical model for holographic grating formation in an acrylamide-based photopolymer, which consists of partial differential equations derived from physical laws. The model is based on the two-way diffusion theory of \cite{izabela}, which assumes short polymer chains are free to diffuse, and generalizes a similar model presented in \cite{josab} by introducing an immobilization rate governed by chain growth and cross-linking. Numerical simulations were carried out in order to investigate the behaviour of the photopolymer system for short and long exposures, with particular emphasis on the effect of recording parameters (such as illumination frequency and intensity), as …

The Dynamics Of Flat Surface Internal Geophysical Waves With Currents, Alan Compelli, Rossen Ivanov

Articles

A two-dimensional water wave system is examined consisting of two discrete incompressible fluid domains separated by a free common interface. In a geophysical context this is a model of an internal wave, formed at a pycnocline or thermocline in the ocean. The system is considered as being bounded at the bottom and top by a flatbed and wave-free surface respectively. A current profile with depth-dependent currents in each domain is considered. The Hamiltonian of the system is determined and expressed in terms of canonical wave-related variables. Limiting behavior is examined and compared to that of other known models. The linearised …

Rainich-Type Conditions For Perfect Fluid Spacetimes, Dionisios Krongos, Charles G. Torre

Research Vignettes

In this worksheet we describe and illustrate a relatively simple set of new Rainich-type conditions on an n-dimensional spacetime which are necessary and sufficient for it to define a perfect fluid solution of the Einstein field equations. Procedures are provided which implement these Rainich-type conditions and which reconstruct the perfect fluid from the metric. These results provide an example of the idea of geometrization of matter fields in general relativity, which is a purely geometrical characterization of matter fields via the Einstein field equations.

A Soft Condensed Matter Approach Towards Mathematical Modelling Of Mass Transport And Swelling In Food Grains, Michael Chapwanya, N. Misra

Articles

Soft condensed matter (SCM) physics has recently gained importance for a large class of engineering materials. The treatment of food materials from a soft matter perspective, however, is only at the surface and is gaining importance for understanding the complex phenomena and structure of foods. In this work, we present a theoretical treatment of navy beans from a SCM perspective to describe the hydration kinetics. We solve the transport equations within a porous matrix and employ the Flory–Huggin’s equation for polymer–solvent mixture to balance the osmotic pressure. The swelling of the legume seed is modelled as a moving boundary with …

One-Dimensional Weakly Nonlinear Model Equations For Rossby Waves, David Henry, Rossen Ivanov

Articles

In this study we explore several possibilities for modelling weakly nonlinear Rossby waves in fluid of constant depth, which propagate predominantly in one direction. The model equations obtained include the BBM equation, as well as the integrable KdV and Degasperis-Procesi equations.

Analysis Of A Partial Differential Equation Model Of Surface Electromigration, Selahittin Cinar

Masters Theses & Specialist Projects

A Partial Differential Equation (PDE) based model combining surface electromigration and wetting is developed for the analysis of the morphological instability of mono-crystalline metal films in a high temperature environment typical to operational conditions of microelectronic interconnects. The atomic mobility and surface energy of such films are anisotropic, and the model accounts for these material properties. The goal of modeling is to describe and understand the time-evolution of the shape of film surface. I will present the formulation of a nonlinear parabolic PDE problem for the height function h(x,t) of the film in the horizontal …

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. …

How To Find Killing Vectors, Charles G. Torre

How to... in 10 minutes or less

We show how to compute the Lie algebra of Killing vector fields of a metric in Maple using the commands KillingVectors and LieAlgebraData. A Maple worksheet and a PDF version can be found below.

Modelling Two-Dimensional Photopolymer Patterns Produced With Multiple-Beam Holography, Dana Mackey, Tsvetanka Babeva, Izabela Naydenova, Vincent Toal

Conference papers

Periodic structures referred to as photonic crystals attract considerable interest due to their potential applications in areas such as nanotechnology, photonics, plasmonics, etc. Among various techniques used for their fabrication, multiple-beam holography is a promising method enabling defect-free structures to be produced in a single step over large areas.

In this paper we use a mathematical model describing photopolymerisation to simulate two-dimensional structures produced by the interference pattern of three noncoplanar beams. The holographic recording of different lattices is studied by variation of certain parameters such as beam wave vectors, time and intensity of illumination.

Stiefel And Grassmann Manifolds In Quantum Chemistry, Eduardo Chiumiento, Michael Melgaard

Articles

We establish geometric properties of Stiefel and Grassmann manifolds which arise in relation to Slatertype variational spaces in many-particle Hartree-Fock theory and beyond. In particular, we prove thatthey are analytic homogeneous spaces and submanifolds of the space of bounded operators on the single-particle Hilbert space. As a by-product we obtain that they are complete Finsler manifolds. These geometric properties underpin state-of-the-art results on existence of solutions to Hartree-Fock type equations.

Heterogeneous Multiscale Modeling Of Advection-Diffusion Problems, David J. Gardner, Daniel R. Reynolds

Mathematics Research

No abstract provided.

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.