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

Quanling Deng

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

ResearcherUnclaimed source record
math.NANumerical AnalysisComputational Engineering, Finance, and Sciencemath.DSphysics.comp-phphysics.data-anquant-ph

Catalog footprint

What is connected

8works
7topics
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

A boundary-penalized isogeometric analysis for second-order hyperbolic equations

Explicit time-marching schemes are popular for solving time-dependent partial differential equations; one of the biggest challenges these methods suffer is increasing the critical time-marching step size that guarantees numerical stability. In general, there are two ways to increase the critical step size. One is to reduce the stiffness of the spatially discretized system, while the other is to design time-marching schemes with larger stability regions. In this paper, we focus on the recently proposed explicit generalized-$α$ method for second-order hyperbolic equations and increase the critical step size by reducing the stiffness of the isogeometric-discretized system. In particular, we apply boundary penalization to lessen the system's stiffness. For $p$-th order $C^{p-1}$ isogeometric elements, we show numerically that the critical step size increases by a factor of $\sqrt{\frac{p^2-3p+6}{4}}$, which indicates the advantages of using the proposed method, especially for high-order elements. Various examples in one, two, and three dimensions validate the performance of the proposed technique.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Isogeometric Analysis of Bound States of a Quantum Three-Body Problem in 1D

In this paper, we initiate the study of isogeometric analysis (IGA) of a quantum three-body problem that has been well-known to be difficult to solve. In the IGA setting, we represent the wavefunctions by linear combinations of B-spline basis functions and solve the problem as a matrix eigenvalue problem. The eigenvalue gives the eigenstate energy while the eigenvector gives the coefficients of the B-splines that lead to the eigenstate. The major difficulty of isogeometric or other finite-element-method-based analyses lies in the lack of boundary conditions and a large number of degrees of freedom for accuracy. For a typical many-body problem with attractive interaction, there are bound and scattering states where bound states have negative eigenvalues. We focus on bound states and start with the analysis for a two-body problem. We demonstrate through various numerical experiments that IGA provides a promising technique to solve the three-body problem.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Superfloe Parameterization with Physics Constraints for Uncertainty Quantification of Sea Ice Floes

The discrete element method (DEM) is providing a new modeling approach for describing sea ice dynamics. It exploits particle-based methods to characterize the physical quantities of each sea ice floe along its trajectory under Lagrangian coordinates. One major challenge in applying the DEM models is the heavy computational cost when the number of floes becomes large. In this paper, an efficient Lagrangian parameterization algorithm is developed, which aims at reducing the computational cost of simulating the DEM models while preserving the key features of the sea ice. The new parameterization takes advantage of a small number of artificial ice floes, named the superfloes, to effectively approximate a considerable number of the floes, where the parameterization scheme satisfies several important physics constraints. The physics constraints guarantee the superfloe parameterized system will have similar short-term dynamical behavior as the full system. These constraints also allow the superfloe parameterized system to accurately quantify the long-range uncertainty, especially the non-Gaussian statistical features, of the full system. In addition, the superfloe parameterization facilitates a systematic noise inflation strategy that significantly advances an ensemble-based data assimilation algorithm for recovering the unobserved ocean field underneath the sea ice. Such a new noise inflation method avoids ad hoc tunings as in many traditional algorithms and is computationally extremely efficient. Numerical experiments based on an idealized DEM model with multiscale features illustrate the success of the superfloe parameterization in quantifying the uncertainty and assimilating both the sea ice and the associated ocean field.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Higher-order generalized-$α$ methods for parabolic problems

We propose a new class of high-order time-marching schemes with dissipation user-control and unconditional stability for parabolic equations. High-order time integrators can deliver the optimal performance of highly-accurate and robust spatial discretizations such as isogeometric analysis. The generalized-$α$ method delivers unconditional stability and second-order accuracy in time and controls the numerical dissipation in the discrete spectrum's high-frequency region. Our goal is to extend the generalized-$alpha$ methodology to obtain a high-order time marching methods with high accuracy and dissipation in the discrete high-frequency range. Furthermore, we maintain the stability region of the original, second-order generalized-$alpha$ method foe the new higher-order methods. That is, we increase the accuracy of the generalized-$α$ method while keeping the unconditional stability and user-control features on the high-frequency numerical dissipation. The methodology solve $k>1, k\in \mathbb{N}$ matrix problems and updates the system unknowns, which correspond to higher-order terms in Taylor expansions to obtain $(3/2k)^{th}$-order method for even $k$ and $(3/2k+1/2)^{th}$-order for odd $k$. A single parameter $ρ^\infty$ controls the dissipation, and the update procedure follows the formulation of the original second-order method. Additionally, we show that our method is A-stable and setting $ρ^\infty=0$ allows us to obtain an L-stable method. Lastly, we extend this strategy to analyze the accuracy order of a generic method.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the Application of Stable Generalized Finite Element Method for Quasilinear Elliptic Two-Point BVP

In this paper, we discuss the application of the Generalized Finite Element Method (GFEM) to approximate the solutions of quasilinear elliptic equations with multiple interfaces in one dimensional space. The problem is characterized by spatial discontinuity of the elliptic coefficient that depends on the unknown solution. It is known that unless the partition of the domain matches the discontinuity configuration, accuracy of standard finite element techniques significantly deteriorates and standard refinement of the partition may not suffice. The GFEM is a viable alternative to overcome this predicament. It is based on the construction of certain enrichment functions supplied to the standard space that capture effects of the discontinuity. This approach is called stable (SGFEM) if it maintains an optimal rate of convergence and the conditioning of GFEM is not worse than that of the standard FEM. A convergence analysis is derived and performance of the method is illustrated by several numerical examples. Furthermore, it is known that typical global formulations such as FEMs do not enjoy the numerical local conservation property that is crucial in many conservation law-based applications. To remedy this issue, a Lagrange multiplier technique is adopted to enforce the local conservation. A numerical example is given to demonstrate the performance of proposed technique.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Construction of Locally Conservative Fluxes for High Order Continuous Galerkin Finite Element Methods

We propose a simple post-processing technique for linear and high order continuous Galerkin Finite Element Methods (CGFEMs) to obtain locally conservative flux field. The post-processing technique requires solving an auxiliary problem on each element independently which results in solving a linear algebra system whose size is low for any order CGFEM. The post-processing could have been done directly from the finite element solution that results in locally conservative flux on the element. However, the normal flux is not continuous at the elemental boundary. To construct locally conservative flux field whose normal component is also continuous, we propose to do the post-processing on the nodal-centered control volumes which are constructed from the original finite element mesh. We show that the post-processed solution converges in an optimal fashion to the true solution in an H1 semi-norm. We present various numerical examples to demonstrate the performance of the post-processing technique.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Locally Conservative Continuous Galerkin FEM for Pressure Equation in Two-Phase Flow Model in Subsurfaces

A typical two-phase model for subsurface flow couples the Darcy equation for pressure and a transport equation for saturation in a nonlinear manner. In this paper, we study a combined method consisting of continuous Galerkin finite element methods (CGFEMs) followed by a post-processing technique for Darcy equation and finite volume method (FVM) with upwind schemes for the saturation transport equation, in which the coupled nonlinear problem is solved in the framework of operator decomposition. The postprocessing technique is applied to CGFEM solutions to obtain locally conservative fluxes which ensures accuracy and robustness of the FVM solver for the saturation transport equation. We applied both upwind scheme and upwind scheme with slope limiter for FVM on triangular meshes in order to eliminate the non-physical oscillations. Various numerical examples are presented to demonstrate the performance of the overall methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Post-processing Technique for Streamline Upwind/Petrov-Galerkin for Advection Dominated Partial Differential Equations

We consider the construction of locally conservative fluxes by means of a simple post-processing technique obtained from the finite element solutions of advection diffusion equations. It is known that a naive calculation of fluxes from these solutions yields non-conservative fluxes. We consider two finite element methods: the usual continuous Galerkin finite element (CGFEM) for solving non dominating advection diffusion equations and the streamline upwind/Petrov-Galerkin (SUPG) for solving advection dominated problems. We then describe the post-processing technique for constructing conservative fluxes from the numerical solutions of the general variational formulation. The post-processing technique requires solving an auxiliary Neumann boundary value problem on each element independently and it produces a locally conservative flux on a vertex centered dual mesh relative to the finite element mesh. We provide a convergence analysis for the post-processing technique. Performance of the technique and the convergence behavior are demonstrated through numerical examples including a set of test problems for advection diffusion equations, advection dominated equations, and drift-diffusion equations.

0 citations0 reviewsNo code linkedOpen access