Search
Filters
total: 141
filtered: 131
Search results for: AUXILIARY SPACE PRECONDITIONING, FEM, GPU, MAXWELL’S EQUATIONS, MULTILEVEL PRECONDITIONING.
-
Block Conjugate Gradient Method with Multilevel Preconditioning and GPU Acceleration for FEM Problems in Electromagnetics
PublicationIn this paper a GPU-accelerated block conjugate gradient solver with multilevel preconditioning is presented for solving large system of sparse equations with multiple right hand-sides (RHSs) which arise in the finite-element analysis of electromagnetic problems. We demonstrate that blocking reduces the time to solution significantly and allows for better utilization of the computing power of GPUs, especially when the system matrix...
-
Tuning a Hybrid GPU-CPU V-Cycle Multilevel Preconditioner for Solving Large Real and Complex Systems of FEM Equations
PublicationThis letter presents techniques for tuning an accelerated preconditioned conjugate gradient solver with a multilevel preconditioner. The solver is optimized for a fast solution of sparse systems of equations arising in computational electromagnetics in a finite element method using higher-order elements. The goal of the tuning is to increase the throughput while at the same time reducing the memory requirements in order to allow...
-
Preconditioners with Low Memory Requirements for Higher-Order Finite-Element Method Applied to Solving Maxwell’s Equations on Multicore CPUs and GPUs
PublicationThis paper discusses two fast implementations of the conjugate gradient iterative method using a hierarchical multilevel preconditioner to solve the complex-valued, sparse systems obtained using the higher order finite-element method applied to the solution of the time-harmonic Maxwell equations. In the first implementation, denoted PCG-V, a classical V-cycle is applied and the system of equations on the lowest level is solved...
-
GENERAL DYNAMIC PROJECTING OF MAXWELL EQUATIONS
PublicationA complete – system of Maxwell equations is splitting into independent subsystems by means of a special dynamic projecting technique. The technique relies upon a direct link between field components that determine correspondent subspaces. The explicit form of links and corresponding subspace evolution equations are obtained in conditions of certain symmetry, it is illustrated by examples of spherical and quasi-one-dimensional waves.
-
A GPU Solver for Sparse Generalized Eigenvalue Problems with Symmetric Complex-Valued Matrices Obtained Using Higher-Order FEM
PublicationThe paper discusses a fast implementation of the stabilized locally optimal block preconditioned conjugate gradient (sLOBPCG) method, using a hierarchical multilevel preconditioner to solve nonHermitian sparse generalized eigenvalue problems with large symmetric complex-valued matrices obtained using the higher-order finite-element method (FEM), applied to the analysis of a microwave resonator. The resonant frequencies of the low-order...
-
GPU Acceleration of Multilevel Solvers for Analysis of Microwave Components With Finite Element Method
PublicationThe letter discusses a fast implementation of the conjugate gradient iterative method with ${rm E}$-field multilevel preconditioner applied to solving real symmetric and sparse systems obtained with vector finite element method. In order to accelerate computations, a graphics processing unit (GPU) was used and significant speed-up (2.61 fold) was achieved comparing to a central processing unit (CPU) based approach. These results...
-
GPU-Accelerated LOBPCG Method with Inexact Null-Space Filtering for Solving Generalized Eigenvalue Problems in Computational Electromagnetics Analysis with Higher-Order FEM
PublicationThis paper presents a GPU-accelerated implementation of the Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) method with an inexact nullspace filtering approach to find eigenvalues in electromagnetics analysis with higherorder FEM. The performance of the proposed approach is verified using the Kepler (Tesla K40c) graphics accelerator, and is compared to the performance of the implementation based on functions from...
-
A New Three-Dimensional Space Vector Modulation for Multilevel Four-Leg Converters Based on the Shape Functions of Tetrahedral Element
PublicationThe paper proposes a new three-dimensional space vector pulse-width modulation (3D-SVPWM) algorithm for multilevel four-leg converters. The proposed PWM duty cycle calculation is based on the shape functions of the threedimensional tetrahedral finite elements. The algorithm ensures synthesis of accurate and undistorted output voltages even under significant imbalance or ripple in the DC-link voltages. At the same time, the algorithm...
-
MULTILEVEL INTEGRATION OF DESIGN AND THE QUALITY OF LIFE IN URBAN SPACE ISSUE
Publication -
Space Vector Pulsewidth Modulation Strategy for Multilevel Cascaded H-Bridge Inverter With DC-Link Voltage Balancing Ability
PublicationSpace vector pulsewidth modulation (SVPWM) algorithms for cascaded H-bridge multilevel (CHB ML) inverter usually provide the possibility of using several combinations of active voltage vectors to generate the same output voltage vector. For preselected H-bridges, some of them may generate output voltages opposite to the assumed direction. This results in the change of the dc-link voltages of these H-bridges in the opposite direction...
-
Space vector modulation in multilevel inverters of the servo drives of the trajectory measurements telescopes
PublicationUsing the MatLab/Simulink mathematical model of a three-phase three-level voltage inverter, the influence of the space-vector modulation (SVM) algorithm on the pulsations of the current (torque) of an AC motor in the range of low rotation speeds is considered. It is shown that the SVM of the second kind does not provide a pulsations level comparable to the pulsations of a sinusoidal pulse-width modulation (SPWM), both in the static...
-
Multilevel model order reduction with generalized compression of boundaries for 3-d FEM electromagnetic analysis
PublicationThis paper presents a multilevel Model Order Reduction technique for a 3-D electromagnetic Finite Element Method analysis. The reduction process is carried out in a hierarchical way and involves several steps which are repeated at each level. This approach brings about versatility and allows one to efficiently analyze complex electromagnetic structures. In the proposed multilevel reduction the entire computational domain is covered...
-
THIRD-ORDER EXPONENTIAL INTEGRATOR FOR LINEAR KLEIN–GORDON EQUATIONS WITH TIME AND SPACE-DEPENDANT MASS
PublicationAllowing for space- and time-dependance of mass in Klein–Gordon equations re- solves the problem of negative probability density and of violation of Lorenz covariance of interaction in quantum mechanics. Moreover it extends their applicability to the domain of quantum cosmology, where the variation in mass may be accompanied by high oscillations....
-
THIRD-ORDER EXPONENTIAL INTEGRATOR FOR LINEAR KLEIN–GORDON EQUATIONS WITH TIME AND SPACE-DEPENDANT MASS
PublicationAllowing for space- and time-dependance of mass in Klein–Gordon equations re- solves the problem of negative probability density and of violation of Lorenz covariance of interaction in quantum mechanics. Moreover it extends their applicability to the domain of quantum cosmology, where the variation in mass may be accompanied by high oscillations....
-
Jacobi and gauss-seidel preconditioned complex conjugate gradient method with GPU acceleration for finite element method
PublicationIn this paper two implementations of iterative solvers for solving complex symmetric and sparse systems resulting from finite element method applied to wave equation are discussed. The problem under investigation is a dielectric resonator antenna (DRA) discretized by FEM with vector elements of the second order (LT/QN). The solvers use the preconditioned conjugate gradient (pcg) method implemented on Graphics Processing Unit (GPU)...
-
Space-Time Conservation Method applied to numerical solution of water hammer equations
PublicationArtykuł poświęcony jest metodzie czasoprzestrzennych objętości skończonych (STC) zastosowanej do przypadku uderzenia hydraulicznego w stalowym przewodzie pracującym pod ciśnieniem. Metoda STC ze względu na swoje własności numeryczne - m.in. wysoką dokładność - może być interesującą alternatywą dla tradycyjnych metod numerycznych, szczególnie w przypadku, gdy efekty numeryczne mają bardzo silny wpływ na rozwiązanie, tym samym utrudniając...
-
Three dimensional fluid equations from distribution function with discontinuity in velocity space
PublicationUkład równań typu hydrodynamicznego dla rozwarstwionego gazu w polu grawitacyjnym jest wyprowadzony z równania BKG przy użyciu metody kawałków ciągłej funkcji dystrybucji. Otrzymany układ równań uogólnia układ Naviera-Stokesa dla dowolnych liczb Knudsena.
-
Mountain pass type periodic solutions for Euler–Lagrange equations in anisotropic Orlicz–Sobolev space
PublicationUsing the Mountain Pass Theorem, we establish the existence of periodic solution for Euler–Lagrange equation. Lagrangian consists of kinetic part (an anisotropic G-function), potential part and a forcing term. We consider two situations: G satisfying at infinity and globally. We give conditions on the growth of the potential near zero for both situations.
-
GPU-accelerated finite element method
PublicationIn this paper the results of the acceleration of computations involved in analysing electromagnetic problems by means of the finite element method (FEM), obtained with graphics processors (GPU), are presented. A 4.7-fold acceleration was achieved thanks to the massive parallelization of the most time-consuming steps of FEM, namely finite-element matrix-generation and the solution of a sparse system of linear equations with the...
-
GPU-Accelerated 3D Mesh Deformation for Optimization Based on the Finite Element Method
PublicationThis paper discusses a strategy for speeding up the mesh deformation process in the design-byoptimization of high-frequency components involving electromagnetic field simulations using the 3D finite element method (FEM). The mesh deformation is assumed to be described by a linear elasticity model of a rigid body; therefore, each time the shape of the device is changed, an auxiliary elasticity finite-element problem must be solved....
-
Single and Dual-GPU Generalized Sparse Eigenvalue Solvers for Finding a Few Low-Order Resonances of a Microwave Cavity Using the Finite-Element Method
PublicationThis paper presents two fast generalized eigenvalue solvers for sparse symmetric matrices that arise when electromagnetic cavity resonances are investigated using the higher-order finite element method (FEM). To find a few loworder resonances, the locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm with null-space deflation is applied. The computations are expedited by using one or two graphical processing...
-
Validated HPTLC method for determination of ledol and alloaromadendrene in the essential oil fractions of Rhododendron tomentosum plants and in vitro cultures and bioautography for their activity screening
PublicationRhododendron tomentosum (Ledum palustre) is a bog shrub used in traditional medicine for treatment of respiratory and rheumatic diseases. Due to the large variability of the chemical composition of its essential oil, depending on the habitat, the in vitro cultures were established as the alternative source of the volatile fraction. There is a need to monitor a quality of the field grown as well as in vitro plant material, especially...
-
Application of the Monte Carlo algorithm for solving volume integral equation in light scattering simulations
PublicationVarious numerical methods were proposed for analysis of the light scattering phenomenon. Important group of these methods is based on solving the volume integral equation describing the light scattering process. The popular method from this group is the discrete dipole approximation (DDA). DDA uses various numerical algorithms to solve the discretized integral equation. In the recent years, the application of the Monte Carlo (MC)...
-
GPU-Accelerated Finite-Element Matrix Generation for Lossless, Lossy, and Tensor Media [EM Programmer's Notebook]
PublicationThis paper presents an optimization approach for limiting memory requirements and enhancing the performance of GPU-accelerated finite-element matrix generation applied in the implementation of the higher-order finite-element method (FEM). It emphasizes the details of the implementation of the matrix-generation algorithm for the simulation of electromagnetic wave propagation in lossless, lossy, and tensor media. Moreover, the impact...
-
Finite element matrix generation on a GPU
PublicationThis paper presents an efficient technique for fast generation of sparse systems of linear equations arising in computational electromagnetics in a finite element method using higher order elements. The proposed approach employs a graphics processing unit (GPU) for both numerical integration and matrix assembly. The performance results obtained on a test platform consisting of a Fermi GPU (1x Tesla C2075) and a CPU (2x twelve-core...
-
Performance evaluation of the parallel object tracking algorithm employing the particle filter
PublicationAn algorithm based on particle filters is employed to track moving objects in video streams from fixed and non-fixed cameras. Particle weighting is based on color histograms computed in the iHLS color space. Particle computations are parallelized with CUDA framework. The algorithm was tested on various GPU devices: a desktop GPU card, a mobile chipset and two embedded GPU platforms. The processing speed depending on the number...
-
Accuracy, Memory and Speed Strategies in GPU-based Finite-Element Matrix-Generation
PublicationThis paper presents strategies on how to optimize GPU-based finite-element matrix-generation that occurs in the finite-element method (FEM) using higher order curvilinear elements. The goal of the optimization is to increase the speed of evaluation and assembly of large finite-element matrices on a single GPU (Graphics Processing Unit) while maintaining the accuracy of numerical integration at the desired level. For this reason,...
-
Communication and Load Balancing Optimization for Finite Element Electromagnetic Simulations Using Multi-GPU Workstation
PublicationThis paper considers a method for accelerating finite-element simulations of electromagnetic problems on a workstation using graphics processing units (GPUs). The focus is on finite-element formulations using higher order elements and tetrahedral meshes that lead to sparse matrices too large to be dealt with on a typical workstation using direct methods. We discuss the problem of rapid matrix generation and assembly, as well as...
-
Vortex flow caused by periodic and aperiodic sound in a relaxing maxwell fluid
PublicationThis paper concerns the description of vortex flow generated by periodic and aperiodic sound in relaxing Maxwell fluid. The analysis is based on governing equation of vorticity mode, which is a result of decomposition of the hydrodynamic equations for fluid flow with relaxation and thermal conductivity into acoustical and non-acoustical parts. The equation governing vorticity mode uses only instantaneous, not averaged over sound...
-
A review on analytical models of brushless permanent magnet machines
PublicationThis study provides an in-depth investigation of the use of analytical and numerical methods in analyzing electrical machines. Although numerical models such as the finite-element method (FEM) can handle complex geometries and saturation effects, they have significant computational burdens, are time-consuming, and are inflexible when it comes to changing machine geometries or input values. Analytical models based on magnetic equivalent...
-
Efficiency of acoustic heating in the Maxwell fluid
PublicationThe nonlinear effects of sound in a fluid describing by the Maxwell model of the viscous stress tensor is the subject of investigation. Among other, viscoelastic biological media belong to this non-newtonian type of fluids. Generation of heating of the medium caused by nonlinear transfer of acoustic energy, is discussed in details. The governing equation of acoustic heating is derived by means of the special linear combination...
-
Efficiency of acoustic heating in the Maxwell fluid
PublicationThe nonlinear effects of sound in a fluid describing by the Maxwell model of the viscous stress tensor is the subject of investigation. Among other, viscoelastic biological media belong to this non-newtonian type of fluids. Generation of heating of the medium caused by nonlinear transfer of acoustic energy, is discussed in details. The governing equation of acoustic heating is derived by means of the special linear combination...
-
Fractional Spectral and Fractional Finite Element Methods: A Comprehensive Review and Future Prospects
PublicationIn this article, we will discuss the applications of the Spectral element method (SEM) and Finite element Method (FEM) for fractional calculusThe so-called fractional Spectral element method (f-SEM) and fractional Finite element method (f-FEM) are crucial in various branches of science and play a significant role. In this review, we discuss the advantages and adaptability of FEM and SEM, which provide the simulations of fractional...
-
Reduced order models in computational electromagnetics (in memory of Ruediger Vahldieck)
PublicationThis paper reviews research of Ruediger Vahldieck's group and the group at the Gdansk University of Technology in the area of model order reduction techniques for accelerating full-wave simulations. The applications of reduced order models to filter design as well as of local and nested(multilevel) macromodels for solving 3D wave equations and wave-guiding problems using finite difference and finite element methods are discussed.
-
Optimization of Data Assignment for Parallel Processing in a Hybrid Heterogeneous Environment Using Integer Linear Programming
PublicationIn the paper we investigate a practical approach to application of integer linear programming for optimization of data assignment to compute units in a multi-level heterogeneous environment with various compute devices, including CPUs, GPUs and Intel Xeon Phis. The model considers an application that processes a large number of data chunks in parallel on various compute units and takes into account computations, communication including...
-
Macro-elements and Model Order Reduction for Efficient Three-Dimensional FEM Analysis
PublicationAn efficient model order reduction (MOR) methodology for three dimensional vector finite element method (FEM) is developed to accelerate simulations of the structures containing features that cause strong variations of mesh density. As the result of presented algorithm, FEM subsystems of equations corresponding to the selected refined region are converted into a very compact sets of linear equations, called macro-elements.Numerical...
-
Diffusion equations with spatially dependent coefficients and fractal Cauer-type networks
PublicationIn this article, we formulate and solve the representation problem for diffusion equations: giving a discretization of the Laplace transform of a diffusion equation under a space discretization over a space scale determined by an increment h > 0, can we construct a continuous in h family of Cauer ladder networks whose constitutive equations match for all h > 0 the discretization. It is proved that for a finite differences discretization...
-
Tuning matrix-vector multiplication on GPU
PublicationA matrix times vector multiplication (matvec) is a cornerstone operation in iterative methods of solving large sparse systems of equations such as the conjugate gradients method (cg), the minimal residual method (minres), the generalized residual method (gmres) and exerts an influence on overall performance of those methods. An implementation of matvec is particularly demanding when one executes computations on a GPU (Graphics...
-
Efficiency of acoustic heating produced in the thermoviscous flow of a fluid with relaxation
PublicationInstantaneous acoustic heating of a fluid with thermodynamic relaxation is the subject of investigation. Among others, viscoelastic biological media described by the Maxwell model of the viscous stress tensor, belong to this type of fluid. The governing equation of acoustic heating is derived by means of the special linear combination of conservation equations in differential form, allowing the reduction of all acoustic terms in...
-
DC-link voltage balancing in cascaded H-Bridge converters
PublicationIn the paper a DC-link voltage balancing strategy for multilevel Cascaded H-Bridge (CHB) converter is proposed. Presented solution bases on optimal choice of active vector durations in Space-Vector Pulse Width Modulation (SV-PWM). It makes it possible to DC-link voltages control and to properly generate the output voltage vector in the case of DC-link voltage unbalance. Results of simulation and experimental researches...
-
Rotational Design Space Reduction for Cost-Efficient Multi-Objective Antenna Optimization
PublicationCost-efficient multi-objective design of antenna structures is presented. Our approach is based on design space reduction algorithm using auxiliary single-objective optimization runs and coordinate system rotation. The initial set of Pareto-optimal solutions is obtained by optimizing a response surface approximation model established in the reduced space using coarse-discretization EM simulation data. The optimization engine is...
-
Deflated Preconditioned Solvers for Parametrized Local Model Order Reduction
PublicationOne of steps in the design of microwave filters is numerical tuning using full-wave simulators. Typically, it is a time-consuming process as it uses advanced computational methods, e.g. the finite-element method (FEM) and it usually requires multiple optimization steps before the specification goals are met. FEM involves solving a large sparse system of equations at many frequency points and therefore its computational cost is...
-
Geometrically nonlinear FEM analysis of 6-parameter resultant shell theory based on 2-D Cosserat constitutive model
PublicationWe develop the elastic constitutive law for the resultant statically and kinematically exact, nonlinear, 6-parameter shell theory. The Cosserat plane stress equations are integrated through-the- thickness under assumption of the Reissner-Mindlin kinematics. The resulting constitutive equations for stress resultant and couple resultants are expressed in terms of two micropolar constants: the micropolar modulus Gc and the micropolar...
-
Investigation of sour substances by five-channel potentiometric taste sensor containing all-solid-state-electrodes ( ASSEs)
PublicationAn elaboration of taste sensor for discrimination of different food products is of great importance for food industry. Potentiometric taste sensor containing ion selective electrodes with lipid/polymer membranes has already been applied commercially in food industry. However, time-consuming and demanding preconditioning method of ions selective electrodes as well as maintenance of electrodes’ bodies are disadvantages of this taste...
-
INVESTIGATION OF SOUR SUBSTANCES BY FIVE-CHANNEL POTENTIOMETRIC TASTE SENSOR CONTAINING ALL-SOLID-STATE- ELECTRODES (ASSEs)
PublicationAn elaboration of taste sensor for discrimination of different food products is of great importance for food industry. Potentiometric taste sensor containing ion selective electrodes with lipid/polymer membranes has already been applied commercially in food industry. However, time-consuming and demanding preconditioning method of ions selective electrodes as well as maintenance of electrodes’ bodies are disadvantages of this taste...
-
The influence of stainless steels microstructure evolution on the mechanical properties of pressure installation elements
PublicationIn this paper results of FEM analysis performed for pipe and 90o elbow, made of super duplex stainless steel, after precipitation of sigma phase (),were presented. Heat treatment conditions simulate accidental overheating that can occur in nuclear installations after cooling circulation fails. For initial and modified microstructure tensile test was performed and numerical description of work hardening curves was done. Verification...
-
The influence of stainless steels microstructure evolution on the mechanical properties of pressure installation elements
PublicationIn this paper results of FEM analysis performed for pipe and 90o elbow, made of super duplex stainless steel, after precipitation of sigma phase (s), were presented. Heat treatment conditions simulate accidental overheating that can occur in nuclear installations after cooling circulation fails. For initial and modified microstructure tensile test was performed and numerical description of work hardening curves was done. Verification...
-
Smooth Particle Hydrodynamics (SPH) approach in simulating large penetration into soil
PublicationA study of Smooth Particle Hydrodynamics (SPH) approach for predicting large soil deformation is presented. Theoretical basics of SPH method, including the equations governing, discussion of the importance of smoothing function length, contact formulation, boundary treatment and finally utilization in hydrocodes simulations are presented. An application of SPH to a real case of large penetrations (crater creating) into soil caused...
-
Numerical Issues and Approximated Models for the Diagnosis of Transmission Pipelines
PublicationThe chapter concerns numerical issues encountered when the pipeline flow process is modeled as a discrete-time state-space model. In particular, issues related to computational complexity and computability are discussed, i.e., simulation feasibility which is connected to the notions of singularity and stability of the model. These properties are critical if a diagnostic system is based on a discrete mathematical model of the flow...
-
Theoretical analysis of the cable-stayed bridge over brda river in Bydgoszcz
PublicationThe topic of this article is the numerical analysis carried out for the cable-stayed tram (road) bridge over Brda river in Bydgoszcz. The goal of numerical studies was the verification of project assumptions and construction solutions elaborated by the autonomous team. The bridge calculations were conducted on three independent FEM models The internal forces and the elements effort under different loading schemes were examined....