Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new user2D Fourier finite element formulation for magnetostatics in curvilinear coordinates with a symmetry direction
65
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a numerical method for the solution of linear magnetostatic problems in domains with a symmetry direction, including axial and translational symmetry. The approach uses a Fourier series decomposition of the vector potential formulation along the symmetry direction and covers both, zeroth (non-oscillatory) and non-zero (oscillatory) harmonics. For the latter it is possible to eliminate one component of the vector potential resulting in a fully transverse vector potential orthogonal to the transverse magnetic field. In addition to the Poisson-like equation for the longitudinal component of the non-oscillatory problem, a general curl-curl Helmholtz equation results for the transverse problem covering both, non-oscillatory and oscillatory case. The derivation is performed in the covariant formalism for curvilinear coordinates with a tensorial permeability and symmetry restrictions on metric and permeability tensor. The resulting variational forms are treated by the usual nodal finite element method for the longitudinal problem and by a two-dimensional edge element method for the transverse problem. The numerical solution can be computed independently for each harmonic which is favourable with regard to memory usage and parallelisation.
rate research
Read More
A high fidelity flow simulation for complex geometries for high Reynolds number ($Re$) flow is still very challenging, which requires more powerful computational capability of HPC system. However, the development of HPC with traditional CPU architecture suffers bottlenecks due to its high power consumption and technical difficulties. Heterogeneous architecture computation is raised to be a promising solution of difficulties of HPC development. GPU accelerating technology has been utilized in low order scheme CFD solvers on structured grid and high order scheme solvers on unstructured meshes. The high order finite difference methods on structured grid possess many advantages, e.g. high efficiency, robustness and low storage, however, the strong dependence among points for a high order finite difference scheme still limits its application on GPU platform. In present work, we propose a set of hardware-aware technology to optimize the efficiency of data transfer between CPU and GPU, and efficiency of communication between GPUs. An in-house multi-block structured CFD solver with high order finite difference methods on curvilinear coordinates is ported onto GPU platform, and obtain satisfying performance with speedup maximum around 2000x over a single CPU core. This work provides efficient solution to apply GPU computing in CFD simulation with certain high order finite difference methods on current GPU heterogeneous computers. The test shows that significant accelerating effects can been achieved for different GPUs.
We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $fin[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.
Planar supersymmetric quantum mechanical systems with separable spectral problem in curvilinear coordinates are analyzed in full generality. We explicitly construct the supersymmetric extension of the Euler/Pauli Hamiltonian describing the motion of a light particle in the field of two heavy fixed Coulombian centers. We shall also show how the SUSY Kepler/Coulomb problem arises in two different limits of this problem: either, the two centers collapse in one center - a problem separable in polar coordinates -, or, one of the two centers flies to infinity - to meet the Coulomb problem separable in parabolic coordinates.
We develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.
We present a practical approach for constructing meshes of general rough surfaces with given autocorrelation functions based on the unstructured meshes of nominally smooth surfaces. The approach builds on a well-known method to construct correlated random numbers from white noise using a decomposition of the autocorrelation matrix. We discuss important details arising in practical applications to the physicalmodeling of surface roughness and provide a software implementation to enable use of the approach with a broad range of numerical methods in various fields of science and engineering.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
