Algebraic multigrid (AMG) methods for directly solving coupled systems of partial differential equations (PDEs) have been extensively used in various types of numerical simulations in engineering. A necessary conditio...
详细信息
In this article, we present a cost-benefit analysis of the approximation in tensor products of Hilbert spaces of Sobolev-analytic type. The Sobolev part is defined on a finite dimensional domain, whereas the analytica...
详细信息
Workshops are used to explore a specific topic, transfer knowledge, solve identified problems or create something new. In funded research projects and other research endeavours, workshops are the mechanism to gather t...
We present a high-order spatial discretization of a continuum gyrokinetic Vlasov model in axisymmetric tokamak edge plasma geometries. Such models describe the phase space advection of plasma species distribution func...
详细信息
We present a high-order spatial discretization of a continuum gyrokinetic Vlasov model in axisymmetric tokamak edge plasma geometries. Such models describe the phase space advection of plasma species distribution functions in the absence of collisions. The gyrokinetic model is posed in a four-dimensional phase space, upon which a grid is imposed when discretized. To mitigate the computational cost associated with high-dimensional grids, we employ a high-order discretization to reduce the grid size needed to achieve a given level of accuracy relative to lower-order methods. Strong anisotropy induced by the magnetic field motivates the use of mapped coordinate grids aligned with magnetic flux surfaces. The natural partitioning of the edge geometry by the separatrix between the closed and open field line regions leads to the consideration of multiple mapped blocks, in what is known as a mapped multiblock (MMB) approach. We describe the specialization of a more general formalism that we have developed for the construction of high-order, finite-volume discretizations on MMB grids, yielding the accurate evaluation of the gyrokinetic Vlasov operator, the metric factors resulting from the MMB coordinate mappings, and the interaction of blocks at adjacent boundaries. Our conservative formulation of the gyrokinetic Vlasov model incorporates the fact that the phase space velocity has zero divergence, which must be preserved discretely to avoid truncation error accumulation. We describe an approach for the discrete evaluation of the gyrokinetic phase space velocity that preserves the divergence-free property to machine precision. A distinguishing feature of an edge geometry is the X point, where the poloidal field component vanishes. The inability to construct fully flux-surface aligned MMB coordinate systems that are smooth up to and through the X point requires the relaxation of alignment in the vicinity of this point. We therefore describe an approach for the generation of s
We extend our system for radio astronomical monitoring by a cross-validation filter, separating near Earth radio frequency interference (RFI) from deep space signals. The filter searches for similar signals in a nearb...
详细信息
We extend our system for radio astronomical monitoring by a cross-validation filter, separating near Earth radio frequency interference (RFI) from deep space signals. The filter searches for similar signals in a nearby frequency band, coming from a different spatial direction than the tested signal. The filter passes the signals which do not have such duplicates. We apply this technique to a database of SETI Institute (***), containing 1.5 millions of sky observations in a frequency range 0.5-11.2 GHz, where our primary selection identified 28 strong signals possessing extraterrestrial (ET) signature. Crossvalidation allows to filter out 24 of those signals as satellite RFI. We present parameters for the remaining 4 signals and discuss statistical significance of these findings.
Fully implicit black-oil simulations result in huge, often very-ill-conditioned, linear systems of equations for different unknowns (e.g., pressure and saturations). It is well-known that the underlying Jacobian matri...
详细信息
Fully implicit black-oil simulations result in huge, often very-ill-conditioned, linear systems of equations for different unknowns (e.g., pressure and saturations). It is well-known that the underlying Jacobian matrices contain both hyperbolic and nearly elliptic subsystems (corresponding to saturations and pressure, respectively). Because a reservoir simulation is typically driven by the behavior of the pressure, constrained-pressure-residual (CPR)-type two-stage preconditioning methods to solve the coupled linear systems are a natural choice and still belong to the most popular approaches. After a suitable extraction and decoupling, the computationally most costly step in such two-stage methods consists in solving the elliptic subsystems accurately enough. Algebraic multigrid (AMG) provides a technique to solve elliptic linear equations very efficiently. Hence, in recent years, corresponding CPR-AMG approaches have been extensively used in practice. Unfortunately, if applied in a straightforward manner, CPR-AMG does not always work as expected. In this paper, we discuss the reasons for the lack of robustness observed in practice, and present remedies. More precisely, we will propose a preconditioning strategy (based on a suitable combination of left and right preconditioning of the Jacobian matrix) that aims at a compromise between the solvability of the pressure subproblem by AMG and the needs of the outer CPR process. The robustness of this new preconditioning strategy will be demonstrated for several industrial test cases, some of which are very ill-conditioned. Furthermore, we will demonstrate that CPR-AMG can be interpreted in a natural way as a special AMG process applied directly to the coupled Jacobian systems.
An application of Virtual Environment for real-time analysis of radio astronomical data is presented. Serving as an interface to radio astronomical databases or being applied to live data from the radio telescope, the...
详细信息
An application of Virtual Environment for real-time analysis of radio astronomical data is presented. Serving as an interface to radio astronomical databases or being applied to live data from the radio telescope, the application supports various data filters measuring signal-to-noise ratio (SNR), Doppler's drift, degree of signal localization on celestial sphere and other useful tools for signal extraction and classification. The application is created using Avango VE framework (***). The application has been used for analysis of radio astronomical data from SETI Institute (***).
We present StarWatch , our application for real-time analysis of radio astronomical data in Virtual Environment. Serving as an interface to radio astronomical databases or being applied to live data from the radio tel...
详细信息
We present StarWatch , our application for real-time analysis of radio astronomical data in Virtual Environment. Serving as an interface to radio astronomical databases or being applied to live data from the radio telescopes, the application supports various data filters measuring signal-to-noise ratio (SNR), Doppler's drift, degree of signal localization on celestial sphere and other useful tools for signal extraction and classification. Originally designed for the database of narrow band signals from SETI Institute (***), the application has been recently extended for the detection of wide band periodic signals, necessary for the search of pulsars. We will also address the detection of week signals possessing arbitrary waveforms and present several data filters suitable for this purpose.
暂无评论