Search results for: numerical solution

Numerical solution of threshold problems in epidemics and population dynamics
PublicationA new algorithm is proposed for the numerical solution of threshold problems in epidemics and population dynamics. These problems are modeled by the delaydifferential equations, where the delay function is unknown and has to be determined from the threshold conditions. The new algorithm is based on embedded pair of continuous Runge–Kutta method of order p = 4 and discrete Runge–Kutta method of order q = 3 which is used for the...

Analytical and numerical solution of a coupled KdV  MKdV system.
PublicationTransformację Darboux zastosowano do całkowania układów równań KdV  MKdV.Reprezentacja Laxa używa 2x2 macierzowe zagadnienie spektralne drugiego rzędu. Numeryczną metodę wprowadzono razem z dowodem zbieżności.

Analytical and numerical solution of a coupled KdVMKdV system.
PublicationTransformacje Darboux zostały użyte do rozwiązania układu równań KdVMKdV.

Numerical solution of boundary value problems with deviated arguments
PublicationZastosowano metodę różnicową, aby wyznaczyć przybliżone rozwiązanie problemubrzegowego z odchylonymi argumentami. Pokazano, że metoda różnicowa, przyodpowiednich warunkach, jest zbieżna do rozwiązania i podano oszacowaniabłędów.

A COMPUTATIONAL ALGORITHM FOR THE NUMERICAL SOLUTION OF NONLINEAR FRACTIONAL INTEGRAL EQUATIONS
Publication 
Numerical solution of fractional neutron point kinetics in nuclear reactor
PublicationThis paper presents results concerning solutions of the fractional neutron point kinetics model for a nuclear reactor. Proposed model consists of a bilinear system of fractional and ordinary differential equations. Three methods to solve the model are presented and compared. The first one entails application of discrete GrünwaldLetnikov definition of the fractional derivative in the model. Second involves building an analog scheme...

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 initialboundary 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 SturmLiouville problem. The horizontal profile is defined by a coupled KdV system which is...

SpaceTime 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...

Vibration of the bridge under moving singular loads  theoretical formulation and numerical solution
PublicationThe paper presents the results of the numerical analysis of a simple vehicle passing over a simply supported bridge span. The bridge is modelled by a EulerBernoulli beam. The vehicle is modelled as a linear, viscoelastic oscillator, moving at a constant speed. The system is described by a set of differential equations of motion and solved numerically using the RungeKutta algorithm. The results are compared with the solution...

On Nonlinear Bending Study of a PiezoFlexomagnetic Nanobeam Based on an AnalyticalNumerical Solution
PublicationAmong various magnetoelastic phenomena, flexomagnetic (FM) coupling can be defined as a dependence between strain gradient and magnetic polarization and, contrariwise, elastic strain and magnetic field gradient. This feature is a higherorder one than piezomagnetic, which is the magnetic response to strain. At the nanoscale, where large strain gradients are expected, the FM effect is significant and could be even dominant. In...

Numerical solution analysis of fractional point kinetics and heat exchange in nuclear reactor
PublicationThe paper presents the neutron point kinetics and heat exchange models for the nuclear reactor. The models consist of a nonlinear system of fractional ordinary differential and algebraic equations. Two numerical algorithms are used to solve them. The first algorithm is application of discrete GrünwaldLetnikov definition of the fractional derivative in the model. The second involves building an analog scheme in the FOMCON Toolbox...

Fractional neutron point kinetics equations for nuclear reactor dynamics – Numerical solution investigations
PublicationThis paper presents results concerning numerical solutions to a fractional neutron point kinetics model for a nuclear reactor. The paper discusses and expands on results presented in (EspinosaParedes et al., 2011). The fractional neutron point kinetics model with six groups of delayed neutron precursors was developed and a numerical solution using the Edwards’ method was proposed (Edwards et al., 2002). The mathematical model...

Numerical Solution of the TwoDimensional Richards Equation Using Alternate Splitting Methods for Dimensional Decomposition
PublicationResearch on seepage flow in the vadose zone has largely been driven by engineering and environmental problems affecting many fields of geotechnics, hydrology, and agricultural science. Mathematical modeling of the subsurface flow under unsaturated conditions is an essential part of water resource management and planning. In order to determine such subsurface flow, the twodimensional (2D) Richards equation can be used. However,...

A new approach to numerical solution of fixedpoint problems and its application to delay differential equations
PublicationW pracy rozpatruje się pewne aproksymacje punktu stałego ciągłego operatora A odwzorowującego przestrzeń metryczną w siebie. Wspomniany punkt stały przybliża się tzw. epsilon przybliżonym punktem stałym z przestrzeni skończenie wymiarowej. Udowodnione zostało twierdzenie dające warunki konieczne i dostateczne istnienia punktu stałego w ogólnej przestrzeni metrycznej. Warunki te wyrażone są w terminach epsilon przybliżonego punktu...

Impact of Energy Slope Averaging Methods on Numerical Solution of 1D Steady Gradually Varied Flow
PublicationIn this paper, energy slope averaging in the onedimensional steady gradually varied flow model is considered. For this purpose, different methods of averaging the energy slope between crosssections are used. The most popular are arithmetic, geometric, harmonic and hydraulic means. However, from the formal viewpoint, the application of different averaging formulas results in different numerical integration formulas. This study...

Balance error generated by numerical diffusion in the solution of Muskingum equation
PublicationIn the paper the conservative properties of the lumped hydrological models with variable parameters are discussed. It is shown that in the case of the nonlinear Muskingum equation the mass balance is not satisfied. The study indicates that the mass balance errors are caused by the improper form of equation and by the numerical diffusion which is generated in the solution. It has been shown that the classical way of derivation...

Balance errors generated by numerical diffusion in the solution of nonlinear open channel flow equations
PublicationThe paper concerns the untypical aspect of application of the dissipative numerical methods to solve nonlinear hyperbolic partial differential equations used in open channel hydraulics. It is shown that in some cases the numerical diffusion generated by the applied method of solution produces not only inaccurate solution but as well as a balance error. This error may occur even for an equation written in the conservative form not...

Numerical solutions for blood flow in elastic vessels
PublicationWe consider the differential–algebraic system for the blood flow and pressure in the systemic arteries. By the operator splitting method, we transform the system into the hyperbolic one, introduce the bicharacteristics, and perform the time–space nonuniform discretization, obtaining the innovative difference scheme. Our results are illustrated with numerical experiments.

Balance errors in numerical solutions of shallow water equations
PublicationThe analysis of the conservative properties of the shallow water equations is presented in the paper. The work focuses on the consistency of numerical solution of these equations with the conservation laws of mass and momentum. The investigations involve two different conservative forms which are solved by an implicit box scheme. The theoretical analysis supported by numerical experiments is carried out for rectangular channel...

Numerical solutions for large deformation problems in geotechnical engineering
PublicationThe problem of large deformations often occurs in geotechnical engineering. Numerical modeling of such issues is usually complex and tricky. The chosen solution has to implicate soilsoil and soilstructure interactions. In this paper, a review of the most popular numerical methods for large deformation problems is presented. The Coupled EulerianLagrangian (CEL) method, the Arbitrary LagrangianEulerian (ALE) method, the Smoothed...

Multimode systems of nonlinear equations: derivation, integrability, and numerical solutions
PublicationWe consider the propagation of electromagnetic pulses in isotropic media taking a thirdorder nonlinearityinto account. We develop a method for transforming Maxwell's equations based on a complete set ofprojection operators corresponding to wavedispersion branches (in a waveguide or in matter) with thepropagation direction taken into account. The most important result of applying the method is a systemof equations describing the...

Application of nonclassical operational calculus to indicate hazards in numerical solutions of engineering problems
PublicationThe article addresses the application of non classical operational calculus to approximative solutions of engineering problems. The engineeringsound examples show that a continuous–discrete problem transformation from differential unequivocal problem to a differential wildcard problem, triggering a change in solution quality. A number of approximative methods are capable to alter both quantitative and qualitative...

Solution of the dikebreak problem using finite volume method and splitting technique
PublicationIn the paper the finite volume method (FVM) is presented for the solution of twodimensional shallow water equations. These equations are frequently used to simulate the dambreak and dikebreak induced flows. The applied numerical algorithm of FVM is based on the wavepropagation algorithm which ensures a stable solution and simultaneously minimizes the numerical errors. The dimensional decomposition according to the coordinate...

Impact of the Finite Element Mesh Structure on the Solution Accuracy of a TwoDimensional Kinematic Wave Equation
PublicationThe paper presents the influence of the finite element mesh structure on the accuracy of the numerical solution of a twodimensional linear kinematic wave equation. This equation was solved using a twolevel scheme for time integration and a modified finite element method with triangular elements for space discretization. The accuracy analysis of the applied scheme was performed using a modified equation method for three different...

ADAPTIVE METHOD FOR THE SOLUTION OF 1D AND 2D ADVECTIONDIFFUSION EQUATIONS USED IN ENVIRONMENTAL ENGINEERING
PublicationThe paper concerns the numerical solution of onedimensional (1D) and twodimensional (2D) advectiondiffusion equations. For the numerical solution of the 1D advectiondiffusion equation a method, originally proposed for solution of the 1D pure advection equation, has been developed. A modified equation analysis carried out for the proposed method allowed increasing of the resulting solution accuracy and consequently, to reduce...

A NUMERICAL STUDY ON THE DYNAMICS OF DENGUE DISEASE MODEL WITH FRACTIONAL PIECEWISE DERIVATIVE
PublicationThe aim of this paper is to study the dynamics of Dengue disease model using a novel piecewise derivative approach in the sense of singular and nonsingular kernels. The singular kernel operator is in the sense of Caputo, whereas the nonsingular kernel operator is the Atangana–Baleanu Caputo operator. The existence and uniqueness of a solution with piecewise derivative are examined for the aforementioned problem. The suggested...

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 wellknown 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...

Numerical and Analytical Investigation of Aluminium Bracket Strengthening
PublicationThis paper focuses on an analytical and numerical investigation of aluminium brackets used to fasten lightweight curtain walls to building facilities. The authors propose a solution to increase the load capacity of aluminium brackets by means of additional cover plates (straps). This paper also includes a short survey of literature and material properties concerning the EN AW6060 T6 aluminium alloy. This paper suggests an initiation...

Solution of coupled integral equations for quantum scattering in the presence of complex potentials
PublicationIn this paper, we present a method to compute solutions of coupled integral equations for quantum scattering problems in the presence of a complex potential. We show how the elastic and absorption cross sections can be obtained from the numerical solution of these equations in the asymptotic region at large radial distances.

FEMbased wave propagation modelling for SHM: Certain numerical issues in 1D structures
PublicationThe numerical modelling of structural elements is an important aspect of modern diagnostic systems. However, the process of numerical implementation requires advanced levels of consideration of multiple aspects. Important issues of that process are the positive and negative aspects of solution applied methods. Therefore the aim of this article is to familiarise the reader with the most important aspects related to the process of...

Identification of Parameters Influencing the Accuracy of the Solution of the Nonlinear Muskingum Equation
PublicationTwo nonlinear versions of the Muskingum equation are considered. The difference between both equations relates to the exponent parameter. In the first version, commonly used in hydrology, this parameter is considered as free, while in the second version, it takes a value resulting from the kinematic wave theory. Consequently, the first version of the equation is dimensionally inconsistent, whereas the proposed second one is consistent. It...

Numerical and experimental investigation of rotational stiffness of zedpurlins connection with sandwich panels
PublicationA rotational resistant stiffness of the zedpurlins connection with sandwich panels is investigated. A simple finite element method model of the connection is proposed. The numerical analysis of the model performed by ABAQUS software, in physically linear and geometrically nonlinear ranges, leads to the rotational resistant stiffness sought. The numerical solution obtained is verified experimentally. Two variants of distribution...

On the necessity of experimental verification of numerical results in biomedical applications
PublicationPorous structures made of metal or biopolymers with a structure similar in shape and mechanical properties to human bone can be easily produced by stereolitography techniques, e.g. selective laser melting (SLM). Numerical techniques, like finite element method (FEM) have great potential in testing new, even the most sophisticated designs, according to their mechanical properties, i.e. strength or stiffness. However, due to different...

Discussion of “Development of an Accurate Time integration Technique for the Assessment of QBased versus hBased Formulations of the Diffusion Wave Equation for Flow Routing” by K. Hasanvand, M.R. Hashemi and M.J. Abedini
PublicationThe discusser read the original with great interest. It seems, however, that some aspects of the original paper need additional comments. The authors of the original paper discuss the accuracy of a numerical solution of the diffusion wave equation formulated with respect to different state variables. The analysis focuses on nonlinear equations in the form of a single transport equation with the discharge Q (volumetric flow rate)...

Impact of diffusion coefficient averaging on solution accuracy of the 2D nonlinear diffusive wave equation for floodplain inundation
PublicationIn the study, the averaging technique of diffusion coefficients in the twodimensional nonlinear diffusive wave equation applied to the floodplain inundation is presented. As a method of solution, the splitting technique and the modified finite element method with linear shape functions are used. On the stage of spatial integration, it is often assumed that diffusion coefficient is constant over element and equal to its average...

Modelling of FloodWave Propagation with Wetdry Front by Onedimensional Diffusive Wave Equation
PublicationA full dynamic model in the form of the shallow water equations (SWE) is often useful for reproducing the unsteady flow in open channels, as well as over a floodplain. However, most of the numerical algorithms applied to the solution of the SWE fail when flood wave propagation over an initially dry area is simulated. The main problems are related to the very small or negative values of water depths occurring in the vicinity of...

Numerical Characterization of Thresholds for the Focusing 1d Nonlinear Schrödinger Equation
PublicationThe focusing nonlinear Schrödinger equation arises in various physical phenomena and it is therefore of interest to determine mathematical conditions on the initial data that guarantee whether the corresponding solution will blow up in finite time or exist globally in time. We focus on solutions to the mass‐supercritical nonlinear Schrödinger equation (1) in 1D case. In particular, we investigate numerical thresholds between blow...

Analytical SteadyState Model of the Pipeline Flow Process
PublicationThe paper addresses the issue of modeling the flow process in transmission pipelines. A base model used for numerical simulation is introduced. Under certain assumptions concerning steady state analysis, the differential equations describing the process are solved analytically for two cases: zero and nonzero inclination angle α. These equations describe a constant flow rate and a corresponding distribution of the pressure along...

Numerical Modelling of Forced Convection of Nanofluids in Smooth, Round Tubes: A Review
PublicationA comprehensive review of published works dealing with numerical modelling of forced convection heat transfer and hydrodynamics of nanofluids is presented. Due to the extensive literature, the review is limited to straight, smooth, circular tubes, as this is the basic geometry in shellandtube exchangers. Works on numerical modelling of forced convection in tubes are presented chronologically in the first part of the article....

Assessment of tensile forces in Sopot Forest Opera membrane by in situ measurements and iterative numerical strategy for inverse problem
PublicationAssessment of tensile forces in newly build Sopot Forest Opera roofing membrane is presented. The procedure is based on in situ measurements and solution of the inverse problem by an iterative procedure. The goal of the analysis is to determine whether the stress state of the membrane is consistent with the design assumptions. The paper contains the description of measurements, used instruments, applied loadings, numerical investigations...

Numerical Modeling of Cone Penetration Test in Slightly Overconsolidated Clay with Arbitrary LagrangianEulerian Formulation
PublicationIn this paper the results of the cone penetration test (CPT) modeling with the arbitrary LagrangianEulerian (ALE) formulation provided by Abaqus software package have been presented. The study compares the cone resistance and sleeve friction obtained in numerical analysis with values measured in soundings performed in the uniform layer of clayey soil in the Koszalin area. The clay layer was found to be slightly overconsolidated...

LaxWendroff 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 LaxWendroff 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 SaintVenant equations describing free...

Operational Enhancement of Numerical Weather Prediction with Data from Realtime Satellite Images
PublicationNumerical weather prediction (NWP) is a rapidly expanding field of science, which is related to meteorology, remote sensing and computer science. Authors present methods of enhancing WRF EMS (Weather Research and Forecast Environmental Modeling System) weather prediction system using data from satellites equipped with AMSU sensor (Advanced Microwave Sounding Unit). The data is acquired with Department of Geoinformatics’ ground...

Numerical Study of Concrete Mesostructure Effect on Lamb Wave Propagation
PublicationThe article presents the results of the numerical investigation of Lamb wave propagation in concrete plates while taking into account the complex concrete mesostructure. Several concrete models with randomly distributed aggregates were generated with the use of the Monte Carlo method. The influence of aggregate ratio and particle size on dispersion curves representing Lamb wave modes was analyzed. The results obtained for heterogeneous...

Geoengineering computer simulation seems attractive but is it the real world?
PublicationCorrect formulation of the differential equation system for equilibriom conditions of subsoil, especially in terms of controlled numerical calculation, is discussed. The problem of solution stability is also considered. The solution of problems, which are illposed, have no practical value in the majority of cases and is this way the engineering prognosis can lead to real disaster. The object of this paper is quite relevant if...

Numerical and quantitative analysis of HIV/AIDS model with modified AtanganaBaleanu in Caputo sense derivative
PublicationFractional calculus plays an important role in the development of control strategies, the study of the dynamical transmission of diseases, and some other reallife problems nowadays. The timefractional HIV/AIDS model is examined using a novel method in this paper. Based on the Atanganaconcept Baleanu’s of a derivative in the Caputo sense, the current modified fractional derivative operator uses singular and nonlocal kernels....

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 computeraided 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 twodimensional flow in the form of a biharmonic equation for the stream function (converted into...

Experimental and numerical study on polymer element used for reduction of temporary steel grandstand vibrations
PublicationThe purpose of the paper is to analyse the effectiveness of a polymer damper in reduction of a temporary steel grandstand vibrations through the experimental and numerical study. The damper considered in the investigation is constructed out of two Lshape steel members bonded with polymer mass of high damping properties. The element has been installed as a diagonal one at the back part of the structure. The method has been compared...

Verification of algorithms determining wave loads on support structure of wind turbine
PublicationThe offshore wind turbines require determination of wave loads on their support structure. This structure is fixed and, therefore, this problem is reduced to solving only the diffraction problem, which is determined by Laplace equation and conditions on the following boundaries: on the support structure, on the sea free surface and on its bottom, and at infinity on free surface. The linear problem was applied to determine the wave...

Comparison of AHP and Numerical Taxonomy Methods Based on Biogas Plant Location Analysis
PublicationThe paper presents a comparison of the multicriteria Analytic Hierarchy Process (AHP) method and numerical taxonomy in biogas plant location selection. Biogas plants are sources that will significantly contribute to the implementation of the provisions of the energy and climate package for Poland by 2030. Increasing the share of energy produced from renewable sources, e.g. biogas plants, will increase the country’s energy security....