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

Jingrun Chen

Jingrun Chen contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisphysics.comp-phComputational Engineering, Finance, and Sciencemath.AP

Trust snapshot

Quick read

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

11 published item(s)

preprint2022arXiv

Bridging Traditional and Machine Learning-based Algorithms for Solving PDEs: The Random Feature Method

One of the oldest and most studied subject in scientific computing is algorithms for solving partial differential equations (PDEs). A long list of numerical methods have been proposed and successfully used for various applications. In recent years, deep learning methods have shown their superiority for high-dimensional PDEs where traditional methods fail. However, for low dimensional problems, it remains unclear whether these methods have a real advantage over traditional algorithms as a direct solver. In this work, we propose the random feature method (RFM) for solving PDEs, a natural bridge between traditional and machine learning-based algorithms. RFM is based on a combination of well-known ideas: 1. representation of the approximate solution using random feature functions; 2. collocation method to take care of the PDE; 3. the penalty method to treat the boundary conditions, which allows us to treat the boundary condition and the PDE in the same footing. We find it crucial to add several additional components including multi-scale representation and rescaling the weights in the loss function. We demonstrate that the method exhibits spectral accuracy and can compete with traditional solvers in terms of both accuracy and efficiency. In addition, we find that RFM is particularly suited for complex problems with complex geometry, where both traditional and machine learning-based algorithms encounter difficulties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Convergence analysis of an implicit finite difference method for the inertial Landau-Lifshitz-Gilbert equation

The Landau-Lifshitz-Gilbert (LLG) equation is a widely used model for fast magnetization dynamics in ferromagnetic materials. Recently, the inertial LLG equation, which contains an inertial term, has been proposed to capture the ultra-fast magnetization dynamics at the sub-picosecond timescale. Mathematically, this generalized model contains the first temporal derivative and a newly introduced second temporal derivative of magnetization. Consequently, it produces extra difficulties in numerical analysis due to the mixed hyperbolic-parabolic type of this equation with degeneracy. In this work, we propose an implicit finite difference scheme based on the central difference in both time and space. A fixed point iteration method is applied to solve the implicit nonlinear system. With the help of a second order accurate constructed solution, we provide a convergence analysis in $H^1$ for this numerical scheme, in the $\ell^\infty (0, T; H_h^1)$ norm. It is shown that the proposed method is second order accurate in both time and space, with unconditional stability and a natural preservation of the magnetization length. In the hyperbolic regime, significant damping wave behaviors of magnetization at a shorter timescale are observed through numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Homogenization of the Landau-Lifshitz-Gilbert equation with natural boundary condition

The full Landau-Lifshitz-Gilbert equation with periodic material coefficients and natural boundary condition is employed to model the magnetization dynamics in composite ferromagnets. In this work, we establish the convergence between the homogenized solution and the original solution via a Lax equivalence theorem kind of argument. There are a few technical difficulties, including: 1) it is proven the classic choice of corrector to homogenization cannot provide the convergence result in the $H^1$ norm; 2) a boundary layer is induced due to the natural boundary condition; 3) the presence of stray field give rise to a multiscale potential problem. To keep the convergence rates near the boundary, we introduce the Neumann corrector with a high-order modification. Estimates on singular integral for disturbed functions and boundary layer are deduced, to conduct consistency analysis of stray field. Furthermore, inspired by length conservation of magnetization, we choose proper correctors in specific geometric space. These, together with a uniform $W^{1,6}$ estimate on original solution, provide the convergence rates in the $H^1$ sense.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Gauss-Seidel projection method with the minimal number of updates for stray field in micromagnetic simulations

Magnetization dynamics in magnetic materials is often modeled by the Landau-Lifshitz equation, which is solved numerically in general. In micromagnetic simulations, the computational cost relies heavily on the time-marching scheme and the evaluation of stray field. Explicit marching schemes are efficient but suffer from severe stability constraints, while nonlinear systems of equations have to be solved in implicit schemes though they are unconditionally stable. A better compromise between stability and efficiency is the semi-implicit scheme, such as the Gauss-Seidel projection method (GSPM) and the second-order backward differentiation formula scheme (BDF2). At each marching step, GSPM solves several linear systems of equations with constant coefficients and updates the stray field several times, while BDF2 updates the stray field only once but solves a larger linear system of equations with variable coefficients and a nonsymmetric structure. In this work, we propose a new method, dubbed as GSPM-BDF2, by combing the advantages of both GSPM and BDF2. Like GSPM, this method is first-order accurate in time and second-order accurate in space, and is unconditionally stable with respect to the damping parameter. However, GSPM-BDF2 updates the stray field only once per time step, leading to an efficiency improvement of about $60\%$ than the state-of-the-art GSPM for micromagnetic simulations. For Standard Problem \#4 and \#5 from National Institute of Standards and Technology, GSPM-BDF2 reduces the computational time over the popular software OOMMF by $82\%$ and $96\%$, respectively. Thus, the proposed method provides a more efficient choice for micromagnetic simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A second-order numerical method for Landau-Lifshitz-Gilbert equation with large damping parameters

A second order accurate numerical scheme is proposed and implemented for the Landau-Lifshitz-Gilbert equation, which models magnetization dynamics in ferromagnetic materials, with large damping parameters. The main advantages of this method are associated with the following features: (1) It only solves linear systems of equations with constant coefficients where fast solvers are available, so that the numerical efficiency has been greatly improved, in comparison with the existing Gauss-Seidel project method. (2) The second-order accuracy in time is achieved, and it is unconditionally stable for large damping parameters. Moreover, both the second-order accuracy and the great efficiency improvement will be verified by several numerical examples in the 1D and 3D simulations. In the presence of large damping parameters, it is observed that this method is unconditionally stable and finds physically reasonable structures while many existing methods have failed. For the domain wall dynamics, the linear dependence of wall velocity with respect to the damping parameter and the external magnetic field will be obtained through the reported simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A comparison study of deep Galerkin method and deep Ritz method for elliptic problems with different boundary conditions

Recent years have witnessed growing interests in solving partial differential equations by deep neural networks, especially in the high-dimensional case. Unlike classical numerical methods, such as finite difference method and finite element method, the enforcement of boundary conditions in deep neural networks is highly nontrivial. One general strategy is to use the penalty method. In the work, we conduct a comparison study for elliptic problems with four different boundary conditions, i.e., Dirichlet, Neumann, Robin, and periodic boundary conditions, using two representative methods: deep Galerkin method and deep Ritz method. In the former, the PDE residual is minimized in the least-squares sense while the corresponding variational problem is minimized in the latter. Therefore, it is reasonably expected that deep Galerkin method works better for smooth solutions while deep Ritz method works better for low-regularity solutions. However, by a number of examples, we observe that deep Ritz method can outperform deep Galerkin method with a clear dependence of dimensionality even for smooth solutions and deep Galerkin method can also outperform deep Ritz method for low-regularity solutions. Besides, in some cases, when the boundary condition can be implemented in an exact manner, we find that such a strategy not only provides a better approximate solution but also facilitates the training process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Enforcing exact boundary and initial conditions in the deep mixed residual method

In theory, boundary and initial conditions are important for the wellposedness of partial differential equations (PDEs). Numerically, these conditions can be enforced exactly in classical numerical methods, such as finite difference method and finite element method. Recent years have witnessed growing interests in solving PDEs by deep neural networks (DNNs), especially in the high-dimensional case. However, in the generic situation, a careful literature review shows that boundary conditions cannot be enforced exactly for DNNs, which inevitably leads to a modeling error. In this work, based on the recently developed deep mixed residual method (MIM), we demonstrate how to make DNNs satisfy boundary and initial conditions automatically in a systematic manner. As a consequence, the loss function in MIM is free of the penalty term and does not have any modeling error. Using numerous examples, including Dirichlet, Neumann, mixed, Robin, and periodic boundary conditions for elliptic equations, and initial conditions for parabolic and hyperbolic equations, we show that enforcing exact boundary and initial conditions not only provides a better approximate solution but also facilitates the training process.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

MIM: A deep mixed residual method for solving high-order partial differential equations

In recent years, a significant amount of attention has been paid to solve partial differential equations (PDEs) by deep learning. For example, deep Galerkin method (DGM) uses the PDE residual in the least-squares sense as the loss function and a deep neural network (DNN) to approximate the PDE solution. In this work, we propose a deep mixed residual method (MIM) to solve PDEs with high-order derivatives. In MIM, we first rewrite a high-order PDE into a first-order system, very much in the same spirit as local discontinuous Galerkin method and mixed finite element method in classical numerical methods for PDEs. We then use the residual of first-order system in the least-squares sense as the loss function, which is in close connection with least-squares finite element method. For aforementioned classical numerical methods, the choice of trail and test functions is important for stability and accuracy issues in many cases. MIM shares this property when DNNs are employed to approximate unknowns functions in the first-order system. In one case, we use nearly the same DNN to approximate all unknown functions and in the other case, we use totally different DNNs for different unknown functions. In most cases, MIM provides better approximations (not only for high-derivatives of the PDE solution but also for the PDE solution itself) than DGM with nearly the same DNN and the same execution time, sometimes by more than one order of magnitude. When different DNNs are used, in many cases, MIM provides even better approximations than MIM with only one DNN, sometimes by more than one order of magnitude. Therefore, we expect MIM to open up a possibly systematic way to understand and improve deep learning for solving PDEs from the perspective of classical numerical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Numerical methods for antiferromagnetics

Compared with ferromagnetic counterparts, antiferromagnetic materials are considered as the future of spintronic applications since these materials are robust against the magnetic perturbation, produce no stray field, and display ultrafast dynamics. There are (at least) two sets of magnetic moments in antiferromagnets (with magnetization of the same magnitude but antiparallel directions) and ferrimagnets (with magnetization of the different magnitude). The coupled dynamics for the bipartite collinear antiferromagnets is modeled by a coupled system of Landau-Lifshitz-Gilbert equations with an additional term originated from the antiferromagnetic exchange, which leads to femtosecond magnetization dynamics. In this paper, we develop three Gauss-Seidel projection methods for micromagnetics simulation in antiferromagnets and ferrimagnets. They are first-order accurate in time and second-order in space, and only solve linear systems of equations with constant coefficients at each step. Femtosecond dynamics, Néel wall structure, and phase transition in presence of an external magnetic field for antiferromagnets are provided with the femtosecond stepsize.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Second-order semi-implicit projection methods for micromagnetics simulations

Micromagnetics simulations require accurate approximation of the magnetization dynamics described by the Landau-Lifshitz-Gilbert equation, which is nonlinear, nonlocal, and has a non-convex constraint, posing interesting challenges in developing numerical methods. In this paper, we propose two second-order semi-implicit projection methods based on the second-order backward differentiation formula and the second-order interpolation formula using the information at previous two temporal steps. Unconditional unique solvability of both methods is proved, with their second-order accuracy verified through numerical examples in both 1D and 3D. The efficiency of both methods is compared to that of another two popular methods. In addition, we test the robustness of both methods for the first benchmark problem with a ferromagnetic thin film material from National Institute of Standards and Technology.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Two improved Gauss-Seidel projection methods for Landau-Lifshitz-Gilbert equation

In this paper, we present two improved Gauss-Seidel projection methods with unconditional stability. The first method updates the gyromagnetic term and the damping term simultaneously and follows by a projection step. The second method introduces two sets of approximate solutions, where we update the gyromagnetic term and the damping term simultaneously for one set of approximate solutions and apply the projection step to the other set of approximate solutions in an alternating manner. Compared to the original Gauss-Seidel projection method which has to solve heat equations $7$ times at each time step, the improved methods solve heat equations $5$ times and $3$ times, respectively. First-order accuracy in time and second-order accuracy in space are verified by examples in both 1D and 3D. In addition, unconditional stability with respect to both the grid size and the damping parameter is confirmed numerically. Application of both methods to a realistic material is also presented with hysteresis loops and magnetization profiles. Compared with the original method, the recorded running times suggest that savings of both methods are about $2/7$ and $4/7$ for the same accuracy requirement, respectively.

0 citations0 reviewsNo code linkedOpen access