4,592 research outputs found
An adaptive, hanging-node, discontinuous isogeometric analysis method for the first-order form of the neutron transport equation with discrete ordinate (SN) angular discretisation
In this paper a discontinuous, hanging-node, isogeometric analysis (IGA) method is developed and applied to the first-order form of the neutron transport equation with a discrete ordinate (SN) angular discretisation in two-dimensional space. The complexities involved in upwinding across curved element boundaries that contain hanging-nodes have been addressed to ensure that the scheme remains conservative. A robust algorithm for cycle-breaking has also been introduced in order to develop a unique sweep ordering of the elements for each discrete ordinates direction. The convergence rate of the scheme has been verified using the method of manufactured solutions (MMS) with a smooth solution. Heuristic error indicators have been used to drive an adaptive mesh refinement (AMR) algorithm to take advantage of the hanging-node discretisation. The effectiveness of this method is demonstrated for three test cases. The first is a homogeneous square in a vacuum with varying mean free path and a prescribed extraneous unit source. The second test case is a radiation shielding problem and the third is a 3×3 “supercell” featuring a burnable absorber. In the final test case, comparisons are made to the discontinuous Galerkin finite element method (DGFEM) using both straight-sided and curved quadratic finite elements
Goal-based h-adaptivity of the 1-D diamond difference discrete ordinate method.
The quantity of interest (QoI) associated with a solution of a partial differential equation (PDE) is not, in general, the solution itself, but a functional of the solution. Dual weighted residual (DWR) error estimators are one way of providing an estimate of the error in the QoI resulting from the discretisation of the PDE. This paper aims to provide an estimate of the error in the QoI due to the spatial discretisation, where the discretisation scheme being used is the diamond difference (DD) method in space and discrete ordinate (SNSN) method in angle. The QoI are reaction rates in detectors and the value of the eigenvalue (Keff)(Keff) for 1-D fixed source and eigenvalue (KeffKeff criticality) neutron transport problems respectively. Local values of the DWR over individual cells are used as error indicators for goal-based mesh refinement, which aims to give an optimal mesh for a given QoI
JIGSAW-GEO (1.0): locally orthogonal staggered unstructured grid generation for general circulation modelling on the sphere
An algorithm for the generation of non-uniform, locally-orthogonal staggered
unstructured spheroidal grids is described. This technique is designed to
generate very high-quality staggered Voronoi/Delaunay meshes appropriate for
general circulation modelling on the sphere, including applications to
atmospheric simulation, ocean-modelling and numerical weather prediction. Using
a recently developed Frontal-Delaunay refinement technique, a method for the
construction of high-quality unstructured spheroidal Delaunay triangulations is
introduced. A locally-orthogonal polygonal grid, derived from the associated
Voronoi diagram, is computed as the staggered dual. It is shown that use of the
Frontal-Delaunay refinement technique allows for the generation of very
high-quality unstructured triangulations, satisfying a-priori bounds on element
size and shape. Grid-quality is further improved through the application of
hill-climbing type optimisation techniques. Overall, the algorithm is shown to
produce grids with very high element quality and smooth grading
characteristics, while imposing relatively low computational expense. A
selection of uniform and non-uniform spheroidal grids appropriate for
high-resolution, multi-scale general circulation modelling are presented. These
grids are shown to satisfy the geometric constraints associated with
contemporary unstructured C-grid type finite-volume models, including the Model
for Prediction Across Scales (MPAS-O). The use of user-defined mesh-spacing
functions to generate smoothly graded, non-uniform grids for multi-resolution
type studies is discussed in detail.Comment: Final revisions, as per: Engwirda, D.: JIGSAW-GEO (1.0): locally
orthogonal staggered unstructured grid generation for general circulation
modelling on the sphere, Geosci. Model Dev., 10, 2117-2140,
https://doi.org/10.5194/gmd-10-2117-2017, 201
Stability of explicit one-step methods for P1-finite element approximation of linear diffusion equations on anisotropic meshes
We study the stability of explicit one-step integration schemes for the
linear finite element approximation of linear parabolic equations. The derived
bound on the largest permissible time step is tight for any mesh and any
diffusion matrix within a factor of , where is the spatial
dimension. Both full mass matrix and mass lumping are considered. The bound
reveals that the stability condition is affected by two factors. The first one
depends on the number of mesh elements and corresponds to the classic bound for
the Laplace operator on a uniform mesh. The other factor reflects the effects
of the interplay of the mesh geometry and the diffusion matrix. It is shown
that it is not the mesh geometry itself but the mesh geometry in relation to
the diffusion matrix that is crucial to the stability of explicit methods. When
the mesh is uniform in the metric specified by the inverse of the diffusion
matrix, the stability condition is comparable to the situation with the Laplace
operator on a uniform mesh. Numerical results are presented to verify the
theoretical findings.Comment: Revised WIAS Preprin
An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems
Heterogeneous anisotropic diffusion problems arise in the various areas of
science and engineering including plasma physics, petroleum engineering, and
image processing. Standard numerical methods can produce spurious oscillations
when they are used to solve those problems. A common approach to avoid this
difficulty is to design a proper numerical scheme and/or a proper mesh so that
the numerical solution validates the discrete counterpart (DMP) of the maximum
principle satisfied by the continuous solution. A well known mesh condition for
the DMP satisfaction by the linear finite element solution of isotropic
diffusion problems is the non-obtuse angle condition that requires the dihedral
angles of mesh elements to be non-obtuse. In this paper, a generalization of
the condition, the so-called anisotropic non-obtuse angle condition, is
developed for the finite element solution of heterogeneous anisotropic
diffusion problems. The new condition is essentially the same as the existing
one except that the dihedral angles are now measured in a metric depending on
the diffusion matrix of the underlying problem. Several variants of the new
condition are obtained. Based on one of them, two metric tensors for use in
anisotropic mesh generation are developed to account for DMP satisfaction and
the combination of DMP satisfaction and mesh adaptivity. Numerical examples are
given to demonstrate the features of the linear finite element method for
anisotropic meshes generated with the metric tensors.Comment: 34 page
Phase-field boundary conditions for the voxel finite cell method: surface-free stress analysis of CT-based bone structures
The voxel finite cell method employs unfitted finite element meshes and voxel quadrature rules to seamlessly
transfer CT data into patient-specific bone discretizations. The method, however, still requires the explicit
parametrization of boundary surfaces to impose traction and displacement boundary conditions, which
constitutes a potential roadblock to automation. We explore a phase-field based formulation for imposing
traction and displacement constraints in a diffuse sense. Its essential component is a diffuse geometry model
generated from metastable phase-field solutions of the Allen-Cahn problem that assumes the imaging data as
initial condition. Phase-field approximations of the boundary and its gradient are then employed to transfer
all boundary terms in the variational formulation into volumetric terms. We show that in the context of the
voxel finite cell method, diffuse boundary conditions achieve the same accuracy as boundary conditions
defined over explicit sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field,
the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method
by analyzing stresses in a human femur and a vertebral body
The diffuse Nitsche method: Dirichlet constraints on phase-field boundaries
We explore diffuse formulations of Nitsche's method for consistently imposing Dirichlet boundary conditions on phase-field approximations of sharp domains. Leveraging the properties of the phase-field gradient, we derive the variational formulation of the diffuse Nitsche method by transferring all integrals associated with the Dirichlet boundary from a geometrically sharp surface format in the standard Nitsche method to a geometrically diffuse volumetric format. We also derive conditions for the stability of the discrete system and formulate a diffuse local eigenvalue problem, from which the stabilization parameter can be estimated automatically in each element. We advertise metastable phase-field solutions of the Allen-Cahn problem for transferring complex imaging data into diffuse geometric models. In particular, we discuss the use of mixed meshes, that is, an adaptively refined mesh for the phase-field in the diffuse boundary region and a uniform mesh for the representation of the physics-based solution fields. We illustrate accuracy and convergence properties of the diffuse Nitsche method and demonstrate its advantages over diffuse penalty-type methods. In the context of imaging based analysis, we show that the diffuse Nitsche method achieves the same accuracy as the standard Nitsche method with sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field, the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method by analyzing stresses in a human vertebral body
- …