skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Meshless X-ray computed tomography with Cauchy-type integration
This paper presents a new algorithm for X-ray Computerized Tomography (CT) based on Bukhgeim’s theory of analytic maps. The reconstruction relies on a Cauchy-type integral formula, where the integration over the boundary replaces the integration in the back- projection operator used in existing algorithms. From the numerical computation stand point, the proposed method recovers the attenuation coefficient at arbitrarily points by utilizing the boundary integration without internal global meshes. This means that it achieves high-parallel efficiency, and it reduces computational resources. Some numerical examples are presented to show feasibility of the proposed algorithm.  more » « less
Award ID(s):
1907097
PAR ID:
10181269
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Keisan suuri kougaku rombunshuu
Volume:
19
Issue:
12
ISSN:
1348-5245
Page Range / eLocation ID:
01-191201
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The present paper proposes a novel numerical scheme to X-ray Computerized Tomography (CT) from partial measurement data. In order to reduce radiation exposure, it is desirable to irradiate X-ray only around region of interest (ROI), while the conventional reconstruction methods such as filtered back projection (FBP) could not work due to its intrinsic limitation of dependency on whole measurement data. The proposed method gives a direct numerical reconstruction employing a Cauchy type boundary integration in $$A$$-analytic theory and a singular integral equation which maps boundary measurement to interior data. Numerical examples using experimental data are also exhibited to show validity of the proposed numerical procedure. 
    more » « less
  2. This paper is concerned with developing an efficient numerical algorithm for the fast implementation of the sparse grid method for computing the d-dimensional integral of a given function. The new algorithm, called the MDI-SG (multilevel dimension iteration sparse grid) method, implements the sparse grid method based on a dimension iteration/reduction procedure. It does not need to store the integration points, nor does it compute the function values independently at each integration point; instead, it reuses the computation for function evaluations as much as possible by performing the function evaluations at all integration points in a cluster and iteratively along coordinate directions. It is shown numerically that the computational complexity (in terms of CPU time) of the proposed MDI-SG method is of polynomial order O(d3Nb)(b≤2) or better, compared to the exponential order O(N(logN)d−1) for the standard sparse grid method, where N denotes the maximum number of integration points in each coordinate direction. As a result, the proposed MDI-SG method effectively circumvents the curse of dimensionality suffered by the standard sparse grid method for high-dimensional numerical integration. 
    more » « less
  3. Abstract We consider the finite element approximation of a coupled fluid‐structure interaction (FSI) system, which comprises a three‐dimensional (3D) Stokes flow and a two‐dimensional (2D) fourth‐order Euler–Bernoulli or Kirchhoff plate. The interaction of these parabolic and hyperbolic partial differential equations (PDE) occurs at the boundary interface which is assumed to be fixed. The vertical displacement of the plate dynamics evolves on the flat portion of the boundary where the coupling conditions are implemented via the matching velocities of the plate and fluid flow, as well as the Dirichlet boundary trace of the pressure. This pressure term also acts as a coupling agent, since it appears as a forcing term on the flat, elastic plate domain. Our main focus in this work is to generate some numerical results concerning the approximate solutions to the FSI model. For this, we propose a numerical algorithm that sequentially solves the fluid and plate subsystems through an effective decoupling approach. Numerical results of test problems are presented to illustrate the performance of the proposed method. 
    more » « less
  4. Abstract In their article “Coupling at a distance HDG and BEM” , Cockburn, Sayas and Solano proposed an iterative coupling of the hybridizable discontinuous Galerkin method (HDG) and the boundary element method (BEM) to solve an exterior Dirichlet problem. The novelty of the numerical scheme consisted of using a computational domain for the HDG discretization whose boundary did not coincide with the coupling interface. In their article, the authors provided extensive numerical evidence for convergence, but the proof of convergence and the error analysis remained elusive at that time. In this article we fill the gap by proving the convergence of a relaxation of the algorithm and providing a priori error estimates for the numerical solution. 
    more » « less
  5. Integrated sensing and communication (ISAC) is emerging as a key technique for next-generation wireless systems. In order to expedite the practical implementation of ISAC in pervasive mobile networks, it is crucial to have widely deployed base stations with radar sensing capabilities. Thus, the utilization of standardized multiple-input multiple-output (MIMO) orthogonal frequency division multiplexing (OFDM) hardware architectures and waveforms is pivotal for realizing seamless integration of effective communication and sensing functionalities. In this paper, we introduce a novel joint anglerange- velocity estimation algorithm for MIMO-OFDM ISAC systems. This approach exclusively depends on the format of conventional MIMO-OFDM waveforms that are widely adopted in wireless communications. Specifically, the angle-range-velocity information of potential targets is jointly extracted by utilizing all the received echo signals within a coherent processing interval (CPI). The proposed joint estimation algorithm can achieve larger signal-to-noise-ratio (SNR) processing gains and higher resolution by fully exploiting the echo signals and jointly estimating the angle-range-velocity information. A theoretical analysis for maximum unambiguous range, resolution, and SNR processing gains is provided to verify the advantages of the proposed joint estimation algorithm. Finally, the results of extensive numerical experiments are presented to demonstrate that the proposed joint estimation approach can achieve significantly lower rootmean- square-error (RMSE) performance for angle/range/velocity estimation for both single- and multi-target scenarios. 
    more » « less