Search results for: FINITE-DIFFERENCE
-
Finite-difference analysis of a loaded hemispherical resonator
PublicationPrzedstawiono analizę osiowo symetrycznej struktury otwartego rezonatora he-misferycznego zawierającego warstwowe ośrodki dielektryczne. W analizie wy-korzystano metodę różnic skończonych w dziedzinie częstotliwości (FDFD). Wcelu osiągnięcia wysokiej dokładności obliczeń zastosowano zaawansowane me-tody modelowania warunków brzegowych i granicznych.
-
Simple and accurate field interpolation in finite difference methods
PublicationW artykule wprowadzono nową metodę interpolacji pól w metodzie różnic skończonych. Metoda jest prosta i daje wyniki bardziej dokładne niż w przypadku interpolacji liniowej. Proponowana metoda interpolacji może zostać użyta do poprawienia dokładności takich metod, jak subgridding.
-
Efficient finite difference analysis of microstructured optical fibers [on line]
PublicationW pracy przedstawiono nową technikę numerycznej analizy mikrostruktur światłowodowych. Technika ta polega na połączeniu metody różnic skończonych i dyskretnych rozwinięć funkcyjnych. Pozwala to na zmniejszenie rozmiaru macierzy powstających przy analizie problemu, a jednocześnie na bezpośrednią implementację radiacyjnych warunków brzegowych. Testy algorytmu przeprowadzono dla kilku wybranych struktur, uzyskując bardzo dobrą zgodność...
-
The finite-difference simulation of x-rays propagation through a system of lenses
PublicationThe propagation of X-ray waves through an optical system consisting of 33 aluminum X-ray refractive lenses is considered. For solving the problem, a finite-difference method is suggested and investigated. It is shown that very small steps of the difference grid are necessary for reliable computation of propagation of X-ray waves through the system of lenses. It is shown that the wave phase is a function very quickly increasing...
-
On the convergence of a nonlinear finite-difference discretization of the generalized Burgers–Fisher equation
PublicationIn this note, we establish analytically the convergence of a nonlinear finite-difference discretization of the generalized Burgers-Fisher equation. The existence and uniqueness of positive, bounded and monotone solutions for this scheme was recently established in [J. Diff. Eq. Appl. 19, 1907{1920 (2014)]. In the present work, we prove additionally that the method is convergent of order one in time, and of order two in space. Some...
-
Applications of the discrete green's function in the finite-difference time-domain method
PublicationIn this paper, applications of the discrete Green's function (DGF) in the three-dimensional (3-D) finite-difference time-domain (FDTD) method are presented. The FDTD method on disjoint domains was developed employing DGF to couple the subdomains as well as to compute the electromagnetic field outside these subdomains. Hence, source and scatterer are simulated in separate subdomains and updating of vacuum cells, being of little...
-
Finite difference approximations for nonlinear first order partial differential equations
PublicationKlasyczne rozwiązania nieliniowych równań różniczkowych o pochodnych cząst-kowych pierwszego rzędu są aproksymowane za pomocą równań quasiliniowych uk-ładów równań różniczkowych. Dowód zbieżności wykorzystuje metody porównawcze
-
Optimizing the computation of a parallel 3D finite difference algorithm for graphics processing units
PublicationThis paper explores the possibilities of using a graphics processing unit for complex 3D finite difference computation via MUSTA‐FORCE and WENO algorithms. We propose a novel algorithm based on the new properties of CUDA surface memory optimized for 2D spatial locality and compare it with 3D stencil computations carried out via shared memory, which is currently considered to be the best approach. A case study was performed for...
-
Finite-difference time-domain analyses of active cloaking for electrically-large objects
PublicationInvisibility cloaking devices constitute a unique and potentially disruptive technology, but only if they can work over broad bandwidths for electrically-large objects. So far, the only known scheme that allows for broadband scattering cancellation from an electrically-large object is based on an active implementation where electric and magnetic sources are deployed over a surface surrounding the object, but whose ‘switching on’...
-
On the convergence of a finite-difference discretization à la Mickens of the generalized Burgers–Huxley equation
PublicationIn this note, we establish the property of convergence for a finite-difference discretization of a diffusive partial differential equation with generalized Burgers convective law and generalized Hodgkin–Huxley reaction. The numerical method was previously investigated in the literature and, amongst other features of interest, it is a fast and nonlinear technique that is capable of preserving positivity, boundedness and monotonicity....
-
The finite difference methods of computation of X-rays propagation through a system of many lenses
PublicationThe propagation of X-ray waves through an optical system consisting of many beryllium X-ray refrac- tive lenses is considered. In order to calculate the propagation of electromagnetic in the optical sys- tem, two differential equations are considered. First equation for an electric field of a monochromatic wave and the second equation derived for complex phase of the same electric The propagation of X-ray waves through an optical system...
-
A simple finite difference approach using unstructured meshes from FEM mesh generators.
PublicationPrezentujemy metodologię używania dwuwymiarowych różnic skończonych wykorzystujących siatki trójkątne z metody FEM. Określiliśmy wyrażenia opisujące pierwsząi drugą pochodną w oparciu o wartości pól w saśiednich punktach. Metodę przetestowaliśmy na strukturze rezonatora z prostokątnym rdzeniem dielektrycznym.
-
A fast high-resolution 3-D Finite Difference Time-Domain scheme with macromodels.
PublicationZaprezentowano szybki trójwymiarowy wysokorozdzielczy algorytm różnic skończonych w dziedzinie czasu. Zwiększona rozdzielczość osiągnięta została poprzez połączenie standardowego schematu FDTD z makromodelami silnie zagęszczonych objętości stworzonych przez zastosowanie redukcji rzędu modelu. Nowy schemat FDTD z makromodelami jest wyprowadzony. Koszt numeryczny nowej metody jest wyestymowany. Doskonałe wyniki zostały zaobserwowane...
-
Approximate analytical boundary conditions for efficient finite difference frequency domain simulations in cylindrical coordinates
PublicationW artykule zaprezentowano prostą technikę analizy rezonatora otwartego. Algorytm łączy w sobie metodę różnic skończonych i rozwinięć funkcyjnych , umożliwiając implementację warunków brzegowych symulujących otwartą przestrzeń. Metoda testowana była w analizie rezonatorów o różnych wymiarach,a otrzymane wyniki dobrze zgadzały się z rezultatami innych metod.
-
Efficient Multiscale Finite Difference Frequency Domain Analysis Using Multiple Macromodels with Compressed Boundaries
PublicationW niniejszym artykule zaprezentowany został nowatorski pomysł na zmniejszenie numerycznej złożoności metody różnic skończonych przy użyciu wielu makromodeli. Skuteczność techniki macromodelingu zależy od liczby portów modelu. Aby zwiększyć efektywność algorytmu, próbki pola na granicy makromodeli są zastępowane przez amplitudy wielomianów Legendre'a. Przedefiniowanie problemu w taki sposób powoduje znaczną redukcję czasu analizy....
-
A new conformal radiation boundary condition for high accuracy finite difference analysis of open waveguides
PublicationW publikacji zaprezentowano nowe warunki brzegowe przeznaczone dla precyzyjnej analizy otwartych prowadnic falowych metodą różnic skończonych. Praca zawiera również liczne przykłady zastosowania proponowanej techniki oraz porównanie otrzymanych rezultatów numerycznych z wynikami otrzymywanymi przy użyciu innych metod.
-
Projection framework for hybrid methods derived from finite difference operators in time and frequency domain
PublicationW artykule przedstawiono ogólną koncepcję tworzenia algorytmów hybrydowych na bazie operatorowego sformułowania metody różnic skończonych. Wykorzystano koncepcję projekcji ortogonalnej w przestrzeni skończeniewymiarowej w celu modyfilacji pierwotnego sformułowania.
-
Application of hybrid finite-difference mode-matching method to analysis of structures loaded with axially-symmetrical posts
PublicationW artykule przedstawiono nową metodę hybrydową do analizy układów falowodowych zawierających dowolne konfiguracje obiektów osiowo-symetrycznych. Metoda oparta jest na połączeniu metody różnic skończonych, metody dopasowania rodzajów oraz iteracyjnej procedury rozpraszania. W pracy przedstawiono badania zbieżności metody. Uzyskane wyniki numeryczne porównano z wynikami odniesienia. Duża zgodność wyników potwierdziła poprawność opracowanego...
-
Existence and uniqueness of monotone and bounded solutions for a finite-difference discretization a` la Mickens of the generalized Burgers–Huxley equation.
PublicationDeparting from a generalized Burgers–Huxley partial differential equation, we provide a Mickens-type, nonlinear, finite-difference discretization of this model. The continuous system is a nonlinear regime for which the existence of travelling-wave solutions has been established previously in the literature. We prove that the method proposed also preserves many of the relevant characteristics of these solutions, such as the positivity,...
-
Hybrid finite-difference/mode-matching method for analysis of scattering from arbitrary configuration of rotationally-symmetrical posts
PublicationW artykule przedstawiona została metoda hybrydowa pozwalająca na analizę zjawiska rozpraszania od dowolnej konfiguracji obiektów o osiowej symetrii. Proponowana metoda oparta jest na reprezentacji obiektu lub dowolnej konfiguracji obiektów poprzez zastępczy obiekt kulisty ze znanymi warunkami brzegowymi. W analizie pojedynczego obiektu wykorzystana została metoda hybrydowa stanowiąca połączenie metody różnic skończonych z metodą...
-
Application of a hybrid finite difference-mode matching method to the analysis of circular cavities loaded with cylindrical metallic objects
PublicationW artykule przedstawiono badania rezonatorów cylindrycznych zawierających obiekty cylindryczne o dowolnym przekroju poprzecznym z wykorzystaniem metody hybrydowej stanowiącej połączenie metody różnic skończonych z metodą dopasowania rodzajów. W ramach badań określono częstotliwości rezonansowe dla wybranych struktur rezonatorów. Wyniki numeryczne zostały zweryfikowane eksperymentalnie.
-
Analysis of scattering from arbitrary configuration of cylindrical objects using hybrid finite-difference mode-matching method
PublicationW artykule przedstawiono nową metodę hybrydową stanowiącą połączenie metody różnic skończonych z metodą dopasowania rodzajów, pozwalającą na analizę obiektów cylindrycznych o dowolnym przekroju poprzecznym. W prezentowanej metodzie poszczególne obiekty reprezentowane są poprzez wielorodzajową macierz impednacyjną zdefiniowaną na powierzchni bocznej zastępczego cylindra zawierającego badany obiekt. Celem określenia parametrów rozproszenia...
-
Analysis of waveguide junctions containing axially-symetrical posts using hybrid finite-difference mode-matching technique
PublicationArtykuł opisuję hybrydową metodę analizy złącz falowodowych zawierających osiowo-symetryczne obiekty rozpraszające. Metoda jest oparta na metodzie dopasowania rodzajów oraz metodzie różnic skończonych w dziedzinie częstotliwości.
-
Data obtained by numerical simulation for X-ray focusing using a finite difference method
Open Research DataThe propagation of X-ray waves through an optical system consisting of many X-ray refractive lenses is considered. For solving the problem for an electromagnetic wave, a finite-difference method is applied.
-
3-D finite-difference time-domain modelling of ground penetrating radar for identification of rebars in complex reinforced concrete structures
PublicationThis paper presents numerical and experimental investigations to identify reinforcing bars using the ground penetrating radar (GPR) method. A novel element of the paper is the inspection of different arrangements of reinforcement bars. Two particular problems, i.e. detection of few adjacent transverse bars and detection of a longitudinal bar located over or under transverse reinforcement, have been raised. An attention was also...
-
Macromodeling in Finite Differences
PublicationRozdział opisuje technikę wykorzystania procesu makromodelowania (wykorzystanie makromodelu do opisu równań różnicowych) w analizie różnic skończonych w dziedzinie czasu. Zawarte informacje pozwalają na samodzielne zaimplementowanie metody i pokazują możliwości poprawienia uzyskiwanych wyników.
-
Improving the accuracy of subgridding scheme in finite differences method based on Legendre polynomials expansion
PublicationIn this article the Legendre polynomials have been used to interpolate the field at the boundary of the meshes of different densities. The numerical verification of the proposed technique has been carried out in frequency domain. It has been shown that the accuracy of the presented method is very high and stable - the error monotonically decreases as a function of the refinement factor.
-
Optimization of the Hardware Layer for IoT Systems using a Trust Region Method with Adaptive Forward Finite Differences
PublicationTrust-region (TR) algorithms represent a popular class of local optimization methods. Owing to straightforward setup and low computational cost, TR routines based on linear models determined using forward finite differences (FD) are often utilized for performance tuning of microwave and antenna components incorporated within the Internet of Things systems. Despite usefulness for design of complex structures, performance of TR methods...
-
Numerical optimization of planar antenna structures using trust-region algorithm with adaptively adjusted finite differences
Open Research DataThe dataset contains initial designs and optimization results for three planar structures that include quasi-patch antenna for WLAN applications, compact spline-parameterized monopole dedicated for ultra-wideband applications, as well as rectifier for energy harvesting with enhanced bandwidth. The numerical results for the first two structures are also...
-
Finite Differences-Finite Elements-Finite Volumes-Boundary Elements
Conferences -
Piotr Sypek dr inż.
PeoplePiotr Sypek received the M.S.E.E. and Ph.D. degrees (with hons.) in microwave engineering from the Gdańsk University of Technology, Gdańsk, Poland, in 2003 and 2012, respectively. He was involved in the design and implementation of parallel algorithms for the formulation and solution of electromagnetic problems executed on CPUs (workstations and clusters) and GPUs. His current research interests include parallel processing in computational...
-
Database of the illustrative simulations of the nonstandard approximation of the generalized Burgers–Huxley equation
Open Research DataThe presented dataset is a result of numerical analysis of a generalized Burgers–Huxley partial differential equation. An analyzed diffusive partial differential equation consist with nonlinear advection and reaction. The reaction term is a generalized form of the reaction law of the Hodgkin–Huxley model, while the advection is a generalized form of...
-
Simulation of unsteady flow over floodplain using the diffusive wave equation and the modified finite element method
PublicationWe consider solution of 2D nonlinear diffusive wave equation in a domain temporarily covered by a layer of water. A modified finite element method with triangular elements and linear shape functions is used for spatial discretization. The proposed modification refers to the procedure of spatial integration and leads to a more general algorithm involving a weighting parameter. The standard finite element method and the finite difference...
-
Database of the convergence analysis results of the nonstandard approximation of the generalized Burgers–Huxley equation for the solution bounded within [0,1].
Open Research DataThe presented dataset is a result of the convergence analysis of the Mickens-type, nonlinear, finite-difference discretization of a generalized Burgers–Huxley partial differential equation.
-
Database of the convergence analysis results of the nonstandard approximation of the generalized Burgers–Huxley equation for the solution bounded within [0, γ^(1/p)].
Open Research DataPresented dataset is a result of the convergence analysis of the Mickens-type, nonlinear, finite-difference discretization of a generalized Burgers–Huxley partial differential equation. The generalized Burgers–Huxley equation is a diffusive partial differential equation with nonlinear advection and diffusion. The boundary problem for this equation possesses...
-
Asymptotic Expansion Method with Respect to Small Parameter for Ternary Diffusion Models
PublicationTernary diffusion models lead to strongly coupled systems of PDEs. We choose the smallest diffusion coefficient as a small parameter in a power series expansion whose components fulfill relatively simple equations. Although this series is divergent, one can use its finite sums to derive feasible numerical approximations, e.g. finite difference methods (FDMs).
-
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.
-
Resonance Frequency Calculation of a Multilayer and Multipatch Spherical Microstrip Structure Using a Hybrid Technique
PublicationThis communication offers a rigorous analysis of the resonance frequency problem of a spherical microstrip structure mounted on a multilayer, dielectric-coated metallic sphere, with an electrically small radius. The structure consists of single or multiple metallic patches with arbitrary shapes. A full-wave analysis is employed with the use of proposed hybrid approach, combining the finite-difference technique with a spectral domain...
-
FEM and experimental investigations of concrete temperature field in the massive stemwall of the bridge abutment
PublicationThe paper deals with the prediction of early-age concrete temperature of cast-in-place stemwall of the bridge abutment. The considered object is an arch bridge located in Gda´nsk. In the case of massive structures, it is particularly important to not exceed the temperature difference between the core and the concrete surface. Too high temperature gradient generates an increase in thermal stresses, what could be the reason of exceeding...
-
A high-accuracy complex-phase method of simulating X-ray propagation through a multi-lens system
PublicationThe propagation of X-ray waves through an optical system consisting of many X-ray refractive lenses is considered. For solving the problem for an electromagnetic wave, a finite-difference method is applied. The error of simulation is analytically estimated and investigated. It was found that a very detailed difference grid is required for reliable and accurate calculations of the propagation of X-ray waves through a multi-lens...
-
Analysis of radiation and scattering problems with the use of hybrid techniques based on the discrete Green's function formulation of the FDTD method
PublicationIn this contribution, simulation scenarios are presented which take advantage of the hybrid techniques based on the discrete Green's function formulation of the finite-difference time-domain (DGF-FDTD) method. DGF-FDTD solutions are compatible with the finite-difference grid and can be applied for perfect hybridization of the FDTD method. The following techniques are considered: (i) DGF-FDTD for antenna simulations, (ii) DGF-based...
-
Updating Finite Element Model of a Wind Turbine Blade Section Using Experimental Modal Analysis Results
PublicationThis paper presents selected results and aspects of themultidisciplinary and interdisciplinary research oriented for the experimental and numerical study of the structural dynamics of a bend-twist coupled full scale section of awind turbine blade structure.Themain goal of the conducted research is to validate finite elementmodel of themodified wind turbine blade section mounted in the flexible support structure accordingly to the...
-
Electromagnetic plane wave scattering from a cylindrical object with an arbitrary cross section using a hybrid technique
PublicationA hybrid technique combining finite-element and mode-matching methods for the analysis of scattering problems in open and closed areas is presented. The main idea of the analysis is based on the utilization of the finite-element method to calculate the post impedance matrix and combine it with external excitation. The discrete analysis, which is the most time- and memory-consuming, is limited here only to the close proximity of...
-
Analysis of Floodplain Inundation Using 2D Nonlinear Diffusive Wave Equation Solved with Splitting Technique
PublicationIn the paper a solution of two-dimensional (2D) nonlinear diffusive wave equation in a partially dry and wet domain is considered. The splitting technique which allows to reduce 2D problem into the sequence of one-dimensional (1D) problems is applied. The obtained 1D equations with regard to x and y are spatially discretized using the modified finite element method with the linear shape functions. The applied modification referring...
-
Coherent-wave Monte Carlo method for simulating light propagation in tissue
PublicationSimulating propagation and scattering of coherent light in turbid media, such as biological tissues, is a complex problem. Numerical methods for solving Helmholtz or wave equation (e.g. finite-difference or finite-element methods) require large amount of computer memory and long computation time. This makes them impractical for simulating laser beam propagation into deep layers of tissue. Other group of methods, based on radiative...
-
Determination of time delay between ventricles contraction using impedance measurements
PublicationThe paper presents a novel approach to assessment of ventricular dyssynchrony basing on multichannel electrical impedance measurements. Using a proper placement of electrodes, the sensitivity approach allows estimating time difference between chambers contraction from over determined nonlinear system of equations. The theoretical considerations which include Finite Element Method simulations were verified using measurements on...
-
How to render FDTD computations more effective using agraphics accelerator.
PublicationGraphics processing units (GPUs) for years have been dedicated mostly to real time rendering. Recently leading GPU manufactures have extended their research area and decided to support also graphics computing. In this paper, we describe an impact of new GPU features on development process of an efficient finite difference time domain (FDTD) implementation.
-
GPR simulations for diagnostics of a reinforced concrete beam
PublicationThe most popular technique for modelling of an electromagnetic field, the finite difference time domain (FDTD) method, has recently become a popular technique as an interpretation tool for ground penetrating radar (GPR) measurements. The aim of this study is to detect the size and the position of damage in a reinforced concrete beam using GPR maps. Numerical simulations were carried out using the finite differ-ence time domain...
-
Using GPUs for Parallel Stencil Computations in Relativistic Hydrodynamic Simulation
PublicationThis paper explores the possibilities of using a GPU for complex 3D finite difference computation. We propose a new approach to this topic using surface memory and compare it with 3D stencil computations carried out via shared memory, which is currently considered to be the best approach. The case study was performed for the extensive computation of collisions between heavy nuclei in terms of relativistic hydrodynamics.
-
The modelling method of discrete-continuous systems
PublicationThe paper introduces a method of discrete-continuous systems modelling. In the proposed method a three-dimensional system is divided into finite elements in only two directions, with the third direction remaining continuous. The thus obtained discrete-continuous model is described by a set of partial differential equations. General difference equations of discrete system are obtained using the rigid finite element method. The limit...
-
Finite element models used in diagnostics of transverse cracks in bridge approach pavement
Open Research DataTransverse cracks in the asphalt pavement were observed on bridge structures next to single-module expansion joints with a 5 meter approach slab set at the depth of 1 m. The finite element (FE) models of the approach pavement were created to investigate the reasons of premature cracking and crack initiation mechanism over the back edge of the abutment...
-
OpenGL accelerated method of the material matrix generation for FDTD simulations
PublicationThis paper presents the accelerated technique of the material matrix generation from CAD models utilized by the finite-difference time-domain (FDTD) simulators. To achieve high performance of these computations, the parallel-processing power of a graphics processing unit was employed with the use of the OpenGL library. The method was integrated with the developed FDTD solver, providing approximately five-fold speedup of the material...
-
Hybrid technique for the analysis of circular waveguide junctions loaded with ferrite posts
PublicationThis study presents a hybrid technique for the analysis of circular waveguide junctions loaded with axially symmetrical ferrite posts of irregular shape. The method is based on a combination of the finite-difference frequency- domain technique with a mode-matching technique. The proposed approach is validated by comparing the presented results with numerical ones obtained from commercial software. The application of a cylindrical...
-
A high-accuracy method of computation of x-ray waves propagation through an optical system consisting of many lenses
PublicationThe propagation of X-ray waves through an optical system consisting of many X-ray refractive lenses is considered. Two differential equations are contemplated for solving the problem for electromagnetic wave propagation: first – an equation for the electric field, second – an equation derived for a complex phase of an electric field. Both equations are solved by the use of a finite-difference method. The simulation error is estimated...
-
Modelowanie przepływu ustalonego niejednostajnego w sieciach kanałów otwartych z uwzględnieniem obiektów hydrotechnicznych
PublicationW pracy sformułowano zagadnienie brzegowe dla równania energii opisującego przepływ ustalony niejednostajny i przedstawiono sposób jego rozwiązania przy pomocy metody różnicowej. Zaproponowana metoda obliczeń nadaje się do analizy przepływu w dendrycznych i pierścieniowych sieciach kanałów otwartych. Ponadto na przykładzie przelewu prostokątnego zaproponowano metodę uwzględnienia w obliczeniach zabudowy hydrotechnicznej. Słowa...
-
Numerical Methods for Partial Differential Equations
e-Learning CoursesCourse description: This course focuses on modern numerical techniques for linear and nonlinear elliptic, parabolic and hyperbolic partial differential equations (PDEs), and integral equations fundamental to a large variety of applications in science and engineering. Topics include: formulations of problems in terms of initial and boundary value problems; finite difference and finite element discretizations; boundary element approach;...
-
Recurrence scheme for FDTD-compatible discrete Green's function derived based on properties of Gauss hypergeometric function
PublicationIn this paper, the formulation of one-dimensional FDTD (Finite-difference time-domain)-compatible discrete Green's function (DGF) is derived based on the Gauss hypergeometric function (GHF). The properties of GHF make it possible to derive the recurrence scheme only in the time domain for the DGF generation. Furthermore, this recurrence scheme is valid for any stable time-step size and can be implemented using standard numerical...
-
Resonance Frequency Calculation of Spherical Microstrip Structure Using Hybrid Technique
PublicationIn this paper the spherical microstrip structure is considered. The structure is composed of a metallic patch with an arbitrary shape placed on a dielectric coated metallic sphere. In the analysis the hybrid technique is utilized. In this approach the finite-difference technique is applied in a cavity model to determine the current basis functions on the patch. Next, using method of moments, the resonance frequency of the structure...
-
Towards an efficient multi-stage Riemann solver for nuclear physics simulations
PublicationRelativistic numerical hydrodynamics is an important tool in high energy nuclear science. However, such simulations are extremely demanding in terms of computing power. This paper focuses on improving the speed of solving the Riemann problem with the MUSTA-FORCE algorithm by employing the CUDA parallel programming model. We also propose a new approach to 3D finite difference algorithms, which employ a GPU that uses surface memory....
-
Numerical FDM modelling of wave propagation in concrete structure
PublicationThe article presents application of finite difference method to damage detection and its size evaluation in concrete structure by elastic wave propagation method. The simulations of wave propagation in concrete structure were performed for six different damage scenarios. Damages were modelled as areas with changed material properties. Investigation focused on the influence of damage size on the energy of wave reflection. Presented...
-
A New Expression for the 3-D Dyadic FDTD-Compatible Green's Function Based on Multidimensional Z-Transform
PublicationIn this letter, a new analytic expression for the time-domain discrete Green's function (DGF) is derived for the 3-D finite-difference time-domain (FDTD) grid. The derivation employs the multidimensional Z-transform and the impulse response of the discretized scalar wave equation (i.e., scalar DGF). The derived DGF expression involves elementary functions only and requires the implementation of a single function in the multiple-precision...
-
A Fortran-95 algorithm to solve the three-dimensional Higgs boson equation in the de Sitter space-time
Open Research DataA numerically efficient finite-difference technique for the solution of a fractional extension of the Higgs boson equation in the de Sitter space-time is designed. The model under investigation is a multidimensional equation with Riesz fractional derivatives of orders in (0,1)U(1,2], which considers a generalized potential and a time-dependent diffusion...
-
Approximate solution for Euler equations of stratified water via numerical solution of coupled KdV system
PublicationWe consider Euler equations with stratified background state that is valid for internal water waves. The solution of the initial-boundary problem for Boussinesq approximation in the waveguide mode is presented in terms of the stream function. The orthogonal eigenfunctions describe a vertical shape of the internal wave modes and satisfy a Sturm-Liouville problem. The horizontal profile is defined by a coupled KdV system which is...
-
Dyskretno-ciągła metoda modelowania układów dynamicznych
PublicationW artykule przedstawiono oryginalną metodę modelowania układów dyskretno-ciągłych. Metoda polega na dyskretyzowaniu układu trójwymiarowego jedynie w dwóch wybranych kierunkach. W trzecim z kierunków układ pozostaje ciągły. Otrzymany w ten sposób model jest modelem dyskretno-ciągłym. Opisany jest za pomocą równań różniczkowych cząstkowych. Ogólne równania różnicowe układu dyskretnego otrzymano, wykorzystując metodę sztywnych elementów...
-
Exact modal absorbing boundary condition for waveguide simulations - discrete Green's function approach
PublicationA modal absorbing boundary condition (ABC) based on the discrete Green's function (DGF) is introduced and applied for termination of waveguides simulated by means of the finite-difference time-domain (FDTD) method. The differences between the developed approach and implementations already demonstrated in the literature are presented. By applying DGF, a consistent theoretical approach to modal ABC in the FDTD method is obtained....
-
Numerical simulation of hardening of concrete plate
PublicationThe paper presents a theoretical formulation of concrete curing in order to predict temperature evolution and strength development. The model of heat flow is based on a well-known Fourier equation. The numerical solution is implemented by means of the Finite Difference Method. In order to verify the model, the in situ temperature measurements at the top plate of a road bridge were carried out. A high agreement between numerical...
-
Coupled Urban Areas Inundation Model with Interaction Between Storm Water System and Surface Flow - Case Study of Sea Level Impact on Seaside Areas Flooding
PublicationInundations are becoming more frequent than ever. What is connected with increasing area of impervious surface in cities. This makes predicting urban flooding and its scale especially important. At the seaside we observe additional conditions such as sea level that makes accurate numerical modelling of issue even harder. With complex approach to the matter which is simultaneous calculation of storm water conduit flow and overland...
-
Lax-Wendroff and McCormack Schemes for Numerical Simulation of Unsteady Gradually and Rapidly Varied Open Channel Flow
PublicationTwo explicit schemes of the finite difference method are presented and analyzed in the paper. The applicability of the Lax-Wendroff and McCormack schemes for modeling unsteady rapidly and gradually varied open channel flow is investigated. For simulation of the transcritical flow the original and improved McCormack scheme is used. The schemes are used for numerical solution of one dimensional Saint-Venant equations describing free...
-
Comparative analysis of numerical with optical soliton solutions of stochastic Gross–Pitaevskii equation in dispersive media
PublicationThis article deals with the stochastic Gross–Pitaevskii equation (SGPE) perturbed with multiplicative time noise. The numerical solutions of the governing model are carried out with the proposed stochastic non-standard finite difference (SNSFD) scheme. The stability of the scheme is proved by using the Von-Neumann criteria and the consistency is shown in the mean square sense. To seek exact solutions, we applied the Sardar subequation...
-
Numerical modeling of GPR field in damage detection of a reinforced concrete footbridge
PublicationThe paper presents a study on the use of the ground penetrating radar (GPR) method in diagnostics of a footbridge. It contains experimental investigations and numerical analyses of the electromagnetic field propagation using the finite difference time domain method (FDTD). The object of research was a reinforced concrete footbridge over a railway line. The calculations of the GPR field propagation were performed on a selected cross-section...
-
On the crack front curvature in bonded joints
PublicationStandard tests of adhesively bonded specimens are likely to produce heterogeneous stress distribution along the crack front and its vicinity. High separation rate mode I dominated fracture test is performed.Observation of post mortem fractured surfaces with an optical microscope reveals characteristic features of mixed mode I/III fracture near the sides of the specimen but not in the middle. At first, finite elements calculations...
-
Simulating coherent light propagation in a random scattering materials using the perturbation expansion
PublicationMultiple scattering of a coherent light plays important role in the optical metrology. Probably the most important phenomenon caused by multiple scattering are the speckle patterns present in every optical imaging method based on coherent or partially coherent light illumination. In many cases the speckle patterns are considered as an undesired noise. However, they were found useful in various subsurface imaging methods such as...
-
Fast implementation of FDTD-compatible green's function on multicore processor
PublicationIn this letter, numerically efficient implementation of the finite-difference time domain (FDTD)-compatible Green's function on a multicore processor is presented. Recently, closed-form expression of this discrete Green's function (DGF) was derived, which simplifies its application in the FDTD simulations of radiation and scattering problems. Unfortunately, the new DGF expression involves binomial coefficients, whose computations...
-
Bending analysis of functionally graded nanoplates based on a higher-order shear deformation theory using dynamic relaxation method
PublicationIn this paper, bending analysis of rectangular functionally graded (FG) nanoplates under a uniform transverse load has been considered based on the modified couple stress theory. Using Hamilton’s principle, governing equations are derived based on a higher-order shear deformation theory (HSDT). The set of coupled equations are solved using the dynamic relaxation (DR) method combined with finite difference (FD) discretization technique...
-
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...
-
Alternative cogeneration thermodynamic cycles for domestic ORC
PublicationThe Organic Flash Cycle (OFC) is suggested as a vapor power cycle that could potentially improve the efficiency of utilization of the heat source. Low and medium temperature finite thermal sources are considered in the cycle. Additionally the OFC’s aim is to reduce temperature difference during heat addition. The study examines 2 different fluids. Comparisons are drawn between the OFC and an optimized basic Organic Rankine Cycle...
-
Acceleration of the DGF-FDTD method on GPU using the CUDA technology
PublicationWe present a parallel implementation of the discrete Green's function formulation of the finite-difference time-domain (DGF-FDTD) method on a graphics processing unit (GPU). The compute unified device architecture (CUDA) parallel computing platform is applied in the developed implementation. For the sake of example, arrays of Yagi-Uda antennas were simulated with the use of DGF-FDTD on GPU. The efficiency of parallel computations...
-
Parallel implementation of the DGF-FDTD method on GPU Using the CUDA technology
PublicationThe discrete Green's function (DGF) formulation of the finite-difference time-domain method (FDTD) is accelerated on a graphics processing unit (GPU) by means of the Compute Unified Device Architecture (CUDA) technology. In the developed implementation of the DGF-FDTD method, a new analytic expression for dyadic DGF derived based on scalar DGF is employed in computations. The DGF-FDTD method on GPU returns solutions that are compatible...
-
Hybridization of the FDTD method with use of the discrete Green's function
PublicationIn this contribution, a hybrid technique is presented which combines the finite-difference time-domain (FDTD) method and the discrete Green's function (DGF) formulation of this method. FDTD is a powerful technique for the analysis of complex penetrable objects but its application is not efficient when the computational domain includes many free-space cells. Therefore, the hybrid method was developed which is applicable to complex...
-
Implementation of FDTD-Compatible Green's Function on Graphics Processing Unit
PublicationIn this letter, implementation of the finite-difference time domain (FDTD)-compatible Green's function on a graphics processing unit (GPU) is presented. Recently, closed-form expression for this discrete Green's function (DGF) was derived, which facilitates its applications in the FDTD simulations of radiation and scattering problems. Unfortunately, implementation of the new DGF formula in software requires a multiple precision...
-
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...
-
Analytical Expression for the Time-Domain Discrete Green's Function of a Plane Wave Propagating in the 2-D FDTD Grid
PublicationIn this letter, a new closed-form expression for the time-domain discrete Green's function (DGF) of a plane wave propagating in the 2-D finite-difference time-domain (FDTD) grid is derived. For the sake of its verification, the time-domain implementation of the analytic field propagator (AFP) technique was developed for the plane wave injection in 2-D total-field/scattered-field (TFSF) FDTD simulations. Such an implementation of...
-
MEMORY EFFECT ANALYSIS USING PIECEWISE CUBIC B-SPLINE OF TIME FRACTIONAL DIFFUSION EQUATION
PublicationThe purpose of this work is to study the memory effect analysis of Caputo–Fabrizio time fractional diffusion equation by means of cubic B-spline functions. The Caputo–Fabrizio interpretation of fractional derivative involves a non-singular kernel that permits to describe some class of material heterogeneities and the effect of memory more effectively. The proposed numerical technique relies on finite difference approach and cubic...
-
Analytical Expression for the Time-Domain Green's Function of a Discrete Plane Wave Propagating in the 3-D FDTD Grid
PublicationIn this paper, a closed-form expression for the time-domain dyadic Green’s function of a discrete plane wave (DPW) propagating in a 3-D finite-difference time-domain (FDTD) grid is derived. In order to verify our findings, the time-domain implementation of the DPW-injection technique is developed with the use of the derived expression for 3-D total-field/scattered-field (TFSF) FDTD simulations. This implementation requires computations...
-
A DISCRETE-CONTINUOUS METHOD OF MECHANICAL SYSTEM MODELLING
PublicationThe paper describes a discrete-continuous method of dynamic system modelling. The presented approach is hybrid in its nature, as it combines the advantages of spatial discretization methods with those of continuous system modelling methods. In the proposed method, a three-dimensional system is discretised in two directions only, with the third direction remaining continuous. The thus obtained discrete-continuous model is described...
-
Global sensitivity analysis of membrane model of abdominal wall with surgical mesh
PublicationThe paper addresses the issue of ventral hernia repair. Finite Element simulations can be helpful in the optimization of hernia parameters. A membrane abdominal wall model is proposed in two variants: a healthy one and including hernia defect repaired by implant. The models include many uncertainties, e.g. due to variability of abdominal wall, intraabdominal pressure value etc. Measuring mechanical properties with high accuracy...
-
Diagnostics of pillars in St. Mary’s Church (Gdańsk, Poland) using the GPR method
PublicationThe main goal of this study was non-destructive evaluation of pillars in the St. Mary’s Church (Gdańsk, Poland) using the ground penetrating radar (GPR) technique. The GPR inspection was conducted on four brick masonry pillars and five pillars strengthened by reinforced concrete jacketing. Data were acquired with a 2 GHz antenna along longitudinal and transverse profiles. The study involved the estimation of the electromagnetic...
-
A novel heterogeneous model of concrete for numerical modelling of ground penetrating radar
PublicationThe ground penetrating radar (GPR) method has increasingly been applied in the non-destructive testing of reinforced concrete structures. The most common approach to the modelling of radar waves is to consider concrete as a homogeneous material. This paper proposes a novel, heterogeneous, numerical model of concrete for exhaustive interpretation of GPR data. An algorithm for determining the substitute values of the material constants...
-
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...
-
Accuracy of the Discrete Green's Function Formulation of the FDTD Method
PublicationThis paper reports an evaluation of the accuracy of the discrete Greens function (DGF) formulation of the finite-difference time-domain (FDTD) method. Recently, the closed-form expression for the DGF and its efficient numerical implementation were presented, which facilitates applications of the DGF in FDTD simulations of radiation and scattering problems. So far, the accuracy of the DGF formulation of the FDTD method has been...
-
Hybrid Technique Combining the FDTD Method and Its Convolution Formulation Based on the Discrete Green's Function
PublicationIn this letter, a technique combining the finite-difference time-domain (FDTD) method and its formulation based on the discrete Green's function (DGF) is presented. The hybrid method is applicable to inhomogeneous dielectric structures that are mutually coupled with wire antennas. The method employs the surface equivalence theorem in the discrete domain to separate the problem into a dielectric domain simulated using the FDTD method...
-
Zastosowanie symulacji metodą różnic skończonych w dziedzinie czasu (FDTD) w kształceniu inżynierów w zakresie optyki i elektrodynamiki
PublicationZrozumienie zjawisk związanych z propagacją fal elektromagnetycznych stanowi kluczowy etap kształcenia inżynierów w dziedzinach związanych z optyką, elektroniką oraz telekomunikacją. Oprócz opanowania aparatu matematycznego oraz metod projektowych istotne jest intuicyjne zrozumienie treści przekazywanych podczas kursów optyki i elektrodynamiki. W realizacji tego celu praktyczną pomoc dydaktyczną stanowić mogą wizualizacje i symulacje...
-
Modelling and simulations in time-fractional electrodynamics based on control engineering methods
PublicationIn this paper, control engineering methods are presented with regard to modelling and simulations of signal propagation in time-fractional (TF) electrodynamics. That is, signal propagation is simulated in electromagnetic media described by Maxwell’s equations with fractional-order constitutive relations in the time domain. We demonstrate that such equations in TF electrodynamics can be considered as a continuous-time system of...
-
Application of the discrete Green's function-based antenna simulations for excitation of the total-field/scattered-field interface in the FDTD method
PublicationIn this article, the discrete Green's function formulation of the finite-difference time-domain (DGF-FDTD) method is proposed for simulation of wire antennas irradiating inhomogeneous dielectric scatterers. Surface equivalence theorem in the discrete domain is used to separate the problem into an inhomogeneous domain and a wire antenna that are simulated with the use of FDTD and DGF-FDTD, respectively. Then, the excitation of the...
-
FDTD Method for Electromagnetic Simulations in Media Described by Time-Fractional Constitutive Relations
PublicationIn this paper, the finite-difference time-domain (FDTD) method is derived for electromagnetic simulations in media described by the time-fractional (TF) constitutive relations. TF Maxwell’s equations are derived based on these constitutive relations and the Grünwald–Letnikov definition of a fractional derivative. Then the FDTD algorithm, which includes memory effects and energy dissipation of the considered media, is introduced....
-
Au–Si plasmonic platforms: synthesis, structure and FDTD simulations
PublicationPlasmonic platforms based on Au nanostructures have been successfully synthesized by directional solidification of a eutectic from Au and the substrate. In order to determine homogeneous shape and space distribution, the influence of annealing conditions and the initial thickness of the Au film on the nanostructures was analyzed. For the surface morphology studies, SEM and AFM measurements were performed. The structure of platforms...
-
Parallel Implementation of the Discrete Green's Function Formulation of the FDTD Method on a Multicore Central Processing Unit
PublicationParallel implementation of the discrete Green's function formulation of the finite-difference time-domain (DGF-FDTD) method was developed on a multicore central processing unit. DGF-FDTD avoids computations of the electromagnetic field in free-space cells and does not require domain termination by absorbing boundary conditions. Computed DGF-FDTD solutions are compatible with the FDTD grid enabling the perfect hybridization of FDTD...
-
On Computational Aspects of Greedy Partitioning of Graphs
PublicationIn this paper we consider a problem of graph P-coloring consisting in partitioning the vertex set of a graph such that each of the resulting sets induces a graph in a given additive, hereditary class of graphs P. We focus on partitions generated by the greedy algorithm. In particular, we show that given a graph G and an integer k deciding if the greedy algorithm outputs a P-coloring with a least k colors is NP-complete for an infinite...
-
Numerical Simulations and Tracer Studies as a Tool to Support Water Circulation Modeling in Breeding Reservoirs
PublicationThe article presents a proposal of a method for computer-aided design and analysis of breeding reservoirs in zoos and aquariums. The method applied involves the use of computer simulations of water circulation in breeding pools. A mathematical model of a pool was developed, and a tracer study was carried out. A simplified model of two-dimensional flow in the form of a biharmonic equation for the stream function (converted into...
-
FDTD Simulations on Disjoint Domains with the Use of Discrete Green's Function Diakoptics
PublicationA discrete Green's function (DGF) approach to couple disjoint domains in the finite-difference time-domain (FDTD) grid is developed. In this method, total-field/scattered-field (TFSF) FDTD domains are associated with simulated objects whereas the interaction between them is modeled with the use of the DGF propagator. Hence, source and scatterer are simulated in separate domains and updating of vacuum cells, being of little interest,...