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

Weizhu Bao

Weizhu Bao contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-scimath.NANumerical Analysisphysics.comp-phcond-mat.softphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
9works
0followers
6topics
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

9 published item(s)

preprint2022arXiv

Diffuse-interface approach to competition between viscous flow and diffusion in pinch-off dynamics

The pinch-off dynamics of a liquid thread has been studied through numerical simulations and theoretical analysis. Occurring at small length scales, the pinch-off dynamics admits similarity solutions that can be classified into the Stokes regime and the diffusion-dominated regime, with the latter being recently experimentally observed in aqueous two-phase systems [Phys. Rev. Lett. 123, 134501 (2019)]. Derived by applying Onsager's variational principle, the Cahn-Hilliard-Navier-Stokes model is employed as a minimal model capable of describing the interfacial motion driven by not only advection but also diffusion. By analyzing the free energy dissipation mechanisms in the model, a characteristic length scale is introduced to measure the competition between diffusion and viscous flow in interfacial motion. This length scale is typically of nanometer scale for systems far from the critical point, but can approach micrometer scale for aqueous two-phase systems close to the critical point. The Cahn-Hilliard-Navier-Stokes model is solved by using an accurate and efficient spectral method in a cylindrical domain with axisymmetry. Ample numerical examples are presented to show the pinch-off processes in the Stokes regime and the diffusion-dominated regime respectively. In particular, the crossover between these two regimes is investigated numerically and analytically to reveal how the scaling behaviors of similarity solutions are to be qualitatively changed as the characteristic length scale is inevitably accessed by the pinching neck of the interface. Discussions are also provided for numerical examples that are performed for the breakup of long liquid filaments and show qualitatively different phenomena in different scaling regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved uniform error bounds of the time-splitting methods for the long-time (nonlinear) Schrödinger equation

We establish improved uniform error bounds for the time-splitting methods for the long-time dynamics of the Schrödinger equation with small potential and the nonlinear Schrödinger equation (NLSE) with weak nonlinearity. For the Schrödinger equation with small potential characterized by a dimensionless parameter $\varepsilon \in (0, 1]$ representing the amplitude of the potential, we employ the unitary flow property of the (second-order) time-splitting Fourier pseudospectral (TSFP) method in $L^2$-norm to prove a uniform error bound at $C(T)(h^m +τ^2)$ up to the long time $T_\varepsilon= T/\varepsilon$ for any $T>0$ and uniformly for $0<\varepsilon\le1$, while $h$ is the mesh size, $τ$ is the time step, $m \ge 2$ depends on the regularity of the exact solution, and $C(T) =C_0+C_1T$ grows at most linearly with respect to $T$ with $C_0$ and $C_1$ two positive constants independent of $T$, $\varepsilon$, $h$ and $τ$. Then by introducing a new technique of {\sl regularity compensation oscillation} (RCO) in which the high frequency modes are controlled by regularity and the low frequency modes are analyzed by phase cancellation and energy method, an improved uniform error bound at $O(h^{m-1} + \varepsilon τ^2)$ is established in $H^1$-norm for the long-time dynamics up to the time at $O(1/\varepsilon)$ of the Schrödinger equation with $O(\varepsilon)$-potential with $m \geq 3$, which is uniformly for $\varepsilon\in(0,1]$. Moreover, the RCO technique is extended to prove an improved uniform error bound at $O(h^{m-1} + \varepsilon^2τ^2)$ in $H^1$-norm for the long-time dynamics up to the time at $O(1/\varepsilon^2)$ of the cubic NLSE with $O(\varepsilon^2)$-nonlinearity strength, uniformly for $\varepsilon \in (0, 1]$. Extensions to the first-order and fourth-order time-splitting methods are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved uniform error bounds on time-splitting methods for the long-time dynamics of the weakly nonlinear Dirac equation

Improved uniform error bounds on time-splitting methods are rigorously proven for the long-time dynamics of the weakly nonlinear Dirac equation (NLDE), where the nonlinearity strength is characterized by a dimensionless parameter $\varepsilon \in (0, 1]$ . We adopt a second order Strang splitting method to discretize the NLDE in time and combine the Fourier pseudospectral method in space for the full-discretization. By employing the {\sl regularity compensation oscillation} (RCO) technique where the high frequency modes are controlled by the regularity of the exact solution and the low frequency modes are analyzed by phase cancellation and energy method, we establish improved uniform error bounds at $O(\varepsilon^2τ^2)$ and $O(h^{m-1}+ \varepsilon^2τ^2)$ for the second-order Strang splitting semi-discretizaion and full-discretization up to the long-time $T_{\varepsilon} = T/\varepsilon^2$ with $T>0$ fixed, respectively. Furthermore, the numerical scheme and error estimates are extended to an oscillatory NLDE which propagates waves with $O(\varepsilon^2)$ wavelength in time. Finally, numerical examples verifying our analytical results are given.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Volume-preserving parametric finite element methods for axisymmetric geometric evolution equations

We propose and analyze volume-preserving parametric finite element methods for surface diffusion, conserved mean curvature flow and an intermediate evolution law in an axisymmetric setting. The weak formulations are presented in terms of the generating curves of the axisymmetric surfaces. The proposed numerical methods are based on piecewise linear parametric finite elements. The constructed fully practical schemes satisfy the conservation of the enclosed volume. In addition, we prove the unconditional stability and consider the distribution of vertices for the discretized schemes. The introduced methods are implicit and the resulting nonlinear systems of equations can be solved very efficiently and accurately via the Newton&#39;s iterative method. Numerical results are presented to show the accuracy and efficiency of the introduced schemes for computing the considered axisymmetric geometric flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An energy-stable parametric finite element method for simulating solid-state dewetting

We propose an energy-stable parametric finite element method (ES-PFEM) for simulating solid-state dewetting of thin films in two dimensions via a sharp-interface model, which is governed by surface diffusion and contact line (point) migration together with proper boundary conditions. By reformulating the relaxed contact angle condition into a Robin-type boundary condition and then treating it as a natural boundary condition, we obtain a new variational formulation for the problem, in which the interface curve and its contact points are evolved simultaneously. Then, the variational problem is discretized in space by using piecewise linear elements. A full discretization is presented by adopting the backward Euler method in time, and the well-posedness and energy dissipation of the full discretization are established. The numerical method is semi-implicit (i.e., a linear system to be solved at each time step and thus efficient), unconditionally energy-stable with respect to the time step, and second-order in space measured by a manifold distance between two curves. In addition, it demonstrates equal mesh distribution when the solution reaches its equilibrium, i.e., long-time dynamics. Numerical results are reported to show accuracy and efficiency as well as some good properties of the proposed numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Power-law scaling for solid-state dewetting of thin films: an Onsager variational approach

We examine the kinetics of surface diffusion-controlled, solid-state dewetting by consideration of the retraction of the contact in a semi-infinite solid thin film on a flat rigid substrate. The analysis is performed within the framework of the Onsager variational principle applied to surface diffusion-controlled morphology evolution. Based on this approach, we derive a simple, reduced-order model to quantitatively analyse the power-law scaling of the dewetting process. Using asymptotic analysis and numerical simulations for the reduced-order model, we find that the retraction distance grows as the $2/5$ power of time and the height of the ridge, adjacent to the contact, grows as the $1/5$ power of time for late time. While the asymptotic analysis focuses on late time and a relatively simple geometric model, the Onsager approach is applicable to all times and descriptions of the morphology of arbitrary complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sharp-interface model for simulating solid-state dewetting in three dimensions

The problem of simulating solid-state dewetting of thin films in three dimensions (3D) by using a sharp-interface approach is considered in this paper. Based on the thermodynamic variation, a speed method is used for calculating the first variation to the total surface energy functional. The speed method shares more advantages than the traditional use of parameterized curves (or surfaces), e.g., it is more intrinsic and its variational structure (related with Cahn-Hoffman $\boldsymbolξ$-vector) is clearer and more direct. By making use of the first variation, necessary conditions for the equilibrium shape of the solid-state dewetting problem is given, and a kinetic sharp-interface model which includes the surface energy anisotropy is also proposed. This sharp-interface model describes the interface evolution in 3D which occurs through surface diffusion and contact line migration. By solving the proposed model, we perform lots of numerical simulations to investigate the evolution of patterned films, e.g., the evolution of a short cuboid and pinch-off of a long cuboid. Numerical simulations in 3D demonstrate the accuracy and efficacy of the sharp-interface approach to capture many of the complexities observed in solid-state dewetting experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A parametric finite element method for solid-state dewetting problems in three dimensions

We propose a parametric finite element method (PFEM) for efficiently solving the morphological evolution of solid-state dewetting of thin films on a flat rigid substrate in three dimensions (3D). The interface evolution of the dewetting problem in 3D is described by a sharp-interface model, which includes surface diffusion coupled with contact line migration. A variational formulation of the sharp-interface model is presented, and a PFEM is proposed for spatial discretization. For temporal discretization, at each time step, we first update the position of the contact line according to the relaxed contact angle condition; then, by using the position of the new contact line as the boundary condition, we solve a linear algebra system resulted from the discretization of PFEM to obtain the new interface surface for the next step. The well-posedness of the solution of the PFEM is also established. Extensive numerical results are reported to demonstrate the accuracy and efficiency of the proposed PFEM and to show the complexities of the dewetting morphology evolution observed in solid-state dewetting experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Computing Ground States of Spin-1 Bose-Einstein Condensates by the Normalized Gradient Flow

In this paper, we propose an efficient and accurate numerical method for computing the ground state of spin-1 Bose-Einstein condensates (BEC) by using the normalized gradient flow or imaginary time method. The key idea is to find a third projection or normalization condition based on the relation between the chemical potentials so that the three projection parameters used in the projection step of the normalized gradient flow are uniquely determined by this condition as well as the other two physical conditions given by the conservation of total mass and total magnetization. This allows us to successfully extend the most popular and powerful normalized gradient flow or imaginary time method for computing the ground state of single component BEC to compute the ground state of spin-1 BEC. An efficient and accurate discretization scheme, the backward-forward Euler sine-pseudospectral method (BFSP), is proposed to discretize the normalized gradient flow. Extensive numerical results on ground states of spin-1 BEC with ferromagnetic/antiferromagnetic interaction and harmonic/optical lattice potential in one/three dimensions are reported to demonstrate the efficiency of our new numerical method.

0 citations0 reviewsNo code linkedOpen access