1,032 research outputs found

    Locally optimal unstructured finite element meshes in 3 dimensions

    Get PDF
    This paper investigates the adaptive finite element solution of a general class of variational problems in three dimensions using a combination of node movement, edge swapping, face swapping and node insertion. The adaptive strategy proposed is a generalization of previous work in two dimensions and is based upon the construction of a hierarchy of locally optimal meshes. Results presented, both for a single equation and a system of coupled equations, suggest that this approach is able to produce better meshes of tetrahedra than those obtained by more conventional adaptive strategies and in a relatively efficient manner

    An inexact Newton method for systems arising from the finite element method

    Get PDF
    In this paper, we introduce an efficient and robust technique for approximating the Jacobian matrix for a nonlinear system of algebraic equations which arises from the finite element discretization of a system of nonlinear partial differential equations. It is demonstrated that when an iterative solver, such as preconditioned GMRES, is used to solve the linear systems of equations that result from the application of Newton's method, this approach is generally more efficient than using matrix-free techniques: the price paid being the extra memory requirement for storing the sparse Jacobian. The advantages of this approach over attempting to calculate the Jacobian exactly or of using other approximations are also discussed. A numerical example is included which is based upon the solution of a 2-d compressible viscous flow problem

    Parallel performance prediction for multigrid codes on distributed memory architectures

    Get PDF
    We propose a model for describing the parallel performance of multigrid software on distributed memory architectures. The goal of the model is to allow reliable predictions to be made as to the execution time of a given code on a large number of processors, of a given parallel system, by only benchmarking the code on small numbers of processors. This has potential applications for the scheduling of jobs in a Grid computing environment where reliable predictions as to execution times on different systems will be valuable. The model is tested for two different multigrid codes running on two different parallel architectures and the results obtained are discussed

    A weakly overlapping parallel domain decomposition preconditioner for the finite element solution of convection-dominated problems in three dimensions

    Get PDF
    In this paper we describe the parallel application of a novel two level additive Schwarz preconditioner to the stable finite element solution of convection-dominated problems in three dimensions. This is a generalization of earlier work, [2,6], in 2-d and 3-d respectively. An algebraic formulation of the preconditioner is presented and the key issues associated with its parallel implementation are discussed. Some computational results are also included which demonstrate empirically the optimality of the preconditioner and its potential for parallel implementation

    A fully implicit, fully adaptive time and space discretisation method for phase-field simulation of binary alloy solidification

    Get PDF
    A fully-implicit numerical method based upon adaptively refined meshes for the simulation of binary alloy solidification in 2D is presented. In addition we combine a second-order fully-implicit time discretisation scheme with variable steps size control to obtain an adaptive time and space discretisation method. The superiority of this method, compared to widely used fully-explicit methods, with respect to CPU time and accuracy, is shown. Due to the high non-linearity of the governing equations a robust and fast solver for systems of nonlinear algebraic equations is needed to solve the intermediate approximations per time step. We use a nonlinear multigrid solver which shows almost h-independent convergence behaviour

    An adaptive, fully implicit multigrid phase-field model for the quantitative simulation of non-isothermal binary alloy solidification

    Get PDF
    Using state-of-the-art numerical techniques, such as mesh adaptivity, implicit time-stepping and a non-linear multi-grid solver, the phase-field equations for the non-isothermal solidification of a dilute binary alloy have been solved. Using the quantitative, thin-interface formulation of the problem we have found that at high Lewis number a minimum in the dendrite tip radius is predicted with increasing undercooling, as predicted by marginal stability theory. Over the dimensionless undercooling range 0.2–0.8 the radius selection parameter, σ*, was observed to vary by over a factor of 2 and in a non-monotonic fashion, despite the anisotropy strength being constant

    Anisotropic adaptivity for the finite element solutions of three-dimensional convection-dominated problems

    Get PDF
    Convection-dominated problems are typified by the presence of strongly directional features such as shock waves or boundary layers. Resolution of numerical solutions using an isotropic mesh can lead to unnecessary refinement in directions parallel to such features. This is particularly important in three dimensions where the grid size increases rapidly during conventional isotropic refinement procedures. In this work, we investigate the use of adaptive finite element methods using anisotropic mesh refinement strategies for convection-dominated problems. The strategies considered here aim to resolve directional features without excessive resolution in other directions, and hence achieve accurate solutions more efficiently. Two such strategies are described here: the first based on minimization of the least-squares residual; the second based on minimizing a finite element error estimate. These are incorporated into an hr-adaptive finite element method and tested on a simple model problem

    Advanced numerical methods for the simulation of alloy solidification with high Lewis number

    Get PDF
    A fully-implicit numerical method based upon adaptively refined meshes for the thermal-solutal simulation of alloy solidification in 2D is presented. In addition we combine an unconditional stable second-order fully-implicit time discretisation scheme with variable step size control to obtain an adaptive time and space discretisation method, where a robust and fast multigrid solver for systems of non-linear algebraic equations is used to solve the intermediate approximations per time step. For the isothermal case, the superiority of this method, compared to widely used fully-explicit methods, with respect to CPU time and accuracy, has been demonstrated and published previously. Here, the new proposed method has been applied to the thermalsolutal case with high Lewis number, where stability issues and time step restrictions have been major constraints in previous research
    • …
    corecore