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

Bjorn Engquist

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

ResearcherUnclaimed source record
math.NAmath-phmath.MPNumerical Analysiscond-mat.mtrl-scimath.APmath.FAphysics.geo-ph

Catalog footprint

What is connected

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

11 published item(s)

preprint2019arXiv

Fast algorithm for computing nonlocal operators with finite interaction distance

Developments of nonlocal operators for modeling processes that traditionally have been described by local differential operators have been increasingly active during the last few years. One example is peridynamics for brittle materials and another is nonstandard diffusion including the use of fractional derivatives. A major obstacle for application of these methods is the high computational cost from the numerical implementation of the nonlocal operators. It is natural to consider fast methods of fast multipole or hierarchical matrix type to overcome this challenge. Unfortunately the relevant kernels do not satisfy the standard necessary conditions. In this work a new class of fast algorithms is developed and analyzed, which is some cases reduces the computational complexity of applying nonlocal operators to essentially the same order of magnitude as the complexity of standard local numerical methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A mathematical and Physical Model Improves Accuracy in Simulating Solid Material Relaxation Modulus and Viscoelastic Responses

We propose a new material viscoelastic model and mathematical solution to simulate relaxation modulus and viscoelastic response. The model formula of relaxation modulus is extended from sigmoidal function considering nonlinear strain hardening and softening. Its physical mechanism can be interpreted by a spring network viscous medium model with only five parameters in a simpler format than the molecular-chain based polymer models to represent general materials. We also developed a three-dimensional finite-element method and robust numerical algorithms to implement this model for solving partial differential equations. We validate the model through both experimental data and numerical simulations on a broad range of materials including bitumen, shape-memory polymer, spider-inspired silk, hydrogel, biomaterials and bone. By satisfying the 2nd law of thermodynamics in the form of Calusius-Duhem inequality, the model is able to simulate creep and sinusoidal deformation, and energy dissipation. As compared to Prony series, the most general model being used often with a large number of model parameters, the proposed model has improved accuracy in fitting experimental data and predicting modulus outside of the experimental range, and the latter one is especially useful for material design. The new model also has higher numerical accuracy while competitive numerical stability and computation speed for convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Homogenization Model for Aberrant Crypt Foci

Several explanations can be found in the literature about the origin of colorectal cancer. There is however some agreement on the fact that the carcinogenic process is a result of several genetic mutations of normal cells. The colon epithelium is characterized by millions of invaginations, very small cavities, called crypts, where most of the cellular activity occurs. It is consensual in the medical community, that a potential first manifestation of the carcinogenic process, observed in conventional colonoscopy images, is the appearance of Aberrant Crypt Foci (ACF). These are clusters of abnormal crypts, morphologically characterized by an atypical behavior of the cells that populate the crypts. In this work an homogenization model is proposed, for representing the cellular dynamics in the colon epithelium. The goal is to simulate and predict, in silico, the spread and evolution of ACF, as it can be observed in colonoscopy images. By assuming that the colon is an heterogeneous media, exhibiting a periodic distribution of crypts, we start this work by describing a periodic model, that represents the ACF cell-dynamics in a two-dimensional setting. Then, homogenization techniques are applied to this periodic model, to find a simpler model, whose solution symbolizes the averaged behavior of ACF at the tissue level. Some theoretical results concerning the existence of solution of the homogenized model are proven, applying a fixed point theorem. Numerical results showing the convergence of the periodic model to the homogenized model are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Numerical methods for multiscale inverse problems

We consider the inverse problem of determining the highly oscillatory coefficient $a^ε$ in partial differential equations of the form $-\nabla\cdot (a^ε\nabla u^ε)+bu^ε= f$ from given measurements of the solutions. Here, $ε$ indicates the smallest characteristic wavelength in the problem ($0<ε\ll1$). In addition to the general difficulty of finding an inverse, the oscillatory nature of the forward problem creates an additional challenge of multiscale modeling, which is hard even for forward computations. The inverse problem in its full generality is typically ill-posed and one common approach is to replace the original problem with an effective parameter estimation problem. We will here include microscale features directly in the inverse problem and avoid ill-posedness by assuming that the microscale can be accurately represented by a low-dimensional parametrization. The basis for our inversion will be a coupling of the parametrization to analytic homogenization or a coupling to efficient multiscale numerical methods when analytic homogenization is not available. We will analyze the reduced problem, $b = 0$, by proving uniqueness of the inverse in certain problem classes and by numerical examples and also include numerical model examples for medical imaging, $b > 0$, and exploration seismology, $b < 0$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optimal Transport for Seismic Full Waveform Inversion

Full waveform inversion is a successful procedure for determining properties of the earth from surface measurements in seismology. This inverse problem is solved by a PDE constrained optimization where unknown coefficients in a computed wavefield are adjusted to minimize the mismatch with the measured data. We propose using the Wasserstein metric, which is related to optimal transport, for measuring this mismatch. Several advantageous properties are proved with regards to convexity of the objective function and robustness with respect to noise. The Wasserstein metric is computed by solving a Monge-Ampere equation. We describe an algorithm for computing its Frechet gradient for use in the optimization. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast integrators for dynamical systems with several temporal scales

We propose a fast integrator to a class of dynamical systems with several temporal scales. The proposed method is developed as an extension of the variable step size Heterogeneous Multiscale Method (VSHMM), which is a two-scale integrator developed by the authors. While iterated applications of multiscale integrators for two different scales increase the computational complexity exponentially as the number of different scales increases, the proposed method, on the other hand, has computational complexity linearly proportional to the number of different scales. This efficiency is achieved by solving different scale components of the vector fields with variable time steps. It is shown that variable time stepping of different force components has an effect of fast integration for the effective force of the slow dynamics. The proposed fast integrator is numerically tested on problems with several different scales which are dissipative and highly oscillatory including multiscale partial differential equations with sparsity in the solution space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale numerical methods for passive advection-diffusion in incompressible turbulent flow fields

We propose a seamless multiscale method which approximates the macroscopic behavior of the passive advection-diffusion equations with steady incompressible velocity fields with multi-spatial scales. The method uses decompositions of the velocity fields in the Fourier space, which are similar to the decomposition in large eddy simulations. It also uses a hierarchy of local domains with different resolutions as in multigrid methods. The effective diffusivity from finer scale is used for the next coarser level computation and this process is repeated up to the coarsest scale of interest. The grids are only in local domains whose sizes decrease depending on the resolution level so that the overall computational complexity increases linearly as the number of different resolution grids increases. The method captures interactions between finer and coarser scales but has to sacrifice some of the interaction between the fine scales. The proposed method is numerically tested with 2D examples including a successful approximation to a continuous spectrum flow.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fast Sweeping Methods for Hyperbolic Systems of Conservation Laws at Steady State II

The idea of using fast sweeping methods for solving stationary systems of conservation laws has previously been proposed for efficiently computing solutions with sharp shocks. We further develop these methods to allow for a more challenging class of problems including problems with sonic points, shocks originating in the interior of the domain, rarefaction waves, and two-dimensional systems. We show that fast sweeping methods can produce higher-order accuracy. Computational results validate the claims of accuracy, sharp shock curves, and optimal computational efficiency.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Application of the Wasserstein metric to seismic signals

Seismic signals are typically compared using travel time difference or $L_2$ difference. We propose the Wasserstein metric as an alternative measure of fidelity or misfit in seismology. It exhibits properties from both of the traditional measures mentioned above. The numerical computation is based on the recent development of fast numerical methods for the Monge-Ampere equation and optimal transport. Applications to waveform inversion and registration are discussed and simple numerical examples are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast Sweeping Methods for Hyperbolic Systems of Conservation Laws at Steady State

Fast sweeping methods have become a useful tool for computing the solutions of static Hamilton-Jacobi equations. By adapting the main idea behind these methods, we describe a new approach for computing steady state solutions to systems of conservation laws. By exploiting the flow of information along characteristics, these fast sweeping methods can compute solutions very efficiently. Furthermore, the methods capture shocks sharply by directly imposing the Rankine-Hugoniot shock conditions. We present convergence analysis and numerics for several one- and two-dimensional examples to illustrate the use and advantages of this approach.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variable Step Size Multiscale Methods for Stiff and Highly Oscillatory Dynamical Systems

We present a new numerical multiscale integrator for stiff and highly oscillatory dynamical systems. The new algorithm can be seen as an improved version of the seamless Heterogeneous Multiscale Method by E, Ren, and Vanden-Eijnden and the method FLAVORS by Tao, Owhadi, and Marsden. It approximates slowly changing quantities in the solution with higher accuracy than these other methods while maintaining the same computational complexity. To achieve higher accuracy, it uses variable mesoscopic time steps which are determined by a special function satisfying moment and regularity conditions. Detailed analytical and numerical comparison between the different methods are given.

0 citations0 reviewsNo code linkedOpen access