Digital Commons Network^™

Preferential Stiffness And The Crack-Tip Fields Of An Elastic Porous Solid Based On The Density-Dependent Moduli Model, Hyun C. Yoon, S. M. Mallikarjunaiah, Dambaru Bhatta

School of Mathematical and Statistical Sciences Faculty Publications and Presentations

In this paper, we study the preferential stiffness and the crack-tip fields for an elastic porous solid of which material properties are dependent upon the density. Such a description is necessary to describe the failure that can be caused by damaged pores in many porous bodies such as ceramics, concrete and human bones. To that end, we revisit a new class of implicit constitutive relations under the assumption of small deformation. Although the constitutive relationship \textit{appears linear} in both the Cauchy stress and linearized strain, the governing equation bestowed from the balance of linear momentum results in a quasi-linear partial …

Second-Order, Fully Decoupled, Linearized, And Unconditionally Stable Scalar Auxiliary Variable Schemes For Cahn–Hilliard–Darcy System, Yali Gao, Xiaoming He, Yufeng Nie

Mathematics and Statistics Faculty Research & Creative Works

In this paper, we establish the fully decoupled numerical methods by utilizing scalar auxiliary variable approach for solving Cahn–Hilliard–Darcy system. We exploit the operator splitting technique to decouple the coupled system and Galerkin finite element method in space to construct the fully discrete formulation. The developed numerical methods have the features of second order accuracy, totally decoupling, linearization, and unconditional energy stability. The unconditionally stability of the two proposed decoupled numerical schemes are rigorously proved. Abundant numerical results are reported to verify the accuracy and effectiveness of proposed numerical methods.

Numerical Analysis Of A Second Order Ensemble Method For Evolutionary Magnetohydrodynamics Equations At Small Magnetic Reynolds Number, John Carter, Nan Jiang

Mathematics and Statistics Faculty Research & Creative Works

We study a second order ensemble method for fast computation of an ensemble of magnetohydrodynamics flows at small magnetic Reynolds number. Computing an ensemble of flow equations with different input parameters is a common procedure for uncertainty quantification in many engineering applications, for which the computational cost can be prohibitively expensive for nonlinear complex systems. We propose an ensemble algorithm that requires only solving one linear system with multiple right-hands instead of solving multiple different linear systems, which significantly reduces the computational cost and simulation time. Comprehensive stability and error analyses are presented proving conditional stability and second order in …

Implementation Of A Least Squares Method To A Navier-Stokes Solver, Jada P. Lytch, Taylor Boatwright, Ja'nya Breeden

Rose-Hulman Undergraduate Mathematics Journal

The Navier-Stokes equations are used to model fluid flow. Examples include fluid structure interactions in the heart, climate and weather modeling, and flow simulations in computer gaming and entertainment. The equations date back to the 1800s, but research and development of numerical approximation algorithms continues to be an active area. To numerically solve the Navier-Stokes equations we implement a least squares finite element algorithm based on work by Roland Glowinski and colleagues. We use the deal.II academic library , the C++ language, and the Linux operating system to implement the solver. We investigate convergence rates and apply the least squares …

A Weak Fractional Calculus Theory And Numerical Methods For Fractional Differential Equations, Mitchell D. Sutton

Doctoral Dissertations

This dissertation is comprised of four integral parts. The first part comprises a self-contained new theory of weak fractional differential calculus in one-dimension. The crux of this new theory is the introduction of a weak fractional derivative notion which is a natural generalization of integer order weak derivatives; it also helps to unify multiple existing fractional derivative definitions.

The second part of this work presents three new families of fractional Sobolev spaces and their accompanying theory in one-dimension. The new construction and theory are based on a newly developed notion of weak fractional derivatives, which are natural generalizations of the …