Matching Items (6)
Filtering by
- All Subjects: MATLAB
- All Subjects: Spectral Methods
- Creators: Platte, Rodrigo
Description
The dynamics of a fluid flow inside 2D square and 3D cubic cavities
under various configurations were simulated and analyzed using a
spectral code I developed.
This code was validated against known studies in the 3D lid-driven
cavity. It was then used to explore the various dynamical behaviors
close to the onset of instability of the steady-state flow, and explain
in the process the mechanism underlying an intermittent bursting
previously observed. A fairly complete bifurcation picture emerged,
using a combination of computational tools such as selective
frequency damping, edge-state tracking and subspace restriction.
The code was then used to investigate the flow in a 2D square cavity
under stable temperature stratification, an idealized version of a lake
with warmer water at the surface compared to the bottom. The governing
equations are the Navier-Stokes equations under the Boussinesq approximation.
Simulations were done over a wide range of parameters of the problem quantifying
the driving velocity at the top (e.g. wind) and the strength of the stratification.
Particular attention was paid to the mechanisms associated with the onset of
instability of the base steady state, and the complex nontrivial dynamics
occurring beyond onset, where the presence of multiple states leads to a
rich spectrum of states, including homoclinic and heteroclinic chaos.
A third configuration investigates the flow dynamics of a fluid in a rapidly
rotating cube subjected to small amplitude modulations. The responses were
quantified by the global helicity and energy measures, and various peak
responses associated to resonances with intrinsic eigenmodes of the cavity
and/or internal retracing beams were clearly identified for the first time.
A novel approach to compute the eigenmodes is also described, making accessible
a whole catalog of these with various properties and dynamics. When the small
amplitude modulation does not align with the rotation axis (precession) we show
that a new set of eigenmodes are primarily excited as the angular velocity
increases, while triadic resonances may occur once the nonlinear regime kicks in.
under various configurations were simulated and analyzed using a
spectral code I developed.
This code was validated against known studies in the 3D lid-driven
cavity. It was then used to explore the various dynamical behaviors
close to the onset of instability of the steady-state flow, and explain
in the process the mechanism underlying an intermittent bursting
previously observed. A fairly complete bifurcation picture emerged,
using a combination of computational tools such as selective
frequency damping, edge-state tracking and subspace restriction.
The code was then used to investigate the flow in a 2D square cavity
under stable temperature stratification, an idealized version of a lake
with warmer water at the surface compared to the bottom. The governing
equations are the Navier-Stokes equations under the Boussinesq approximation.
Simulations were done over a wide range of parameters of the problem quantifying
the driving velocity at the top (e.g. wind) and the strength of the stratification.
Particular attention was paid to the mechanisms associated with the onset of
instability of the base steady state, and the complex nontrivial dynamics
occurring beyond onset, where the presence of multiple states leads to a
rich spectrum of states, including homoclinic and heteroclinic chaos.
A third configuration investigates the flow dynamics of a fluid in a rapidly
rotating cube subjected to small amplitude modulations. The responses were
quantified by the global helicity and energy measures, and various peak
responses associated to resonances with intrinsic eigenmodes of the cavity
and/or internal retracing beams were clearly identified for the first time.
A novel approach to compute the eigenmodes is also described, making accessible
a whole catalog of these with various properties and dynamics. When the small
amplitude modulation does not align with the rotation axis (precession) we show
that a new set of eigenmodes are primarily excited as the angular velocity
increases, while triadic resonances may occur once the nonlinear regime kicks in.
ContributorsWu, Ke (Author) / Lopez, Juan (Thesis advisor) / Welfert, Bruno (Thesis advisor) / Tang, Wenbo (Committee member) / Platte, Rodrigo (Committee member) / Herrmann, Marcus (Committee member) / Arizona State University (Publisher)
Created2019
Description
Chebfun is a collection of algorithms and an open-source software system in object-oriented Matlab that extends familiar powerful methods of numerical computation involving numbers to continuous or piecewise-continuous functions. The success of this strategy is based on the mathematical fact that smooth functions can be represented very efficiently by polynomial interpolation at Chebyshev points or by trigonometric interpolation at equispaced points for periodic functions. More recently, the system has been extended to handle bivariate functions and vector fields. These two new classes of objects are called Chebfun2 and Chebfun2v, respectively. We will show that Chebfun2 and Chebfun2v, and can be used to accurately and efficiently perform various computations on parametric surfaces in two or three dimensions, including path trajectories and mean and Gaussian curvatures. More advanced surface computations such as mean curvature flows are also explored. This is also the first work to use the newly implemented trigonometric representation, namely Trigfun, for computations on surfaces.
ContributorsPage-Bottorff, Courtney Michelle (Author) / Platte, Rodrigo (Thesis director) / Kostelich, Eric (Committee member) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
Description
Using object-oriented programming in MATLAB, a collection of functions, named Fourfun, has been created to allow quick and accurate approximations of periodic functions with Fourier expansions. To increase efficiency and reduce the number of computations of the Fourier transform, Fourfun automatically determines the number of nodes necessary for representations that are accurate to close to machine precision. Common MATLAB functions have been overloaded to keep the syntax of the Fourfun class as consistent as possible with the general MATLAB syntax. We show that the system can be used to efficiently solve several differential equations. Comparisons with Chebfun, a similar system based on Chebyshev polynomial approximations, are provided.
ContributorsMcleod, Kristyn Noelle (Author) / Platte, Rodrigo (Thesis director) / Gelb, Anne (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / School of International Letters and Cultures (Contributor)
Created2014-05
Description
High-order methods are known for their accuracy and computational performance when applied to solving partial differential equations and have widespread use
in representing images compactly. Nonetheless, high-order methods have difficulty representing functions containing discontinuities or functions having slow spectral decay in the chosen basis. Certain sensing techniques such as MRI and SAR provide data in terms of Fourier coefficients, and thus prescribe a natural high-order basis. The field of compressed sensing has introduced a set of techniques based on $\ell^1$ regularization that promote sparsity and facilitate working with functions having discontinuities. In this dissertation, high-order methods and $\ell^1$ regularization are used to address three problems: reconstructing piecewise smooth functions from sparse and and noisy Fourier data, recovering edge locations in piecewise smooth functions from sparse and noisy Fourier data, and reducing time-stepping constraints when numerically solving certain time-dependent hyperbolic partial differential equations.
in representing images compactly. Nonetheless, high-order methods have difficulty representing functions containing discontinuities or functions having slow spectral decay in the chosen basis. Certain sensing techniques such as MRI and SAR provide data in terms of Fourier coefficients, and thus prescribe a natural high-order basis. The field of compressed sensing has introduced a set of techniques based on $\ell^1$ regularization that promote sparsity and facilitate working with functions having discontinuities. In this dissertation, high-order methods and $\ell^1$ regularization are used to address three problems: reconstructing piecewise smooth functions from sparse and and noisy Fourier data, recovering edge locations in piecewise smooth functions from sparse and noisy Fourier data, and reducing time-stepping constraints when numerically solving certain time-dependent hyperbolic partial differential equations.
ContributorsDenker, Dennis (Author) / Gelb, Anne (Thesis advisor) / Archibald, Richard (Committee member) / Armbruster, Dieter (Committee member) / Boggess, Albert (Committee member) / Platte, Rodrigo (Committee member) / Saders, Toby (Committee member) / Arizona State University (Publisher)
Created2016
Description
This Creative Project was carried out in coordination with the capstone project, Around the Corner Imaging with Terahertz Waves. This capstone project deals with a system designed to implement Around the Corner, or Non Line-of-Sight (NLoS) Imaging. This document discusses the creation of a GUI using MATLAB to control the Terahertz Imaging system. The GUI was developed in response to a need for synchronization, ease of operation, easy parameter modification, and data management. Along the way, many design decisions were made ranging from choosing a software platform to determining how variables should be passed. These decisions and considerations are discussed in this document. The resulting GUI has measured up to the design criteria and will be able to be used by anyone wishing to use the Terahertz Imaging System for further research in the field of Around the Corner or NLoS Imaging.
ContributorsWood, Jacob Cannon (Author) / Trichopoulos, Georgios (Thesis director) / Aberle, James (Committee member) / Electrical Engineering Program (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2019-05
Description
The purpose of this project is to analyze the MIT OpenCourseWare coffee can radar design and modify it to be better suited for drone based synthetic aperture radar (SAR) applications while maintaining the low-cost aspect of the original design. The MIT coffee can radar can function as a ranged radar, a Doppler radar, or as SAR. Through simulations and research, the suggestions for how to modify the radar resulted in swapping the coffee can monopole antennas for patch antenna arrays or helical ordinary end-fire antennas, adding an Arduino for automatic recording of output pulses, and switching from a breadboard construction to a PCB to shrink form factor and keep costs and construction time low.
ContributorsRivera, Danielle (Author) / Trichopoulos, Georgios (Thesis director) / Aberle, James (Committee member) / Department of Information Systems (Contributor) / Electrical Engineering Program (Contributor) / Barrett, The Honors College (Contributor)
Created2020-12