Search

SUMMARY In this paper, we develop a thermodynamically consistent four-species model of tumor growth on the basis of the continuum theory of mixtures. Unique to this model is the incorporation of nutrient within the mixture as opposed to being modeled with an auxiliary reaction-diffusion equation. The formulation involves systems of highly nonlinear partial differential equations of surface effects through diffuse-interface models. A mixed finite element spatial discretization is developed and implemented to provide numerical results demonstrating the range of solutions this model can produce. A time-stepping algorithm is then presented for this system, which is shown to be first order accurate and energy gradient stable. The results of an array of numerical experiments are presented, which demonstrate a wide range of solutions produced by various choices of model parameters.Copyright © 2011 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

The travel time functional measures the time taken for a particle trajectory to travel from a given initial position to the boundary of the domain. Such evaluation is paramount in the postclosure safety assessment of deep geological storage facilities for radioactive waste where leaked, nonsorbing solutes can be transported to the surface of the site by the surrounding groundwater. The accurate simulation of this transport can be attained using standard dual-weighted-residual techniques to derive goal-oriented a posteriori error bounds. This work provides a key aspect in obtaining a suitable error estimate for the travel time functional: the evaluation of its Gateaux derivative. A mixed finite element method is implemented to approximate Darcy's equations, and numerical experiments are presented to test the performance of the proposed error estimator. In particular, we consider a test case inspired by the Sellafield site located in Cumbria, UK. [ABSTRACT FROM AUTHOR]

Recent developments in the context of minimum residual finite element methods are paving the way for designing quasi-optimal discretization methods in non-standard function spaces, such as q-type Sobolev spaces. For q → 1, these methods have demonstrated huge potential in avoiding the notorious Gibbs phenomena, i.e., the occurrence of spurious non-physical oscillations near thin layers and jump discontinuities. In this work we provide theoretical results that explain some of these numerical observations. In particular, we investigate the Gibbs phenomena for q-best approximations of discontinuities in finite element spaces with 1 ≤ q < ∞. We prove sufficient conditions on meshes in one and two dimensions such that over- and undershoots vanish in the limit q → 1. Moreover, we include examples of meshes such that Gibbs phenomena remain present even for q = 1 and demonstrate that our results can be used to design meshes so as to eliminate the Gibbs phenomenon. [ABSTRACT FROM AUTHOR]

We introduce the concept of machine-learning minimal-residual (ML-MRes) finite element discretizations of partial differential equations (PDEs), which resolve quantities of interest with striking accuracy, regardless of the underlying mesh size. The methods are obtained within a machine-learning framework during which the parameters defining the method are tuned against available training data. In particular, we use a provably stable parametric Petrov–Galerkin method that is equivalent to a minimal-residual formulation using a weighted norm. While the trial space is a standard finite element space, the test space has parameters that are tuned in an off-line stage. Finding the optimal test space therefore amounts to obtaining a goal-oriented discretization that is completely tailored towards the quantity of interest. We use an artificial neural network to define the parametric family of test spaces. Using numerical examples for the Laplacian and advection equation in one and two dimensions, we demonstrate that the ML-MRes finite element method has superior approximation of quantities of interest even on very coarse meshes. [ABSTRACT FROM AUTHOR]

We present an adaptive methodology for the solution of (linear and) nonlinear time dependent problems that is especially tailored for massively parallel computations. The basic concept is to solve for large blocks of space-time unknowns instead of marching sequentially in time. The methodology is a combination of a computationally efficient implementation of a parallel-in-space-time finite element solver coupled with a posteriori space-time error estimates and a parallel mesh generator. While we focus on spatial adaptivity in this work, the methodology enables simultaneous adaptivity in both space and time domains. We explore this basic concept in the context of a variety of time steppers including Θ-schemes and backward difference formulas. We specifically illustrate this framework with applications involving time dependent linear, quasi-linear, and semilinear diffusion equations. We focus on investigating how the coupled space-time refinement indicators for this class of problems affect spatial adaptivity. Finally, we show good scaling behavior up to 150,000 processors on the NCSA Blue Waters machine. This conceptually simple methodology enables scaling on next generation multicore machines by simultaneously solving for a large number of timesteps, and reducing computational overhead by locally refining spatial blocks that can track localized features. This methodology also opens up the possibility of efficiently incorporating adjoint equations for error estimators and inverse design problems, since blocks of space-time are simultaneously solved and stored in memory. [ABSTRACT FROM AUTHOR]

There is tremendous potential in using neural networks to optimize numerical methods. In this paper, we introduce and analyze a framework for the neural optimization of discrete weak formulations , suitable for finite element methods. The main idea of the framework is to include a neural-network function acting as a control variable in the weak form. Finding the neural control that (quasi-) minimizes a suitable cost (or loss) functional, then yields a numerical approximation with desirable attributes. In particular, the framework allows in a natural way the incorporation of known data of the exact solution, or the incorporation of stabilization mechanisms (e.g., to remove spurious oscillations). The main result of our analysis pertains to the well-posedness and convergence of the associated constrained-optimization problem. In particular, we prove under certain conditions, that the discrete weak forms are stable, and that quasi-minimizing neural controls exist, which converge quasi-optimally. We specialize the analysis results to Galerkin, least squares and minimal-residual formulations, where the neural-network dependence appears in the form of suitable weights. Elementary numerical experiments support our findings and demonstrate the potential of the framework. [ABSTRACT FROM AUTHOR]

We provide an abstract framework for optimal goal-oriented adaptivity for finite element methods and boundary element methods in the spirit of [C. Carstensen et al., Comput. Math. Appl., 67 (2014), pp. 1195-1253]. We prove that this framework covers standard discretizations of general second-order linear elliptic PDEs and hence generalizes available results [R. Becker, E. Estecahandy, and D. Trujillo, SIAM J. Numer. Anal., 49 (2011), pp. 2451-2469, M. S. Mommer and R. Stevenson, SIAM J. Numer. Anal., 47 (2009), pp. 861-886] beyond the Poisson equation. [ABSTRACT FROM AUTHOR]

Abstract: The numerical solution of fluid–structure-interaction problems poses a paradox in that most of the computational resources are consumed by the subsystem of least practical interest, viz., the fluid. Goal-oriented adaptive discretization methods provide a paradigm to bypass this paradox. Based on the solution of a dual problem, the contribution of local residuals to the error in a specific goal functional is estimated, and only the regions that yield a dominant contribution are refined. In the present work, we address a fundamental complication in the application of goal-oriented adaptivity to fluid–structure-interaction problems, namely, that the treatment of the interface conditions has nontrivial consequences for the properties of the dual problem. In the context of a linearized model problem, we consider two equivalent discretizations differing only on the formulation of the interface coupling terms. By means of an adjoint consistency analysis, we show that only one of these discretizations is adjoint consistent. Numerical experiments convey that the two discretizations behave very differently for the dual problem, and that the adjoint-consistent discretization yields more reliable error estimates. Based on the adjoint-consistent discretization, we finally present some h- and hp-adaptive results, confirming that tremendous savings in computational cost can be realized through the use of goal-oriented refinement strategies. The numerical experiments illustrate that the goal-oriented approach effectively equilibrates the error contributions of the fluid and structure subsystems, which is imperative for efficiently resolving the coupled fluid–structure-interaction problem, and which cannot be accomplished by uniform or residual-based refinement strategies. [ABSTRACT FROM AUTHOR]

Abstract Goal-oriented adaptivity is a powerful tool to accurately approximate physically relevant solution features for partial differential equations. In time dependent problems, we seek to represent the error in the quantity of interest as an integral over the whole space–time domain. A full space–time variational formulation allows such representation. Most authors employ implicit time marching schemes to perform goal-oriented adaptivity as it is known that they can be reinterpreted as Galerkin methods. In this work, we consider variational forms for explicit methods in time. We derive an appropriate error representation and propose a goal-oriented adaptive algorithm in space. For that, we derive the forward Euler method in time employing a discontinuous-in-time Petrov–Galerkin formulation. In terms of time domain adaptivity, we impose the Courant–Friedrichs–Lewy condition to ensure the stability of the method. We provide some numerical results in 1D space + time for the diffusion and advection–diffusion equations to show the performance of the proposed explicit-in-time goal-oriented adaptive algorithm. [ABSTRACT FROM AUTHOR]