Matching Items (6)
Filtering by
- All Subjects: Applied Mathematics
![189270-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2023-08/189270-Thumbnail%20Image.png?versionId=zeEDzs8d9nl_7oEBJn8D5VnxAPk3b6tu&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240618/us-west-2/s3/aws4_request&X-Amz-Date=20240618T191639Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=d35dbdf46407c5fc01446479f32dc0749350080e999f5ad6c84f5f55c69fb060&itok=rvcqiz1p)
Description
The variable projection method has been developed as a powerful tool for solvingseparable nonlinear least squares problems. It has proven effective in cases where
the underlying model consists of a linear combination of nonlinear functions, such as
exponential functions. In this thesis, a modified version of the variable projection
method to address a challenging semi-blind deconvolution problem involving mixed
Gaussian kernels is employed. The aim is to recover the original signal accurately
while estimating the mixed Gaussian kernel utilized during the convolution process.
The numerical results obtained through the implementation of the proposed algo-
rithm are presented. These results highlight the method’s ability to approximate the
true signal successfully. However, accurately estimating the mixed Gaussian kernel
remains a challenging task. The implementation details, specifically focusing on con-
structing a simplified Jacobian for the Gauss-Newton method, are explored. This
contribution enhances the understanding and practicality of the approach.
ContributorsDworaczyk, Jordan Taylor (Author) / Espanol, Malena (Thesis advisor) / Welfert, Bruno (Committee member) / Platte, Rodrigo (Committee member) / Arizona State University (Publisher)
Created2023
![157240-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/157240-Thumbnail%20Image.png?versionId=MXEUT_fSxAIVod2l.QJzJyR3hA2uL78h&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240618/us-west-2/s3/aws4_request&X-Amz-Date=20240618T191639Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=7630a9c376f19e7eb2800251284559e148d78ff4c0b37a76d9d1078379995bd9&itok=VNuTCsYa)
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
![154804-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/154804-Thumbnail%20Image.png?versionId=Md1XKPA3eVGiv6B.17FWK87.tFZQxKBR&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240618/us-west-2/s3/aws4_request&X-Amz-Date=20240618T191639Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=c166453e932f60284b9afbaaa2ff941be964208e8fa63db0bcaa100947e5b88d&itok=YAZCwYIE)
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
![155789-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/155789-Thumbnail%20Image.png?versionId=jqf8qC0gHxWpZR8ZPEUp00xzvasd0otx&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240618/us-west-2/s3/aws4_request&X-Amz-Date=20240618T191639Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=f1f05a97862763801ec2d63f813fae9c6fce66230080f751e5f531165b9d4fa3&itok=KBomkXNt)
Description
Functional magnetic resonance imaging (fMRI) is one of the popular tools to study human brain functions. High-quality experimental designs are crucial to the success of fMRI experiments as they allow the collection of informative data for making precise and valid inference with minimum cost. The primary goal of this study is on identifying the best sequence of mental stimuli (i.e. fMRI design) with respect to some statistically meaningful optimality criteria. This work focuses on two related topics in this research field. The first topic is on finding optimal designs for fMRI when the design matrix is uncertain. This challenging design issue occurs in many modern fMRI experiments, in which the design matrix of the statistical model depends on both the selected design and the experimental subject's uncertain behavior during the experiment. As a result, the design matrix cannot be fully determined at the design stage that makes it difficult to select a good design. For the commonly used linear model with autoregressive errors, this study proposes a very efficient approach for obtaining high-quality fMRI designs for such experiments. The proposed approach is built upon an analytical result, and an efficient computer algorithm. It is shown through case studies that our proposed approach can outperform the existing method in terms of computing time, and the quality of the obtained designs. The second topic of the research is to find optimal designs for fMRI when a wavelet-based technique is considered in the fMRI data analysis. An efficient computer algorithm to search for optimal fMRI designs for such cases is developed. This algorithm is inspired by simulated annealing and a recently proposed algorithm by Saleh et al. (2017). As demonstrated in the case studies, the proposed approach makes it possible to efficiently obtain high-quality designs for fMRI studies, and is practically useful.
ContributorsZhou, Lin (Author) / Kao, Ming-Hung (Thesis advisor) / Welfert, Bruno (Thesis advisor) / Jackiewicz, Zdzislaw (Committee member) / Reiser, Mark R. (Committee member) / Stufken, John (Committee member) / Taylor, Jesse Earl (Committee member) / Arizona State University (Publisher)
Created2017
![156420-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-09/156420-Thumbnail%20Image.png?versionId=MHNyO1Ca_BPTHVygm3USdkV_QzM5Hv4Z&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240618/us-west-2/s3/aws4_request&X-Amz-Date=20240618T191639Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=c9b3961c113b8dce35e48a268abf462a93afbc9d936566e62066c11575eeaf04&itok=Y-OLOaVu)
Description
The Kuramoto model is an archetypal model for studying synchronization in groups
of nonidentical oscillators where oscillators are imbued with their own frequency and
coupled with other oscillators though a network of interactions. As the coupling
strength increases, there is a bifurcation to complete synchronization where all oscillators
move with the same frequency and show a collective rhythm. Kuramoto-like
dynamics are considered a relevant model for instabilities of the AC-power grid which
operates in synchrony under standard conditions but exhibits, in a state of failure,
segmentation of the grid into desynchronized clusters.
In this dissertation the minimum coupling strength required to ensure total frequency
synchronization in a Kuramoto system, called the critical coupling, is investigated.
For coupling strength below the critical coupling, clusters of oscillators form
where oscillators within a cluster are on average oscillating with the same long-term
frequency. A unified order parameter based approach is developed to create approximations
of the critical coupling. Some of the new approximations provide strict lower
bounds for the critical coupling. In addition, these approximations allow for predictions
of the partially synchronized clusters that emerge in the bifurcation from the
synchronized state.
Merging the order parameter approach with graph theoretical concepts leads to a
characterization of this bifurcation as a weighted graph partitioning problem on an
arbitrary networks which then leads to an optimization problem that can efficiently
estimate the partially synchronized clusters. Numerical experiments on random Kuramoto
systems show the high accuracy of these methods. An interpretation of the
methods in the context of power systems is provided.
of nonidentical oscillators where oscillators are imbued with their own frequency and
coupled with other oscillators though a network of interactions. As the coupling
strength increases, there is a bifurcation to complete synchronization where all oscillators
move with the same frequency and show a collective rhythm. Kuramoto-like
dynamics are considered a relevant model for instabilities of the AC-power grid which
operates in synchrony under standard conditions but exhibits, in a state of failure,
segmentation of the grid into desynchronized clusters.
In this dissertation the minimum coupling strength required to ensure total frequency
synchronization in a Kuramoto system, called the critical coupling, is investigated.
For coupling strength below the critical coupling, clusters of oscillators form
where oscillators within a cluster are on average oscillating with the same long-term
frequency. A unified order parameter based approach is developed to create approximations
of the critical coupling. Some of the new approximations provide strict lower
bounds for the critical coupling. In addition, these approximations allow for predictions
of the partially synchronized clusters that emerge in the bifurcation from the
synchronized state.
Merging the order parameter approach with graph theoretical concepts leads to a
characterization of this bifurcation as a weighted graph partitioning problem on an
arbitrary networks which then leads to an optimization problem that can efficiently
estimate the partially synchronized clusters. Numerical experiments on random Kuramoto
systems show the high accuracy of these methods. An interpretation of the
methods in the context of power systems is provided.
ContributorsGilg, Brady (Author) / Armbruster, Dieter (Thesis advisor) / Mittelmann, Hans (Committee member) / Scaglione, Anna (Committee member) / Strogatz, Steven (Committee member) / Welfert, Bruno (Committee member) / Arizona State University (Publisher)
Created2018
![155323-Thumbnail Image.png](https://d1rbsgppyrdqq4.cloudfront.net/s3fs-public/styles/width_400/public/2021-08/155323-Thumbnail%20Image.png?versionId=LUy__sn120b4klL1h_yg.8qhz0xnMxIT&X-Amz-Content-Sha256=UNSIGNED-PAYLOAD&X-Amz-Algorithm=AWS4-HMAC-SHA256&X-Amz-Credential=AKIASBVQ3ZQ42ZLA5CUJ/20240618/us-west-2/s3/aws4_request&X-Amz-Date=20240618T191639Z&X-Amz-SignedHeaders=host&X-Amz-Expires=120&X-Amz-Signature=05904d2e5ca531334442026df8ac91603910a20a2a3e3095f09ea1e33652277c&itok=tLwnGe7x)
Description
The three-dimensional flow contained in a rapidly rotating circular
split cylinder is studied numerically solving the Navier--Stokes
equations. The cylinder is completely filled with fluid
and is split at the midplane. Three different types of boundary
conditions were imposed, leading to a variety of instabilities and
complex flow dynamics.
The first configuration has a strong background rotation and a small
differential rotation between the two halves. The axisymmetric flow
was first studied identifying boundary layer instabilities which
produce inertial waves under some conditions. Limit cycle states and
quasiperiodic states were found, including some period doubling
bifurcations. Then, a three-dimensional study was conducted
identifying low and high azimuthal wavenumber rotating waves due to
G’ortler and Tollmien–-Schlichting type instabilities. Over most of
the parameter space considered, quasiperiodic states were found where
both types of instabilities were present.
In the second configuration, both cylinder halves are in exact
counter-rotation, producing an O(2) symmetry in the system. The basic state flow dynamic
is dominated by the shear layer created
in the midplane. By changing the speed rotation and the aspect ratio
of the cylinder, the flow loses symmetries in a variety of ways
creating static waves, rotating waves, direction reversing waves and
slow-fast pulsing waves. The bifurcations, including infinite-period
bifurcations, were characterized and the flow dynamics was elucidated.
Additionally, preliminary experimental results for this case are
presented.
In the third set up, with oscillatory boundary conditions, inertial
wave beams were forced imposing a range of frequencies. These beams
emanate from the corner of the cylinder and from the split at the
midplane, leading to destructive/constructive interactions which
produce peaks in vorticity for some specific frequencies. These
frequencies are shown to be associated with the resonant Kelvin
modes. Furthermore, a study of the influence of imposing a phase
difference between the oscillations of the two halves of the cylinder
led to the interesting result that different Kelvin
modes can be excited depending on the phase difference.
split cylinder is studied numerically solving the Navier--Stokes
equations. The cylinder is completely filled with fluid
and is split at the midplane. Three different types of boundary
conditions were imposed, leading to a variety of instabilities and
complex flow dynamics.
The first configuration has a strong background rotation and a small
differential rotation between the two halves. The axisymmetric flow
was first studied identifying boundary layer instabilities which
produce inertial waves under some conditions. Limit cycle states and
quasiperiodic states were found, including some period doubling
bifurcations. Then, a three-dimensional study was conducted
identifying low and high azimuthal wavenumber rotating waves due to
G’ortler and Tollmien–-Schlichting type instabilities. Over most of
the parameter space considered, quasiperiodic states were found where
both types of instabilities were present.
In the second configuration, both cylinder halves are in exact
counter-rotation, producing an O(2) symmetry in the system. The basic state flow dynamic
is dominated by the shear layer created
in the midplane. By changing the speed rotation and the aspect ratio
of the cylinder, the flow loses symmetries in a variety of ways
creating static waves, rotating waves, direction reversing waves and
slow-fast pulsing waves. The bifurcations, including infinite-period
bifurcations, were characterized and the flow dynamics was elucidated.
Additionally, preliminary experimental results for this case are
presented.
In the third set up, with oscillatory boundary conditions, inertial
wave beams were forced imposing a range of frequencies. These beams
emanate from the corner of the cylinder and from the split at the
midplane, leading to destructive/constructive interactions which
produce peaks in vorticity for some specific frequencies. These
frequencies are shown to be associated with the resonant Kelvin
modes. Furthermore, a study of the influence of imposing a phase
difference between the oscillations of the two halves of the cylinder
led to the interesting result that different Kelvin
modes can be excited depending on the phase difference.
ContributorsGutierrez Castillo, Paloma (Author) / Lopez, Juan M. (Thesis advisor) / Herrmann, Marcus (Committee member) / Platte, Rodrigo (Committee member) / Welfert, Bruno (Committee member) / Tang, Wenbo (Committee member) / Arizona State University (Publisher)
Created2017