Numerical Analysis and Computation 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. …

Controlled Manipulation And Transport By Microswimmers In Stokes Flows, Jake Buzhardt

All Dissertations

Remotely actuated microscale swimming robots have the potential to revolutionize many aspects of biomedicine. However, for the longterm goals of this field of research to be achievable, it is necessary to develop modelling, simulation, and control strategies which effectively and efficiently account for not only the motion of individual swimmers, but also the complex interactions of such swimmers with their environment including other nearby swimmers, boundaries, other cargo and passive particles, and the fluid medium itself. The aim of this thesis is to study these problems in simulation from the perspective of controls and dynamical systems, with a particular focus …

Modeling Single And Multiple Pacemaker Interaction In Jellyfish Locomotion, Alexander Hoover

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Boundary Integral Equation Methods For Superhydrophobic Flow And Integrated Photonics, Kosuke Sugita

Dissertations

This dissertation presents fast integral equation methods (FIEMs) for solving two important problems encountered in practical engineering applications.

The first problem involves the mixed boundary value problem in two-dimensional Stokes flow, which appears commonly in computational fluid mechanics. This problem is particularly relevant to the design of microfluidic devices, especially those involving superhydrophobic (SH) flows over surfaces made of composite solid materials with alternating solid portions, grooves, or air pockets, leading to enhanced slip.

The second problem addresses waveguide devices in two dimensions, governed by the Helmholtz equation with Dirichlet conditions imposed on the boundary. This problem serves as a …

Hydrodynamic And Physicochemical Interactions Between An Active Janus Particle And An Inactive Particle, Jessica S. Rosenberg

Dissertations, Theses, and Capstone Projects

Active matter is an area of soft matter science in which units consume energy and turn it into autonomous motion. Groups of these units – whether flocks of birds, bacterial colonies, or even collections of synthetically-made active particles – may exhibit complex behavior on large scales. While the large-scale picture is of great importance, so is the microscopic scale. Studying the individual particles that make up active matter will allow us to understand how they move, and whether and under what circumstances their activity can be controlled.

Here we delve into the world of active matter by studying colloidal-sized (100 …

The Magnetic Field Of Protostar-Disk-Outflow Systems, Mahmoud Sharkawi

Electronic Thesis and Dissertation Repository

Recent observations of protostellar cores reveal complex magnetic field configurations that are distorted in the innermost disk region. Unlike the prestellar phase, where the magnetic field geometry is simpler with an hourglass configuration, magnetic fields in the protostellar phase are sculpted by the formation of outflows and rapid rotation. This gives rise to a significant azimuthal (or toroidal) component that has not yet been analytically modelled in the literature. Moreover, the onset of outflows, which act as angular momentum transport mechanisms, have received considerable attention in the past few decades. Two mechanisms: 1) the driving by the gradient of a …

Practical Implementation Of The Immersed Interface Method With Triangular Meshes For 3d Rigid Solids In A Fluid Flow, Norah Hakami

Mathematics Theses and Dissertations

When employing the immersed interface method (IIM) to simulate a fluid flow around a moving rigid object, the immersed object can be replaced by a virtual fluid enclosed by singular forces on the interface between the real and virtual fluids. These forces represent the impact of the rigid motion on the fluid flow and cause jump discontinuities across the interface in the whole flow field. Then, the IIM resolves the fluid flow on a fixed computational domain by directly incorporating the jump conditions across the interface into numerical schemes. Previous development of the method is limited to simple smooth boundaries. …

Impact Of Spallation And Internal Radiation On Fibrous Ablative Materials, Raghava Sai Chaitanya Davuluri

Theses and Dissertations--Mechanical Engineering

Space vehicles are equipped with Thermal Protection Systems (TPS) that encounter high heat rates and protect the payload while entering a planetary atmosphere. For most missions that interest NASA, ablative materials are used as TPS. These materials undergo several mass and energy transfer mechanisms to absorb intense heat. The size and construction of the TPS are based on the composition of the planetary atmosphere and the impact of various ablative mechanisms on the flow field and the material. Therefore, it is essential to quantify the rates of different ablative phenomena to model TPS accurately. In this work, the impact of …

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

Faculty Publications and Other Works -- Mathematics

In this report we review some preliminary work on the numerical solution of 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 accurate simulation. This requires sophisticated numerical algorithms and efficient code to realize predictions over desired space and time scales. …

An Examination Of Fontan Circulation Using Differential Equation Models And Numerical Methods, Vanessa Maybruck

Honors Student Research

Certain congenital heart defects can lead to the development of only a single pumping chamber, or ventricle, in the heart instead of the usual two ventricles. Individuals with this defect undergo a corrective, three-part surgery, the third step of which is the Fontan procedure, but as the patients age, their cardiovascular health will likely deteriorate. Using computational fluid dynamics and differential equations, Fontan circulation can be modeled to investigate why the procedure fails and how Fontan failure can be maximally prevented. Borrowing from well-established literature on RC circuits, the differential equation models simulate systemic blood flow in a piecewise, switch-like …

Numerical Reconstruction Of Spalled Particle Trajectories In An Arc-Jet Environment, Raghava S. C. Davuluri, Sean C. C. Bailey, Kaveh A. Tagavi, Alexandre Martin

Mechanical Engineering Faculty Publications

To evaluate the effects of spallation on ablative material, it is necessary to evaluate the mass loss. To do so, a Lagrangian particle trajectory code is used to reconstruct trajectories that match the experimental data for all kinematic parameters. The results from spallation experiments conducted at the NASA HYMETS facility over a wedge sample were used. A data-driven adaptive methodology was used to adapts the ejection parameters until the numerical trajectory matches the experimental data. The preliminary reconstruction results show that the size of the particles seemed to be correlated with the location of the ejection event. The size of …

Improving The Temporal Accuracy Of Turbulence Models And Resolving The Implementation Issues Of Fluid Flow Modeling, Kyle J. Schwiebert

Dissertations, Master's Theses and Master's Reports

A sizeable proportion of the work in this thesis focuses on a new turbulence model, dubbed ADC (the approximate deconvolution model with defect correction). The ADC is improved upon using spectral deferred correction, a means of constructing a higher order ODE solver. Since both the ADC and SDC are based on a predictor-corrector approach, SDC is incorporated with essentially no additional computational cost. We will show theoretically and using numerical tests that the new scheme is indeed higher order in time than the original, and that the benefits of defect correction, on which the ADC is based, are preserved.

The …

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

Articles

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

Accurate Spectral Algorithms For Solving Variable-Order Fractional Percolation Equations, M. A. Abdelkawy

Applications and Applied Mathematics: An International Journal (AAM)

A high accurate spectral algorithm for one-dimensional variable-order fractional percolation equations (VO-FPEs) is considered.We propose a shifted Legendre Gauss-Lobatto collocation (SL-GLC) method in conjunction with shifted Chebyshev Gauss-Radau collocation (SC-GR-C) method to solve the proposed problem. Firstly, the solution and its space fractional derivatives are expanded as shifted Legendre polynomials series. Then, we determine the expansion coefficients by reducing the VO-FPEs and its conditions to a system of ordinary differential equations (SODEs) in time. The numerical approximation of SODEs is achieved by means of the SC-GR-C method. The under-study’s problem subjected to the Dirichlet or non-local boundary conditions is presented …

Numerical Simulations Of Nonlinear Waves And Their Stability: Stokes Waves And Nonlinear Schroedinger Equation, Anastassiya Semenova

Mathematics & Statistics ETDs

The present work offers an investigation of dynamics and stability of nonlinear waves in Hamiltonian systems. The first part of the manuscript discusses the classical problem of water waves on the surface of an ideal fluid in 2D. We demonstrate how to construct the Stokes waves, and how to apply a continuation method to find waves in close vicinity to the limiting Stokes wave. We provide new insight into the stability of the Stokes waves by identifying previously inaccessible branches of instability in the equations of motion for the fluid. We provide numerical evidence that pairs of unstable eigenvalues of …

A Study Of The Design Of Adaptive Camber Winglets, Justin J. Rosescu

Master's Theses

A numerical study was conducted to determine the effect of changing the camber of a winglet on the efficiency of a wing in two distinct flight conditions. Camber was altered via a simple plain flap deflection in the winglet, which produced a constant camber change over the winglet span. Hinge points were located at 20%, 50% and 80% of the chord and the trailing edge was deflected between -5° and +5°. Analysis was performed using a combination of three-dimensional vortex lattice method and two-dimensional panel method to obtain aerodynamic forces for the entire wing, based on different winglet camber configurations. …

Local Non-Similar Solution Of Powell-Eyring Fluid Flow Over A Vertical Flat Plate, Hemangini Shukla, Hema C. Surati, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

Our objective is to obtain the non-similarity solution of non-Newtonian fluid for Powell-Eyring model by a local non-similarity method. Here, free stream velocity is considered in power-law form (𝑈=𝑥^m). The governing equations are transformed using non-similar transformations and derived equations are treated as ordinary differential equations. Non-similar solutions are obtained for different values of power-law index 𝑚 and stream-wise location 𝜉. Influence of various parameters on velocity and temperature field are presented graphically using MATLAB bvp4c solver.

Meshless Modeling Of Flow Dispersion And Progressive Piping In Poroelastic Levees, Anthony Khoury, Eduardo Divo, Alain J. Kassab, Sai Kakuturu, Lakshmi Reddi

Publications

Performance data on earth dams and levees continue to indicate that piping is one of the major causes of failure. Current criteria for prevention of piping in earth dams and levees have remained largely empirical. This paper aims at developing a mechanistic understanding of the conditions necessary to prevent piping and to enhance the likelihood of self-healing of cracks in levees subjected to hydrodynamic loading from astronomical and meteorological (including hurricane storm surge-induced) forces. Systematic experimental investigations are performed to evaluate erosion in finite-length cracks as a result of transient hydrodynamic loading. Here, a novel application of the localized collocation …

Numerical Solution Of 3rd Order Ode Using Fdm: On A Moving Surface In Mhd Flow Of Sisko Fluid, Manisha Patel, Jayshri Patel, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

A Similarity group theoretical technique is used to transform the governing nonlinear partial differential equations of two dimensional MHD boundary layer flow of Sisko fluid into nonlinear ordinary differential equations. Then the resulting third order nonlinear ordinary differential equation with corresponding boundary conditions is linearised by Quasi linearization method. Numerical solution of the linearised third order ODE is obtained using Finite Difference method (FDM). Graphical presentation of the solution is given.

Solitary And Periodic Exact Solutions Of The Viscosity-Capillarity Van Der Waals Gas Equations, Emad A. Az-Zo’Bi

Applications and Applied Mathematics: An International Journal (AAM)

Periodic and soliton solutions are derived for the (1+1)-dimensional van der Waals gas system in the viscosity-capillarity regularization form. The system is handled via the e-^φ(ξ) -expansion method. The obtained solutions have been articulated by the hyperbolic, trigonometric, exponential and rational functions with arbitrary constants. Mathematical analysis and numerical graphs are provided for some solitons, periodic and kink solitary wave solutions to visualize the dynamics of equations. Obtained results reveal that the method is very influential and effective tool for solving nonlinear partial differential equations in applied mathematics.

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 …

Call For Abstracts - Resrb 2019, July 8-9, Wrocław, Poland, Wojciech M. Budzianowski

Wojciech Budzianowski

No abstract provided.

Numerical Treatment For The Flow Of Casson Fluid And Heat Transfer Model Over An Unsteady Stretching Surface In The Presence Of Internal Heat Generation/Absorption And Thermal Radiation, Mohammed M. Babatin

Applications and Applied Mathematics: An International Journal (AAM)

Several important industrial and engineering problems are very difficult to solve analytically since they are high nonlinear. The Chebyshev spectral collocation method possesses an ability to predict the solution behavior for a system of high nonlinear ordinary differential equations. This method which is based on differentiated Chebyshev polynomials is introduced to obtain an approximate solution to the system of ordinary differential equations which physically describe the flow and heat transfer problem of an unsteady Casson fluid model taking into consideration both heat generation and radiation effects in the temperature equation. Based on the spectral collocation method, the obtained solution is …

Quasi-Linearization Method With Rational Legendre Collocation Method For Solving Mhd Flow Over A Stretching Sheet With Variable Thickness And Slip Velocity Which Embedded In A Porous Medium, M. M. Khader

Applications and Applied Mathematics: An International Journal (AAM)

The quasi-linearization method (QLM) and the rational Legendre functions are introduced here to present the numerical solution for the Newtonian fluid flow past an impermeable stretching sheet which embedded in a porous medium with a power-law surface velocity, variable thickness and slip velocity. Firstly, due to the high nonlinearity which yielded from the ordinary differential equation which describes the proposed physical problem, we construct a sequence of linear ODEs by using the QLM, hence the resulted equations become a system of linear algebraic equations. The comparison with the available results in the literature review proves that the obtained results via …

Numerical Studies For Mhd Flow And Gradient Heat Transport Past A Stretching Sheet With Radiation And Heat Production Via Dtm, Khadijah M. Abualnaja

Applications and Applied Mathematics: An International Journal (AAM)

This paper presents a numerically study for the effect of the internal heat generation, magnetic field and thermal radiation effects on the flow and gradient heat transfer of a Newtonian fluid over a stretching sheet. By using a similarity transformation, the governing PDEs can be transformed into a coupled non-linear system of ODEs with variable coefficients. Numerical solutions for these equations subject to appropriate boundary conditions are obtained by using the differential transformation method (DTM). The effects of various physical parameters such as viscosity parameter, the suction parameter, the radiation parameter, internal heat generation or absorption parameter and the Prandtl …

Patient-Specific Multiscale Computational Fluid Dynamics Assessment Of Embolization Rates In The Hybrid Norwood: Effects Of Size And Placement Of The Reverse Blalock–Taussig Shunt, Ray Prather, John Seligson, Marcus Ni, Eduardo Divo, Alain J. Kassab, William Decampli

Publications

The hybrid Norwood operation is performed to treat hypoplastic left heart syndrome. Distal arch obstruction may compromise flow to the brain. In a variant of this procedure, a synthetic graft (reverse Blalock–Taussig shunt) is placed between the pulmonary trunk and innominate artery to improve upper torso blood flow. Thrombi originating in the graft may embolize to the brain. In this study, we used computational fluid dynamics and particle tracking to investigate the patterns of particle embolization as a function of the anatomic position of the reverse Blalock–Taussig shunt. The degree of distal arch obstruction and position of particle origin influence …

Risk Assessment Of Dropped Cylindrical Objects In Offshore Operations, Adelina Steven

University of New Orleans Theses and Dissertations

Dropped object are defined as any object that fall under its own weight from a previously static position or fell due to an applied force from equipment or a moving object. It is among the top ten causes of injuries and fatality in oil and gas industry. To solve this problem, several in-house tools and guidelines is developed over time to assess the risk of dropped objects on the sub-sea structures. This thesis focuses on compiling and comparing those methods in hope to improve the recommended practices available in the market. A simple modification is done on the in-house tools …

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 …

High Accuracy Methods And Regularization Techniques For Fluid Flows And Fluid-Fluid Interaction, Mustafa Aggul

Dissertations, Master's Theses and Master's Reports

This dissertation contains several approaches to resolve irregularity issues of CFD problems, including a decoupling of non-linearly coupled fluid-fluid interaction, due to high Reynolds number. New models present not only regularize the linear systems but also produce high accurate solutions both in space and time. To achieve this goal, methods solve a computationally attractive artificial viscosity approximation of the target problem, and then utilize a correction approach to make it high order accurate. This way, they all allow the usage of legacy code | a frequent requirement in the simulation of fluid flows in complex geometries. In addition, they all …

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 …