Source author record

Chiu-Yen Kao

Chiu-Yen Kao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA math.OC math.SP math.AP math.DG nlin.PS Computational Geometry math.PR nlin.AO

Catalog footprint

What is connected

8works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Flat tori with large Laplacian eigenvalues in dimensions up to eight

We consider the optimization problem of maximizing the $k$-th Laplacian eigenvalue, $λ_{k}$, over flat $d$-dimensional tori of fixed volume. For $k=1$, this problem is equivalent to the densest lattice sphere packing problem. For larger $k$, this is equivalent to the NP-hard problem of finding the $d$-dimensional (dual) lattice with longest $k$-th shortest lattice vector. As a result of extensive computations, for $d \leq 8$, we obtain a sequence of flat tori, $T_{k,d}$, each of volume one, such that the $k$-th Laplacian eigenvalue of $T_{k,d}$ is very large; for each (finite) $k$ the $k$-th eigenvalue exceeds the value in (the $k\to \infty$ asymptotic) Weyl's law by a factor between 1.54 and 2.01, depending on the dimension. Stationarity conditions are derived and numerically verified for $T_{k,d}$ and we describe the degeneration of the tori as $k \to \infty$.

preprint2020arXiv

Computation of free boundary minimal surfaces via extremal Steklov eigenvalue problems

Recently Fraser and Schoen showed that the solution of a certain extremal Steklov eigenvalue problem on a compact surface with boundary can be used to generate a free boundary minimal surface, i.e., a surface contained in the ball that has (i) zero mean curvature and (ii) meets the boundary of the ball orthogonally (doi:10.1007/s00222-015-0604-x). In this paper, we develop numerical methods that use this connection to realize free boundary minimal surfaces. Namely, on a compact surface, $Σ$, with genus $γ$ and $b$ boundary components, we maximize $σ_j(Σ,g) \ L(\partial Σ, g)$ over a class of smooth metrics, $g$, where $σ_j(Σ,g)$ is the $j$-th nonzero Steklov eigenvalue and $L(\partial Σ, g)$ is the length of $\partial Σ$. Our numerical method involves (i) using conformal uniformization of multiply connected domains to avoid explicit parameterization for the class of metrics, (ii) accurately solving a boundary-weighted Steklov eigenvalue problem in multi-connected domains, and (iii) developing gradient-based optimization methods for this non-smooth eigenvalue optimization problem. For genus $γ=0$ and $b=2,\dots, 9, 12, 15, 20$ boundary components, we numerically solve the extremal Steklov problem for the first eigenvalue. The corresponding eigenfunctions generate a free boundary minimal surface, which we display in striking images. For higher eigenvalues, numerical evidence suggests that the maximizers are degenerate, but we compute local maximizers for the second and third eigenvalues with $b=2$ boundary components and for the third and fifth eigenvalues with $b=3$ boundary components.

preprint2016arXiv

A New Algorithm to Simulate the First Exit Times of a Vector of Brownian Motions, with an Application to Finance

We provide a new methodology to simulate the first exit times of a vector of Brownian motions from an orthant. This new approach can be used to simulate the first exit times of dimension higher than two. When at least one Brownian motion has non-zero drift, the joint density function of the first exit times in N dimensions needs to be known, or approximated. However, when the drifts are all zero, a simpler simulation algorithm is obtained without using the joint density function.

preprint2016arXiv

Maximization of Laplace-Beltrami eigenvalues on closed Riemannian surfaces

Let $(M,g)$ be a connected, closed, orientable Riemannian surface and denote by $λ_k(M,g)$ the $k$-th eigenvalue of the Laplace-Beltrami operator on $(M,g)$. In this paper, we consider the mapping $(M, g)\mapsto λ_k(M,g)$. We propose a computational method for finding the conformal spectrum $Λ^c_k(M,[g_0])$, which is defined by the eigenvalue optimization problem of maximizing $λ_k(M,g)$ for $k$ fixed as $g$ varies within a conformal class $[g_0]$ of fixed volume $textrm{vol}(M,g) = 1$. We also propose a computational method for the problem where $M$ is additionally allowed to vary over surfaces with fixed genus, $γ$. This is known as the topological spectrum for genus $γ$ and denoted by $Λ^t_k(γ)$. Our computations support a conjecture of N. Nadirashvili (2002) that $Λ^t_k(0) = 8 πk$, attained by a sequence of surfaces degenerating to a union of $k$ identical round spheres. Furthermore, based on our computations, we conjecture that $Λ^t_k(1) = \frac{8π^2}{\sqrt{3}} + 8π(k-1)$, attained by a sequence of surfaces degenerating into a union of an equilateral flat torus and $k-1$ identical round spheres. The values are compared to several surfaces where the Laplace-Beltrami eigenvalues are well-known, including spheres, flat tori, and embedded tori. In particular, we show that among flat tori of volume one, the $k$-th Laplace-Beltrami eigenvalue has a local maximum with value $λ_k = 4π^2 \left\lceil \frac{k}{2} \right\rceil^2 \left( \left\lceil \frac{k}{2} \right\rceil^2 - \frac{1}{4}\right)^{-\frac{1}{2}}$. Several properties are also studied computationally, including uniqueness, symmetry, and eigenvalue multiplicity.

preprint2013arXiv

Convergent finite difference methods for one-dimensional fully nonlinear second order partial differential equations

This paper develops a new framework for designing and analyzing convergent finite difference methods for approximating both classical and viscosity solutions of second order fully nonlinear partial differential equations (PDEs) in 1-D. The goal of the paper is to extend the successful framework of monotone, consistent, and stable finite difference methods for first order fully nonlinear Hamilton-Jacobi equations to second order fully nonlinear PDEs such as Monge-Ampère and Bellman type equations. New concepts of consistency, generalized monotonicity, and stability are introduced; among them, the generalized monotonicity and consistency, which are easier to verify in practice, are natural extensions of the corresponding notions of finite difference methods for first order fully nonlinear Hamilton-Jacobi equations. The main component of the proposed framework is the concept of "numerical operator", and the main idea used to design consistent, monotone and stable finite difference methods is the concept of "numerical moment". These two new concepts play the same roles as the "numerical Hamiltonian" and the "numerical viscosity" play in the finite difference framework for first order fully nonlinear Hamilton-Jacobi equations. In the paper, two classes of consistent and monotone finite difference methods are proposed for second order fully nonlinear PDEs. The first class contains Lax-Friedrichs-like methods which also are proved to be stable and the second class contains Godunov-like methods. Numerical results are also presented to gauge the performance of the proposed finite difference methods and to validate the theoretical results of the paper.

preprint2011arXiv

Bounded domain problem for the modified Buckley-Leverett equation

The focus of the present study is the modified Buckley-Leverett (MBL) equation describing two-phase flow in porous media. The MBL equation differs from the classical Buckley-Leverett (BL) equation by including a balanced diffusive-dispersive combination. The dispersive term is a third order mixed derivatives term, which models the dynamic effects in the pressure difference between the two phases. The classical BL equation gives a monotone water saturation profile for any Riemann problem; on the contrast, when the dispersive parameter is large enough, the MBL equation delivers non-monotone water saturation profile for certain Riemann problems as suggested by the experimental observations. In this paper, we first show that the solution of the finite interval [0,L] boundary value problem converges to that of the half-line [0,+\infty) boundary value problem for the MBL equation as L-> +\infty. This result provides a justification for the use of the finite interval boundary value problem in numerical studies for the half line problem. Furthermore, we extend the classical central schemes for the hyperbolic conservation laws to solve the MBL equation which is of pseudo-parabolic type. Numerical results confirm the existence of non-monotone water saturation profiles consisting of constant states separated by shocks.

preprint2010arXiv

Numerical study of the KP equation for non-periodic waves

The Kadomtsev-Petviashvili (KP) equation describes weakly dispersive and small amplitude waves propagating in a quasi-two dimensional situation. Recently a large variety of exact soliton solutions of the KP equation has been found and classified. Those soliton solutions are localized along certain lines in a two-dimensional plane and decay exponentially everywhere else, and they are called line-soliton solutions in this paper. The classification is based on the far-field patterns of the solutions which consist of a finite number of line-solitons. In this paper, we study the initial value problem of the KP equation with V- and X-shape initial waves consisting of two distinct line-solitons by means of the direct numerical simulation. We then show that the solution converges asymptotically to some of those exact soliton solutions. The convergence is in a locally defined $L^2$-sense. The initial wave patterns considered in this paper are related to the rogue waves generated by nonlinear wave interactions in shallow water wave problem.

preprint2009arXiv

A moving boundary model motivated by electric breakdown: II. Initial value problem

An interfacial approximation of the streamer stage in the evolution of sparks and lightning can be formulated as a Laplacian growth model regularized by a 'kinetic undercooling' boundary condition. Using this model we study both the linearized and the full nonlinear evolution of small perturbations of a uniformly translating circle. Within the linear approximation analytical and numerical results show that perturbations are advected to the back of the circle, where they decay. An initially analytic interface stays analytic for all finite times, but singularities from outside the physical region approach the interface for $t\to\infty$, which results in some anomalous relaxation at the back of the circle. For the nonlinear evolution numerical results indicate that the circle is the asymptotic attractor for small perturbations, but larger perturbations may lead to branching. We also present results for more general initial shapes, which demonstrate that regularization by kinetic undercooling cannot guarantee smooth interfaces globally in time.