Numerical Analysis and Computation Commons^™

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

A mathematical model for the evolution of pulsed laser-irradiated, molten metallic films has been developed using the lubrication theory. The heat transfer problem that incorporates the absorbed heat from a single laser beam or the interfering laser beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the reflectivity, the peak laser beam …

Random Walks With Elastic And Reflective Lower Boundaries, Lucas Clay Devore

Masters Theses & Specialist Projects

No abstract provided.

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

A mathematical model for the evolution of pulsed laser-irradiated, molten metallic films has been developed using the lubrication theory. The heat transfer problem that incorporates the absorbed heat from a single laser beam or the interfering laser beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the reflectivity, the peak laser beam …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

A mathematical model for the evolution of pulsed laser-irradiated, molten metallic films has been developed using the lubrication theory. The heat transfer problem that incorporates the absorbed heat from a single laser beam or the interfering laser beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the reflectivity, the peak laser beam …

Higher Order Difference Schemes For Heat Equation, Jianzhong Wang

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we construct the explicit difference schemes for the heat equation with arbitrary high orders. We also show the validity of the new schemes by numerical simulations.

Targeted Maximum Likelihood Estimation: A Gentle Introduction, Susan Gruber, Mark J. Van Der Laan

U.C. Berkeley Division of Biostatistics Working Paper Series

This paper provides a concise introduction to targeted maximum likelihood estimation (TMLE) of causal effect parameters. The interested analyst should gain sufficient understanding of TMLE from this introductory tutorial to be able to apply the method in practice. A program written in R is provided. This program implements a basic version of TMLE that can be used to estimate the effect of a binary point treatment on a continuous or binary outcome.

Waves In Inhomogeneous Solids, Arkadi Berezovski, Mihhail Berezovski, Juri Engelbrecht

Publications

The paper aims at presenting a numerical technique used in simulating the propagation of waves in inhomogeneous elastic solids. The basic governing equations are solved by means of a finite-volume scheme that is faithful, accurate, and conservative. Furthermore, this scheme is compatible with thermodynamics through the identification of the notions of numerical fluxes (a notion from numerics) and of excess quantities (a notion from irreversible thermodynamics). A selection of one-dimensional wave propagation problems is presented, the simulation of which exploits the designed numerical scheme. This selection of exemplary problems includes (i) waves in periodic media for weakly nonlinear waves with …

Shrinkage Estimation Of Expression Fold Change As An Alternative To Testing Hypotheses Of Equivalent Expression, Zahra Montazeri, Corey M. Yanofsky, David R. Bickel

COBRA Preprint Series

Research on analyzing microarray data has focused on the problem of identifying differentially expressed genes to the neglect of the problem of how to integrate evidence that a gene is differentially expressed with information on the extent of its differential expression. Consequently, researchers currently prioritize genes for further study either on the basis of volcano plots or, more commonly, according to simple estimates of the fold change after filtering the genes with an arbitrary statistical significance threshold. While the subjective and informal nature of the former practice precludes quantification of its reliability, the latter practice is equivalent to using a …

Variational Analysis In Semi-Infinite And Infinite Programming, Ii: Necessary Optimality Conditions, M J. Cánovas, M A. Lopez, Boris S. Mordukhovich, J Parra

Mathematics Research Reports

This paper concerns applications of advanced techniques of variational analysis and generalized differentiation to problems of semi-infinite and infinite programming with feasible solution sets defined by parameterized systems of infinitely many linear inequalities of the type intensively studied in the preceding development [5] from our viewpoint of robust Lipschitzian stability. We present meaningful interpretations and practical examples of such models. The main results establish necessary optimality conditions for a broad class of semi-infinite and infinite programs, where objectives are generally described by nonsmooth and nonconvex functions on Banach spaces and where infinite constraint inequality systems are indexed by arbitrary sets. …

Caching And Visualizing Statistical Analyses, Roger D. Peng, Duncan Temple Lang

Johns Hopkins University, Dept. of Biostatistics Working Papers

We present the cacher and CodeDepends packages for R, which provide tools for (1) caching and analyzing the code for statistical analyses and (2) distributing these analyses to others in an efficient manner over the web. The cacher package takes objects created by evaluating R expressions and stores them in key-value databases. These databases of cached objects can subsequently be assembled into “cache packages” for distribution over the web. The cacher package also provides tools to help readers examine the data and code in a statistical analysis and reproduce, modify, or improve upon the results. In addition, readers can easily …

Variational Analysis In Semi-Infinite And Infinite Programming, I: Stability Of Linear Inequality Systems Of Feasible Solutions, M J. Cánovas, M A. Lopez, Boris S. Mordukhovich, J Parra

Mathematics Research Reports

This paper concerns applications of advanced techniques of variational analysis and generalized differentiation to parametric problems of semi-infinite and infinite programming, where decision variables run over finite-dimensional and infinite-dimensional spaces, respectively. Part I is primarily devoted to the study of robust Lipschitzian stability of feasible solutions maps for such problems described by parameterized systems of infinitely many linear inequalities in Banach spaces of decision variables indexed by an arbitrary set T. The parameter space of admissible perturbations under consideration is formed by all bounded functions on T equipped with the standard supremum norm. Unless the index set is finite, this …

Variational Iteration Method For Solving Telegraph Equations, Syed T. Mohyud-Din, Muhammad A. Noor, Khalida I. Noor

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we apply the variational iteration method (VIM) for solving telegraph equations, which arise in the propagation of electrical signals along a telegraph line. The suggested algorithm is more efficient and easier to handle as compare to the decomposition method. Numerical results show the efficiency and accuracy of the proposed VIM.

Circular Nonlinear Subdivision Schemes For Curve Design, Jian-Ao Lian, Yonghui Wang, Yonggao Yang

Applications and Applied Mathematics: An International Journal (AAM)

Two new families of nonlinear 3-point subdivision schemes for curve design are introduced. The first family is ternary interpolatory and the second family is binary approximation. All these new schemes are circular-invariant, meaning that new vertices are generated from local circles formed by three consecutive old vertices. As consequences of the nonlinear schemes, two new families of linear subdivision schemes for curve design are established. The 3-point linear binary schemes, which are corner-cutting depending on the choices of the tension parameter, are natural extensions of the Lane-Riesenfeld schemes. The four families of both nonlinear and linear subdivision schemes are implemented …

Generalized Probabilistic Bowling Distributions, Jennifer Lynn Hohn

Masters Theses & Specialist Projects

Have you ever wondered if you are better than the average bowler? If so, there are a variety of ways to compute the average score of a bowling game, including methods that account for a bowler’s skill level. In this thesis, we discuss several different ways to generate bowling scores randomly. For each distribution, we give results for the expected value and standard deviation of each frame's score, the expected value of the game’s final score, and the correlation coefficient between the score of the first and second roll of a single frame. Furthermore, we shall generalize the results in …

Tight Lower Bound For The Sparse Travelling Salesman Problem, Fredrick Mtenzi

Conference papers

The Sparse Travelling Salesman Problem (Sparse TSP) which is a variant of the classical Travelling Salesman Problem (TSP) is the problem of finding the shortest route of the salesman when visiting cities in a region making sure that each city is visited at least once and returning home at the end. In the Sparse TSP, the distance between cities may not obey the triangle inequality; this makes the use of algorithms and formulations designed for the TSP to require modifications in order to produce near-optimal results. A lower bound for optmisation problems gives us the quality guarantee of the near-optimal …

Noisy Signal Recovery Via Iterative Reweighted L1-Minimization, Deanna Needell

CMC Faculty Publications and Research

Compressed sensing has shown that it is possible to reconstruct sparse high dimensional signals from few linear measurements. In many cases, the solution can be obtained by solving an L1-minimization problem, and this method is accurate even in the presence of noise. Recent a modified version of this method, reweighted L1-minimization, has been suggested. Although no provable results have yet been attained, empirical studies have suggested the reweighted version outperforms the standard method. Here we analyze the reweighted L1-minimization method in the noisy case, and provide provable results showing an improvement in the error bound over the standard bounds.

Fractional Calculus: Definitions And Applications, Joseph M. Kimeu

Masters Theses & Specialist Projects

No abstract provided.

Determining The Orbit Locations Of Turkish Airborne Early Warning And Control Aircraft Over The Turkish Air Space, Nebi Sarikaya

Theses and Dissertations

The technology improvement affects the military needs of individual countries. The new doctrine of defense for many countries emphasizes detecting threats as far away as you can from your homeland. Today, the military uses both ground RADAR and Airborne Early Warning and Control (AEW&C) Aircraft. AEW&C aircraft has become vital to detect low altitude threats that a ground RADAR cannot detect because of obstacles on the earth. Turkey has ordered four AEW&C aircraft for her air defense system because of the lack of complete coverage by ground RADAR. This research provides optimal orbit locations that can be updated according to …

Efficient Evaluation Of Ranking Procedures When The Number Of Units Is Large With Application To Snp Identification, Thomas A. Louis, Ingo Ruczinski

Johns Hopkins University, Dept. of Biostatistics Working Papers

Simulation-based assessment is a popular and frequently necessary approach to evaluation of statistical procedures. Sometimes overlooked is the ability to take advantage of underlying mathematical relations and we focus on this aspect. We show how to take advantage of large-sample theory when conducting a simulation using the analysis of genomic data as a motivating example. The approach uses convergence results to provide an approximation to smaller-sample results, results that are available only by simulation. We consider evaluating and comparing a variety of ranking-based methods for identifying the most highly associated SNPs in a genome-wide association study, derive integral equation representations …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Agegnehu Atena, Mikhail Khenner

Mathematics Faculty Publications

In this paper the lubrication-type dynamical model is developed of a molten, pulsed laser-irradiated metallic film. The heat transfer problem that incorporates the absorbed heat from a single beam or interfering beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the peak laser beam intensity, the film optical thickness, the Biot and …

Thermocapillary Effects In Driven Dewetting And Self-Assembly Of Pulsed Laser-Irradiated Metallic Films, Agegnehu Atena, Mikhail Khenner

Mathematics Faculty Publications

In this paper the lubrication-type dynamical model is developed of a molten, pulsed laser-irradiated metallic film. The heat transfer problem that incorporates the absorbed heat from a single beam or interfering beams is solved analytically. Using this temperature field, we derive the 3D long-wave evolution PDE for the film height. To get insights into dynamics of dewetting, we study the 2D version of the evolution equation by means of a linear stability analysis and by numerical simulations. The stabilizing and destabilizing effects of various system parameters, such as the peak laser beam intensity, the film optical thickness, the Biot and …

The Beauty Of Mathematics And The Mathematics Of Beauty: Continued Fractions And The Golden Ratio, Jessica Tush

Inquiry: The University of Arkansas Undergraduate Research Journal

This project begins with a look at the history of simple continued fractions and how we have arrived where we are today. We then move through a study of simple continued fractions, beginning first with rational numbers and moving to irrational numbers. Continuing further in the pursuit of joining mathematics and art, we define the specific continued fraction that gives rise to the Fibonacci sequence and the Golden Ratio~ (phi, pronounced 'Jai"). These two notions form a direct link to art and the properties that we hope to examine. I have taken an analytic approach to showing that the Golden …

Polynomial Extension Operators. Part Ii, Leszek Demkowicz, Jay Gopalakrishnan, Joachim Schöberl

Mathematics and Statistics Faculty Publications and Presentations

Consider the tangential trace of a vector polynomial on the surface of a tetrahedron. We construct an extension operator that extends such a trace function into a polynomial on the tetrahedron. This operator can be continuously extended to the trace space of H(curl ). Furthermore, it satisfies a commutativity property with an extension operator we constructed in Part I of this series. Such extensions are a fundamental ingredient of high order finite element analysis.

A Covert Encryption Method For Applications In Electronic Data Interchange, Jonathan Blackledge, Dmitry Dubovitskiy

Articles

A principal weakness of all encryption systems is that the output data can be ‘seen’ to be encrypted. In other words, encrypted data provides a ‘flag’ on the potential value of the information that has been encrypted. In this paper, we provide a new approach to ‘hiding’ encrypted data in a digital image.

In conventional (symmetric) encryption, the plaintext is usually represented as a binary stream and encrypted using an XOR type operation with a binary cipher. The algorithm used is ideally designed to: (i) generate a maximum entropy cipher so that there is no bias with regard to any …

Singular Superposition/Boundary Element Method For Reconstruction Of Multi-Dimensional Heat Flux Distributions With Application To Film Cooling Holes, Mahmood Silieti, Eduardo Divo, Alain J. Kassab

Publications

A hybrid singularity superposition/boundary element-based inverse problem method for the reconstruction of multi-dimensional heat flux distributions is developed. Cauchy conditions are imposed at exposed surfaces that are readily reached for measurements while convective boundary conditions are unknown at surfaces that are not amenable to measurements such as the walls of the cooling holes. The purpose of the inverse analysis is to determine the heat flux distribution along cooling hole surfaces. This is accomplished in an iterative process by distributing a set of singularities (sinks) inside the physical boundaries of the cooling hole (usually along cooling hole centerline) with a given …

An H-Adaptive Finite-Element Technique For Constructing 3d Wind Fields, Darrell Pepper, Xiuling Wang

Mechanical Engineering Faculty Research

An h-adaptive, mass-consistent finite-element model (FEM) has been developed for constructing 3D wind fields over irregular terrain utilizing sparse meteorological tower data. The element size in the computational domain is dynamically controlled by an a posteriori error estimator based on the L2 norm. In the h-adaptive FEM algorithm, large element sizes are typically associated with smooth flow regions and small errors; small element sizes are attributed to fast-changing flow regions and large errors. The adaptive procedure employed in this model uses mesh refinement–unrefinement to satisfy error criteria. Results are presented for wind fields using sparse data obtained from two regions …