Matching Items (42)
Description
Earth-system models describe the interacting components of the climate system and
technological systems that affect society, such as communication infrastructures. Data
assimilation addresses the challenge of state specification by incorporating system
observations into the model estimates. In this research, a particular data
assimilation technique called the Local Ensemble Transform Kalman Filter (LETKF) is
applied to the ionosphere, which is a domain of practical interest due to its effects
on infrastructures that depend on satellite communication and remote sensing. This
dissertation consists of three main studies that propose strategies to improve space-
weather specification during ionospheric extreme events, but are generally applicable
to Earth-system models:
Topic I applies the LETKF to estimate ion density with an idealized model of
the ionosphere, given noisy synthetic observations of varying sparsity. Results show
that the LETKF yields accurate estimates of the ion density field and unobserved
components of neutral winds even when the observation density is spatially sparse
(2% of grid points) and there is large levels (40%) of Gaussian observation noise.
Topic II proposes a targeted observing strategy for data assimilation, which uses
the influence matrix diagnostic to target errors in chosen state variables. This
strategy is applied in observing system experiments, in which synthetic electron density
observations are assimilated with the LETKF into the Thermosphere-Ionosphere-
Electrodynamics Global Circulation Model (TIEGCM) during a geomagnetic storm.
Results show that assimilating targeted electron density observations yields on
average about 60%–80% reduction in electron density error within a 600 km radius of
the observed location, compared to 15% reduction obtained with randomly placed
vertical profiles.
Topic III proposes a methodology to account for systematic model bias arising
ifrom errors in parametrized solar and magnetospheric inputs. This strategy is ap-
plied with the TIEGCM during a geomagnetic storm, and is used to estimate the
spatiotemporal variations of bias in electron density predictions during the
transitionary phases of the geomagnetic storm. Results show that this strategy reduces
error in 1-hour predictions of electron density by about 35% and 30% in polar regions
during the main and relaxation phases of the geomagnetic storm, respectively.
technological systems that affect society, such as communication infrastructures. Data
assimilation addresses the challenge of state specification by incorporating system
observations into the model estimates. In this research, a particular data
assimilation technique called the Local Ensemble Transform Kalman Filter (LETKF) is
applied to the ionosphere, which is a domain of practical interest due to its effects
on infrastructures that depend on satellite communication and remote sensing. This
dissertation consists of three main studies that propose strategies to improve space-
weather specification during ionospheric extreme events, but are generally applicable
to Earth-system models:
Topic I applies the LETKF to estimate ion density with an idealized model of
the ionosphere, given noisy synthetic observations of varying sparsity. Results show
that the LETKF yields accurate estimates of the ion density field and unobserved
components of neutral winds even when the observation density is spatially sparse
(2% of grid points) and there is large levels (40%) of Gaussian observation noise.
Topic II proposes a targeted observing strategy for data assimilation, which uses
the influence matrix diagnostic to target errors in chosen state variables. This
strategy is applied in observing system experiments, in which synthetic electron density
observations are assimilated with the LETKF into the Thermosphere-Ionosphere-
Electrodynamics Global Circulation Model (TIEGCM) during a geomagnetic storm.
Results show that assimilating targeted electron density observations yields on
average about 60%–80% reduction in electron density error within a 600 km radius of
the observed location, compared to 15% reduction obtained with randomly placed
vertical profiles.
Topic III proposes a methodology to account for systematic model bias arising
ifrom errors in parametrized solar and magnetospheric inputs. This strategy is ap-
plied with the TIEGCM during a geomagnetic storm, and is used to estimate the
spatiotemporal variations of bias in electron density predictions during the
transitionary phases of the geomagnetic storm. Results show that this strategy reduces
error in 1-hour predictions of electron density by about 35% and 30% in polar regions
during the main and relaxation phases of the geomagnetic storm, respectively.
ContributorsDurazo, Juan, Ph.D (Author) / Kostelich, Eric J. (Thesis advisor) / Mahalov, Alex (Thesis advisor) / Tang, Wenbo (Committee member) / Moustaoui, Mohamed (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2018
Description
A continuously and stably stratified fluid contained in a square cavity subjected to harmonic body forcing is studied numerically by solving the Navier-Stokes equations under the Boussinesq approximation. Complex dynamics are observed near the onset of instability of the basic state, which is a flow configuration that is always an exact analytical solution of the governing equations. The instability of the basic state to perturbations is first studied with linear stability analysis (Floquet analysis), revealing a multitude of intersecting synchronous and subharmonic resonance tongues in parameter space. A modal reduction method for determining the locus of basic state instability is also shown, greatly simplifying the computational overhead normally required by a Floquet study. Then, a study of the nonlinear governing equations determines the criticality of the basic state's instability, and ultimately characterizes the dynamics of the lowest order spatial mode by the three discovered codimension-two bifurcation points within the resonance tongue. The rich dynamics include a homoclinic doubling cascade that resembles the logistic map and a multitude of gluing bifurcations.
The numerical techniques and methodologies are first demonstrated on a homogeneous fluid contained within a three-dimensional lid-driven cavity. The edge state technique and linear stability analysis through Arnoldi iteration are used to resolve the complex dynamics of the canonical shear-driven benchmark problem. The techniques here lead to a dynamical description of an instability mechanism, and the work serves as a basis for the remainder of the dissertation.
The numerical techniques and methodologies are first demonstrated on a homogeneous fluid contained within a three-dimensional lid-driven cavity. The edge state technique and linear stability analysis through Arnoldi iteration are used to resolve the complex dynamics of the canonical shear-driven benchmark problem. The techniques here lead to a dynamical description of an instability mechanism, and the work serves as a basis for the remainder of the dissertation.
ContributorsYalim, Jason (Author) / Welfert, Bruno D. (Thesis advisor) / Lopez, Juan M. (Thesis advisor) / Jones, Donald (Committee member) / Tang, Wenbo (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2019
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
The main goal of this project is to study approximations of functions on circular and spherical domains using the cubed sphere discretization. On each subdomain, the function is approximated by windowed Fourier expansions. Of particular interest is the dependence of accuracy on the different choices of windows and the size of the overlapping regions. We use Matlab to manipulate each of the variables involved in these computations as well as the overall error, thus enabling us to decide which specific values produce the most accurate results. This work is motivated by problems arising in atmospheric research.
ContributorsSopa, Megan Grace (Author) / Platte, Rodrigo (Thesis director) / Kostelich, Eric (Committee member) / Department of Information Systems (Contributor) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2018-05
Description
The detection and characterization of transients in signals is important in many wide-ranging applications from computer vision to audio processing. Edge detection on images is typically realized using small, local, discrete convolution kernels, but this is not possible when samples are measured directly in the frequency domain. The concentration factor edge detection method was therefore developed to realize an edge detector directly from spectral data. This thesis explores the possibilities of detecting edges from the phase of the spectral data, that is, without the magnitude of the sampled spectral data. Prior work has demonstrated that the spectral phase contains particularly important information about underlying features in a signal. Furthermore, the concentration factor method yields some insight into the detection of edges in spectral phase data. An iterative design approach was taken to realize an edge detector using only the spectral phase data, also allowing for the design of an edge detector when phase data are intermittent or corrupted. Problem formulations showing the power of the design approach are given throughout. A post-processing scheme relying on the difference of multiple edge approximations yields a strong edge detector which is shown to be resilient under noisy, intermittent phase data. Lastly, a thresholding technique is applied to give an explicit enhanced edge detector ready to be used. Examples throughout are demonstrate both on signals and images.
ContributorsReynolds, Alexander Bryce (Author) / Gelb, Anne (Thesis director) / Cochran, Douglas (Committee member) / Viswanathan, Adityavikram (Committee member) / School of Mathematical and Statistical Sciences (Contributor) / Barrett, The Honors College (Contributor)
Created2016-05
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
Physical limitations of Magnetic Resonance Imaging (MRI) introduce different errors in the image reconstruction process. The discretization and truncation of data under discrete Fourier transform causes oscillations near jump discontinuities, a phenomenon known as the Gibbs effect. Using Gaussian-based approximations rather than the discrete Fourier transform to reconstruct images serves to diminish the Gibbs effect slightly, especially when coupled with filtering. Additionally, a simplifying assumption is made that, during signal collection, the amount of transverse magnetization decay at a point does not depend on that point's position in space. Though this methodology significantly reduces operational run-time, it nonetheless introduces geometric error, which can be mitigated using Single-Shot (SS) Parse.
ContributorsNeufer, Ian Douglas (Author) / Platte, Rodrigo (Thesis director) / Gelb, Anne (Committee member) / Barrett, The Honors College (Contributor) / School of Mathematical and Statistical Sciences (Contributor)
Created2015-05
Description
The reconstruction of piecewise smooth functions from non-uniform Fourier data arises in sensing applications such as magnetic resonance imaging (MRI). This thesis presents a new polynomial based resampling method (PRM) for 1-dimensional problems which uses edge information to recover the Fourier transform at its integer coefficients, thereby enabling the use of the inverse fast Fourier transform algorithm. By minimizing the error of the PRM approximation at the sampled Fourier modes, the PRM can also be used to improve on initial edge location estimates. Numerical examples show that using the PRM to improve on initial edge location estimates and then taking of the PRM approximation of the integer frequency Fourier coefficients is a viable way to reconstruct the underlying function in one dimension. In particular, the PRM is shown to converge more quickly and to be more robust than current resampling techniques used in MRI, and is particularly amenable to highly irregular sampling patterns.
ContributorsGutierrez, Alexander Jay (Author) / Platte, Rodrigo (Thesis director) / Gelb, Anne (Committee member) / Viswanathan, Adityavikram (Committee member) / Barrett, The Honors College (Contributor) / School of International Letters and Cultures (Contributor) / School of Mathematical and Statistical Sciences (Contributor)
Created2013-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
Divergence-free vector field interpolants properties are explored on uniform and scattered nodes, and also their application to fluid flow problems. These interpolants may be applied to physical problems that require the approximant to have zero divergence, such as the velocity field in the incompressible Navier-Stokes equations and the magnetic and electric fields in the Maxwell's equations. In addition, the methods studied here are meshfree, and are suitable for problems defined on complex domains, where mesh generation is computationally expensive or inaccurate, or for problems where the data is only available at scattered locations.
The contributions of this work include a detailed comparison between standard and divergence-free radial basis approximations, a study of the Lebesgue constants for divergence-free approximations and their dependence on node placement, and an investigation of the flat limit of divergence-free interpolants. Finally, numerical solvers for the incompressible Navier-Stokes equations in primitive variables are implemented using discretizations based on traditional and divergence-free kernels. The numerical results are compared to reference solutions obtained with a spectral
method.
The contributions of this work include a detailed comparison between standard and divergence-free radial basis approximations, a study of the Lebesgue constants for divergence-free approximations and their dependence on node placement, and an investigation of the flat limit of divergence-free interpolants. Finally, numerical solvers for the incompressible Navier-Stokes equations in primitive variables are implemented using discretizations based on traditional and divergence-free kernels. The numerical results are compared to reference solutions obtained with a spectral
method.
ContributorsAraujo Mitrano, Arthur (Author) / Platte, Rodrigo (Thesis advisor) / Wright, Grady (Committee member) / Welfert, Bruno (Committee member) / Gelb, Anne (Committee member) / Renaut, Rosemary (Committee member) / Arizona State University (Publisher)
Created2016