Numerical Analysis and Computation Commons^™

Integrability, Recursion Operators And Soliton Interactions, Boyka Aneva, Georgi Grahovski, Rossen Ivanov, Dimitar Mladenov

Book chapter/book

This volume contains selected papers based on the talks,presentedat the Conference Integrability, Recursion Operators and Soliton Interactions, held in Sofia, Bulgaria (29-31 August 2012) at the Institute for Nuclear Research and Nuclear Energy of the Bulgarian Academy of Sciences. Included are also invited papers presenting new research developments in the thematic area. The Conference was dedicated to the 65-th birthday of our esteemed colleague and friend Vladimir Gerdjikov. The event brought together more than 30 scientists, from 6 European countries to celebrate Vladimir's scientific achievements. All participants enjoyed a variety of excellent talks in a friendly and stimulating atmosphere. …

A Posteriori Error Estimates For Surface Finite Element Methods, Fernando F. Camacho

Theses and Dissertations--Mathematics

Problems involving the solution of partial differential equations over surfaces appear in many engineering and scientific applications. Some of those applications include crystal growth, fluid mechanics and computer graphics. Many times analytic solutions to such problems are not available. Numerical algorithms, such as Finite Element Methods, are used in practice to find approximate solutions in those cases.

In this work we present L2 and pointwise a posteriori error estimates for Adaptive Surface Finite Elements solving the Laplace-Beltrami equation −△_Γ u = f . The two sources of errors for Surface Finite Elements are a Galerkin error, and a …

Dispersive Waves In Microstructured Solids, A. Berezovski, J. Engelbrecht, A. Salupere, K. Tamm, T. Peets, Mihhail Berezovski

Publications

The wave motion in micromorphic microstructured solids is studied. The mathematical model is based on ideas of Mindlin and governing equations are derived by making use of the Euler–Lagrange formalism. The same result is obtained by means of the internal variables approach. Actually such a model describes internal fields in microstructured solids under external loading and the interaction of these fields results in various physical effects. The emphasis of the paper is on dispersion analysis and wave profiles generated by initial or boundary conditions in a one-dimensional case.

An Exponential Matrix Method For Numerical Solutions Of Hantavirus Infection Model, Şuayip Yüzbaşi, Mehmet Sezer

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, a new matrix method based on exponential polynomials and collocation points is proposed to obtain approximate solutions of Hantavirus infection model corresponding to a class of systems of nonlinear ordinary differential equations. The method converts the model problem into a system of nonlinear algebraic equations by means of the matrix operations and the collocation points. The reliability and efficiency of the proposed scheme is demonstrated by the numerical applications and all numerical computations have been made by using a computer program written in Maple.

Projected Surface Finite Elements For Elliptic Equations, Necibe Tuncer

Applications and Applied Mathematics: An International Journal (AAM)

In this article, we define a new finite element method for numerically approximating solutions of elliptic partial differential equations defined on “arbitrary” smooth surfaces S in RN+1. By “arbitrary” smooth surfaces, we mean surfaces that can be implicitly represented as level sets of smooth functions. The key idea is to first approximate the surface S by a polyhedral surface Sh, which is a union of planar triangles whose vertices lie on S; then to project Sh onto S. With this method, we can also approximate the eigenvalues and eigenfunctions of th Laplace-Beltrami operator on these “arbitrary” surfaces.

Influence Of Microstructure On Thermoelastic Wave Propagation, Arkadi Berezovski, Mihhail Berezovski

Publications

Numerical simulations of the thermoelastic response of a microstructured material on a thermal loading are performed in the one-dimensional setting to examine the influence of temperature gradient effects at the microstructure level predicted by the thermoelastic description of microstructured solids (Berezovski et al. in J. Therm. Stress. 34:413–430, 2011). The system of equations consisting of a hyperbolic equation of motion, a parabolic macroscopic heat conduction equation, and a hyperbolic evolution equation for the microtemperature is solved by a finite-volume numerical scheme. Effects of microtemperature gradients exhibit themselves on the macrolevel due to the coupling of equations of the macromotion …

Two-Dimensional Hydrodynamic Modeling Of Two-Phase Flow For Understanding Geyser Phenomena In Urban Stormwater System, Zhiyu S. Shao

Theses and Dissertations--Civil Engineering

During intense rain events a stormwater system can fill rapidly and undergo a transition from open channel flow to pressurized flow. This transition can create large discrete pockets of trapped air in the system. These pockets are pressurized in the horizontal reaches of the system and then are released through vertical vents. In extreme cases, the transition and release of air pockets can create a geyser feature.

The current models are inadequate for simulating mixed flows with complicated air-water interactions, such as geysers. Additionally, the simulation of air escaping in the vertical dropshaft is greatly simplified, or completely ignored, in …

Time-Stepping For Laser Ablation, Harihar Khanal, David Autrique, Vasilios Alexiades

Publications

Nanosecond laser ablation is a popular technique, applied in many areas of science and technology such as medicine, archaeology, chemistry, environmental and materials sciences. We outline a computational model for radiative and collisional processes occurring during ns-laser ablation, and compare the performance of various low and high order time-stepping algorithms.

Hydrodynamic Modeling Of Ns-Laser Ablation, David Autrique, Vasilios Alexiades, Harihar Khanal

Publications

Laser ablation is a versatile and widespread technique, applied in an increasing number of medical, industrial and analytical applications. A hydrodynamic multiphase model describing nanosecond-laser ablation (ns- LA) is outlined. The model accounts for target heating and mass removal mechanisms as well as plume expansion and plasma formation. A copper target is placed in an ambient environment consisting of helium and irradiated by a nanosecond-laser pulse. The effect of variable laser settings on the ablation process is explored in 1-D numerical simulations.

Numerical Studies For Solving Fractional Riccati Differential Equation, N. H. Sweilam, M. M. Khader, A. M. S. Mahdy

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, finite difference method (FDM) and Pade'-variational iteration method (Pade'- VIM) are successfully implemented for solving the nonlinear fractional Riccati differential equation. The fractional derivative is described in the Caputo sense. The existence and the uniqueness of the proposed problem are given. The resulting nonlinear system of algebraic equations from FDM is solved by using Newton iteration method; moreover the condition of convergence is verified. The convergence's domain of the solution is improved and enlarged by Pade'-VIM technique. The results obtained by using FDM is compared with Pade'-VIM. It should be noted that the Pade'-VIM is preferable because …

Introducing An Efficient Modification Of The Variational Iteration Method By Using Chebyshev Polynomials, M. M. Khader

Applications and Applied Mathematics: An International Journal (AAM)

In this article an efficient modification of the variational iteration method (VIM) is presented using Chebyshev polynomials. Special attention is given to study the convergence of the proposed method. The new modification is tested for some examples to demonstrate reliability and efficiency of the proposed method. A comparison of our numerical results those of the conventional numerical method, the fourth-order Runge-Kutta method (RK4) are given. The comparison shows that the solution using our modification is fast-convergent and is in excellent conformance with the exact solution. Finally, we conclude that the proposed method can be applied to a large class of …

A Duhamel Integral Based Approach To Identify An Unknown Radiation Term In A Heat Equation With Non-Linear Boundary Condition, R. Pourgholi, M. Abtahi, A. Saeedi

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, we consider the determination of an unknown radiation term in the nonlinear boundary condition of a linear heat equation from an overspecified condition. First we study the existence and uniqueness of the solution via an auxiliary problem. Then a numerical method consisting of zeroth-, first-, and second-order Tikhonov regularization method to the matrix form of Duhamel's principle for solving the inverse heat conduction problem (IHCP) using temperature data containing significant noise is presented. The stability and accuracy of the scheme presented is evaluated by comparison with the Singular Value Decomposition (SVD) method. Some numerical experiments confirm the …

Pointwise Schauder Estimates Of Parabolic Equations In Carnot Groups, Heather Arielle Griffin

Graduate Theses and Dissertations

Schauder estimates were a historical stepping stone for establishing uniqueness and smoothness of solutions for certain classes of partial differential equations. Since that time, they have remained an essential tool in the field. Roughly speaking, the estimates state that the Holder continuity of the coefficient functions and inhomogeneous term implies the Holder continuity of the solution and its derivatives. This document establishes pointwise Schauder estimates for second order parabolic equations where the traditional role of derivatives are played by vector fields generated by the first layer of the Lie algebra stratification for a Carnot group. The Schauder estimates are shown …

On The Stability Of A Microstructure Model, Mihhail Berezovski, Arkadi Berezovski

Publications

Abstract

The asymptotic stability of solutions of the Mindlin-type microstructure model for solids is analyzed in the paper. It is shown that short waves are asymptotically stable even in the case of a weakly non-convex free energy dependence on microdeformation.

Research highlights

The Mindlin-type microstructure model cannot describe properly short wave propagation in laminates. A modified Mindlin-type microstructure model with weakly non-convex free energy resolves this discrepancy. It is shown that the improved model with weakly non-convex free energy is asymptotically stable for short waves.

Controlling Nanoparticles Formation In Molten Metallic Bilayers By Pulsed-Laser Interference Heating, Mikhail Khenner, Sagar Yadavali, Ramki Kalyanaraman

Mathematics Faculty Publications

The impacts of the two-beam interference heating on the number of core-shell and embedded nanoparticles and on nanostructure coarsening are studied numerically based on the non-linear dynamical model for dewetting of the pulsed-laser irradiated, thin (< 20 nm) metallic bilayers. The model incorporates thermocapillary forces and disjoining pressures, and assumes dewetting from the optically transparent substrate atop of the reflective support layer, which results in the complicated dependence of light reflectivity and absorption on the thicknesses of the layers. Stabilizing thermocapillary effect is due to the local thickness-dependent, steady- state temperature profile in the liquid, which is derived based on the mean substrate temperature estimated from the elaborate thermal model of transient heating and melting/freezing. Linear stability analysis of the model equations set for Ag/Co bilayer predicts the dewetting length scales in the qualitative agreement with experiment.

Singular Solutions Of Coss-Coupled Epdiff Equations: Waltzing Peakons And Compacton Pairs, Colin Cotter, Darryl Holm, Rossen Ivanov, James Percival

Conference papers

We introduce EPDiff equations as Euler-Poincare´ equations related to Lagrangian provided by a metric, invariant under the Lie Group Diff(Rⁿ). Then we proceed with a particular form of EPDiff equations, a cross coupled two-component system of Camassa-Holm type. The system has a new type of peakon solutions, 'waltzing' peakons and compacton pairs.

Wave Propagation And Dispersion In Microstructured Solids, Arkadi Berezovski, Juri Engelbrecht, Mihhail Berezovski

Publications

A series of numerical simulations is carried on in order to understand the accuracy of dispersive wave models for microstructured solids. The computations are performed by means of the finite-volume numerical scheme, which belongs to the class of wave-propagation algorithms. The dispersion effects are analyzed in materials with different internal structures: microstructure described by micromorphic theory, regular laminates, laminates with substructures, etc., for a large range of material parameters and wavelengths.

Formation Of Organized Nanostructures From Unstable Bilayers Of Thin Metallic Liquids, Mikhail Khenner, Sagar Yadavali, Ramki Kalyanaraman

Mathematics Faculty Publications

Dewetting of pulsed-laser irradiated, thin (< 20 nm), optically reflective metallic bilayers on an optically transparent substrate with a reflective support layer is studied within the lubrication equations model. A steady-state bilayer film thickness (h) dependent temperature profile is derived based on the mean substrate temperature estimated from the elaborate thermal model of transient heating and melting/freezing. Large thermocapillary forces are observed along the plane of the liquid-liquid and liquid-gas interfaces due to this h-dependent temperature, which, in turn, is strongly influenced by the h-dependent laser light reflection and absorption. Consequently the dewetting is a result of the competition between thermocapillary and intermolecular forces. A linear analysis of the dewetting length scales established that the non-isothermal calculations better predict the experimental results as compared to the isothermal case within the bounding Hamaker coefficients. Subsequently, a computational non-linear dynamics study of the dewetting pathway was performed for Ag/Co and Co/Ag bilayer systems to predict the morphology evolution. We found that the systems evolve towards formation of different morphologies, including core-shell, embedded, or stacked nanostructure morphologies.

A Dynamical Study Of The Evolution Of Pressure Waves Propagating Through A Semi-Infinite Region Of Homogeneous Gas Combustion Subject To A Time-Harmonic Signal At The Boundary, John Eslick

LSU New Orleans Theses and Dissertations

In this dissertation, the evolution of a pressure wave driven by a harmonic signal on the boundary during gas combustion is studied. The problem is modeled by a nonlinear, hyperbolic partial differential equation. Steady-state behavior is investigated using the perturbation method to ensure that enough time has passed for any transient effects to have dissipated. The zeroth, first and second-order perturbation solutions are obtained and their moduli are plotted against frequency. It is seen that the first and second-order corrections have unique maxima that shift to the right as the frequency decreases and to the left as the frequency increases. …

Two-Scale Microstructure Dynamics, Arkadi Berezovski, Mihhail Berezovski, Juri Engelbrecht

Publications

Wave propagation in materials with embedded two different microstructures is considered. Each microstructure is characterized by its own length scale. The dual internal variables approach is adopted yielding in a Mindlin-type model including both microstructures. Equations of motion for microstructures are coupled with the balance of linear momentum for the macromotion, but not coupled with each other. Corresponding dispersion curves are provided and scale separation is pointed out.

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 …

Waves In Microstructured Solids: A Unified Viewpoint Of Modelling, Arkadi Berezovski, Juri Engelbrecht, Mihhail Berezovski

Publications

The basic ideas for describing the dispersive wave motion in microstructured solids are discussed in the one-dimensional setting because then the differences between various microstructure models are clearly visible. An overview of models demonstrates a variety of approaches, but the consistent structure of the theory is best considered from the unified viewpoint of internal variables. It is shown that the unification of microstructure models can be achieved using the concept of dual internal variables.

On The Stability Of A Microstructure Model, Mihhail Berezovski, Arkadi Berezovski

Publications

The asymptotic stability of solutions of the Mindlin-type microstructure model for solids is analyzed in the paper. It is shown that short waves are asymptotically stable even in the case of a weakly non-convex free energy dependence on microdeformation.

Modeling Diverse Physics Of Nanoparticle Self-Assembly In Pulsed Laser-Irradiated Metallic Films, Mikhail Khenner

Mathematics Faculty Publications

Presents physics behind dewetting of thin liquid films and mathematical/computational modeling tools (Educational/Research presentation for senior physics majors).

Stability Of A Strongly Anisotropic Thin Epitaxial Film In A Wetting Interaction With Elastic Substrate, Mikhail Khenner, Wondimu T. Tekalign, Margo S. Levine

Mathematics Faculty Publications

The linear dispersion relation for longwave surface perturbations, as derived by Levine et al. Phys. Rev. B 75, 205312 (2007) is extended to include a smooth surface energy anisotropy function with a variable anisotropy strength (from weak to strong, such that sharp corners and slightly curved facets occur on the corresponding Wulff shape). Through detailed parametric studies it is shown that a combination of a wetting interaction and strong anisotropy, and even a wetting interaction alone results in complicated linear stability characteristics of strained and unstrained films.

Stability Of A Strongly Anisotropic Thin Epitaxial Film In A Wetting Interaction With Elastic Substrate, Mikhail Khenner, Wondimu Tekalign, Margo Levine

Mathematics Faculty Publications

The linear dispersion relation for longwave surface perturbations, as derived by Levine et al. Phys. Rev. B 75, 205312 (2007) is extended to include a smooth surface energy anisotropy function with a variable anisotropy strength (from weak to strong, such that sharp corners and slightly curved facets occur on the corresponding Wulff shape). Through detailed parametric studies it is shown that a combination of a wetting interaction and strong anisotropy, and even a wetting interaction alone results in complicated linear stability characteristics of strained and unstrained films.

Dispersive Wave Equations For Solids With Microstructure, A. Berezovski, Juri Engelbrecht, Mihhail Berezovski

Publications

The dispersive wave motion in solids with microstructure is considered in the one-dimensional setting in order to understand better the mechanism of dispersion. It is shown that the variety of dispersive wave propagation models derived by homogenization, continualisation, and generalization of continuum mechanics can be unified in the framework of dual internal variables theory.

Application Of Differential Transform Method To The Generalized Burgers–Huxley Equation, J. Biazar, F. Mohammadi

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, the differential transform method (DTM) will be applied to the generalized Burgers-Huxley equation, and some special cases of the equation, say, Huxley equation and Fitzhugh-Nagoma equation. The DTM produces an approximate solution for the equation, with few and easy computations. Numerical comparison between differential transform method, Adomian decomposition method and Variational iteration method for Burgers-Huxley, Huxley equation and Fitzhugh-Nagoma equation reveal that differential transform method is simple, accurate and efficient.

Approximating Solutions For Ginzburg – Landau Equation By Hpm And Adm, J. Biazar, M. Partovi, Z. Ayati

Applications and Applied Mathematics: An International Journal (AAM)

In this paper, an analytical approximation to the solution of Ginzburg-Landauis discussed. A Homotopy perturbation method introduced by He is employed to derive the analytic approximation solution and results compared with those of the Adomian decomposition method. Two examples are presented to show the capability of the methods. The results reveal that the methods are almost equally effective and promising.

Application Of Homotopy Perturbation Method To Biological Population Model, Pradip Roul

Applications and Applied Mathematics: An International Journal (AAM)

In this article, a well-known analytical approximation method, so-called the Homotopy perturbation method (HPM) is adopted for solving the nonlinear partial differential equations arising in the spatial diffusion of biological populations. The resulting solutions are compared with those of the existing solutions obtained by employing the Adomian’s decomposition method. The comparison reveals that our approximate solutions are in very good agreement with the solutions by Adomian’s method. Moreover, the results show that the proposed method is a more reliable, efficient and convenient one for solving the non-linear differential equations.