BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Thomas Y. Hou

Thomas Y. Hou contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.APmath.NANumerical AnalysisInformation Theorymath.ITmath.OCphysics.flu-dynBiological PhysicsCell BehaviorMachine Learning

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

10 published item(s)

preprint2026arXiv

$L^2$-based stability of blowup with log correction for semilinear heat equation

We propose an alternative proof of the classical result of Type-I blowup with log correction for the semilinear heat equation. Compared with previous proofs, we use a novel idea of enforcing stable normalizations for perturbations around the approximate profile and we establish a weighted $H^k$ stability, thereby avoiding the use of a topological argument and the analysis of a linearized spectrum. Consequently, this approach can be adopted even if we only have a numerical profile and do not have explicit information on the spectrum of its linearized operator. This result generalizes the $L^2$-based stability framework beyond exactly self-similar blowup and can be adapted to higher dimensions. Numerical results corroborate the effectiveness of our normalization, even in the large perturbation regime beyond our theoretical setting.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Asymptotic Escape of Spurious Critical Points on the Low-rank Matrix Manifold

We show that on the manifold of fixed-rank and symmetric positive semi-definite matrices, the Riemannian gradient descent algorithm almost surely escapes some spurious critical points on the boundary of the manifold. Our result is the first to partially overcome the incompleteness of the low-rank matrix manifold without changing the vanilla Riemannian gradient descent algorithm. The spurious critical points are some rank-deficient matrices that capture only part of the eigen components of the ground truth. Unlike classical strict saddle points, they exhibit very singular behavior. We show that using the dynamical low-rank approximation and a rescaled gradient flow, some of the spurious critical points can be converted to classical strict saddle points in the parameterized domain, which leads to the desired result. Numerical experiments are provided to support our theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Finite time blowup of 2D Boussinesq and 3D Euler equations with $C^{1,α}$ velocity and boundary

Inspired by the numerical evidence of a potential 3D Euler singularity by Luo-Hou [30,31] and the recent breakthrough by Elgindi [11] on the singularity formation of the 3D Euler equation without swirl with $C^{1,α}$ initial velocity, we prove the finite time singularity for the 2D Boussinesq and the 3D axisymmetric Euler equations with boundary and $C^{1,α}$ initial data for the velocity (and density in the case of Boussinesq equations). Our finite time blowup solution for the 3D Euler equations and the singular solution considered in [30,31] share many essential features, including the symmetry properties of the solution, the flow structure, and the sign of the solution in each quadrant, except that we use $C^{1,α}$ initial velocity. We use a dynamic rescaling formulation and follow the general framework of analysis developed in [11]. We also use some strategy proposed in our recent joint work with Huang in [7] and adopt several methods of analysis in [11] to establish the linear and nonlinear stability of an approximate self-similar profile. The nonlinear stability enables us to prove that the solution of the 3D Euler or the 2D Boussinesq equations with $C^{1,α}$ initial data will develop a finite time singularity. Moreover, the velocity field has finite energy before the singularity time. In the previous version of this paper, we proved the blowup results for the 3D axisymmetric Euler equations with initial data $(u_0^θ)^2, u_0^r, u_0^z \in C^{1,α}$. Though the velocity $u^r, u^z$ in the axisymmetric setting is $C^{1,α}$, our interpretation that the velocity is $C^{1,α}$ is not correct since the velocity in 3D also depends on $u^θ$, which is not $C^{1,α}$. This oversight can be fixed easily with a minor change in the construction of the approximate steady state and a minor modification to localize the approximate steady state.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On stability and instability of $C^{1,α}$ singular solutions to the 3D Euler and 2D Boussinesq equations

Singularity formation of the 3D incompressible Euler equations is known to be extremely challenging. In [18], Elgindi proved that the 3D axisymmetric Euler equations with no swirl and $C^{1,α}$ initial velocity develops a finite time singularity. Inspired by Elgindi's work, we proved that the 3D axisymmetric Euler and 2D Boussinesq equations with $C^{1,α}$ initial velocity and boundary develop a stable asymptotically (or approximately) self-similar finite time singularity [8]. On the other hand, the authors of [35,52] recently showed that blowup solutions to the 3D Euler equations are hydrodynamically unstable. The instability results obtained in [35,52] require some strong regularity assumption on the initial data, which is not satisfied by the $C^{1,α}$ velocity field. In this paper, we generalize the analysis of [8,18,35,52] to show that the blowup solutions of the 3D Euler and 2D Boussinesq equations with $C^{1,α}$ velocity are unstable under the notion of stability introduced in [35,52]. These two seemingly contradictory results reflect the difference of the two approaches in studying the stability of 3D Euler blowup solutions. The stability analysis of the blowup solution obtained in [8,18] is based on the stability of a dynamically rescaled blowup profile in space and time, which is nonlinear in nature. The linear stability analysis in [35,52] is performed by directly linearizing the 3D Euler equations around a blowup solution in the original variables. It does not take into account the changes in the blowup time, the dynamic changes of the rescaling rate of the perturbed blowup profile and the blowup exponent of the original 3D Euler equations using a perturbed initial condition when there is an approximate self-similar blowup profile. Such information has been used in an essential way in establishing the nonlinear stability of the blowup profile in [8,18,19].

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Potential singularity formation of incompressible axisymmetric Euler equations with degenerate viscosity coefficients

In this paper, we present strong numerical evidences that the incompressible axisymmetric Euler equations with degenerate viscosity coefficients and smooth initial data of finite energy develop a potential finite-time locally self-similar singularity at the origin. An important feature of this potential singularity is that the solution develops a two-scale traveling wave that travels towards the origin. The two-scale feature is characterized by the scaling property that the center of the traveling wave is located at a ring of radius $O((T-t)^{1/2})$ surrounding the symmetry axis while the thickness of the ring collapses at a rate $O(T-t)$. The driving mechanism for this potential singularity is due to an antisymmetric vortex dipole that generates a strong shearing layer in both the radial and axial velocity fields. Without the viscous regularization, the $3$D Euler equations develop a sharp front and some shearing instability in the far field. On the other hand, the Navier-Stokes equations with a constant viscosity coefficient regularize the two-scale solution structure and do not develop a finite-time singularity for the same initial data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Potential singularity of the 3D Euler equations in the interior domain

Whether the 3D incompressible Euler equations can develop a finite time singularity from smooth initial data is one of the most challenging problems in nonlinear PDEs. In this paper, we present some new numerical evidence that the 3D axisymmetric incompressible Euler equations with smooth initial data of finite energy develop a potential finite time singularity at the origin. This potential singularity is different from the blow-up scenario revealed by Luo-Hou in \cite{luo2014potentially,luo2014toward}, which occurs on the boundary. Our initial condition has a simple form and shares several attractive features of a more sophisticated initial condition constructed by Hou-Huang in \cite{Hou-Huang-2021,Hou-Huang-2022}. One important difference between these two blow-up scenarios is that the solution for our initial data has a one-scale structure instead of a two-scale structure reported in \cite{Hou-Huang-2021,Hou-Huang-2022}. More importantly, the solution seems to develop nearly self-similar scaling properties that are compatible with those of the 3D Navier-Stokes equations. We will present numerical evidence that the 3D Euler equations seem to develop a potential finite time singularity. Moreover, the nearly self-similar profile seems to be very stable to the small perturbation of the initial data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The potentially singular behavior of the 3D Navier-Stokes equations

Whether the 3D incompressible Navier-Stokes equations can develop a finite time singularity from smooth initial data is one of the most challenging problems in nonlinear PDEs. In this paper, we present some new numerical evidence that the incompressible axisymmetric Navier-Stokes equations with smooth initial data of finite energy seem to develop potentially singular behavior at the origin. This potentially singular behavior is induced by a potential finite time singularity of the 3D Euler equations that we reported in the companion paper (arXiv:2107.05870). We present numerical evidence that the 3D Navier--Stokes equations develop nearly self-similar singular scaling properties with maximum vorticity increased by a factor of 10^7. We have applied several blow-up criteria to study the potentially singular behavior of the Navier--Stokes equations. The Beale-Kato-Majda blow-up criterion and the blow-up criteria based on the growth of enstrophy and negative pressure seem to imply that the Navier--Stokes equations using our initial data develop a potential finite time singularity. We have also examined the Ladyzhenskaya-Prodi-Serrin regularity criteria. Our numerical results for the cases of (p,q) = (4,8), (6,4), (9,3) and (p,q)=(\infty,2) provide strong evidence for the potentially singular behavior of the Navier--Stokes equations. Our numerical study shows that while the global L^3 norm of the velocity grows very slowly, the localized version of the L^3 norm of the velocity experiences rapid dynamic growth relative to the localized L^3 norm of the initial velocity. This provides further evidence for the potentially singular behavior of the NavieStokes equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Analysis of Asymptotic Escape of Strict Saddle Sets in Manifold Optimization

In this paper, we provide some analysis on the asymptotic escape of strict saddles in manifold optimization using the projected gradient descent (PGD) algorithm. One of our main contributions is that we extend the current analysis to include non-isolated and possibly continuous saddle sets with complicated geometry. We prove that the PGD is able to escape strict critical submanifolds under certain conditions on the geometry and the distribution of the saddle point sets. We also show that the PGD may fail to escape strict saddles under weaker assumptions even if the saddle point set has zero measure and there is a uniform escape direction. We provide a counterexample to illustrate this important point. We apply this saddle analysis to the phase retrieval problem on the low-rank matrix manifold, prove that there are only a finite number of saddles, and they are strict saddles with high probability. We also show the potential application of our analysis for a broader range of manifold optimization problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Function Approximation via The Subsampled Poincar\' e Inequality

Function approximation and recovery via some sampled data have long been studied in a wide array of applied mathematics and statistics fields. Analytic tools, such as the Poincaré inequality, have been handy for estimating the approximation errors in different scales. The purpose of this paper is to study a generalized Poincar\' e inequality, where the measurement function is of subsampled type, with a small but non-zero lengthscale that will be made precise. Our analysis identifies this inequality as a basic tool for function recovery problems. We discuss and demonstrate the optimality of the inequality concerning the subsampled lengthscale, connecting it to existing results in the literature. In application to function approximation problems, the approximation accuracy using different basis functions and under different regularity assumptions is established by using the subsampled Poincaré inequality. We observe that the error bound blows up as the subsampled lengthscale approaches zero, due to the fact that the underlying function is not regular enough to have well-defined pointwise values. A weighted version of the Poincar\' e inequality is proposed to address this problem; its optimality is also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

A Minimal Mechanosensing Model Predicts Keratocyte Evolution on Flexible Substrates

A mathematical model is proposed for shape evolution and locomotion of fish epidermal keratocytes on elastic substrates. The model is based on mechanosensing concepts: cells apply contractile forces onto the elastic substrate, while cell shape evolution depends locally on the substrate stress generated by themselves or external mechanical stimuli acting on the substrate. We use the level set method to study the behavior of the model numerically, and predict a number of distinct phenomena observed in experiments, such as (i) symmetry breaking from the stationary centrosymmetric to the well-known steadily propagating crescent shape, (ii) asymmetric bipedal oscillations and traveling waves in the lamellipodium leading edge (iii) response to mechanical stress externally applied to the substrate (tensotaxis), (iv) changing direction of motion towards an interface with a rigid substrate (durotaxis) and (v) the configuration of substrate wrinkles induced by contractile forces applied by the keratocyte.

0 citations0 reviewsNo code linkedOpen access