Search
Search results for: FINITE DIFFERENCE
-
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...
-
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...
-
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...
-
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;...
-
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...
-
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...
-
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...
-
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...
-
Mixed, quantum-classical description of electron density transfer in the collision process
PublicationIn this work, we investigate an ion-atom model describing the time-dependent evolution of electron density during the collision. For a S3+- H system, numerical simulations are based on classical trajectory calculations, and the electron density behaviour is described with the time-dependent Schrödinger equation. We apply the finite difference method to obtain quantitative insights into the charge transfer dynamics, providing detailed...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
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...
-
Crank–Nicolson FDTD Method in Media Described by Time-Fractional Constitutive Relations
PublicationIn this contribution, we present the Crank-Nicolson finite-difference time-domain (CN-FDTD) method, implemented for simulations of wave propagation in media described by time-fractional (TF) constitutive relations. That is, the considered constitutive relations involve fractional-order (FO) derivatives based on the Grünwald-Letnikov definition, allowing for description of hereditary properties and memory effects of media and processes....
-
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...
-
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...
-
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,...