Numerical Analysis and Computation Commons^™

Shear Driven Micro-Fluidic Pump For Cardiovascular Applications, Nihad E. Daidzic

Aviation Department Publications

A valveless shear-driven micro-fluidic pump design (SDMFP) for hemodynamic applications is presented in this work. One of the possible medical and biomedical applications is in-vivo hemodynamic (human blood circulation) support/assist. One or more SDMFPs can be inserted/implanted into vascular lumens in a form of a stent/duct in series and/or in parallel (bypass duct) to support blood circulation in-vivo. A comprehensive review of various micro-pump designs up to about mid 2000’s is given in [1,2]. Many of micropump designs considered are not suitable for in-vivo or even in-vitro medical/biomedical applications.

Operating principles, design, and SDMFP features are given in [3]. A …

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 …

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 …

Application Of An Rbf Blending Interpolation Method To Problems With Shocks, Michael Harris, Eduardo Divo, Alain J. Kassab

Publications

Radial basis functions (RBF) have become an area of research in recent years, especially in the use of solving partial differential equations (PDE). Radial basis functions have an impressive capability in interpolating scattered data, even for data with discontinuities. Although, for infinitely smooth radial basis functions such as the multi-quadrics and inverse multi-quadrics, 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 depending on the field, such as in locations of sharp gradients or shocks. Typically, the shape parameter is chosen to maintain a high …

Kinetic Monte Carlo Methods For Computing First Capture Time Distributions In Models Of Diffusive Absorption, Daniel Schmidt

HMC Senior Theses

In this paper, we consider the capture dynamics of a particle undergoing a random walk above a sheet of absorbing traps. In particular, we seek to characterize the distribution in time from when the particle is released to when it is absorbed. This problem is motivated by the study of lymphocytes in the human blood stream; for a particle near the surface of a lymphocyte, how long will it take for the particle to be captured? We model this problem as a diffusive process with a mixture of reflecting and absorbing boundary conditions. The model is analyzed from two approaches. …

A Numerical Study Of Construction Of Honey Bee Comb, Pamela Guerrero, Pamela C. Guerrero

Murray State Theses and Dissertations

We use finite difference methods in the treatment of an existing system of partial differential equations that captures the dynamics of parallel honeycomb construction in a bee hive. We conduct an uncertainty analysis by calculating the partial rank correlation coefficient for the parameters to find which are most important to the outcomes of the model. We then use an eFAST method to determine both the individual and total sensitivity index for the parameters. Afterwards we examine our numerical model under varying initial conditions and parameter values, and compare ratios found from local data with the golden mean by fitting images …

A Physics-Based Approach To Modeling Wildland Fire Spread Through Porous Fuel Beds, Tingting Tang

Theses and Dissertations--Mechanical and Aerospace Engineering

Wildfires are becoming increasingly erratic nowadays at least in part because of climate change. CFD (computational fluid dynamics)-based models with the potential of simulating extreme behaviors are gaining increasing attention as a means to predict such behavior in order to aid firefighting efforts. This dissertation describes a wildfire model based on the current understanding of wildfire physics. The model includes physics of turbulence, inhomogeneous porous fuel beds, heat release, ignition, and firebrands. A discrete dynamical system for flow in porous media is derived and incorporated into the subgrid-scale model for synthetic-velocity large-eddy simulation (LES), and a general porosity-permeability model is …

Multi-Scale Cardiovascular Flow Analysis By An Integrated Meshless-Lumped Parameter Model, Leonardo A. Bueno, Eduardo A. Divo, Alain J. Kassab

Publications

A computational tool that integrates a Radial basis function (RBF)-based Meshless solver with a Lumped Parameter model (LPM) is developed to analyze the multi-scale and multi-physics interaction between the cardiovascular flow hemodynamics, the cardiac function, and the peripheral circulation. The Meshless solver is based on localized RBF collocations at scattered data points which allows for automation of the model generation via CAD integration. The time-accurate incompressible flow hemodynamics are addressed via a pressure-velocity correction scheme where the ensuing Poisson equations are accurately and efficiently solved at each time step by a Dual-Reciprocity Boundary Element method (DRBEM) formulation that takes advantage …

Direct Sampling Methods For Inverse Scattering Problems, Ala Mahmood Nahar Al Zaalig

Dissertations, Master's Theses and Master's Reports

Recently, direct sampling methods became popular for solving inverse scattering problems to estimate the shape of the scattering object. They provide a simple tool to directly reconstruct the shape of the unknown scatterer. These methods are based on choosing an appropriate indicator function f on R^d, d=2 or 3, such that f(z) decides whether z lies inside or outside the scatterer. Consequently, we can determine the location and the shape of the unknown scatterer.

In this thesis, we first present some sampling methods for shape reconstruction in inverse scattering problems. These methods, which are described in Chapter 1, …

Computational Fluid Dynamics In A Terminal Alveolated Bronchiole Duct With Expanding Walls: Proof-Of-Concept In Openfoam, Jeremy Myers

Theses and Dissertations

Mathematical Biology has found recent success applying Computational Fluid Dynamics (CFD) to model airflow in the human lung. Detailed modeling of flow patterns in the alveoli, where the oxygen-carbon dioxide gas exchange occurs, has provided data that is useful in treating illnesses and designing drug-delivery systems. Unfortunately, many CFD software packages have high licensing fees that are out of reach for independent researchers. This thesis uses three open-source software packages, Gmsh, OpenFOAM, and ParaView, to design a mesh, create a simulation, and visualize the results of an idealized terminal alveolar sac model. This model successfully demonstrates that OpenFOAM can be …

A Numerical Scheme For Generalized Fractional Optimal Control Problems, N. Singha, C. Nahak

Applications and Applied Mathematics: An International Journal (AAM)

This paper introduces a generalization of the Fractional Optimal Control Problem (GFOCP). Proposed generalizations differ in terms of explaining the constraint involved in the dynamical system of the control problem. We assume the constraint as an arbitrary function of fractional derivatives and fractional integrals. By this assumption the restriction on constraint, to be of some prescribed function of fractional operators, is removed. Deduction of necessary optimality conditions followed by particular cases and examples has been provided. Additionally, we construct a solution scheme for the suggested class of (GFOCP)’s. The formulation of this scheme is done by implementing the Adomian decomposition …

Fast Method Of Particular Solutions For Solving Partial Differential Equations, Anup Raja Lamichhane

Dissertations

Method of particular solutions (MPS) has been implemented in many science and engineering problems but obtaining the closed-form particular solutions, the selection of the good shape parameter for various radial basis functions (RBFs) and simulation of the large-scale problems are some of the challenges which need to overcome. In this dissertation, we have used several techniques to overcome such challenges.

The closed-form particular solutions for the Matérn and Gaussian RBFs were not known yet. With the help of the symbolic computational tools, we have derived the closed-form particular solutions of the Matérn and Gaussian RBFs for the Laplace and biharmonic …

How Steep Is Steep? Learning Curves In Training Undergraduates To Do Fluid-Structure Interaction Modeling, Nicholas A. Battista

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Krylov Subspace Spectral Method With Multigrid For A Time-Dependent, Variable-Coefficient Partial Differential Equation, Haley Renee Dozier

Master's Theses

Krylov Subspace Spectral (KSS) methods are traditionally used to solve time-dependent, variable-coefficient PDEs. They are high-order accurate, component-wise methods that are efficient with variable input sizes.

This thesis will demonstrate how one can make KSS methods even more efficient by using a Multigrid-like approach for low-frequency components. The essential ingredients of Multigrid, such as restriction, residual correction, and prolongation, are adapted to the timedependent case. Then a comparison of KSS, KSS with Multigrid, KSS-EPI and standard Krylov projection methods will be demonstrated.

An Algorithm For The Machine Calculation Of Minimal Paths, Robert Whitinger

Electronic Theses and Dissertations

Problems involving the minimization of functionals date back to antiquity. The mathematics of the calculus of variations has provided a framework for the analytical solution of a limited class of such problems. This paper describes a numerical approximation technique for obtaining machine solutions to minimal path problems. It is shown that this technique is applicable not only to the common case of finding geodesics on parameterized surfaces in R³, but also to the general case of finding minimal functionals on hypersurfaces in Rⁿ associated with an arbitrary metric.

A Robust Uniform B-Spline Collocation Method For Solving The Generalized Phi-Four Equation, W. K. Zahra, W. A. Ouf, M. S. El-Azab

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we develop a numerical solution based on cubic B-spline collocation method. By applying Von-Neumann stability analysis, the proposed technique is shown to be unconditionally stable. The accuracy of the presented method is demonstrated by a test problem. The numerical results are found to be in good agreement with the exact solution.

Explicitly Separating Growth And Motility In A Glioblastoma Tumor Model, Tracy Stepien, Erica Rutter, Meng Fan, Yang Kuang

Biology and Medicine Through Mathematics Conference

No abstract provided.

Thermoelastic Waves In Microstructured Solids, Arkadi Berezovski, Mihhail Berezovski

Publications

Thermoelastic wave propagation suggests a coupling between elastic deformation and heat conduction in a body. Microstructure of the body influences the both processes. Since energy is conserved in elastic deformation and heat conduction is always dissipative, the generalization of classical elasticity theory and classical heat conduction is performed differently. It is shown in the paper that a hyperbolic evolution equation for microtemperature can be obtained in the framework of the dual internal variables approach keeping the parabolic equation for the macrotemperature. The microtemperature is considered as a macrotemperature fluctuation. Numerical simulations demonstrate the formation and propagation of thermoelastic waves in …

Stereographic Visualization Of Bose-Einstein Condensate Clouds To Measure The Gravitational Constant, Ed Wesley Wells

College of Graduate Studies: Theses & Dissertations

This thesis describes a set of tools that can be used for the rapid design of atom interferometer schemes suitable for measuring Newton's Universal Gravitation constant also known as "Big G". This tool set is especially applicable to Bose--Einstein--condensed systems present in NASA's Cold Atom Laboratory experiment to be deployed to the International Space Station in 2017. These tools include a method of approximating the solutions of the nonlinear Schrödinger or Gross--Pitaevskii equation (GPE) using the Lagrangian Variational Method. They also include a set of software tools for translating the approximate solutions of the GPE into images of the optical …

A High Accuracy Minimally Invasive Regularization Technique For Navier-Stokes Equations At High Reynolds Number, Mustafa Aggul

Dissertations, Master's Theses and Master's Reports

A method is presented, that combines the defect and deferred correction approaches to approximate solutions of Navier-Stokes equations at high Reynolds number. The method is of high accuracy in both space and time, and it allows for the usage of legacy codes (a frequent requirement in the simulation of turbulent flows in complex geometries). The two-step method is considered here; in order to obtain a regularization that is second order accurate in space and time, the method computes a low-order accurate, stable and computationally inexpensive approximation (Backward Euler with artificial viscosity) twice. The results are readily extendable to the higher …

Defect-Deferred Correction Method For The Two-Domain Convection-Dominated Convection-Diffusion Problem, Dilek Erkmen

Dissertations, Master's Theses and Master's Reports

We present a method for solving a fluid-fluid interaction problem (two convection-dominated convection-diusion problems adjoined by an interface), which is a simplifed version of the atmosphere ocean coupling problem. The method resolves some of the issues that can be crucial to the fluid-fluid interaction problems: it is a partitioned time stepping method, yet it is of high order accuracy in both space and time (the two-step algorithm considered in this report provides second order accuracy); it allows for the usage of the legacy codes (which is a common requirement when resolving flows in complex geometries), yet it can be applied …

Pattern Formation Of Elastic Waves And Energy Localization Due To Elastic Gratings, A. Berezovski, J. Engelbrecht, Mihhail Berezovski

Publications

Elastic wave propagation through diffraction gratings is studied numerically in the plane strain setting. The interaction of the waves with periodically ordered elastic inclusions leads to a self-imaging Talbot effect for the wavelength equal or close to the grating size. The energy localization is observed at the vicinity of inclusions in the case of elastic gratings. Such a localization is absent in the case of rigid gratings.

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.

On The Selection Of A Good Shape Parameter For Rbf Approximation And Its Application For Solving Pdes, Lei-Hsin Kuo

Dissertations

Meshless methods utilizing Radial Basis Functions~(RBFs) are a numerical method that require no mesh connections within the computational domain. They are useful for solving numerous real-world engineering problems. Over the past decades, after the 1970s, several RBFs have been developed and successfully applied to recover unknown functions and to solve Partial Differential Equations (PDEs).
However, some RBFs, such as Multiquadratic (MQ), Gaussian (GA), and Matern functions, contain a free variable, the shape parameter, c. Because c exerts a strong influence on the accuracy of numerical solutions, much effort has been devoted to developing methods for determining shape parameters which provide …

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 …

Time Integration Methods Of Fundamental Solutions And Approximate Fundamental Solutions For Nonlinear Elliptic Partial Differential Equations, Corey Leon Jones

Dissertations

A time-dependent method is coupled with the Method of Approximate Particular Solutions (MAPS) of Delta-shaped basis functions, the Method of Fundamental Solutions (MFS), and the Method of Approximate Fundamental Solutions (MAFS) to solve a second order nonlinear elliptic partial differential equation (PDE) on regular and irregular shaped domains. The nonlinear PDE boundary value problem is first transformed into a time-dependent quasilinear problem by introducing a fictitious time. Forward Euler integration is then used to ultimately convert the problem into a sequence of time-dependent linear nonhomogeneous modified Helmholtz boundary value problems on which the superposition principle is applied to split the …

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 …

Numerical Solution For The Systems Of Variable-Coefficient Coupled Burgers’ Equation By Two-Dimensional Legendre Wavelets Method, Hossein Aminikhah, Sakineh Moradian

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, a numerical method for solving the systems of variable-coefficient coupled Burgers’ equation is proposed. The method is based on two-dimensional Legendre wavelets. Two-dimensional operational matrices of integration are introduced and then employed to find a solution to the systems of variable-coefficient coupled Burgers’ equation. Two examples are presented to illustrate the capability of the method. It is shown that the numerical results are in good agreement with the exact solutions for each problem.

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 …

An Applied Functional And Numerical Analysis Of A 3-D Fluid-Structure Interactive Pde, Thomas J. Clark

Department of Mathematics: Dissertations, Theses, and Student Research

We will present qualitative and numerical results on a partial differential equation (PDE) system which models a certain fluid-structure dynamics. In Chapter \ref{ChWellposedness}, the wellposedness of this PDE model is established by means of constructing for it a nonstandard semigroup generator representation; this representation is essentially accomplished by an appropriate elimination of the pressure. This coupled PDE model involves the Stokes system which evolves on a three dimensional domain $\mathcal{O}$ being coupled to a fourth order plate equation, possibly with rotational inertia parameter $\rho >0$, which evolves on a flat portion $\Omega$ of the boundary of $\mathcal{O}$. The coupling on …