Applied Mathematics Commons^™

Wavelet-Based Analysis Of Neutron-Induced Photon Spectral Data, Bruce Kessler, Alexander Barzilov, Phillip Womble

Mathematics Faculty Publications

Neutron-based methods of non-destructive inter- rogation of objects for the purpose of their characterization are well-established techniques, employed in the field of bulk material analysis, contraband detection, unexploded ordnance, etc. The characteristic gamma rays produced in nuclear reactions initiated by neutrons in the volume of the irradiated object (inelastic neutron scattering, thermal neutron capture, and activation) are used for the elemental identification. In many real-world applications, an automated spectral analysis is needed, and many algorithms are used for that purpose. The Applied Physics Institute at Western Kentucky University has recently started to employ a mathematical spectrum analysis technique based on …

A Nonnegative Analysis Of Politics, Tim Chartier, Charles D. Wessell

Math Faculty Publications

The article investigates how linear algebra can recover mathematical information from the electronic messages using the Enron Email Sets in Pennsylvania. It states that a term-by-email matrix has been created to cluster algorithm, which allows one to mine through data and discover meaning. Moreover, nonnegative matrix factorization (NMF) enables one to interpret the resulting factorization in terms of the original problem.

A Numerical Method For Solving The Elliptic And Elasticity Interface Problems, Liqun Wang

Doctoral Dissertations

Interface problems arise when dealing with physical problems composed of different materials or of the same material at different states. Because of the irregularity along interfaces, many common numerical methods do not work, or work poorly, for interface problems. Matrix-coefficient elliptic and elasticity equations with oscillatory solutions and sharp-edged interfaces are especially complicated and challenging for most existing methods. An accurate and efficient method is desired.

In 1999, the boundary condition capturing method was proposed to deal with Poisson equations with interfaces whose variable coefficients and solutions may be discontinuous. In 2003, a weak formulation was derived. Built on previous …

Shape Reconstruction And Classification Using The Response Matrix, Wei Wang

Doctoral Dissertations

This dissertation presents a novel method for the inverse scattering problem for extended target. The acoustic or electromagnetic wave is scattered by the target and received by all the transducers around the target. The scattered field on all the transducers forms the response matrix which contains the information of the geometry of the target. The objective of the inverse scattering problem is to reconstruct the shape of the scatter using the Response Matrix.

There are two types of numerical methods for solving the inverse problem: the direct imaging method and the iterative method. Two direct imaging methods, MUSIC method and …

Results In Lattices, Ortholattices, And Graphs, Jianning Su

Doctoral Dissertations

This dissertation contains two parts: lattice theory and graph theory. In the lattice theory part, we have two main subjects. First, the class of all distributive lattices is one of the most familiar classes of lattices. We introduce "π-versions" of five familiar equivalent conditions for distributivity by applying the various conditions to 3-element antichains only. We prove that they are inequivalent concepts, and characterize them via exclusion systems. A lattice L satisfies D0π, if a ✶ (b ✶ c) ≤ (a ✶ b) ✶ c for all 3-element antichains { a, b, c}. We consider …

Quantitative Stability Of Linear Infinite Inequality Systems Under Block Perturbations With Applications To Convex Systems, M J. Cánovas, M A. Lopez, Boris S. Mordukhovich, J Parra

Mathematics Research Reports

The original motivation for this paper was to provide an efficient quantitative analysis of convex infinite (or semi-infinite) inequality systems whose decision variables run over general infinite-dimensional (resp. finite-dimensional) Banach spaces and that are indexed by an arbitrary fixed set J. Parameter perturbations on the right-hand side of the inequalities are required to be merely bounded, and thus the natural parameter space is loo(J). Our basic strategy consists of linearizing the parameterized convex system via splitting convex inequalities into linear ones by using the Fenchel-Legendre conjugate. This approach yields that arbitrary bounded right-hand side perturbations of the convex system turn …

Preliminary Analysis Of An Agent-Based Model For A Tick-Borne Disease, Holly Gaff

Biological Sciences Faculty Publications

Ticks have a unique life history including a distinct set of life stages and a single blood meal per life stage. This makes tick-host interactions more complex from a mathematical perspective. In addition, any model of these interactions must involve a significant degree of stochasticity on the individual tick level. In an attempt to quantify these relationships, I have developed an individual-based model of the interactions between ticks and their hosts as well as the transmission of tick-borne disease between the two populations. The results from this model are compared with those from previously published differential equation based population models. …

An Extensible Mathematical Model Of Glucose Metabolism, Caleb L. Adams

Mathematics & Statistics Theses & Dissertations

The American Diabetes Association reports that diabetes is the fifth leading cause of death by disease in the United States. An estimated 23.6 million individuals, or seven percent of the population, have diabetes. Nearly one-third are unaware that they have the disease. The total of the direct and indirect medical costs associated with diabetes in 2007 was projected to be $174 billion, or approximately one out of every ten health care dollars.

One must understand the glucose regulatory system of the healthy body to understand diabetes. Blood glucose concentration returns to a constant level after eating and is maintained during …

A Least Squares Closure Approximation For Liquid Crystalline Polymers, Traci Ann Sievenpiper

Mathematics & Statistics Theses & Dissertations

An introduction to existing closure schemes for the Doi-Hess kinetic theory of liquid crystalline polymers is provided. A new closure scheme is devised based on a least squares fit of a linear combination of the Doi, Tsuji-Rey, Hinch-Leal I, and Hinch-Leal II closure schemes. The orientation tensor and rate-of-strain tensor are fit separately using data generated from the kinetic solution of the Smoluchowski equation. The known behavior of the kinetic solution and existing closure schemes at equilibrium is compared with that of the new closure scheme. The performance of the proposed closure scheme in simple shear flow for a variety …

Set-Based Corral Control In Stochastic Dynamical Systems: Making Almost Invariant Sets More Invariant, Eric Forgoston, Lora Billings, Philip Yecko, Ira Schwartz

Department of Mathematics Faculty Scholarship and Creative Works

We consider the problem of stochastic prediction and control in a time-dependent stochastic environment, such as the ocean, where escape from an almost invariant region occurs due to random fluctuations. We determine high-probability control-actuation sets by computing regions of uncertainty, almost invariant sets, and Lagrangian coherent structures. The combination of geometric and probabilistic methods allows us to design regions of control, which provide an increase in loitering time while minimizing the amount of control actuation. We show how the loitering time in almost invariant sets scales exponentially with respect to the control actuation, causing an exponential increase in loitering times …

A “Peak” At The Algorithm Behind “Peaklet Analysis” Software, Bruce Kessler

Mathematics Faculty Publications

In response to a problem posed by faculty at the Applied Physics Institute at Western Kentucky University, the speaker has developed an algorithm for providing an automated analysis of spectrum data for the purpose of determining the elemental composition of the item generating the data. A full, non-provisional patent application has been filed on the idea, and a full marketing campaign has started to license software implementing the algorithm. This presentation will give a brief explanation of the mathematics in use in the algorithm, and will give some examples of the software in action.

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 Of Bacillus Spores: Inactivation And Outgrowth, Alexis X. Hurst

Theses and Dissertations

This research models and analyzes the thermochemical damage produced in Bacillus spores by short, high-temperature exposures as well the repair process within damaged Bacillus spores. Thermochemical damage in spores is significantly due to reaction with water, hydrolysis reactions. Applying heat to the spore causes absorbed and chemically bound water molecules become mobile within the spore. These mobile water molecules react by hydrolysis reactions to degrade DNA and enzyme molecules in the spore. In order to survive the thermal inactivation, the spore must repair the damaged DNA during spore germination. The DNA repair process, as well as other germination functions, is …

Quantitative Interpretation Of A Genetic Model Of Carcinogenesis Using Computer Simulations, Donghai Dai, Brandon Beck, Xiaofang Wang, Cory Howk, Yi Li

Mathematics and Statistics Faculty Publications

The genetic model of tumorigenesis by Vogelstein et al. (V theory) and the molecular definition of cancer hallmarks by Hanahan and Weinberg (W theory) represent two of the most comprehensive and systemic understandings of cancer. Here, we develop a mathematical model that quantitatively interprets these seminal cancer theories, starting from a set of equations describing the short life cycle of an individual cell in uterine epithelium during tissue regeneration. The process of malignant transformation of an individual cell is followed and the tissue (or tumor) is described as a composite of individual cells in order to quantitatively account for intra-tumor …

A Primer On Chaos And Fractals, Bruce Kessler

Mathematics Faculty Publications

This is a prelude to a performance of the play "Arcadia" exclusively for the science and math majors of Lipscomb University. One of the main characters of the story is a mathematical genius, and has realized the power and limitations of iterations in generating mathematical models and structures, although she is living in the early 1800's. This talk gives an introduction to the ideas of chaos theory, fractals, and randomness.

Maximal Sensitive Dependence And The Optimal Path To Epidemic Extinction, Eric Forgoston, Simone Bianco, Leah B. Shaw, Ira B. Schwartz

Department of Applied Mathematics and Statistics Faculty Scholarship and Creative Works

Extinction of an epidemic or a species is a rare event that occurs due to a large, rare stochastic fluctuation. Although the extinction process is dynamically unstable, it follows an optimal path that maximizes the probability of extinction. We show that the optimal path is also directly related to the finite-time Lyapunov exponents of the underlying dynamical system in that the optimal path displays maximum sensitivity to initial conditions. We consider several stochastic epidemic models, and examine the extinction process in a dynamical systems framework. Using the dynamics of the finite-time Lyapunov exponents as a constructive tool, we demonstrate that …

Complete Characterizations Of Local Weak Sharp Minima With Applications To Semi-Infinite Optimization And Complementarity, Boris S. Mordukhovich, Naihua Xiu, Jinchuan Zhou

Mathematics Research Reports

In this paper we identify a favorable class of nonsmooth functions for which local weak sharp minima can be completely characterized in terms of normal cones and subdifferentials, or tangent cones and subderivatives, or their mixture in finite-dimensional spaces. The results obtained not only significantly extend previous ones in the literature, but also allow us to provide new types of criteria for local weak sharpness. Applications of the developed theory are given to semi-infinite programming and to semi-infinite complementarity problems.

Global Well-Posedness And Asymptotic Behavior Of A Class Of Initial-Boundary-Value Problems Of The Kdv Equation On A Finite Domain, Ivonne Rivas, Muhammad Usman, Bingyu Zhang

Mathematics Faculty Publications

In this paper, we study a class of initial boundary value problem (IBVP) of the Korteweg- de Vries equation posed on a ?nite interval with nonhomogeneous boundary conditions. The IBVP is known to be locally well-posed, but its global L2 a priori estimate is not available and therefore it is not clear whether its solutions exist globally or blow up in finite time. It is shown in this paper that the solutions exist globally as long as their initial value and the associated boundary data are small, and moreover, those solutions decay exponentially if their boundary data decay exponentially.

Projective-Planar Graphs With No K3,4-Minor, John Maharry, Dan Slilaty

Mathematics and Statistics Faculty Publications

An exact structure is described to classify the projective‐planar graphs that do not contain a K3, 4‐minor.

Sharp Asymptotics Of The Lp Approximation Error For Interpolation On Block Partitions, Yuliya Babenko, Tatyana Leskevich, Jean-Marie Mirebeau

Faculty Articles

Adaptive approximation (or interpolation) takes into account local variations in the behavior of the given function, adjusts the approximant depending on it, and hence yields the smaller error of approximation. The question of constructing optimal approximating spline for each function proved to be very hard. In fact, no polynomial time algorithm of adaptive spline approximation can be designed and no exact formula for the optimal error of approximation can be given. Therefore, the next natural question would be to study the asymptotic behavior of the error and construct asymptotically optimal sequences of partitions. In this paper we provide sharp asymptotic …