BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Suchuan Dong

Suchuan Dong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynmath.NAphysics.comp-phMachine LearningNumerical AnalysisComputational Engineering, Finance, and Science

Catalog footprint

What is connected

8works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

8 published item(s)

preprint2022arXiv

Numerical Approximation of Partial Differential Equations by a Variable Projection Method with Artificial Neural Networks

We present a method for solving linear and nonlinear PDEs based on the variable projection (VarPro) framework and artificial neural networks (ANN). For linear PDEs, enforcing the boundary/initial value problem on the collocation points leads to a separable nonlinear least squares problem about the network coefficients. We reformulate this problem by the VarPro approach to eliminate the linear output-layer coefficients, leading to a reduced problem about the hidden-layer coefficients only. The reduced problem is solved first by the nonlinear least squares method to determine the hidden-layer coefficients, and then the output-layer coefficients are computed by the linear least squares method. For nonlinear PDEs, enforcing the boundary/initial value problem on the collocation points leads to a nonlinear least squares problem that is not separable, which precludes the VarPro strategy for such problems. To enable the VarPro approach for nonlinear PDEs, we first linearize the problem with a Newton iteration, using a particular form of linearization. The linearized system is solved by the VarPro framework together with ANNs. Upon convergence of the Newton iteration, the network coefficients provide the representation of the solution field to the original nonlinear problem. We present ample numerical examples with linear and nonlinear PDEs to demonstrate the performance of the method herein. For smooth field solutions, the errors of the current method decrease exponentially as the number of collocation points or the number of output-layer coefficients increases. We compare the current method with the ELM method from a previous work. Under identical conditions and network configurations, the current method exhibits an accuracy significantly superior to the ELM method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical Computation of Partial Differential Equations by Hidden-Layer Concatenated Extreme Learning Machine

The extreme learning machine (ELM) method can yield highly accurate solutions to linear/nonlinear partial differential equations (PDEs), but requires the last hidden layer of the neural network to be wide to achieve a high accuracy. If the last hidden layer is narrow, the accuracy of the existing ELM method will be poor, irrespective of the rest of the network configuration. In this paper we present a modified ELM method, termed HLConcELM (hidden-layer concatenated ELM), to overcome the above drawback of the conventional ELM method. The HLConcELM method can produce highly accurate solutions to linear/nonlinear PDEs when the last hidden layer of the network is narrow and when it is wide. The new method is based on a type of modified feedforward neural networks (FNN), termed HLConcFNN (hidden-layer concatenated FNN), which incorporates a logical concatenation of the hidden layers in the network and exposes all the hidden nodes to the output-layer nodes. HLConcFNNs have the interesting property that, given a network architecture, when additional hidden layers are appended to the network or when extra nodes are added to the existing hidden layers the representation capacity of the HLConcFNN associated with the new architecture is guaranteed to be not smaller than that of the original network architecture. Here representation capacity refers to the set of all functions that can be exactly represented by the neural network of a given architecture. We present ample benchmark tests with linear/nonlinear PDEs to demonstrate the computational accuracy and performance of the HLConcELM method and the superiority of this method to the conventional ELM from previous works.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Marangoni Convection-Driven Laser Fountains and Waves on Free Surfaces of Liquids

It is well accepted that an outward Marangoni convection from a low surface tension region will make the surface depressed. Here, we report that this established perception is only valid for thin liquid films. Using surface laser heating, we show that in deep liquids a laser beam actually pulls up the fluid above the free surface generating fountains with different shapes. Whereas with decreasing liquid depth a transition from fountain to indentation with fountain in-indentation is observed. Further, high-speed imaging reveals a transient surface process before steady elevation is formed, and this dynamic deformation is subsequently utilized to resonantly excite giant surface waves by a modulated laser beam. Computational fluid dynamics models reveal the underlying flow patterns and quantify the depth-dependent and time-resolved surface deformations. Our discoveries and techniques have upended the century-old perception and opened up a new regime of interdisciplinary research and applications of Marangoni-induced interface phenomena and optocapillary fluidic surfaces-the control of fluids with light.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On Computing the Hyperparameter of Extreme Learning Machines: Algorithm and Application to Computational PDEs, and Comparison with Classical and High-Order Finite Elements

We consider the use of extreme learning machines (ELM) for computational partial differential equations (PDE). In ELM the hidden-layer coefficients in the neural network are assigned to random values generated on $[-R_m,R_m]$ and fixed, where $R_m$ is a user-provided constant, and the output-layer coefficients are trained by a linear or nonlinear least squares computation. We present a method for computing the optimal value of $R_m$ based on the differential evolution algorithm. The presented method enables us to illuminate the characteristics of the optimal $R_m$ for two types of ELM configurations: (i) Single-Rm-ELM, in which a single $R_m$ is used for generating the random coefficients in all the hidden layers, and (ii) Multi-Rm-ELM, in which multiple $R_m$ constants are involved with each used for generating the random coefficients of a different hidden layer. We adopt the optimal $R_m$ from this method and also incorporate other improvements into the ELM implementation. In particular, here we compute all the differential operators involving the output fields of the last hidden layer by a forward-mode auto-differentiation, as opposed to the reverse-mode auto-differentiation in a previous work. These improvements significantly reduce the network training time and enhance the ELM performance. We systematically compare the computational performance of the current improved ELM with that of the finite element method (FEM), both the classical second-order FEM and the high-order FEM with Lagrange elements of higher degrees, for solving a number of linear and nonlinear PDEs. It is shown that the current improved ELM far outperforms the classical FEM. Its computational performance is comparable to that of the high-order FEM for smaller problem sizes, and for larger problem sizes the ELM markedly outperforms the high-order FEM.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Roadmap for Discretely Energy-Stable Schemes for Dissipative Systems Based on a Generalized Auxiliary Variable with Guaranteed Positivity

We present a framework for devising discretely energy-stable schemes for general dissipative systems based on a generalized auxiliary variable. The auxiliary variable, a scalar number, can be defined in terms of the energy functional by a general class of functions, not limited to the square root function adopted in previous approaches. The current method has another remarkable property: the computed values for the generalized auxiliary variable are guaranteed to be positive on the discrete level, regardless of the time step sizes or the external forces. This property of guaranteed positivity is not available in previous approaches. A unified procedure for treating the dissipative governing equations and the generalized auxiliary variable on the discrete level has been presented. The discrete energy stability of the proposed numerical scheme and the positivity of the computed auxiliary variable have been proved for general dissipative systems. The current method, termed gPAV (generalized Positive Auxiliary Variable), requires only the solution of linear algebraic equations within a time step. With appropriate choice of the operator in the algorithm, the resultant linear algebraic systems upon discretization involve only constant and time-independent coefficient matrices, which only need to be computed once and can be pre-computed. Several specific dissipative systems are studied in relative detail using the gPAV framework. Ample numerical experiments are presented to demonstrate the performance of the method, and the robustness of the scheme at large time step sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Convective-like Energy-Stable Open Boundary Condition for Simulations of Incompressible Flows

We present a new energy-stable open boundary condition, and an associated numerical algorithm, for simulating incompressible flows with outflow/open boundaries. This open boundary condition ensures the energy stability of the system, even when strong vortices or backflows occur at the outflow boundary. Under certain situations it can be reduced to a form that can be analogized to the usual convective boundary condition. One prominent feature of this boundary condition is that it provides a control over the velocity on the outflow/open boundary. This is not available with the other energy-stable open boundary conditions from previous works. Our numerical algorithm treats the proposed open boundary condition based on a rotational velocity-correction type strategy. It gives rise to a Robin-type condition for the discrete pressure and a Robin-type condition for the discrete velocity on the outflow/open boundary, respectively at the pressure and the velocity sub-steps. We present extensive numerical experiments on a canonical wake flow and a jet flow in open domain to test the effectiveness and performance of the method developed herein. Simulation results are compared with the experimental data as well as with other previous simulations to demonstrate the accuracy of the current method. Long-time simulations are performed for a range of Reynolds numbers, at which strong vortices and backflows occur at the outflow/open boundaries. The results show that our method is effective in overcoming the backflow instability, and that it allows for the vortices to discharge from the domain in a fairly natural fashion even at high Reynolds numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A pressure correction scheme for generalized form of energy-stable open boundary conditions for incompressible flows

We present a generalized form of open boundary conditions, and an associated numerical algorithm, for simulating incompressible flows involving open or outflow boundaries. The generalized form represents a family of open boundary conditions, which all ensure the energy stability of the system, even in situations where strong vortices or backflows occur at the open/outflow boundaries. Our numerical algorithm for treating these open boundary conditions is based on a rotational pressure correction-type strategy, with a formulation suitable for $C^0$ spectral-element spatial discretizations. We have introduced a discrete equation and associated boundary conditions for an auxiliary variable. The algorithm contains constructions that prevent a numerical locking at the open/outflow boundary. In addition, we have also developed a scheme with a provable unconditional stability for a sub-class of the open boundary conditions. Extensive numerical experiments have been presented to demonstrate the performance of our method for several flow problems involving open/outflow boundaries. We compare simulation results with the experimental data to demonstrate the accuracy of our algorithm. Long-time simulations have been performed for a range of Reynolds numbers at which strong vortices or backflows occur at the open/outflow boundaries. We show that the open boundary conditions and the numerical algorithm developed herein produce stable simulations in such situations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Physical Formulation and Numerical Algorithm for Simulating N Immiscible Incompressible Fluids Involving General Order Parameters

We present a physical formulation, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of $N$ ($N\geqslant 2$) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The introduction of general order parameters leads to a more strongly coupled system of phase field equations, in contrast to that with certain special choice of the order parameters. However, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. From the simulation perspective, the increased complexity in the form of the phase field equations with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the ($N-1$) strongly-coupled phase field equations for general order parameters into $2(N-1)$ Helmholtz-type equations that are completely de-coupled from one another, leading to a computational complexity essentially the same as that of the simpler phase field equations associated with special choice of order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir-de Gennes theory of floating liquid lenses we show that the method produces physically accurate results for multiple fluid phases.

0 citations0 reviewsNo code linkedOpen access