Numerical Analysis and Computation Commons^™

Laminar Boundary Layer Flow Of Sisko Fluid, Manisha Patel, Jayshri Patel, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

The problem of steady two dimensional laminar boundary layer flow of non-Newtonian fluid is analyzed in the present paper. Sisko fluid model, one of the various fluid models of non- Newtonian fluid, is considered for stress-strain relationship. Similarity and numerical solutions obtained for the defined flow problem.

Numerical Solutions Of Generalized Burgers' Equations For Some Incompressible Non-Newtonian Fluids, Yupeng Shu

LSU New Orleans Theses and Dissertations

The author presents some generalized Burgers' equations for incompressible and isothermal flow of viscous non-Newtonian fluids based on the Cross model, the Carreau model, and the Power-Law model and some simple assumptions on the flows. The author numerically solves the traveling wave equations for the Cross model, the Carreau model, the Power-Law model by using industrial data. The author proves existence and uniqueness of solutions to the traveling wave equations of each of the three models. The author also provides numerical estimates of the shock thickness as well as maximum strain $\varepsilon_{11}$ for each of the fluids.

Global Optimized Isothermal And Nonlinear Models Of Earth’S Standard Atmosphere, Nihad E. Daidzic, Ph.D.,

International Journal of Aviation, Aeronautics, and Aerospace

Both, a global isothermal temperature model and a nonlinear quadratic temperature model of the ISA was developed and presented here. Constrained optimization techniques in conjunction with the least-square-root approximations were used to design best-fit isothermal models for ISA pressure and density changes up to 47 geopotential km for NLPAM, and 86 orthometric km for ISOAM respectively. The mass of the dry atmosphere and the relevant fractional-mass scale heights have been computed utilizing the very accurate eight-point Gauss-Legendre numerical quadrature for both ISOAM and NLPAM. Both, the ISOAM and the NLPAM represent viable alternatives to ISA in many practical applications and …

Implicit-Explicit Higher-Order Time Integration Schemes For Computations Of Structural Dynamics With Fluid-Structure Interaction, José C. Pedro, Mapundi K. Banda, Precious Sibanda

Applications and Applied Mathematics: An International Journal (AAM)

In this paper higher order implicit Runge-Kutta schemes are applied to fluid-structure interaction (FSI) simulations. A staggered approach with a structural predictor is applied to an FSI problem. The equations governing the dynamics of the structure are integrated in time by the Explicit Single Diagonal Implicit Runge-Kutta (ESDIRK) schemes and the arbitrary high order finite volume scheme is taken as the fluid solver. The performance of the ESDIRK scheme of order of convergence three to five is tested. Comparative studies with other time integration schemes which have been successfully applied to FSI problems are undertaken. Comparisons to test the performance …

Transition Orbits Of Walking Droplets, Joshua Parker

Physics

It was recently discovered that millimeter-sized droplets bouncing on the surface of an oscillating bath of the same fluid can couple with the surface waves it produces and begin walking across the fluid bath. These walkers have been shown to behave similarly to quantum particles; a few examples include single-particle diffraction, tunneling, and quantized orbits. Such behavior occurs because the drop and surface waves depend on each other to exist, making this the first and only known macroscopic pilot-wave system. In this paper, the quantized orbits between two identical drops are explored. By sending a perturbation to a pair of …

Efficient General Computational Method For Estimation Of Standard Atmosphere Parameters, Nihad E. Daidzic Ph.D., Sc.D.

International Journal of Aviation, Aeronautics, and Aerospace

Knowledge of standard air temperature, pressure, density, speed of sound, and viscosity as a function of altitude is essential information in aircraft design, performance testing, pressure altimeter calibration, and several other aeronautical engineering and aviation science applications. A new efficient computational method for rapid calculations of standard atmospheric parameters up to 86 orthometric km is presented. Additionally, mass and weight of each standard atmospheric layer were calculated using a numerical integration method. The sum of all fractional masses and weights represents the total mass and weight of Earth’s atmosphere. The results obtained here agree well with measurements and models of …

A Contribution Toward Better Understanding Of Overbanking Tendency In Fixed-Wing Aircraft, Nihad E. Daidzic

Journal of Aviation Technology and Engineering

The phenomenon of overbanking tendency for a rigid-body, fixed-wing aircraft is investigated. Overbanking tendency is defined as a spontaneous, unbalanced rolling moment that keeps increasing an airplane’s bank angle in steep turns and must be arrested by opposite aileron action. As stated by the Federal Aviation Administration, the overbanking tendency may lead to a loss of control, especially in instrument meteorological conditions. It was found in this study that the speed differential over wing halves in horizontal turns indeed creates a rolling moment that achieves maximum values for bank angles between 45 and 55 degrees. However, this induced rolling moment …

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 …

Low Mach Number Fluctuating Hydrodynamics Of Diffusively Mixing Fluids, Aleksandar Donev, Andy J. Nonaka, Yifei Sun, Thomas Fai, Alejandro Garcia, John B. Bell

Faculty Publications

We formulate low Mach number fluctuating hydrodynamic equations appropriate for modeling diffusive mixing in isothermal mixtures of fluids with different density and transport coefficients. These equations eliminate the fast isentropic fluctuations in pressure associated with the propagation of sound waves by replacing the equation of state with a local thermodynamic constraint. We demonstrate that the low Mach number model preserves the spatio-temporal spectrum of the slower diffusive fluctuations. We develop a strictly conservative finite-volume spatial discretization of the low Mach number fluctuating equations in both two and three dimensions. We construct several explicit Runge-Kutta temporal integrators that strictly maintain the …

Evolution Of Perturbations In Flow Field Mechanics, Samantha R. Bell, David Forliti, Nils Sedano, Kriss Vanderhyde

STAR Program Research Presentations

This project explores the stability analysis of a given flow field. Specifically, where the peak disturbance occurs in a flow as this is the disturbance that is most likely to occur. In rocket combustion, it is important to understand where the maximum disturbance occurs so that the mixing of fuel can be stabilized. The instabilities are the results of frequencies in the area surrounding the flow field. The linear stability governing equations are employed to better understand the disturbance. The governing equations for continuity and momentum in the x and y directions are used to form an equation for the …

Energy-Driven Pattern Formation In Planar Dipole-Dipole Systems, Jaron P. Kent-Dobias

HMC Senior Theses

A variety of two-dimensional fluid systems, known as dipole-mediated systems, exhibit a dipole-dipole interaction between their fluid constituents. The com- petition of this repulsive dipolar force with the cohesive fluid forces cause these systems to form intricate and patterned structures in their boundaries. In this thesis, we show that the microscopic details of any such system are irrelevant in the macroscopic limit and contribute only to a constant offset in the system’s energy. A numeric model is developed, and some important stable domain morphologies are characterized. Previously unresolved bifurcating branches are explored. Finally, by applying a random energy background to …

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

Relative Equilibria Of Isosceles Triatomic Molecules In Classical Approximation, Damaris Miriam Mckinley

Theses and Dissertations (Comprehensive)

In this thesis we study relative equilibria of di-atomic and isosceles tri-atomic molecules in classical approximations with repulsive-attractive interaction. For di-atomic systems we retrieve well-known results. The main contribution consists of the study of the existence and stability of relative equilibria in a three-atom system formed by two identical atoms of mass $m$ and a third of mass $m_3$, constrained in an isosceles configuration at all times.

Given the shape of the binary potential only, we discuss the existence of equilibria and relative equilibria. We represent the results in the form of energy-momentum diagrams. We find that fixing the masses …

Characterization Of The Drilling Via The Vibration Augmenter Of Rotary-Drills And Sound Signal Processing Of Impacted Pipe As A Potential Water Height Assessment Tool, Nicholas Morris

STAR Program Research Presentations

The focus of the internship has been on two topics: a) Characterize the drilling performance of a novel percussive augmenter – this drill was developed by the JPL’s Advanced Technologies Group and its performance was characterized; and b) Examine the feasibility of striking a pipe as a means of assessing the water height inside the pipe. The purpose of this investigation is to examine the possibility of using a simple method of applying impacts to a pipe wall and determining the water height from the sonic characteristic differences including damping, resonance frequencies, etc. Due to multiple variables that are relevant …

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 …

Preconditioning Visco-Resistive Mhd For Tokamak Plasmas, Daniel R. Reynolds, Ravi Samtaney, Hilari C. Tiedeman

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.

Numerical Methods For Fluid-Structure Interaction - A Review, Gene Hou, Jin Wang, Anita Layton

Mechanical & Aerospace Engineering Faculty Publications

The interactions between incompressible fluid flows and immersed structures are nonlinear multi-physics phenomena that have applications to a wide range of scientific and engineering disciplines. In this article, we review representative numerical-methods based on conforming and non-conforming meshes that are currently available for computing fluid-structure interaction problems, with an emphasis on some of the recent developments in the field. A goal is to categorize the selected methods and assess their accuracy and efficiency. We discuss challenges faced by researchers in this field, and we emphasize the importance of interdisciplinary effort for advancing the study in fluid-structure interactions

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.

Paleomagnetism And Investigation Of 40 Ma Lavas, Liverpool Range, New South Whales, Australia, Nathan M. Padilla

Physics

The main focus of this project is the continued study of a reversal of the earth’s magnetic field recorded from lavas in the Liverpool Range of New South Whales, Australia. This reverse-to-normal transition, recently dated at ~40 Ma, was first reported in Nature in 1986. [2] In March 2011 some 200+ cores were drilled from several sections about the volcanic range—Jemmy’s Creek, Bald Hill, Rock Creek, Yarraman, and Coolah Tops Road. Here we focus on paleomagnetic findings from samples drilled from the most extensive section, that being along the trail near Jemmy’s Creek. Results from alternating field demagnetization show the …

Quantum Mechanical Model Of The Hydrogen Atom, Andres Ortiz^, Snehanshu Saha*

COURI Symposium Abstracts, Spring 2011

The quantum mechanical wave function of an object under a certain potential is determined by the Schrödinger equation. It is known that the latter equation can be solved analytically only for a finite number of distinct potentials, being a Coulomb potential one of them. For an electron in a hydrogen atom, we analytically find the angular parts of its wave function by the usual method. Under the appropriate change of variable, the radial part of the equation is turned into a Legendre differential equation whose solutions are Legendre polynomials. Using finite difference method and shooting method taking into account physically …

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.

Finite Element Analysis In Porous Media For Incompressible Flow Of Contamination From Nuclear Waste, Abbas Al-Bayati, Saad A. Manaa, Ekhlass S. Ahmed

Applications and Applied Mathematics: An International Journal (AAM)

A non-linear parabolic system is used to describe incompressible nuclear waste disposal contamination in porous media, in which both molecular diffusion and dispersion are considered. The Galerkin method is applied for the pressure equation. For the brine, radionuclide and heat, a kind of partial upwind finite element scheme is constructed. Examples are included to demonstrate certain aspects of the theory and illustrate the capabilities of the kind of partial upwind finite element approach.

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