• Nuclear Engineering and Technology
• /
• v.55 no.2
• /
• pp.769-779
• /
• 2023

This paper describes an hp-angular adaptivity algorithm in the discrete ordinates method for Boltzmann transport applications with strong angular effects. This adaptivity uses discontinuous finite element quadrature sets with different degrees, which updates both angular mesh and the degree of the underlying discontinuous finite element basis functions, allowing different angular local refinement to be applied in space. The regular and goal-based error metrics are considered in this algorithm to locate some regions to be refined. A mapping algorithm derived by moment conservation is developed to pass the angular solution between spatial regions with different quadrature sets. The proposed method is applied to some test problems that demonstrate the ability of this hp-angular adaptivity to resolve complex fluxes with relatively few angular unknowns. Results illustrate that a reduction to approximately 1/50 in quadrature ordinates for a given accuracy compared with uniform angular discretization. This method therefore offers a highly efficient angular adaptivity for investigating difficult particle transport problems.

• Geophysics and Geophysical Exploration
• /
• v.13 no.1
• /
• pp.118-127
• /
• 2010

Resistivity logging instruments are designed to measure the electrical resistivity of a formation, and this can be directly interpreted to provide a water-saturation profile. However, resistivity logs are sensitive to borehole and shoulder-bed effects, which often result in misinterpretation of the results. These effects are emphasised more in the presence of tool eccentricity. For precise interpretation of short- and long-normal logging measurements in the presence of tool eccentricity, we simulate and analyse eccentricity effects by combining the use of a Fourier series expansion in a new system of coordinates with a 2D goal-oriented high-order self-adaptive hp finite-element refinement strategy, where h denotes the element size and p the polynomial order of approximation within each element. The algorithm automatically performs local mesh refinement to construct an optimal grid for the problem under consideration. In addition, the proper combination of h and p refinements produces highly accurate simulations even in the presence of high electrical resistivity contrasts. Numerical results demonstrate that our algorithm provides highly accurate and reliable simulation results. Eccentricity effects are more noticeable when the borehole is large or resistive, or when the formation is highly conductive.