Numerical Analysis and Computation Commons^™

Modeling, Simulation And Control Of Microrobots For The Microfactory., Zhong Yang

Electronic Theses and Dissertations

Future assembly technologies will involve higher levels of automation in order to satisfy increased microscale or nanoscale precision requirements. Traditionally, assembly using a top-down robotic approach has been well-studied and applied to the microelectronics and MEMS industries, but less so in nanotechnology. With the boom of nanotechnology since the 1990s, newly designed products with new materials, coatings, and nanoparticles are gradually entering everyone’s lives, while the industry has grown into a billion-dollar volume worldwide. Traditionally, nanotechnology products are assembled using bottom-up methods, such as self-assembly, rather than top-down robotic assembly. This is due to considerations of volume handling of large …

Mathematical Modeling Of Egain Droplets Sliding Down An Inclined Plane, Shawn Koohy, Kathryn Massey, Luis Schneegans, Megan Vezzetti, Keith Hillaire, Souradip Chattopadhyay, Jessie Chen, Carmen Lee, Karen Daniels

Undergraduate Research Symposium

Eutectic Gallium-Indium (EGaIn) is a room-temperature liquid metal alloy that dramatically changes its surface tension and dynamics under an applied electric field. EGaIn has been used heavily in soft electronics engineering due to its high conductivity, malleability, and safety. However, the absence of mathematical modeling in the current literature makes its behavior difficult to understand and predict. In this study, we present a one-dimensional lubrication model for the dynamics of an EGaIn droplet moving along an inclined plane. Our model incorporates essential physical effects and parameters including oxidation, capillary action, diffusion, gravity, and Marangoni effects. In particular, we incorporate effects …

Practical Implementation Of The Immersed Interface Method With Triangular Meshes For 3d Rigid Solids In A Fluid Flow, Norah Hakami

Mathematics Theses and Dissertations

When employing the immersed interface method (IIM) to simulate a fluid flow around a moving rigid object, the immersed object can be replaced by a virtual fluid enclosed by singular forces on the interface between the real and virtual fluids. These forces represent the impact of the rigid motion on the fluid flow and cause jump discontinuities across the interface in the whole flow field. Then, the IIM resolves the fluid flow on a fixed computational domain by directly incorporating the jump conditions across the interface into numerical schemes. Previous development of the method is limited to simple smooth boundaries. …

Impact Of Spallation And Internal Radiation On Fibrous Ablative Materials, Raghava Sai Chaitanya Davuluri

Theses and Dissertations--Mechanical and Aerospace Engineering

Space vehicles are equipped with Thermal Protection Systems (TPS) that encounter high heat rates and protect the payload while entering a planetary atmosphere. For most missions that interest NASA, ablative materials are used as TPS. These materials undergo several mass and energy transfer mechanisms to absorb intense heat. The size and construction of the TPS are based on the composition of the planetary atmosphere and the impact of various ablative mechanisms on the flow field and the material. Therefore, it is essential to quantify the rates of different ablative phenomena to model TPS accurately. In this work, the impact of …

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 …

Modelling Spherical Aberration Detection In An Analog Holographic Wavefront Sensor, Emma Branigan, Suzanne Martin, Matthew Sheehan, Kevin Murphy

Conference Papers

The analog holographic wavefront sensor (AHWFS) is a simple and robust solution to wavefront sensing in turbulent environments. Here, the ability of a photopolymer based AHWFS to detect refractively generated spherical aberration is modelled and verified.

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 …

Modeling Cherenkov Light Detection Timing For The Very Energetic Radiation Imaging Telescope Array System, Keilan Finn Ramirez

Physics

The Very Energetic Radiation Imaging Telescope Array System (VERITAS) is an array of four 12-meter telescopes which use the Imaging Atmospheric Cherenkov Technique to conduct high-energy gamma-ray astronomy. VERITAS detects magnitude and location information associated with Cherenkov light, and uses this information to indirectly observe gamma-rays through a software reconstruction process. VERITAS also records timing information corresponding to Cherenkov light detection, and this additional information could theoretically be incorporated into the reconstruction process to improve the accuracy of gamma-ray observations. The first step to including timing information is to understand when Cherenkov light detection would be expected from a known …

The Exact Factorization Equations For One- And Two-Level Systems, Bart Rosenzweig

Theses and Dissertations

Exact Factorization is a framework for studying quantum many-body problems. This decomposes the wavefunctions of such systems into conditional and marginal components. We derive corresponding evolution equations for molecular systems whose conditional electronic subsystems are described by one or two Born-Oppenheimer levels and develop a program for their mathematical study.

An Examination Of Fontan Circulation Using Differential Equation Models And Numerical Methods, Vanessa Maybruck

Honors Student Research

Certain congenital heart defects can lead to the development of only a single pumping chamber, or ventricle, in the heart instead of the usual two ventricles. Individuals with this defect undergo a corrective, three-part surgery, the third step of which is the Fontan procedure, but as the patients age, their cardiovascular health will likely deteriorate. Using computational fluid dynamics and differential equations, Fontan circulation can be modeled to investigate why the procedure fails and how Fontan failure can be maximally prevented. Borrowing from well-established literature on RC circuits, the differential equation models simulate systemic blood flow in a piecewise, switch-like …

Computational Modelling Enables Robust Multidimensional Nanoscopy, Matthew D. Lew

Electrical & Systems Engineering Publications and Presentations

The following sections are included:

• Present State of Computational Modelling in Fluorescence Nanoscopy

• Recent Contributions to Computational Modelling in Fluorescence Nanoscopy

• Outlook on Computational Modelling in Fluorescence Nanoscopy

• Acknowledgments

• References

Numerical Reconstruction Of Spalled Particle Trajectories In An Arc-Jet Environment, Raghava S. C. Davuluri, Sean C. C. Bailey, Kaveh A. Tagavi, Alexandre Martin

Mechanical Engineering Faculty Publications

To evaluate the effects of spallation on ablative material, it is necessary to evaluate the mass loss. To do so, a Lagrangian particle trajectory code is used to reconstruct trajectories that match the experimental data for all kinematic parameters. The results from spallation experiments conducted at the NASA HYMETS facility over a wedge sample were used. A data-driven adaptive methodology was used to adapts the ejection parameters until the numerical trajectory matches the experimental data. The preliminary reconstruction results show that the size of the particles seemed to be correlated with the location of the ejection event. The size of …

Five-Wave Resonances In Deep Water Gravity Waves: Integrability, Numerical Simulations And Experiments, Dan Lucas, Marc Perlin, Dian-Yong Liu, Shane Walsh, Rossen Ivanov, Miguel D. Bustamante

Articles

In this work we consider the problem of finding the simplest arrangement of resonant deep water gravity waves in one-dimensional propagation, from three perspectives: Theoretical, numerical and experimental. Theoretically this requires using a normal-form Hamiltonian that focuses on 5-wave resonances. The simplest arrangement is based on a triad of wave vectors K₁ + K₂ = K₃ (satisfying specific ratios) along with their negatives, corresponding to a scenario of encountering wave packets, amenable to experiments and numerical simulations. The normal-form equations for these encountering waves in resonance are shown to be non-integrable, but they admit an integrable reduction …

Development Of Synthetic Coal Char Simulant For Microwave Conversion Studies: A Computationally-Driven Approach, Kevin A. Hager

Graduate Theses, Dissertations, and Problem Reports (ETD)

Recent experimental demonstration of new reaction windows for coal char/methane reactions that are less energy-intensive, provides innovation for modular reactors. However, the correlation of the exact mechanism for the enhancement of these reaction windows is not certain. This study investigates the simplification of these experimental studies by developing a well-characterized coal char simulant. The approach involves using a computational approach to screen macroscopic composition to replicate the dielectric and compositional response of actual char. This study is focused on PRB coal char. A discrete element method (DEM) technique was used to simulate the packing of coal chars to give the …

Improving The Temporal Accuracy Of Turbulence Models And Resolving The Implementation Issues Of Fluid Flow Modeling, Kyle J. Schwiebert

Dissertations, Master's Theses and Master's Reports

A sizeable proportion of the work in this thesis focuses on a new turbulence model, dubbed ADC (the approximate deconvolution model with defect correction). The ADC is improved upon using spectral deferred correction, a means of constructing a higher order ODE solver. Since both the ADC and SDC are based on a predictor-corrector approach, SDC is incorporated with essentially no additional computational cost. We will show theoretically and using numerical tests that the new scheme is indeed higher order in time than the original, and that the benefits of defect correction, on which the ADC is based, are preserved.

The …

Accurate Spectral Algorithms For Solving Variable-Order Fractional Percolation Equations, M. A. Abdelkawy

Applications and Applied Mathematics: An International Journal (AAM)

A high accurate spectral algorithm for one-dimensional variable-order fractional percolation equations (VO-FPEs) is considered.We propose a shifted Legendre Gauss-Lobatto collocation (SL-GLC) method in conjunction with shifted Chebyshev Gauss-Radau collocation (SC-GR-C) method to solve the proposed problem. Firstly, the solution and its space fractional derivatives are expanded as shifted Legendre polynomials series. Then, we determine the expansion coefficients by reducing the VO-FPEs and its conditions to a system of ordinary differential equations (SODEs) in time. The numerical approximation of SODEs is achieved by means of the SC-GR-C method. The under-study’s problem subjected to the Dirichlet or non-local boundary conditions is presented …

From Wave Propagation To Spin Dynamics: Mathematical And Computational Aspects, Oleksii Beznosov

Mathematics & Statistics ETDs

In this work we concentrate on two separate topics which pose certain numerical challenges. The first topic is the spin dynamics of electrons in high-energy circular accelerators. We introduce a stochastic differential equation framework to study spin depolarization and spin equilibrium. This framework allows the mathematical study of known equations and new equations modelling the spin distribution of an electron bunch. A spin distribution is governed by a so-called Bloch equation, which is a linear Fokker-Planck type PDE, in general posed in six dimensions. We propose three approaches to approximate solutions, using analytical and modern numerical techniques. We also present …

Numerical Simulations Of Nonlinear Waves And Their Stability: Stokes Waves And Nonlinear Schroedinger Equation, Anastassiya Semenova

Mathematics & Statistics ETDs

The present work offers an investigation of dynamics and stability of nonlinear waves in Hamiltonian systems. The first part of the manuscript discusses the classical problem of water waves on the surface of an ideal fluid in 2D. We demonstrate how to construct the Stokes waves, and how to apply a continuation method to find waves in close vicinity to the limiting Stokes wave. We provide new insight into the stability of the Stokes waves by identifying previously inaccessible branches of instability in the equations of motion for the fluid. We provide numerical evidence that pairs of unstable eigenvalues of …

A Study Of The Design Of Adaptive Camber Winglets, Justin J. Rosescu

Master's Theses

A numerical study was conducted to determine the effect of changing the camber of a winglet on the efficiency of a wing in two distinct flight conditions. Camber was altered via a simple plain flap deflection in the winglet, which produced a constant camber change over the winglet span. Hinge points were located at 20%, 50% and 80% of the chord and the trailing edge was deflected between -5° and +5°. Analysis was performed using a combination of three-dimensional vortex lattice method and two-dimensional panel method to obtain aerodynamic forces for the entire wing, based on different winglet camber configurations. …

A Computationally-Efficient Bound For The Variance Of Measuring The Orientation Of Single Molecules, Tingting Wu, Tianben Ding, Hesam Mazidi, Oumeng Zhang, Matthew D. Lew

Electrical & Systems Engineering Publications and Presentations

Modulating the polarization of excitation light, resolving the polarization of emitted fluorescence, and point spread function (PSF) engineering have been widely leveraged for measuring the orientation of single molecules. Typically, the performance of these techniques is optimized and quantified using the Cramér-Rao bound (CRB), which describes the best possible measurement variance of an unbiased estimator. However, CRB is a local measure and requires exhaustive sampling across the measurement space to fully characterize measurement precision. We develop a global variance upper bound (VUB) for fast quantification and comparison of orientation measurement techniques. Our VUB tightly bounds the diagonal elements of the …

Effects Of Aperiodicity And Frustration On The Magnetic Properties Of Artificial Quasicrystals, Barry Farmer

Theses and Dissertations--Physics and Astronomy

Quasicrystals have been shown to exhibit physical properties that are dramatically different from their periodic counterparts. A limited number of magnetic quasicrystals have been fabricated and measured, and they do not exhibit long-range magnetic order, which is in direct conflict with simulations that indicate such a state should be accessible. This dissertation adopts a metamaterials approach in which artificial quasicrystals are fabricated and studied with the specific goal of identifying how aperiodicity affects magnetic long-range order. Electron beam lithography techniques were used to pattern magnetic thin films into two types of aperiodic tilings, the Penrose P2, and Ammann-Beenker tilings. SQUID …

Interpretations Of Bicoherence In Space & Lab Plasma Dynamics, Gregory Allen Riggs

Graduate Theses, Dissertations, and Problem Reports (ETD)

The application of bicoherence analysis to plasma research, particularly in non-linear, coupled-wave regimes, has thus far been significantly belied by poor resolution in time, and/or outright destruction of frequency information. Though the typical power spectrum cloaks the phase-coherency between frequencies, Fourier transforms of higher-order convolutions provide an n-dimensional spectrum which is adept at elucidating n-wave phase coherence. As such, this investigation focuses on the utility of the normalized bispectrum for detection of wave-wave coupling in general, with emphasis on distinct implications within the scope of non-linear plasma physics. Interpretations of bicoherent features are given for time series from …

Local Non-Similar Solution Of Powell-Eyring Fluid Flow Over A Vertical Flat Plate, Hemangini Shukla, Hema C. Surati, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

Our objective is to obtain the non-similarity solution of non-Newtonian fluid for Powell-Eyring model by a local non-similarity method. Here, free stream velocity is considered in power-law form (𝑈=𝑥^m). The governing equations are transformed using non-similar transformations and derived equations are treated as ordinary differential equations. Non-similar solutions are obtained for different values of power-law index 𝑚 and stream-wise location 𝜉. Influence of various parameters on velocity and temperature field are presented graphically using MATLAB bvp4c solver.

Derivation Of Direct Explicit Integrators Of Rk Type For Solving Class Of Seventh-Order Ordinary Differential Equations, Mohammed S. Mechee, Jawad K. Mshachal

Karbala International Journal of Modern Science

The main contribution of this work is the development of direct explicit methods of Runge-Kutta (RK) type for solving class of seventh-order ordinary differential equations (ODEs) to improve computational efficiency. For this purpose, we have generalized RK, RKN, RKD, RKT, RKFD and RKM methods for solving class of first-, second-, third-, fourth-, and fifth-order ODEs. Using Taylor expansion approach, we have derived the algebraic equations of the order conditions for the proposed RKM integrators up to the tenth-order. Based on these order conditions, two RKM methods of fifth- and sixth-order with four- and five-stage are derived. The zero stability of …

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 …

Solitary And Periodic Exact Solutions Of The Viscosity-Capillarity Van Der Waals Gas Equations, Emad A. Az-Zo’Bi

Applications and Applied Mathematics: An International Journal (AAM)

Periodic and soliton solutions are derived for the (1+1)-dimensional van der Waals gas system in the viscosity-capillarity regularization form. The system is handled via the e-^φ(ξ) -expansion method. The obtained solutions have been articulated by the hyperbolic, trigonometric, exponential and rational functions with arbitrary constants. Mathematical analysis and numerical graphs are provided for some solitons, periodic and kink solitary wave solutions to visualize the dynamics of equations. Obtained results reveal that the method is very influential and effective tool for solving nonlinear partial differential equations in applied mathematics.

Numerical Solution Of 3rd Order Ode Using Fdm: On A Moving Surface In Mhd Flow Of Sisko Fluid, Manisha Patel, Jayshri Patel, M. G. Timol

Applications and Applied Mathematics: An International Journal (AAM)

A Similarity group theoretical technique is used to transform the governing nonlinear partial differential equations of two dimensional MHD boundary layer flow of Sisko fluid into nonlinear ordinary differential equations. Then the resulting third order nonlinear ordinary differential equation with corresponding boundary conditions is linearised by Quasi linearization method. Numerical solution of the linearised third order ODE is obtained using Finite Difference method (FDM). Graphical presentation of the solution is given.

Surface Waves Over Currents And Uneven Bottom, Alan Compelli, Rossen Ivanov, Calin I. Martin, Michail D. Todorov

Articles

The propagation of surface water waves interacting with a current and an uneven bottom is studied. Such a situation is typical for ocean waves where the winds generate currents in the top layer of the ocean. The role of the bottom topography is taken into account since it also influences the local wave and current patterns. Specific scaling of the variables is selected which leads to approximations of Boussinesq and KdV types. The arising KdV equation with variable coefficients, dependent on the bottom topography, is studied numerically when the initial condition is in the form of the one soliton solution …

A Detection And Data Acquisition System For Precision Beta Decay Spectroscopy, Aaron P. Jezghani

Theses and Dissertations--Physics and Astronomy

Free neutron and nuclear beta decay spectroscopy serves as a robust laboratory for investigations of the Standard Model of Particle Physics. Observables such as decay product angular correlations and energy spectra overconstrain the Standard Model and serve as a sensitive probe for Beyond the Standard Model physics. Improved measurement of these quantities is necessary to complement the TeV scale physics being conducted at the Large Hadron Collider. The UCNB, ⁴⁵Ca, and Nab experiments aim to improve upon existing measurements of free neutron decay angular correlations and set new limits in the search for exotic couplings in beta decay. To …

Quasi-Linearization Method With Rational Legendre Collocation Method For Solving Mhd Flow Over A Stretching Sheet With Variable Thickness And Slip Velocity Which Embedded In A Porous Medium, M. M. Khader

Applications and Applied Mathematics: An International Journal (AAM)

The quasi-linearization method (QLM) and the rational Legendre functions are introduced here to present the numerical solution for the Newtonian fluid flow past an impermeable stretching sheet which embedded in a porous medium with a power-law surface velocity, variable thickness and slip velocity. Firstly, due to the high nonlinearity which yielded from the ordinary differential equation which describes the proposed physical problem, we construct a sequence of linear ODEs by using the QLM, hence the resulted equations become a system of linear algebraic equations. The comparison with the available results in the literature review proves that the obtained results via …