Physical Sciences and Mathematics Commons^™

Rigid Body Constrained Motion Optimization And Control On Lie Groups And Their Tangent Bundles, Brennan S. Mccann

Doctoral Dissertations and Master's Theses

Rigid body motion requires formulations where rotational and translational motion are accounted for appropriately. Two Lie groups, the special orthogonal group SO(3) and the space of quaternions H, are commonly used to represent attitude. When considering rigid body pose, that is spacecraft position and attitude, the special Euclidean group SE(3) and the space of dual quaternions DH are frequently utilized. All these groups are Lie groups and Riemannian manifolds, and these identifications have profound implications for dynamics and controls. The trajectory optimization and optimal control problem on Riemannian manifolds presents significant opportunities for theoretical development. Riemannian optimization is an attractive …

Response Of Planetary Waves And Tides To The 2019 Southern Hemisphere Ssw And Q2dw Enhancement In Jan-Feb 2020 Observed By Condor Meteor Radar In Chile And Adelaide Meteor Radar In Australia, Alan Liu, Zishun Qiao, Iain Reid, Javier Fuentes, Chris Adami

Publications

A new multi-static meteor radar (CONDOR) has recently been installed in northern Chile. This CONDOR meteor radar (30.3°S, 70.7°W) and the Adelaide meteor radar (35°S, 138°E) have provided longitudinally spaced observations of the mean winds, tides and planetary waves of the PW-tides interaction cases we present here. We have observed a Quasi-6-Day Wave (Q6DW) enhancement in MLT winds at the middle latitudes (30.3°S, 35°S) during the unusual minor South Hemisphere SSW 2019 by the ground-based meteor radars. Tidal analysis also indicates modulation of the Q6DW w/ amplitude ~15 [m/s] and diurnal tides w/ amplitude ~60 [m/s]. Another case we present …

Manufacturability And Analysis Of Topologically Optimized Continuous Fiber Reinforced Composites, Jesus A. Ferrand

Doctoral Dissertations and Master's Theses

Researchers are unlocking the potential of Continuous Fiber Reinforced Composites for producing components with greater strength-to-weight ratios than state of the art metal alloys and unidirectional composites. The key is the emerging technology of topology optimization and advances in additive manufacturing. Topology optimization can fine tune component geometry and fiber placement all while satisfying stress constraints. However, the technology cannot yet robustly guarantee manufacturability. For this reason, substantial post-processing of an optimized design consisting of manual fiber replacement and subsequent Finite Element Analysis (FEA) is still required.

To automate this post-processing in two dimensions, two (2) algorithms were developed. The …

A Data Driven Modeling Approach For Store Distributed Load And Trajectory Prediction, Nicholas Peters

Doctoral Dissertations and Master's Theses

The task of achieving successful store separation from aircraft and spacecraft has historically been and continues to be, a critical issue for the aerospace industry. Whether it be from store-on-store wake interactions, store-parent body interactions or free stream turbulence, a failed case of store separation poses a serious risk to aircraft operators. Cases of failed store separation do not simply imply missing an intended target, but also bring the risk of collision with, and destruction of, the parent body vehicle. Given this risk, numerous well-tested procedures have been developed to help analyze store separation within the safe confines of wind …

Computational Models To Detect Radiation In Urban Environments: An Application Of Signal Processing Techniques And Neural Networks To Radiation Data Analysis, Jose Nicolas Gachancipa

Beyond: Undergraduate Research Journal

Radioactive sources, such as uranium-235, are nuclides that emit ionizing radiation, and which can be used to build nuclear weapons. In public areas, the presence of a radioactive nuclide can present a risk to the population, and therefore, it is imperative that threats are identified by radiological search and response teams in a timely and effective manner. In urban environments, such as densely populated cities, radioactive sources may be more difficult to detect, since background radiation produced by surrounding objects and structures (e.g., buildings, cars) can hinder the effective detection of unnatural radioactive material. This article presents a computational model …

Statistical Characteristics Of High-Frequency Gravity Waves Observed By An Airglow Imager At Andes Lidar Observatory, Alan Z. Liu, Bing Cao

Publications

The long-term statistical characteristics of high-frequency quasi-monochromatic gravity waves are presented using multi-year airglow images observed at Andes Lidar Observatory (ALO, 30.3° S, 70.7° W) in northern Chile. The distribution of primary gravity wave parameters including horizontal wavelength, vertical wavelength, intrinsic wave speed, and intrinsic wave period are obtained and are in the ranges of 20–30 km, 15–25 km, 50–100 m s−1, and 5–10 min, respectively. The duration of persistent gravity wave events captured by the imager approximately follows an exponential distribution with an average duration of 7–9 min. The waves tend to propagate against the local background winds and …

A Meshless Approach To Computational Pharmacokinetics, Anthony Matthew Khoury

Doctoral Dissertations and Master's Theses

The meshless method is an incredibly powerful technique for solving a variety of problems with unparalleled accuracy and efficiency. The pharmacokinetic problem of transdermal drug delivery (TDDD) is one such topic and is of significant complexity. The locally collocated meshless method (LCMM) is developed in solution to this topic. First, the meshless method is formulated to model this transport phenomenon and is then validated against an analytical solution of a pharmacokinetic problem set, to demonstrate this accuracy and efficiency. The analytical solution provides a locus by which convergence behavior are evaluated, demonstrating the super convergence of the locally collocated meshless …

A Mathematical Analysis Of The Wind Triangle Problem And An Inquiry Of True Airspeed Calculations In Supersonic Flight, Leonard T. Huang, Lisa I. Cummings

International Journal of Aviation, Aeronautics, and Aerospace

In the first half of this paper, we present a fresh perspective toward the Wind Triangle Problem in aerial navigation by deriving necessary and sufficient conditions, which we call "go/no-go conditions", for the existence/non-existence of a solution of the problem. Although our derivation is based on simple trigonometry and basic properties of quadratic functions, it is mathematically rigorous. We also offer examples to demonstrate how easy it is to check these conditions graphically. In the second half of this paper, we use function theory to re-examine another problem in aerial navigation, namely, that of computing true airspeed — even in …

Dynamics Of Discontinuities In Elastic Solids, Arkadi Berezovski, Mihhail Berezovski

Publications

The paper is devoted to evolving discontinuities in elastic solids. A discontinuity is represented as a singular set of material points. Evolution of a discontinuity is driven by the configurational force acting at such a set. The main attention is paid to the determination of the velocity of a propagating discontinuity. Martensitic phase transition fronts and brittle cracks are considered as representative examples.

Meshless Modeling Of Flow Dispersion And Progressive Piping In Poroelastic Levees, Anthony Khoury, Eduardo Divo, Alain J. Kassab, Sai Kakuturu, Lakshmi Reddi

Publications

Performance data on earth dams and levees continue to indicate that piping is one of the major causes of failure. Current criteria for prevention of piping in earth dams and levees have remained largely empirical. This paper aims at developing a mechanistic understanding of the conditions necessary to prevent piping and to enhance the likelihood of self-healing of cracks in levees subjected to hydrodynamic loading from astronomical and meteorological (including hurricane storm surge-induced) forces. Systematic experimental investigations are performed to evaluate erosion in finite-length cracks as a result of transient hydrodynamic loading. Here, a novel application of the localized collocation …

Full Field Computing For Elastic Pulse Dispersion In Inhomogeneous Bars, A. Berezovski, R. Kolman, M. Berezovski, D. Gabriel, V. Adamek

Publications

In the paper, the finite element method and the finite volume method are used in parallel for the simulation of a pulse propagation in periodically layered composites beyond the validity of homogenization methods. The direct numerical integration of a pulse propagation demonstrates dispersion effects and dynamic stress redistribution in physical space on example of a one-dimensional layered bar. Results of numerical simulations are compared with analytical solution constructed specifically for the considered problem. Analytical solution as well as numerical computations show the strong influence of the composition of constituents on the dispersion of a pulse in a heterogeneous bar and …

Patient-Specific Multiscale Computational Fluid Dynamics Assessment Of Embolization Rates In The Hybrid Norwood: Effects Of Size And Placement Of The Reverse Blalock–Taussig Shunt, Ray Prather, John Seligson, Marcus Ni, Eduardo Divo, Alain J. Kassab, William Decampli

Publications

The hybrid Norwood operation is performed to treat hypoplastic left heart syndrome. Distal arch obstruction may compromise flow to the brain. In a variant of this procedure, a synthetic graft (reverse Blalock–Taussig shunt) is placed between the pulmonary trunk and innominate artery to improve upper torso blood flow. Thrombi originating in the graft may embolize to the brain. In this study, we used computational fluid dynamics and particle tracking to investigate the patterns of particle embolization as a function of the anatomic position of the reverse Blalock–Taussig shunt. The degree of distal arch obstruction and position of particle origin influence …

Godunov-Type Upwind Flux Schemes Of The Two-Dimensional Finite Volume Discrete Boltzmann Method, Leitao Chen, Laura Schaefer

Publications

A simple unified Godunov-type upwind approach that does not need Riemann solvers for the flux calculation is developed for the finite volume discrete Boltzmann method (FVDBM) on an unstructured cell-centered triangular mesh. With piecewise-constant (PC), piecewise-linear (PL) and piecewise-parabolic (PP) reconstructions, three Godunov-type upwind flux schemes with different orders of accuracy are subsequently derived. After developing both a semi-implicit time marching scheme tailored for the developed flux schemes, and a versatile boundary treatment that is compatible with all of the flux schemes presented in this paper, numerical tests are conducted on spatial accuracy for several single-phase flow problems. Four major …

Numerical Simulation Of Energy Localization In Dynamic Materials, Arkadi Berezovski, Mihhail Berezovski

Publications

Dynamic materials are artificially constructed in such a way that they may vary their characteristic properties in space or in time, or both, by an appropriate arrangement or control. These controlled changes in time can be provided by the application of an external (non-mechanical) field, or through a phase transition. In principle, all materials change their properties with time, but very slowly and smoothly. Changes in properties of dynamic materials should be realized in a short or quasi-nil time lapse and over a sufficiently large material region. Wave propagation is a characteristic feature for dynamic materials because it is also …

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 …

A Coupled Localized Rbf Meshless/Drbem Formulation For Accurate Modeling Of Incompressible Fluid Flows, Leonardo Bueno, Eduardo Divo, Alain J. Kassab

Publications

Velocity-pressure coupling schemes for the solution of incompressible fluid flow problems in Computational Fluid Dynamics (CFD) rely on the formulation of Poisson-like equations through projection methods. The solution of these Poisson-like equations represent the pressure correction and the velocity correction to ensure proper satisfaction of the conservation of mass equation at each step of a time-marching scheme or at each level of an iteration process. Inaccurate solutions of these Poisson-like equations result in meaningless instantaneous or intermediate approximations that do not represent the proper time-accurate behavior of the flow. The fact that these equations must be solved to convergence at …

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.

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. …

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 …

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 …

Burglary Crime Analysis Using Logistic Regression, Daniel Antolos, Dahai Liu, Andrei Ludu, Dennis Vincenzi

Andrei Ludu

This study used a logistic regression model to investigate the relationship between several predicting factors and burglary occurrence probability with regard to the epicenter. These factors include day of the week, time of the day, repeated victimization, connectors and barriers. Data was collected from a local police report on 2010 burglary incidents. Results showed the model has various degrees of significance in terms of predicting the occurrence within difference ranges from the epicenter. Follow-up refined multiple comparisons of different sizes were observed to further discover the pattern of prediction strength of these factors. Results are discussed and further research directions …

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 …