Open Access. Powered by Scholars. Published by Universities.®

Numerical Analysis and Computation Commons™

Open Access. Powered by Scholars. Published by Universities.®

Partial Differential Equations

Institution
Keyword
Publication Year
Publication
Publication Type
File Type

Articles 31 - 60 of 265

Full-Text Articles in Numerical Analysis and Computation

Domain Decomposition Methods For Fluid-Structure Interaction Problems Involving Elastic, Porous, Or Poroelastic Structures, Hemanta Kunwar May 2024

Domain Decomposition Methods For Fluid-Structure Interaction Problems Involving Elastic, Porous, Or Poroelastic Structures, Hemanta Kunwar

All Dissertations

We introduce two global-in-time domain decomposition methods, namely the Steklov-Poincare method and Schwarz waveform relaxation (SWR) method using Robin transmission conditions (or the Robin method), for solving fluid-structure interaction systems involving elastic, porous, or poroelastic structure. These methods allow us to formulate the coupled system as a space-time interface problem and apply iterative algorithms directly to the evolutionary problem. Each time-dependent fluid and the structure subdomain problem is solved independently, which enables the use of different time discretization schemes and time step sizes in the subsystems. This leads to an efficient way of simulating time-dependent multiphysics phenomena. For the fluid-porous …

Go to article

Tools For Biomolecular Modeling And Simulation, Xin Yang Apr 2024

Tools For Biomolecular Modeling And Simulation, Xin Yang

Mathematics Theses and Dissertations

Electrostatic interactions play a pivotal role in understanding biomolecular systems, influencing their structural stability and functional dynamics. The Poisson-Boltzmann (PB) equation, a prevalent implicit solvent model that treats the solvent as a continuum while describes the mobile ions using the Boltzmann distribution, has become a standard tool for detailed investigations into biomolecular electrostatics. There are two primary methodologies: grid-based finite difference or finite element methods and body-fitted boundary element methods. This dissertation focuses on developing fast and accurate PB solvers, leveraging both methodologies, to meet diverse scientific needs and overcome various obstacles in the field.

Go to article

Advancements In Elliptic Mesh Redistribution And Automated Overset Grid Generation For Computational Fluid Dynamics, Chase Ashby Jan 2024

Advancements In Elliptic Mesh Redistribution And Automated Overset Grid Generation For Computational Fluid Dynamics, Chase Ashby

Theses and Dissertations--Mathematics

Automated and robust structured curvilinear overset grid generation remains a significant challenge for the application of computational fluid dynamics. This thesis intro- duces a specialized anisotropic Mach cone aligned adaptation algorithm for low-boom simulations and an automated procedure for generating structured curvilinear overset grids for immersed boundary simulations. Results demonstrate that anisotropic mesh adaptation can reduce computational resource usage by over 50% while maintaining sonic boom prediction accuracy across the entire sonic boom domain compared to user-constructed Mach cone aligned grids. The proposed immersed overset grid gen- eration procedure eliminates the need for complex surface generation algorithms, producing grids suitable …

Go to article

Multiscale Modelling Of Brain Networks And The Analysis Of Dynamic Processes In Neurodegenerative Disorders, Hina Shaheen Jan 2024

Multiscale Modelling Of Brain Networks And The Analysis Of Dynamic Processes In Neurodegenerative Disorders, Hina Shaheen

Theses and Dissertations (Comprehensive)

The complex nature of the human brain, with its intricate organic structure and multiscale spatio-temporal characteristics ranging from synapses to the entire brain, presents a major obstacle in brain modelling. Capturing this complexity poses a significant challenge for researchers. The complex interplay of coupled multiphysics and biochemical activities within this intricate system shapes the brain's capacity, functioning within a structure-function relationship that necessitates a specific mathematical framework. Advanced mathematical modelling approaches that incorporate the coupling of brain networks and the analysis of dynamic processes are essential for advancing therapeutic strategies aimed at treating neurodegenerative diseases (NDDs), which afflict millions of …

Go to article

A Class Of Game-Theoretic And Fokker-Planck Optimal Control Frameworks In Colon And Esophageal Cancer, Mesfer Alajmi Phd Jan 2024

A Class Of Game-Theoretic And Fokker-Planck Optimal Control Frameworks In Colon And Esophageal Cancer, Mesfer Alajmi Phd

Mathematics Dissertations - Archive

In this dissertation, we first present a new stochastic framework for parameter estimation and uncertainty quantification in colon cancer-induced immune responses. A stochastic process that captures the system's inherent randomness determines the dynamics of colon cancer. The stochastic framework is based on the Fokker-Planck equation, which represents the evolution of the probability density function corresponding to the stochastic process. We formulate an optimization problem that takes individual patient data with randomness present and solves it to obtain the unknown parameters corresponding to the individual tumor characteristics. Furthermore, we perform a sensitivity analysis of the optimal parameter set to identify the …

Go to article

Discontinuous Galerkin Methods For Compressible Miscible Displacements And Applications In Reservoir Simulation, Yue Kang Jan 2024

Discontinuous Galerkin Methods For Compressible Miscible Displacements And Applications In Reservoir Simulation, Yue Kang

Dissertations, Master's Theses and Master's Reports

This dissertation contains research on discontinuous Galerkin (DG) methods applied to the system of compressible miscible displacements, which is widely adopted to model surfactant flooding in enhanced oil recovery (EOR) techniques. In most scenarios, DG methods can effectively simulate problems in miscible displacements.
However, if the problem setting is complex, the oscillations in the numerical results can be detrimental, with severe overshoots leading to nonphysical numerical approximations. The first way to address this issue is to apply the bound-preserving
technique. Therefore, we adopt a bound-preserving Discontinuous Galerkin method
with a Second-order Implicit Pressure Explicit Concentration (SIPEC) time marching
method to …

Go to article

Mathematical Modeling Of Coupled Heat And Mass Transfer In Metal-Hydride Hydrogen Storage Systems, Muhammad Hasnain Jan 2024

Mathematical Modeling Of Coupled Heat And Mass Transfer In Metal-Hydride Hydrogen Storage Systems, Muhammad Hasnain

College of Graduate Studies: Theses & Dissertations

As a promising clean energy carrier hydrogen has recently gained significant interest, but its efficient and safe storage is a major challenge. Compared to the gaseous state and liquid state, metal hydrides (MH) offer a potentially more effective storage approach for hydrogen. However, the main challenge in this approach is the low thermal conductivity of the MH bed that leads to low heat transfer and ultimately to higher charging and discharging times. The purpose of this work is to develop an in-house comprehensive heat and mass transfer model for hydrogen sorption in MH reactors to simulate the dynamic behavior of …

Go to article

Numerical Modeling Of Thermal Runaway In Lithium-Ion Batteries Using Decomposition Kinetics And Inter-Cell Contact Resistance, Shehzad Khan Jan 2024

Numerical Modeling Of Thermal Runaway In Lithium-Ion Batteries Using Decomposition Kinetics And Inter-Cell Contact Resistance, Shehzad Khan

College of Graduate Studies: Theses & Dissertations

Lithium-ion batteries (LIBs) are central in numerous high-demand applications due to their high energy density and prolonged cycle life. Despite these advantages, their susceptibility to thermal runaway (TR) poses a significant safety risk, with the potential for catastrophic failures. This study focuses on the thermal behavior of prismatic lithium-ion cells, using a finite volume-based partial differential equation (PDE) solver developed in MATLAB and JULIA to model TR behavior. This solver accurately simulates transient behaviors, convection, diffusion, and source terms across various coordinate systems. By engaging in a series of increasing complex case studies, this research aims to identify the critical …

Go to article

Simulation Of Wave Propagation In Granular Particles Using A Discrete Element Model, Syed Tahmid Hussan Jan 2024

Simulation Of Wave Propagation In Granular Particles Using A Discrete Element Model, Syed Tahmid Hussan

College of Graduate Studies: Theses & Dissertations

The understanding of Bender Element mechanism and utilization of Particle Flow Code (PFC) to simulate the seismic wave behavior is important to test the dynamic behavior of soil particles. Both discrete and finite element methods can be used to simulate wave behavior. However, Discrete Element Method (DEM) is mostly suitable, as the micro scaled soil particle cannot be fully considered as continuous specimen like a piece of rod or aluminum. Recently DEM has been widely used to study mechanical properties of soils at particle level considering the particles as balls. This study represents a comparative analysis of Voigt and Best …

Go to article

Computational Modeling Using A Novel Continuum Approach Coupled With Pathway-Informed Neural Networks To Optimize Dynein-Mediated Centrosome Positioning In Polarized Cells, Arkaprovo Ghosal, Padmanabhan Seshaiyar Dr., Adriana Dawes Dr., General Genomics Inc. Nov 2023

Computational Modeling Using A Novel Continuum Approach Coupled With Pathway-Informed Neural Networks To Optimize Dynein-Mediated Centrosome Positioning In Polarized Cells, Arkaprovo Ghosal, Padmanabhan Seshaiyar Dr., Adriana Dawes Dr., General Genomics Inc.

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Go to article

Lnksc Method On Pde-Constrained Optimization For Mcf-7 Breast Cancer Cell Growth Predictions And Treatment Response With Gold Nanoparticles, Widodo Samyono, Shakhawat Bhuiyan Nov 2023

Lnksc Method On Pde-Constrained Optimization For Mcf-7 Breast Cancer Cell Growth Predictions And Treatment Response With Gold Nanoparticles, Widodo Samyono, Shakhawat Bhuiyan

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Go to article

Helices In Fluids And Applications To Modeling In Biology, Eva M. Strawbridge Nov 2023

Helices In Fluids And Applications To Modeling In Biology, Eva M. Strawbridge

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Go to article

A Tale Of Two Diagonalizations: Methods To Diagonalize A 1-D Piecewise Constant Indefinite Schrödinger Operator, Sarah Wright Aug 2023

A Tale Of Two Diagonalizations: Methods To Diagonalize A 1-D Piecewise Constant Indefinite Schrödinger Operator, Sarah Wright

Dissertations

We present two numerical methods for computing the solution of a partial differential equation (PDE) for modeling acoustic pressure, known as an extra-wide angle parabolic equation, that features the square root of a differential operator. The differential operator is the negative of an indefinite Schrödinger operator with a piecewise constant potential. This work primarily deals with the 3-piece case; however, a generalization is made to the case of an arbitrary number of pieces. In the first method, the Rayleigh-Secant Method, through restriction to a judiciously chosen lower-dimensional subspace, approximate eigenfunctions are used to obtain estimates for the eigenvalues of the …

Go to article

A Robust Model Reduction For The Optimal Boundary Feedback Stabilization Of Magnetizable Piezoelectric Beams, Md Rafi As Sadeq Ibn Emran Aug 2023

A Robust Model Reduction For The Optimal Boundary Feedback Stabilization Of Magnetizable Piezoelectric Beams, Md Rafi As Sadeq Ibn Emran

Masters Theses & Specialist Projects

The one-dimensional Partial Differential Equation (PDE) model of the wave equation, which describes wave dynamics of a wide-range of controlled mechanical systems, with a state feedback controller at its boundary is known to have exponentially stable solutions. Moreover, it is also reported that the several model reductions of the wave equation by the standard Finite Differences and Finite Elements approximations suffer from the lack of exponential stability (and exact observability without a state feedback controller) uniformly as the discretization parameter tends to zero. This is due to the loss of a uniform gap among the high-frequency eigenvalues as the discretization …

Go to article

Null Space Removal In Finite Element Discretizations, Pengfei Jia Aug 2023

Null Space Removal In Finite Element Discretizations, Pengfei Jia

All Theses

Partial differential equations are frequently utilized in the mathematical formulation of physical problems. Boundary conditions need to be applied in order to obtain the unique solution to such problems. However, some types of boundary conditions do not lead to unique solutions because the continuous problem has a null space. In this thesis, we will discuss how to solve such problems effectively. We first review the foundation of all three problems and prove that Laplace problem, linear elasticity problem and Stokes problem can be well posed if we restrict the test and trial space in the continuous and discrete finite element …

Go to article

Solving The Cable Equation, A Second-Order Time Dependent Pde For Non-Ideal Cables With Action Potentials In The Mammalian Brain Using Kss Methods, Nirmohi Charbe Jun 2023

Solving The Cable Equation, A Second-Order Time Dependent Pde For Non-Ideal Cables With Action Potentials In The Mammalian Brain Using Kss Methods, Nirmohi Charbe

Master's Theses

In this thesis we shall perform the comparisons of a Krylov Subspace Spectral method with Forward Euler, Backward Euler and Crank-Nicolson to solve the Cable Equation. The Cable Equation measures action potentials in axons in a mammalian brain treated as an ideal cable in the first part of the study. We shall subject this problem to the further assumption of a non-ideal cable. Assume a non-uniform cross section area along the longitudinal axis. At the present time, the effects of torsion, curvature and material capacitance are ignored. There is particular interest to generalize the application of the PDEs including and …

Go to article

Advancements In Fluid Simulation Through Enhanced Conservation Schemes, Sean Ingimarson May 2023

Advancements In Fluid Simulation Through Enhanced Conservation Schemes, Sean Ingimarson

All Dissertations

To better understand and solve problems involving the natural phenomenon of fluid and air flows, one must understand the Navier-Stokes equations. Branching several different fields including engineering, chemistry, physics, etc., these are among the most important equations in mathematics. However, these equations do not have analytic solutions save for trivial solutions. Hence researchers have striven to make advancements in varieties of numerical models and simulations. With many variations of numerical models of the Navier-Stokes equations, many lose important physical meaningfulness. In particular, many finite element schemes do not conserve energy, momentum, or angular momentum. In this thesis, we will study …

Go to article

Effects Of Slip On Highly Viscous Thin-Film Flows Inside Vertical Tubes (Constant Radius, Constricted And Flexible), Mark S. Schwitzerlett Jan 2023

Effects Of Slip On Highly Viscous Thin-Film Flows Inside Vertical Tubes (Constant Radius, Constricted And Flexible), Mark S. Schwitzerlett

Theses and Dissertations

Viscous liquid film flows in a tube arise in numerous industrial and biological applications, including the transport of mucus in human airways. Previous modeling studies have typically used no-slip boundary conditions, but in some applications the effects of slip at the boundary may not be negligible. We derive a long-wave model based on lubrication theory which allows for slippage along the boundary. Linear stability analysis verifies the impact of slip-length on the speed, growth rate, and wavelength of the most unstable mode. Nonlinear simulations demonstrate the impact of slip-length on plug formation and wave dynamics. These simulations are conducted for …

Go to article

Using Modflow To Assess Groundwater Storage Enhancement Via A Floodplain Infiltration Basin, Lindsay Henning Jan 2023

Using Modflow To Assess Groundwater Storage Enhancement Via A Floodplain Infiltration Basin, Lindsay Henning

All Master's Theses

Delaying groundwater discharge into rivers until it is critically needed during baseflow conditions provides promise for lowering elevated stream temperatures and improving habitat for aquatic species. Increasing groundwater storage may accomplish this in locations where excess spring runoff can be captured and allowed to infiltrate into the subsurface for later beneficial use, a process known as Managed Aquifer Recharge (MAR). Here, MAR via an infiltration basin is considered at a site along the Teanaway River in central Washington State. The effects of simulated ephemeral ponds of sizes varying from 554 m3 to 2430 m3 (0.449 acre-feet to 1.97 …

Go to article

Novel Approaches To Compute Manifold Operators With The Radial Basis Functions Method, Jacob James Blazejewski Jan 2023

Novel Approaches To Compute Manifold Operators With The Radial Basis Functions Method, Jacob James Blazejewski

Dissertations, Master's Theses and Master's Reports

The bulk of this dissertation is mainly composed of four chapters, which are organized as follows: Chapter 1 provides an introduction to the Radial Basis Functions (RBF) method by briefly outlining its historical developments and reviewing the RBF interpolation and the RBF-Finite Difference (FD) methodologies, and their advantages/disadvantages. Chapter 2 describes the Orthogonal Gradients (OGr) method and the Fast OGr method and how these can be used to compute differential operators restricted to hypersurfaces and space curves ($\Gamma$) embedded in R3. We will highlight a challenge of pairing Fast OGr with RBF-FD on nearly flat local clusters and how to …

Go to article

(R1894) Invariant Solution For Two-Dimensional And Axisymmetric Jet Of Power-Law Fluids, Bhavixa Bhagat, M. G. Timol Dec 2022

(R1894) Invariant Solution For Two-Dimensional And Axisymmetric Jet Of Power-Law Fluids, Bhavixa Bhagat, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

An invariant solution is derived using the Lie symmetry technique for steady laminar two-dimensional and axisymmetric boundary layer jet flow of incompressible power-law fluids with appropriate boundary conditions. Using symmetry, the nonlinear partial differential equation of the jet flow problem is transformed into a nonlinear ordinary differential equation. The resultant nonlinear ordinary differential equation with boundary conditions is converted to an initial value problem using the Lie symmetry technique. A numerical solution for the resulting initial value problem is derived using Fehlberg’s fourth-fifth order Runge-Kutta method through Maple software. The graphical representation of the characteristics of the velocity field for …

Go to article

On Analysis Of Effectiveness Controlling Covid-19 With Quarantine And Vaccination Compartments In Indonesia, Prihantini Prihantini Nov 2022

On Analysis Of Effectiveness Controlling Covid-19 With Quarantine And Vaccination Compartments In Indonesia, Prihantini Prihantini

Annual Symposium on Biomathematics and Ecology Education and Research

No abstract provided.

Go to article

(Si10-067) Numerical Study Of The Time Fractional Burgers’ Equation By Using Explicit And Implicit Schemes, Swapnali Doley, A. Vanav Kumar, L. Jino Oct 2022

(Si10-067) Numerical Study Of The Time Fractional Burgers’ Equation By Using Explicit And Implicit Schemes, Swapnali Doley, A. Vanav Kumar, L. Jino

Applications and Applied Mathematics: An International Journal (AAM)

The study discusses the numerical solution for a time fractional Burgers’ equation using explicit (scheme 1) and implicit scheme (scheme 2), respectively. The approximation of the differential equation is discretized using the finite difference method (FDM). A non-linear term present in the Burgers’ equation is approximated using the time-averaged values. The Von-Neumann analysis shows that the Scheme 1 is conditionally stable and Scheme 2 is unconditionally stable. The numerical solutions are compared with the exact solutions and are good in agreement. Also, the error is estimated between exact and numerical solutions.

Go to article

Navier-Stokes Equations In One And Two Dimensions, Jon Nerdal Jul 2022

Navier-Stokes Equations In One And Two Dimensions, Jon Nerdal

LSU Master's Theses

The Navier-Stokes equations are an important tool in understanding and describing fluid flow. We investigate different formulations of the incompressible Navier-Stokes equations in the one-dimensional case along an axis and in the two-dimensional case in a circular pipe without swirl. For the one-dimensional case we show that the velocity approximations are remarkably accurate and we suggest that understanding this simple axial behaviour is an important starting point for further exploration in higher dimensions. The complexity of the boundary is then increased with the two-dimensional case of fluid flow through the cross section of a circular pipe, where we investigate two …

Go to article

A Bidirectional Formulation For Walk On Spheres, Yang Qi Jun 2022

A Bidirectional Formulation For Walk On Spheres, Yang Qi

Dartmouth College Master’s Theses

Poisson’s equations and Laplace’s equations are important linear partial differential equations (PDEs)
widely used in many applications. Conventional methods for solving PDEs numerically often need to
discretize the space first, making them less efficient for complex shapes. The random walk on spheres
method (WoS) is a grid-free Monte-Carlo method for solving PDEs that does not need to discrete the
space. We draw analogies between WoS and classical rendering algorithms, and find that the WoS
algorithm is conceptually identical to forward path tracing.
We show that solving the Poisson’s equation is equivalent to solving the Green’s function for every
pair of …

Go to article

Mathematical Modeling Of Brain Cancer Growth Using A Level-Set Method, Gbocho M. Terasaki May 2022

Mathematical Modeling Of Brain Cancer Growth Using A Level-Set Method, Gbocho M. Terasaki

Biology and Medicine Through Mathematics Conference

No abstract provided.

Go to article

A Conservative Numerical Scheme For The Multilayer Shallow Water Equations, Evan Butterworth May 2022

A Conservative Numerical Scheme For The Multilayer Shallow Water Equations, Evan Butterworth

All Theses

An energy-conserving numerical scheme is developed for the multilayer shallow water equations (SWE’s). The scheme is derived through the Hamiltonian formulation of the inviscid shallow water flows related to the vorticity-divergence variables. Through the employment of the skew-symmetric Poisson bracket, the continuous system for the multilayer SWE’s is shown to preserve an infinite number of quantities, most notably the energy and enstrophy. An energy-preserving numerical scheme is then developed through the careful discretization of the Hamiltonian and the Poisson bracket, ensuring the skew-symmetry of the latter. This serves as the groundwork for developing additional schemes that preserve other conservation properties …

Go to article

On The Consistency Of Alternative Finite Difference Schemes For The Heat Equation, Tran April Apr 2022

On The Consistency Of Alternative Finite Difference Schemes For The Heat Equation, Tran April

Rose-Hulman Undergraduate Mathematics Journal

While the well-researched Finite Difference Method (FDM) discretizes every independent variable into algebraic equations, Method of Lines discretizes all but one dimension, leaving an Ordinary Differential Equation (ODE) in the remaining dimension. That way, ODE's numerical methods can be applied to solve Partial Differential Equations (PDEs). In this project, Linear Multistep Methods and Method of Lines are used to numerically solve the heat equation. Specifically, the explicit Adams-Bashforth method and the implicit Backward Differentiation Formulas are implemented as Alternative Finite Difference Schemes. We also examine the consistency of these schemes.

Go to article

A Spatially And Temporally Second Order Method For Solving Parabolic Interface Problems, Kumudu Gamage, Yan Peng Apr 2022

A Spatially And Temporally Second Order Method For Solving Parabolic Interface Problems, Kumudu Gamage, Yan Peng

College of Sciences Posters

Parabolic interface problems have many applications in physics and biology, such as hyperthermia treatment of cancer, underground water flow, and food engineering. Here we present an algorithm for solving two-dimensional parabolic interface problems where the coefficient and the forcing term have a discontinuity across the interface. The Crank-Nicolson scheme is used for time discretization, and the direct immersed interface method is used for spatial discretization. The proposed method is second order in both space and time for both solution and gradients in maximum norm.

Go to article

A Meshless Approach To Computational Pharmacokinetics, Anthony Matthew Khoury Apr 2022

A Meshless Approach To Computational Pharmacokinetics, Anthony Matthew Khoury

Doctoral Dissertations and Master's Theses

The meshless method is an incredibly powerful technique for solving a variety of problems with unparalleled accuracy and efficiency. The pharmacokinetic problem of transdermal drug delivery (TDDD) is one such topic and is of significant complexity. The locally collocated meshless method (LCMM) is developed in solution to this topic. First, the meshless method is formulated to model this transport phenomenon and is then validated against an analytical solution of a pharmacokinetic problem set, to demonstrate this accuracy and efficiency. The analytical solution provides a locus by which convergence behavior are evaluated, demonstrating the super convergence of the locally collocated meshless …

Go to article

First Prev Next Last