We discuss the application of the finite-element (FE) method to ab initio solid-state: electronic-structure calculations. In this method, the basis functions are strictly local, piecewise polynomials. Because the basi...
详细信息
We discuss the application of the finite-element (FE) method to ab initio solid-state: electronic-structure calculations. In this method, the basis functions are strictly local, piecewise polynomials. Because the basis is composed of polynomials, the method is completely general and its convergence can be controlled systematically. Because the basis functions are strictly local in real space, the method allows for variable resolution in real space;produces sparse, structured matrices, enabling the effective use of iterative solution methods;and is well suited to parallel implementation. The method thus combines the significant advantages of both real-space-grid and basis-oriented approaches and so promises to be particularly well suited for large, accurate ab initio calculations. We discuss the construction and properties of the required FE bases and develop in detail their use in the solution of the Schrodinger and Poisson equations subject to boundary conditions appropriate for a periodic solid. We present results for the Schrodinger equation illustrating the rapid, variational convergence of the method in electronic band-structure calculations. We present results for the Poisson equation illustrating the rapid convergence of the method, both pointwise and in the L-2 norm, and its linear scaling with system size in the context of a model charge-density and Si pseudo-charge-density. Finally, we discuss the application of the method to large-scale ab initio positron distribution and lifetime calculations in solids and present results for a host of systems within the range of a conventional LMTO based approach for comparison, as well as results for systems well beyond the range of the conventional approach. The largest such calculation, involving a unit cell of 4092 atoms, was shown to be well within the range of the FE approach on existing computational platforms. (C) 2001 Elsevier Science B.V. All rights reserved.
finite-element methods for the approximation of the solution of streamfunction-vorticity equations are considered. Among the issues dealt with are multiply connected domains, the use of low-order elements, the incorpo...
详细信息
finite-element methods for the approximation of the solution of streamfunction-vorticity equations are considered. Among the issues dealt with are multiply connected domains, the use of low-order elements, the incorporation of a variety of boundary conditions into the methodology, error estimates, and the recovery of the primitive variables. Various numerical examples are also provided.
The main goal of this work is the generation of ground-truth data for the validation of atrophy measurement techniques, commonly used in the study of neurodegenerative diseases such as dementia. Several techniques hav...
详细信息
The main goal of this work is the generation of ground-truth data for the validation of atrophy measurement techniques, commonly used in the study of neurodegenerative diseases such as dementia. Several techniques have been used to measure atrophy in cross-sectional and longitudinal studies, but it is extremely difficult to compare their performance since they have been applied to different patient populations. Furthermore, assessment of performance based on phantom measurements or simple scaled images overestimates these techniques' ability to capture the complexity of neurodegeneration of the human brain. We propose a method for atrophy simulation in structural magnetic resonance (MR) images based on finite-element methods. The method produces cohorts of brain images with known change that is physically and clinically plausible, providing data for objective evaluation of atrophy measurement techniques. Atrophy is simulated in different tissue compartments or in different neuroanatomical structures with a phenomenological model. This model of diffuse global and regional atrophy is based on volumetric measurements such as the brain or the hippocampus, from patients with known disease and guided by clinical knowledge of the relative pathological involvement of regions and tissues. The consequent biomechanical readjustment of structures is modelled using conventional physics-based techniques based on biomechanical tissue properties and simulating plausible tissue deformations with finite-element methods. A thermoelastic model of tissue deformation is employed, controlling the rate of progression of atrophy by means of a set of thermal coefficients, each one corresponding to a different type of tissue. Tissue characterization is performed by means of the meshing of a labelled brain atlas, creating a reference volumetric mesh that will be introduced to a finite-element solver to create the simulated deformations. Preliminary work on the simulation of acquisition artefac
The investigation and reconstruction are fundamental to enable designers to improve the subway vehicle's crashworthiness. In this paper, a crashworthiness study of a subway vehicle collision accident in Beijing wa...
详细信息
The investigation and reconstruction are fundamental to enable designers to improve the subway vehicle's crashworthiness. In this paper, a crashworthiness study of a subway vehicle collision accident in Beijing was conducted. Based on field investigation, the full-scale detail finite-element model, has established to calculate the dynamic behaviour of the vehicle in the accident. Then, the numerical simulation assists with the field investigation to identifying the structural weakness, which leading the injury of the occupant. Finally, an enhancement design for the structural weak point was considered to improve the vehicle crashworthiness in the same accident scenario.
Automatically optimizing the design of a microwave device can be prohibitively time-consuming when a numerical electromagnetic-field analysis is necessary at each iteration. However, the time taken for the field analy...
详细信息
Automatically optimizing the design of a microwave device can be prohibitively time-consuming when a numerical electromagnetic-field analysis is necessary at each iteration. However, the time taken for the field analysis depends on the accuracy required, and in the early stage of the optimization relatively inaccurate solutions are adequate. This idea is exploited in a scheme that combines a quasi-Newton constrained optimizer with a two-dimensional p-adaptive finite-element method for finding scattering parameters. The scheme has been tested on three H-plane rectangular waveguide devices: a T-junction, a miter bend with a dielectric column, and a two-cavity iris-coupled filter. Time savings of more than an order of magnitude were obtained, compared to the standard approach of requiring equally high accuracy throughout the optimization.
We propose a novel finite-element method for polygonal meshes. The resulting scheme is hp-adaptive, where h and p are a measure of, respectively, the size and the number of degrees of freedom of each polygon. Moreover...
详细信息
We propose a novel finite-element method for polygonal meshes. The resulting scheme is hp-adaptive, where h and p are a measure of, respectively, the size and the number of degrees of freedom of each polygon. Moreover, it is locally meshfree, since it is possible to arbitrarily choose the locations of the degrees of freedom inside each polygon. Our construction is based on nodal kernel functions, whose support consists of all polygons that contain a given node. This ensures a significantly higher sparsity compared to standard meshfree approximations. In this work, we choose axis-aligned quadrilaterals as polygonal primitives and maximum entropy approximants as kernels. However, any other convex approximation scheme and convex polygons can be employed. We study the optimal placement of nodes for regular elements, ie, those that are not intersected by the boundary, and propose a method to generate a suitable mesh. Finally, we show via numerical experiments that the proposed approach provides good accuracy without undermining the sparsity of the resulting matrices.
In this paper, various time-domain finite-element methods for the simulation of transient electromagnetic wave phenomena are discussed, Detailed descriptions of test/trial spaces, explicit and implicit formulations, n...
详细信息
In this paper, various time-domain finite-element methods for the simulation of transient electromagnetic wave phenomena are discussed, Detailed descriptions of test/trial spaces, explicit and implicit formulations, nodal and edge/facet element basis functions are given, along with the numerical stability properties of the different methods, The advantages and disadvantages of mass lumping are examined, Finally, the various formulations are compared on the basis of their numerical dispersion performance.
The perfectly matched layer (PML) concept introduced by Berenger is implemented for nodal-based finite-element frequency-domain methods. Starting from a scalar/vector potential framework, anisotropic media-equivalent ...
详细信息
The perfectly matched layer (PML) concept introduced by Berenger is implemented for nodal-based finite-element frequency-domain methods. Starting from a scalar/vector potential framework, anisotropic media-equivalent gauge conditions are developed for both coupled and uncoupled (i.e., direct field) scalar/vector field formulations. The resulting discrete system of equations are shown to be identical for both the anisotropic and stretched coordinate viewpoints of PML mesh termination on node-based finiteelements, Reaching this equivalency requires that special attention be paid to the basis/weighting functions used within the PML region, specifically, a material dependency is found to be essential. The alternative but identical stretched coordinate approach provides the perspective needed to realize a scheme for generalizing the PML to non-Cartesian mesh terminations which are more natural in the finite-element context, Several benchmark problems and associated numerical results are presented to demonstrate the performance of the PML on node-based finiteelements.
For bounded smooth domains, we study how the solution of the Stokes problem is bounded in terms of the data when the domain changes. We show that standard bounds for a fixed domain hold with the same constant for wide...
详细信息
For bounded smooth domains, we study how the solution of the Stokes problem is bounded in terms of the data when the domain changes. We show that standard bounds for a fixed domain hold with the same constant for wide classes of domains. This is done by first reviewing the original results of Cattabriga and then, in terms of geometric properties of the domains, by specifying when to apply Cattabriga's intermediate results. We do the same with standard strong extension and trace theorems. We apply these results to an elegant technique of analysis due to Wahlbin to overcome the disparity between a curved domain and the domains where finite-element computations are carried out in practice without resorting to numerical quadrature.
暂无评论