Numerical Analysis and Computation Commons^™

An Rbf Interpolation Blending Scheme For Effective Shock-Capturing, M. Harris, Eduardo Divo, Alain J. Kassab

Publications

In recent years significant focus has been given to the study of Radial basis functions (RBF), especially in their use on solving partial differential equations (PDE). RBF have an impressive capability of inter- polating scattered data, even when this data presents localized discontinuities. However, for infinitely smooth RBF such as the Multiquadrics, inverse Multiquadrics, and Gaussian, the shape parameter must be chosen properly to obtain accurate approximations while avoiding ill-conditioning of the interpolating matrices. The optimum shape parameter can vary significantly depending on the field, particularly in locations of steep gradients, shocks, or discontinuities. Typically, the shape parameter is chosen …

Application Of An Rbf Blending Interpolation Method To Problems With Shocks, Michael Harris, Eduardo Divo, Alain J. Kassab

Publications

Radial basis functions (RBF) have become an area of research in recent years, especially in the use of solving partial differential equations (PDE). Radial basis functions have an impressive capability in interpolating scattered data, even for data with discontinuities. Although, for infinitely smooth radial basis functions such as the multi-quadrics and inverse multi-quadrics, the shape parameter must be chosen properly to obtain accurate approximations while avoiding ill-conditioning of the interpolating matrices. The optimum shape parameter can vary depending on the field, such as in locations of sharp gradients or shocks. Typically, the shape parameter is chosen to maintain a high …

Numerical Simulation Of Acoustic Emission During Crack Growth In 3-Point Bending Test, Mihhail Berezovski, Arkadi Berezovski

Publications

Numerical simulation of acoustic emission by crack propagation in 3-point bending tests is performed to investigate how the interaction of elastic waves generates a detectable signal. It is shown that the use of a kinetic relation for the crack tip velocity combined with a simple crack growth criterion provides the formation of waveforms similar to those observed in experiments.

Multi-Scale Cardiovascular Flow Analysis By An Integrated Meshless-Lumped Parameter Model, Leonardo A. Bueno, Eduardo A. Divo, Alain J. Kassab

Publications

A computational tool that integrates a Radial basis function (RBF)-based Meshless solver with a Lumped Parameter model (LPM) is developed to analyze the multi-scale and multi-physics interaction between the cardiovascular flow hemodynamics, the cardiac function, and the peripheral circulation. The Meshless solver is based on localized RBF collocations at scattered data points which allows for automation of the model generation via CAD integration. The time-accurate incompressible flow hemodynamics are addressed via a pressure-velocity correction scheme where the ensuing Poisson equations are accurately and efficiently solved at each time step by a Dual-Reciprocity Boundary Element method (DRBEM) formulation that takes advantage …

Long And Short-Range Air Navigation On Spherical Earth, Nihad E. Daidzic

International Journal of Aviation, Aeronautics, and Aerospace

Global range air navigation implies non-stop flight between any two airports on Earth. Such effort would require airplanes with the operational air range of at least 12,500 NM which is about 40-60% longer than anything existing in commercial air transport today. Air transportation economy requires flying shortest distance, which in the case of spherical Earth are Orthodrome arcs. Rhumb-line navigation has little practical use in long-range flights, but has been presented for historical reasons and for comparison. Database of about 50 major international airports from every corner of the world has been designed and used in testing and route validation. …

Breakdown Of Itcz-Like Pv Patterns, Ajay Raghavendra, Thomas A. Guinn

Beyond: Undergraduate Research Journal

The Inter-Tropical Convergence Zone (ITCZ) is a zonal belt of intense convection, responsible for the genesis of over 80% of all tropical cyclones. This region of intense diabatic heating and shear results in a maximum of Ertel's potential vorticity (PV) meeting Rayleigh's necessary condition for barotropic instability. A fundamental issue is understanding the necessary precursor events leading to the breakdown of the ITCZ and subsequent formation of tropical cyclones. Our research examines the non-linear PV dynamics of the breakdown of both finite-length and infinite-length vorticity strips of varying widths and shapes, simulating the ITCZ found near the tropical eastern Pacific …

Thermoelastic Waves In Microstructured Solids, Arkadi Berezovski, Mihhail Berezovski

Publications

Thermoelastic wave propagation suggests a coupling between elastic deformation and heat conduction in a body. Microstructure of the body influences the both processes. Since energy is conserved in elastic deformation and heat conduction is always dissipative, the generalization of classical elasticity theory and classical heat conduction is performed differently. It is shown in the paper that a hyperbolic evolution equation for microtemperature can be obtained in the framework of the dual internal variables approach keeping the parabolic equation for the macrotemperature. The microtemperature is considered as a macrotemperature fluctuation. Numerical simulations demonstrate the formation and propagation of thermoelastic waves in …

Signal Flow Graph Approach To Efficient Dst I-Iv Algorithms, Sirani M. Perera

Publications

In this paper, fast and efficient discrete sine transformation (DST) algorithms are presented based on the factorization of sparse, scaled orthogonal, rotation, rotation-reflection, and butterfly matrices. These algorithms are completely recursive and solely based on DST I-IV. The presented algorithms have low arithmetic cost compared to the known fast DST algorithms. Furthermore, the language of signal flow graph representation of digital structures is used to describe these efficient and recursive DST algorithms having (n�1) points signal flow graph for DST-I and n points signal flow graphs for DST II-IV.

Pattern Formation Of Elastic Waves And Energy Localization Due To Elastic Gratings, A. Berezovski, J. Engelbrecht, Mihhail Berezovski

Publications

Elastic wave propagation through diffraction gratings is studied numerically in the plane strain setting. The interaction of the waves with periodically ordered elastic inclusions leads to a self-imaging Talbot effect for the wavelength equal or close to the grating size. The energy localization is observed at the vicinity of inclusions in the case of elastic gratings. Such a localization is absent in the case of rigid gratings.

Global Optimized Isothermal And Nonlinear Models Of Earth’S Standard Atmosphere, Nihad E. Daidzic, Ph.D.,

International Journal of Aviation, Aeronautics, and Aerospace

Both, a global isothermal temperature model and a nonlinear quadratic temperature model of the ISA was developed and presented here. Constrained optimization techniques in conjunction with the least-square-root approximations were used to design best-fit isothermal models for ISA pressure and density changes up to 47 geopotential km for NLPAM, and 86 orthometric km for ISOAM respectively. The mass of the dry atmosphere and the relevant fractional-mass scale heights have been computed utilizing the very accurate eight-point Gauss-Legendre numerical quadrature for both ISOAM and NLPAM. Both, the ISOAM and the NLPAM represent viable alternatives to ISA in many practical applications and …

Efficient General Computational Method For Estimation Of Standard Atmosphere Parameters, Nihad E. Daidzic Ph.D., Sc.D.

International Journal of Aviation, Aeronautics, and Aerospace

Knowledge of standard air temperature, pressure, density, speed of sound, and viscosity as a function of altitude is essential information in aircraft design, performance testing, pressure altimeter calibration, and several other aeronautical engineering and aviation science applications. A new efficient computational method for rapid calculations of standard atmospheric parameters up to 86 orthometric km is presented. Additionally, mass and weight of each standard atmospheric layer were calculated using a numerical integration method. The sum of all fractional masses and weights represents the total mass and weight of Earth’s atmosphere. The results obtained here agree well with measurements and models of …

A Fast Algorithm For The Inversion Of Quasiseparable Vandermonde-Like Matrices, Sirani M. Perera, Grigory Bonik, Vadim Olshevsky

Publications

The results on Vandermonde-like matrices were introduced as a generalization of polynomial Vandermonde matrices, and the displacement structure of these matrices was used to derive an inversion formula. In this paper we first present a fast Gaussian elimination algorithm for the polynomial Vandermonde-like matrices. Later we use the said algorithm to derive fast inversion algorithms for quasiseparable, semiseparable and well-free Vandermonde-like matrices having O(n2) complexity. To do so we identify structures of displacement operators in terms of generators and the recurrence relations(2-term and 3-term) between the columns of the basis transformation matrices for quasiseparable, semiseparable and well-free polynomials. Finally we …

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.

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 …

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.

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.

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.

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.

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.

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.

Deformation Waves In Microstructured Materials: Theory And Numerics, Juri Engelbrecht, Arkadi Berezovski, Mihhail Berezovski

Publications

A linear model of the microstructured continuum based on Mindlin theory is adopted which can be represented in the framework of the internal variable theory. Fully coupled systems of equations for macro-motion and microstructure evolution are represented in the form of conservation laws. A modification of wave propagation algorithm is used for numerical calculations. Results of direct numerical simulations of wave propagation in periodic medium are compared with similar results for the continuous media with the modelled microstructure. It is shown that the proper choice of material constants should be made to match the results obtained by both approaches

Elements Of Study On Dynamic Materials, Marine Rousseau, Gerard A. Maugin, Mihhail Berezovski

Publications

As a preliminary study to more complex situations of interest in small-scale technology, this paper envisages the elementary propagation properties of elastic waves in one-spatial dimension when some of the properties (mass density, elasticity) may vary suddenly in space or in time, the second case being of course more original. Combination of the two may be of even greater interest. Toward this goal, a critical examination of what happens to solutions at the crossing of pure space-like and time-like material discontinuities is given together with simple solutions for smooth transitions and numerical simulations in the discontinuous case. The effects on …

Temporal Scales For Transport Patterns In The Gulf Of Finland, Bert Viikmae, Tarmo Soomere, Mikk Viidebaum, Mihhail Berezovski

Publications

The basic time scales for current-induced net transport of surface water and associated time scales of reaching the nearshore in the Gulf of Finland, the Baltic Sea, are analysed based on Lagrangian trajectories of water particles reconstructed from three-dimensional velocity fields by the Rossby Centre circulation model for 1987–1991. The number of particles reaching the nearshore exhibits substantial temporal variability whereas the rate of leaving the gulf is almost steady. It is recommended to use an about 3 grid cells wide nearshore area as a substitute to the coastal zone and about 10–15 day long trajectories for calculations of the …

Waves In Materials With Microstructure: Numerical Simulation, Mihhail Berezovski, Arkadi Berezovski, Juri Engelbrecht

Publications

Results of numerical experiments are presented in order to compare direct numerical calculations of wave propagation in a laminate with prescribed properties and corresponding results obtained for an effective medium with the microstructure modelling. These numerical experiments allowed us to analyse the advantages and weaknesses of the microstructure model.

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 …

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 …

Numerical Simulation Of Waves And Fronts In Inhomogeneous Solids, A. Berezovski, M. Berezovski, J. Engelbrecht, G. A. Maugin

Publications

Dynamic response of inhomogeneous materials exhibits new effects, which often do not exist in homogeneous media. It is quite natural that most of studies of wave and front propagation in inhomogeneous materials are associated with numerical simulations. To develop a numerical algorithm and to perform the numerical simulations of moving fronts we need to formulate a kinetic law of progress relating the driving force and the velocity of the discontinuity. The velocity of discontinuity is determined by means of the non-equilibrium jump relations at the front. The obtained numerical method generalizes the wave-propagation algorithm to the case of moving discontinuities …

Multiphoton Response Of Retinal Rod Photoreceptors, Vasilios Alexiades, Harihar Khanal

Publications

Phototransduction is the process by which light is converted into an electrical response in retinal photoreceptors. Rod photoreceptors contain a stack of (about 1000) disc membranes packed with photopigment rhodopsin molecules, which absorb the photons. We present computational experiments which show the profound effect on the response of the distances (how many discs apart) photons happen to be absorbed at. This photon-distribution effect alone can account for much of the observed variability in response.