Fast algorithm for the three-dimensional Poisson equation in infinite domains

2020 ◽  
Author(s):  
Chunxiong Zheng ◽  
Xiang Ma
Keyword(s):  
Poisson Equation ◽  
Fast Algorithm ◽  
Root Function ◽  
Computational Cost ◽  
Three Dimensional ◽  
Chebyshev Approximation ◽  
Laplacian Operator ◽  
Square Root ◽  
Infinite Domain ◽  
Infinite Domains

Abstract This paper is concerned with a fast finite element method for the three-dimensional Poisson equation in infinite domains. Both the exterior problem and the strip-tail problem are considered. Exact Dirichlet-to-Neumann (DtN)-type artificial boundary conditions (ABCs) are derived to reduce the original infinite-domain problems to suitable truncated-domain problems. Based on the best relative Chebyshev approximation for the square-root function, a fast algorithm is developed to approximate exact ABCs. One remarkable advantage is that one need not compute the full eigensystem associated with the surface Laplacian operator on artificial boundaries. In addition, compared with the modal expansion method and the method based on Pad$\acute{\textrm{e}}$ approximation for the square-root function, the computational cost of the DtN mapping is further reduced. An error analysis is performed and numerical examples are presented to demonstrate the efficiency of the proposed method.

3D simulation of emulsion flow using the boundary element method on the heterogeneous computational systems

2012 ◽  
Vol 9 (1) ◽  
pp. 142-146
Author(s):  
O.A. Solnyshkina
Keyword(s):  
Boundary Element Method ◽  
Boundary Element ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Problem Size ◽  
Low Reynolds Numbers ◽  
Infinite Domains ◽  
The Matrix ◽  
Computational Systems ◽  
Element Method

In this work the 3D dynamics of two immiscible liquids in unbounded domain at low Reynolds numbers is considered. The numerical method is based on the boundary element method, which is very efficient for simulation of the three-dimensional problems in infinite domains. To accelerate calculations and increase the problem size, a heterogeneous approach to parallelization of the computations on the central (CPU) and graphics (GPU) processors is applied. To accelerate the iterative solver (GMRES) and overcome the limitations associated with the size of the memory of the computation system, the software component of the matrix-vector product

Research on the application of 3D reconstruction technology in traditional handicraft design

2021 ◽  
pp. 002072092110052
Author(s):  
GuoLong Zhang
Keyword(s):  
Poisson Equation ◽  
Optimization Algorithm ◽  
Three Dimensional ◽  
Reconstruction Algorithm ◽  
Social Production ◽  
Design Constraint ◽  
Reconstruction Process ◽  
Key Steps ◽  
Shielding Design ◽  
Reconstruction Model

The use of computer technology for three-dimensional (3 D) reconstruction is one of the important development directions of social production. The purpose is to find a new method that can be used in traditional handicraft design, and to explore the application of 3 D reconstruction technology in it. Based on the description and analysis of 3 D reconstruction technology, the 3 D reconstruction algorithm based on Poisson equation is analyzed, and the key steps and problems of the method are clarified. Then, by introducing the shielding design constraint, a 3 D reconstruction algorithm based on shielded Poisson equation is proposed. Finally, the performance of two algorithms is compared by reconstructing the 3 D image of rabbit. The results show that: when the depth value of the algorithm is 11, the surface of the rabbit image obtained by the proposed optimization algorithm is smoother, and the details are more delicate and fluent; under different depth values, with the increase of the depth value, the number of vertices and faces of the two algorithms increase, and the optimal depth values of 3 D reconstruction are more than 8. However, the proposed optimization algorithm has more vertices, and performs better in the reconstruction process; the larger the depth value is, the more time and memory are consumed in 3 D reconstruction, so it is necessary to select the appropriate depth value; the shielding parameters of the algorithm have a great impact on the fineness of the reconstruction model. The larger the parameter is, the higher the fineness is. In a word, the proposed 3 D reconstruction algorithm based on shielded Poisson equation has better practicability and superiority.

scholarly journals The ICON Single-Column Mode

Atmosphere ◽  
2021 ◽  
Vol 12 (7) ◽  
pp. 906
Author(s):  
Ivan Bašták Ďurán ◽  
Martin Köhler ◽  
Astrid Eichhorn-Müller ◽  
Vera Maurer ◽  
Juerg Schmidli ◽  
...  
Keyword(s):  
Data Storage ◽  
Model Development ◽  
Computational Cost ◽  
Three Dimensional ◽  
Shallow Convection ◽  
Modeling Framework ◽  
Additional Tool ◽  
Eddy Simulation ◽  
Single Column ◽  
Stratocumulus Clouds

The single-column mode (SCM) of the ICON (ICOsahedral Nonhydrostatic) modeling framework is presented. The primary purpose of the ICON SCM is to use it as a tool for research, model evaluation and development. Thanks to the simplified geometry of the ICON SCM, various aspects of the ICON model, in particular the model physics, can be studied in a well-controlled environment. Additionally, the ICON SCM has a reduced computational cost and a low data storage demand. The ICON SCM can be utilized for idealized cases—several well-established cases are already included—or for semi-realistic cases based on analyses or model forecasts. As the case setup is defined by a single NetCDF file, new cases can be prepared easily by the modification of this file. We demonstrate the usage of the ICON SCM for different idealized cases such as shallow convection, stratocumulus clouds, and radiative transfer. Additionally, the ICON SCM is tested for a semi-realistic case together with an equivalent three-dimensional setup and the large eddy simulation mode of ICON. Such consistent comparisons across the hierarchy of ICON configurations are very helpful for model development. The ICON SCM will be implemented into the operational ICON model and will serve as an additional tool for advancing the development of the ICON model.

scholarly journals Retinal vessel segmentation with constrained-based nonnegative matrix factorization and 3D modified attention U-Net

2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Yang Yu ◽  
Hongqing Zhu
Keyword(s):  
Matrix Factorization ◽  
Nonnegative Matrix Factorization ◽  
Computational Cost ◽  
Three Dimensional ◽  
Nonnegative Matrix ◽  
Retinal Vessel ◽  
Classification Problem ◽  
Image Resolution ◽  
Retinal Vessels ◽  
Segmentation Task

AbstractDue to the complex morphology and characteristic of retinal vessels, it remains challenging for most of the existing algorithms to accurately detect them. This paper proposes a supervised retinal vessels extraction scheme using constrained-based nonnegative matrix factorization (NMF) and three dimensional (3D) modified attention U-Net architecture. The proposed method detects the retinal vessels by three major steps. First, we perform Gaussian filter and gamma correction on the green channel of retinal images to suppress background noise and adjust the contrast of images. Then, the study develops a new within-class and between-class constrained NMF algorithm to extract neighborhood feature information of every pixel and reduce feature data dimension. By using these constraints, the method can effectively gather similar features within-class and discriminate features between-class to improve feature description ability for each pixel. Next, this study formulates segmentation task as a classification problem and solves it with a more contributing 3D modified attention U-Net as a two-label classifier for reducing computational cost. This proposed network contains an upsampling to raise image resolution before encoding and revert image to its original size with a downsampling after three max-pooling layers. Besides, the attention gate (AG) set in these layers contributes to more accurate segmentation by maintaining details while suppressing noises. Finally, the experimental results on three publicly available datasets DRIVE, STARE, and HRF demonstrate better performance than most existing methods.

scholarly journals Data-Informed Decomposition for Localized Uncertainty Quantification of Dynamical Systems

Vibration ◽  
2020 ◽  
Vol 4 (1) ◽  
pp. 49-63
Author(s):  
Waad Subber ◽  
Sayan Ghosh ◽  
Piyush Pandita ◽  
Yiming Zhang ◽  
Liping Wang
Keyword(s):  
Dynamical Systems ◽  
Computational Cost ◽  
Three Dimensional ◽  
Region Of Interest ◽  
System Response ◽  
Computational Domain ◽  
User Interest ◽  
Numerical Resolution ◽  
Multi Scale ◽  
Operation Conditions

Industrial dynamical systems often exhibit multi-scale responses due to material heterogeneity and complex operation conditions. The smallest length-scale of the systems dynamics controls the numerical resolution required to resolve the embedded physics. In practice however, high numerical resolution is only required in a confined region of the domain where fast dynamics or localized material variability is exhibited, whereas a coarser discretization can be sufficient in the rest majority of the domain. Partitioning the complex dynamical system into smaller easier-to-solve problems based on the localized dynamics and material variability can reduce the overall computational cost. The region of interest can be specified based on the localized features of the solution, user interest, and correlation length of the material properties. For problems where a region of interest is not evident, Bayesian inference can provide a feasible solution. In this work, we employ a Bayesian framework to update the prior knowledge of the localized region of interest using measurements of the system response. Once, the region of interest is identified, the localized uncertainty is propagate forward through the computational domain. We demonstrate our framework using numerical experiments on a three-dimensional elastodynamic problem.

Correlation of Electrical, Structural, and Optical Properties of Erbium In Silicon

1993 ◽  
Vol 301 ◽  
Author(s):  
J. L. Benton ◽  
D. J. Eaglesham ◽  
M. Almonte ◽  
P. H. Citrin ◽  
M. A. Marcus ◽  
...  
Keyword(s):  
Root Function ◽  
Free Exciton ◽  
Excitation Power ◽  
Optically Active ◽  
Photoluminescence Intensity ◽  
Square Root ◽  
Structural And Optical Properties ◽  
X Ray ◽  
X Ray Absorption ◽  
Donor Activity

ABSTRACTAn understanding of the electrical, structural, and optical properites of Er in Si is necessary to evaluate this system as an opto-electronic material. Extended x-ray absorption fine structure, EXAFS, measurements of Er-implanted Si show that the optically active impurity complex is Er surrounded by an O cage of 6 atoms. The Er photoluminescence intensity is a square root function of excitation power, while the free exciton intensity increases linearly. The square root dependence of the 1.54μm-intensity is independent of measurement temperature and independent of co-implanted species. Ion-implantation of Er in Si introduces donor activity, but spreading resistance carrier concentration profiles indicate that these donors do not effect the optical activity of the Er.

Adaptivity and Three Dimensional Hierarchical Boundary Elements for Rigid Acoustic Scattering Problems

1995 ◽  
Author(s):  
Steven J. Newhouse ◽  
Ian C. Mathews
Keyword(s):  
Boundary Element ◽  
Acoustic Analysis ◽  
Three Dimensional ◽  
Acoustic Scattering ◽  
Error Estimator ◽  
Infinite Domain ◽  
Scattering Problems ◽  
Effective Form ◽  
Mesh Density ◽  
Acoustic Scattering Problems

Abstract The boundary element method is an established numerical tool for the analysis of acoustic pressure fields in an infinite domain. There is currently no well established method of estimating the surface pressure error distribution for an arbitrary three dimensional body. Hierarchical shape functions have been used as a highly effective form of p refinement in many finite and boundary element applications. Their ability to be used as an error estimator in acoustic analysis has never been fully exploited. This paper studies the influence of mesh density and interpolation order on several acoustic scattering problems. A hierarchical error estimator is implemented and its effectiveness verified against the spherical problem. A coarse cylindrical mesh is then refined using the new error estimator until the solution has converged. The effectiveness of this analysis is shown by comparing the error indicators derived during the analysis to the solution generated from a very fine cylindrical mesh.

scholarly journals Prediction of Residual Stresses in a Multipass Pipe Weld by a Novel 3D Finite Element Approach

2018 ◽  
Author(s):  
Hui Huang ◽  
Jian Chen ◽  
Blair Carlson ◽  
Hui-Ping Wang ◽  
Paul Crooker ◽  
...  
Keyword(s):  
Finite Element ◽  
Residual Stresses ◽  
High Performance ◽  
Large Scale ◽  
Graphics Processing Unit ◽  
Computational Cost ◽  
Three Dimensional ◽  
Processing Unit ◽  
Girth Welds ◽  
Welding Processes

Due to enormous computation cost, current residual stress simulation of multipass girth welds are mostly performed using two-dimensional (2D) axisymmetric models. The 2D model can only provide limited estimation on the residual stresses by assuming its axisymmetric distribution. In this study, a highly efficient thermal-mechanical finite element code for three dimensional (3D) model has been developed based on high performance Graphics Processing Unit (GPU) computers. Our code is further accelerated by considering the unique physics associated with welding processes that are characterized by steep temperature gradient and a moving arc heat source. It is capable of modeling large-scale welding problems that cannot be easily handled by the existing commercial simulation tools. To demonstrate the accuracy and efficiency, our code was compared with a commercial software by simulating a 3D multi-pass girth weld model with over 1 million elements. Our code achieved comparable solution accuracy with respect to the commercial one but with over 100 times saving on computational cost. Moreover, the three-dimensional analysis demonstrated more realistic stress distribution that is not axisymmetric in hoop direction.

Sign in / Sign up

Export Citation Format

Share Document