Applied Mathematics Commons^™

Modelling The Flow Of Aqueous Humor In Schlemm’S Canal In The Eye, Ram Avtar, Swati Srivastava, Rashmi Srivastava

Applications and Applied Mathematics: An International Journal (AAM)

A simple mathematical model for the transient flow of aqueous humor in the canal of Schlemm is developed to investigate the acceleration effects of a sudden elevation in the intraocular pressure on the flow characteristics of the aqueous humor in the canal. The model treats a canal segment as a tube of elliptic cross-section. Exact analytical solution to the model is obtained using separation of variables method. The effects of some important model parameters on the maximum and minimum shear stresses exerted on the Schlemm’s canal epithelial cells (wall) by flowing aqueous humor are investigated for the steady-state flow.

Smarticles: A Method For Identifying And Correcting Instability And Error Caused By Explicit Integration Techniques In Physically Based Simulations, Susan Aileen Marano

Master's Theses

Using an explicit integration method in physically based animations has many advantages including conceptual and computational simplicity, however, it re- quires small time steps to ensure low numerical instability. Simulations with large numbers of individually interacting components such as cloth, hair, and fluid models, are limited by the sections of particles most susceptible to error. This results in the need for smaller time steps than required for the majority of the system. These sections can be diverse and dynamic, quickly changing in size and location based on forces in the system. Identifying and handling these trou- blesome sections could allow …

Viscosity Dependence Of Faraday Wave Formation Thresholds, Lisa Michelle Slaughter

Physics

This experiment uses an electromagnetic shaker to produce standing wave patterns on the surface of a vertically oscillating sample of silicon liquid. These surface waves, known as Faraday waves, form shapes such as squares, lines, and hexagons. They are known to be dependent upon the frequency and amplitude of the forcing as well as on the viscosity and depth of the liquid in the dish. At a depth of 4mm and for various silicon liquids having kinematic viscosities of 10, 20, and 38 cSt, we determined the acceleration at which patterns form for frequencies between 10 and 60 Hz. For …

Dispersion Of A Solute In Hartmann Two-Fluid Flow Between Two Parallel Plates, J. P. Kumar, J. C. Umavathi

Applications and Applied Mathematics: An International Journal (AAM)

The paper presents an analytical solution for the dispersion of a solute in a conducting immiscible fluid flowing between two parallel plates in the presence of a transverse magnetic field. The fluids in both the regions are incompressible, electrically conducting and the transport properties are assumed to be constant. The channel walls are assumed to be electrically insulating. Separate solutions for each fluid are obtained and these solutions are matched at the interface using suitable matching conditions. The results are tabulated for various values of viscosity ratio, pressure gradient and Hartman number on the effective Taylor dispersion coefficient and volumetric …

Hydro-Thermal Convective Solutions For An Aquifer System Heated From Below, Dambaru Bhatta

Applications and Applied Mathematics: An International Journal (AAM)

We investigate the effect of hydro-thermal convection in an aquifer system. It is assumed that the aquifer is bounded below and above by impermeable boundaries and it is heated from below. The solution of the governing system is expressed in terms of the basic steady state solution and perturbed solution. We obtain the critical Rayleigh number and critical wavenumber using Runge-Kutta method in combination of shooting method and present the marginal stability curve. The amplitude equation is derived by introducing the adjoint system. After amplitude is obtained, we compute the linear solutions for super-critical and sub-critical cases. Numerical results for …

Simulations Of Surfactant Driven Thin Film Flow, Shreyas Kumar

HMC Senior Theses

This thesis is intended to fulfill the requirements of the Math and Physics departments at Harvey Mudd College. We begin with a brief introduction to the study of surfactant dynamics followed by some background on the experimental framework our work is related to. We then go through a derivation of the model we use, and explore in depth the nature of the Equation of State (EoS), the relationship between the surface tension on a fluid and the surfactant concentration. We consider the effect of using an empirical equation of state on the results of the simulations and compare the new …

Computational Fluid Dynamics (Cfd) Modeling Of A Laboratory Scale Coal Gasifier, Kiel S. Schultheiss

College of Graduate Studies: Theses & Dissertations

Furthering gasification technology is an essential part of advancing clean coal technologies. In order to seek insight into the appropriate operations for the formation of synthetic gas (syngas) a numerical simulation was performed to predict the phenomena of coal gasification in a laboratory scale entrained-flow coal gasifier. The mesh for the model was developed with ICEM CFD software and the chemical and physical phenomena were modeled using the fluid flow solver ANSYS FLUENT. Mesh independence was verified. The model was validated with experimental data from several studies performed on a laboratory scale gasifier.

Systematic examination of the model was performed …

Mhd Mixed Convective Flow Of Viscoelastic And Viscous Fluids In A Vertical Porous Channel, R. Sivaraj, B. R. Kumar, J. Prakash

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we analyze the problem of steady, mixed convective, laminar flow of two incompressible, electrically conducting and heat absorbing immiscible fluids in a vertical porous channel filled with viscoelastic fluid in one region and viscous fluid in the other region. A uniform magnetic field is applied in the transverse direction, the fluids rise in the channel driven by thermal buoyancy forces associated with thermal radiation. The equations are modeled using the fully developed flow conditions. An exact solution is obtained for the velocity, temperature, skin friction and Nusselt number distributions. The physical interpretation to these expressions is examined …

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.

Simulations Of Surfactant Spreading, Jeffrey Wong

HMC Senior Theses

Thin liquid films driven by surface tension gradients are studied in diverse applications, including the spreading of a droplet and fluid flow in the lung. The nonlinear partial differential equations that govern thin films are difficult to solve analytically, and must be approached through numerical simulations. We describe the development of a numerical solver designed to solve a variety of thin film problems in two dimensions. Validation of the solver includes grid refinement studies and comparison to previous results for thin film problems. In addition, we apply the solver to a model of surfactant spreading and make comparisons with theoretical …

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

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.

Bound For The Complex Growth Rate In Thermosolutal Convection Coupled With Cross-Diffusions, Hari Mohan Sharma

Applications and Applied Mathematics: An International Journal (AAM)

Thermosolutal convection problem of the Veronis’ type coupled with cross–diffusion is considered in the present paper. A semi -circle theorem that prescribes upper limit for the complex growth rate of oscillatory motions of neutral or growing amplitude in such a manner that it naturally culminates in sufficient conditions precluding the non- existence of such motions is derived. Further, results for thermosolutal convection problems with or without the individual consideration of Dufour and Soret effects follow as a consequence.

Effects Of Radiation And Variable Viscosity On Mhd Free Convective Flow And Mass Transfer Over A Stretching Sheet With Chemical Reaction, M. A. Seddeek, A. A. Almushigeh

Applications and Applied Mathematics: An International Journal (AAM)

A similarity solution is proposed for the analysis of steady free convection heat and mass transfer over a stretching sheet. The effect of radiation, chemical reaction and variable viscosity on hydromagnetic heat and mass transfer in the presence of magnetic field are investigated. The governing partial differential equations are transformed to the ordinary differential equations using similarity variables, and then solved numerically by means of the fourth-order Runge– Kutta method with shooting technique. A comparison with exact solution is performed and the results are found to be in excellent. Numerical results for the velocity, temperature and concentration as well as …

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.

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 …

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, Agegnehu Atena, Mikhail Khenner

Mathematics Faculty Publications

In this paper the lubrication-type dynamical model is developed of a molten, pulsed laser-irradiated metallic film. The heat transfer problem that incorporates the absorbed heat from a single beam or interfering 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 peak laser beam intensity, the film optical thickness, the Biot and …

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

Mathematics Faculty Publications

In this paper the lubrication-type dynamical model is developed of a molten, pulsed laser-irradiated metallic film. The heat transfer problem that incorporates the absorbed heat from a single beam or interfering 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 peak laser beam intensity, the film optical thickness, the Biot and …

Stability Of Traveling Waves In Thin Liquid Films Driven By Gravity And Surfactant, Ellen Peterson, Michael Shearer, Thomas P. Witelski, Rachel Levy

All HMC Faculty Publications and Research

A thin layer of fluid flowing down a solid planar surface has a free surface height described by a nonlinear PDE derived via the lubrication approximation from the Navier Stokes equations. For thin films, surface tension plays an important role both in providing a significant driving force and in smoothing the free surface. Surfactant molecules on the free surface tend to reduce surface tension, setting up gradients that modify the shape of the free surface. In earlier work [12, 13J a traveling wave was found in which the free surface undergoes three sharp transitions, or internal layers, and the surfactant …

Converging Flow Between Coaxial Cones, O. Hall, A. D. Gilbert, C. P. Hills

Articles

Fluid flow governed by the Navier-Stokes equation is considered in a domain bounded by two cones with the same axis. In the first, 'non-parallel' case, the two cones have the same apex and different angles θ = α and β in spherical polar coordinates (r, θ, φ). In the second, 'parallel' case, the two cones have the same opening angle α, parallel walls separated by a gap h and apices separated by a distance h/sinα. Flows are driven by a source Q at the origin, the apex of the lower cone in the parallel case. The Stokes solution for the …

Nonaxisymmetric Stokes Flow Between Concentric Cones, O. Hall, C. P. Hills, A. D. Gilbert

Articles

We study the fully three-dimensional Stokes flow within a geometry consisting of two infinite cones with coincident apices. The Stokes approximation is valid near the apex and we consider the dominant flow features as it is approached. The cones are assumed to be stationary and the flow to be driven by an arbitrary far-field disturbance. We express the flow quantities in terms of eigenfunction expansions and allow for the first time for nonaxisymmetric flow regimes through an azimuthal wave number. The eigenvalue problem is solved numerically for successive wave numbers. Both real and complex sequences of eigenvalues are found, their …

Singular Superposition/Boundary Element Method For Reconstruction Of Multi-Dimensional Heat Flux Distributions With Application To Film Cooling Holes, Mahmood Silieti, Eduardo Divo, Alain J. Kassab

Publications

A hybrid singularity superposition/boundary element-based inverse problem method for the reconstruction of multi-dimensional heat flux distributions is developed. Cauchy conditions are imposed at exposed surfaces that are readily reached for measurements while convective boundary conditions are unknown at surfaces that are not amenable to measurements such as the walls of the cooling holes. The purpose of the inverse analysis is to determine the heat flux distribution along cooling hole surfaces. This is accomplished in an iterative process by distributing a set of singularities (sinks) inside the physical boundaries of the cooling hole (usually along cooling hole centerline) with a given …

Peristaltic Pumping Of A Non-Newtonian Fluid, Amit Medhavi

Applications and Applied Mathematics: An International Journal (AAM)

The flow induced by sinusoidal peristaltic motion of the tube wall of a non-Newtonian fluid obeying Herschel-Bulkley equation (a general rheological equation that represents a powerlaw, Bingham and Newtonian fluid for particular choice of parameters) under long wavelength and low Reynolds number approximation is investigated. The results obtained for flow rate, pressure drop and friction force are discussed both qualitatively and quantitatively and compared with other related studies. It is found that the pressure drop increases with the flow rate and yield stress but decreases with the increasing amplitude ratio. The flow behaviour index shows significant impact on the magnitude …