Search
Filters
total: 233
filtered: 192
-
Catalog
Chosen catalog filters
Search results for: BIOHEAT EQUATION, IMPLICIT NUMERICAL SCHEME
-
Thermal ablation modeling via the bioheat equation and its numerical treatment
PublicationThe phenomenon of thermal ablation is described by Pennes’ bioheat equation. This model is based on Newton’s law of cooling. Many approximate methods have been considered because of the importance of this issue. We propose an implicit numerical scheme which has better stability properties than other approaches.
-
Thermal ablation modeling via bioheat equation
PublicationWe consider Pennes’ bioheat equation and discuss an implicit numerical scheme which has better stability properties than other approaches. Our discussion concerns Carthesian geometry problems, however it carries over to spherical geometry models and more complicated shapes.
-
Computationally Effcient Solution of a 2D Diffusive Wave Equation Used for Flood Inundation Problems
PublicationThis paper presents a study dealing with increasing the computational efficiency in modeling floodplain inundation using a two-dimensional diffusive wave equation. To this end, the domain decomposition technique was used. The resulting one-dimensional diffusion equations were approximated in space with the modified finite element scheme, whereas time integration was carried out using the implicit two-level scheme. The proposed...
-
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...
-
Balance errors generated by numerical diffusion in the solution of non-linear 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...
-
PROPERTIES OF ONE DIMENSIONAL OPEN-CHANNEL STEADY FLOW EQUATIONS
PublicationIn this paper properties of discrete forms of one dimensional steady gradually varied flow equations are discussed. Such forms of flow equations are obtained as a result of approximation of their differential forms, which is required to solve them numerically. For such purpose explicit or implicit numerical approximation schemes for ordinary differential equations can be applied. It turns out that dependently on the chosen approximation...
-
Coupled nonlinear Schrödinger equations in optic fibers theory
PublicationIn this paper a detailed derivation and numerical solutions of CoupledNonlinear Schr¨odinger Equations for pulses of polarized electromagnetic wavesin cylindrical fibers has been reviewed. Our recent work has been compared withsome previous ones and the advantage of our new approach over other methods hasbeen assessed. The novelty of our approach lies is an attempt to proceed withoutloss of information within the frame of basic...
-
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 non-linear 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...
-
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...
-
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 Solution of the Two-Dimensional 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 two-dimensional (2D) Richards equation can be used. However,...
-
On convergence and stability of a numerical scheme of Coupled Nonlinear Schrödinger Equations
PublicationRozważamy rozwiązania numeryczne układu sprężynowych równań nieliniowych Schrödingera. Udowodniliśmy stabilność i zbieżność. Testujemy za pomocą rozwiązań solitonowych.
-
Asymptotic numerical solver for the linear Klein–Gordon equation with space- and time-dependent mass
Publication -
Numerical simulation of transient flow in storm sewers using standard and improved mccormack scheme
PublicationProblem przepływu wody w rurach kanalizacji deszczowej wiąże się z pewnymi specyficznymi zjawiskami zachodzącymi w przewodach podczas zjawisk burzowych. Jeśli rury zaczynają być całkowicie wypełnione wodą, można zaobserwować przejście z ruchu ze swobodną powierzchnią do ruchu pod ciśnieniem i odwrotnie. Takie zjawisko można zaobserwować również w kanałach kontrolowanych przez urządzenia sterujące, np. zasuwy. Ponadto ruch szybkozmienny...
-
Analytical-numerical approach to solve the transport equation for steady gradually varied flow in open channel
PublicationW pracy przedstawiono metodę rozwiązania równania transportu adwekcyjno-dyfuzyjnego w przypadku ustalonego niejednostajnego przepływu w kanałach otwartych. Metoda wykorzystuje technikę dekompozycji. Do rozwiązania równania adwekcji-dyfuzji zastosowano analityczne rozwiązanie w postaci odpowiedzi impulsowej liniowego równania adwekcji-dyfuzji. Dokonano adaptacji metody dla przypadku ze zmiennymi parametrami. Do rozwiązania drugiej...
-
Impact of the Finite Element Mesh Structure on the Solution Accuracy of a Two-Dimensional Kinematic Wave Equation
PublicationThe paper presents the influence of the finite element mesh structure on the accuracy of the numerical solution of a two-dimensional linear kinematic wave equation. This equation was solved using a two-level 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...
-
Elastoplastic material law in 6-parameter nonlinear shell theory
PublicationWe develop the elastoplastic constitutive relations for nonlinear exact 6-parameter shell theory. A J2-type theory with strain hardening is formulated that takes into account asymmetric membrane strain measures. The incremental equations are solved using implicit Euler scheme with closest point projection algorithm. The presented test example shows the correctness of the proposed approach. Influence of micropolar material parameters...
-
Inverse Flood Routing Using Simplified Flow Equations
PublicationThe paper considers the problem of inverse flood routing in reservoir operation strategy. The aim of the work is to investigate the possibility of determining the hydrograph at the upstream end based on the hydrograph required at the downstream end using simplified open channel flow models. To accomplish this, the linear kinematic wave equation, the diffusive wave equation and the linear Muskingum equation are considered. To achieve...
-
An isogeometric finite element formulation for boundary and shell viscoelasticity based on a multiplicative surface deformation split
PublicationThis work presents a numerical formulation to model isotropic viscoelastic material behavior for membranes and thin shells. The surface and the shell theory are formulated within a curvilinear coordinate system,which allows the representation of general surfaces and deformations. The kinematics follow from Kirchhoff–Love theory and the discretization makes use of isogeometric shape functions. A multiplicative split of the surface...
-
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...
-
A robust design of a numerically demanding compact rat-race coupler
PublicationA fast and accurate design procedure of a computationally expensive microwave circuit has been presented step-by-step and experimentally validated on the basis of a compact rat-race coupler (RRC) comprising slow-wave resonant structures (SWRSs). The final compact RRC solution has been obtained by means of a sequential optimization scheme exploiting the implicit space mapping (ISM) algorithm. A well-suited surrogate optimization...
-
Straightened characteristics of McKendrick-von Foerster equation
PublicationWe study the McKendrick-von Foerster equation with renewal (that is the age-structured model, with total population dependent coefficient and nonlinearity). By using a change of variables, the model is then transformed to a standard age-structured model in which the total population dependent coefficient of the transport term reduces to a constant 1. We use this transformation to get existence, uniqueness of solutions of the problem...
-
Alternative Approach to Convolution Term of Viscoelasticity in Equations of Unsteady Pipe Flow
PublicationIn the paper the selected aspects concerning description of viscoelastic behavior of pipe walls during unsteady flow are analyzed. The alternative convolution expression of the viscoelastic term is presented and compared with the corresponding term referring to unsteady friction. Both approaches indicate similarities in the forms of impulse response functions and the parameter properties. The flow memory was introduced into convolution...
-
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...
-
Electronically Excited States in Solution via a Smooth Dielectric Model Combined with Equation-of-Motion Coupled Cluster Theory
PublicationWe present a method for computing excitation energies for molecules in solvent, based on the combination of a minimal parameter implicit solvent model and the equation-of-motion coupled-cluster singles and doubles method (EOM-CCSD). In this method, the solvent medium is represented by a smoothly varying dielectric function, constructed directly from the quantum mechanical electronic density using only two tunable parameters. The...
-
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...
-
Discrete and continuous fractional persistence problems – the positivity property and applications
PublicationIn this article, we study the continuous and discrete fractional persistence problem which looks for the persistence of properties of a given classical (α=1) differential equation in the fractional case (here using fractional Caputo’s derivatives) and the numerical scheme which are associated (here with discrete Grünwald–Letnikov derivatives). Our main concerns are positivity, order preserving ,equilibrium points and stability...
-
Some new soliton solutions to the higher dimensional Burger–Huxley and Shallow water waves equation with couple of integration architectonic
PublicationIn this paper, we retrieve some traveling wave, periodic solutions, bell shaped, rational, kink and anti-kink type and Jacobi elliptic functions of Burger’s equation and Shallow water wave equation with the aid of various integration schemes like improved -expansion scheme and Jacobi elliptic function method respectively. We also present our solutions graphically in various dimensions.
-
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...
-
Minimal parameter implicit solvent model for ab initioelectronic-structure calculations
PublicationAbstract - We present an implicit solvent model for ab initio electronic-structure calculations which is fully self-consistent and is based on direct solution of the nonhomogeneous Poisson equation. The solute cavity is naturally defined in terms of an isosurface of the electronic density according to the formula of Fattebert and Gygi (J. Comput. Chem., 23 (2002) 662). While this model depends on only two parameters, we demonstrate...
-
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...
-
Discussion of “Development of an Accurate Time integration Technique for the Assessment of Q-Based versus h-Based 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)...
-
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...
-
Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
PublicationIn this paper, the conservative properties of the Muskingum equation, commonly applied to solve river flood routing, are analysed. The aim of this analysis is to explain the causes ofthe mass balance error, which is observed in the numerical solutions of its non-linear form. The linear Muskingum model has been considered as a semi-discrete form of the kinematic wave equation and therefore it was possible to derive its two non-linear...
-
FPGA Acceleration of Matrix-Assembly Phase of RWG-Based MoM
PublicationIn this letter, the field-programmable-gate-array accelerated implementation of matrix-assembly phase of the method of moments (MoM) is presented. The solution is based on a discretization of the frequency-domain mixed potential integral equation using the Rao-Wilton-Glisson basis functions and their extension to wire-to-surface junctions. To take advantage of the given hardware resources (i.e., Xilinx Alveo U200 accelerator card),...
-
Stability and limit load analysis of a cold-formed channel section column
PublicationThe paper presents stability and limit load analysis of a steel column 1440 mm high of a cold-formed channel section, subjected to a combination of compression and bending. Experimental results were compared to the resistance of a code procedure and to the outcome of numerical non-linear analysis. Comparison was made of numerical solutions by means of static (Riks) and dynamic (Explicit and Implicit) methods. Perfectly elastic...
-
Numerical simulation of cold flow and combustion in a swirl stabilized combustor
PublicationA numerical simulation model was developed to investigate the cold flow and combustion using Ansys FLUENT 2021R1. The governing equations were solved using the pressure-based method, and pressure–velocity coupling was performed using the SIMPLE method. To model the turbulent process, the RSM model was used. Non-premixed model is chosen to solve the chemical kinetics between fuel and oxigen. Radiation heat transfer was calculated...
-
ADAPTIVE METHOD FOR THE SOLUTION OF 1D AND 2D ADVECTION-DIFFUSION EQUATIONS USED IN ENVIRONMENTAL ENGINEERING
PublicationThe paper concerns the numerical solution of one-dimensional (1D) and two-dimensional (2D) advection-diffusion equations. For the numerical solution of the 1D advection-diffusion 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...
-
MODELLING OF TRANSIENT FLOW IN STORM SEWERS
PublicationThe paper focuses on the assessment of second-order explicit numerical scheme for unsteady flows in sewers. In order to simulate the pressurized flow the 'Preissmann slot' concept is implemented. For simulation of the transcritical flow the original and improved McCormack scheme is used. The calculated results are compared with numerical solutions and laboratory measurements published in the technical literature. Moreover, the...
-
Underfrequency Load Shedding: An Innovative Algorithm Based on Fuzzy Logic
PublicationIn contemporary power systems, the load shedding schemes are typically based on disconnecting a pre-specified amount of load after the frequency drops below a predetermined value. The actual conditions at the time of disturbance may largely dier from the assumptions, which can lead to non-optimal or ineective operation of the load shedding scheme. For many years, increasing the eectiveness of the underfrequency load shedding (UFLS)...
-
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...
-
Comparison of Average Energy Slope Estimation Formulas for One-dimensional Steady Gradually Varied Flow
PublicationTo find the steady flow water surface profile, it is possible to use Bernoulli’s equation, which is a discrete form of the differential energy equation. Such an approach requires the average energy slope between cross-sections to be estimated. In the literature, many methods are proposed for estimating the average energy slope in this case, such as the arithmetic mean, resulting in the standard step method, the harmonic mean and...
-
Finite element simulation of cross shaped window panel supports
PublicationThe aim of the work is to verify suitability of cross-shaped window panel supports for mullion-transom wall systems. The Finite Element Method (FEM) is chosen to determine the behaviour of stainless steel elements under loading. The advanced non-linear numerical simulations are carried out using an implicit FEM software package MSC.Marc. This study is proposed to initiate the comprehensive investigation of mechanical properties...
-
A pore-scale thermo–hydro-mechanical model for particulate systems
PublicationA pore scale numerical method dedicated to the simulation of heat transfer and associated thermo–hydro-mechanical couplings in granular media is described. The proposed thermo–hydro-mechanical approach builds on an existing hydromechanical model that employs the discrete element method for simulating the mechanical behavior of dense sphere packings and combines it with the finite volume method for simulating pore space fluid flow...
-
Hyperelastic Microcantilever AFM: Efficient Detection Mechanism Based on Principal Parametric Resonance
PublicationThe impetus of writing this paper is to propose an efficient detection mechanism to scan the surface profile of a micro-sample using cantilever-based atomic force microscopy (AFM), operating in non-contact mode. In order to implement this scheme, the principal parametric resonance characteristics of the resonator are employed, benefiting from the bifurcation-based sensing mechanism. It is assumed that the microcantilever is made...
-
On possible applications of media described by fractional-order models in electromagnetic cloaking
PublicationThe purpose of this paper is to open a scientific discussion on possible applications of media described by fractional-order (FO) models (FOMs) in electromagnetic cloaking. A 2-D cloak based on active sources and the surface equivalence theorem is simulated. It employs a medium described by FOM in communication with sources cancelling the scattered field. A perfect electromagnetic active cloak is thereby demonstrated with the use...
-
Dimensionally Consistent Nonlinear Muskingum Equation
PublicationAlthough the Muskingum equation was proposed nearly 75 years ago, it is still a subject of active research. Despite of its simple form, the real properties of this equation have not been comprehensively explained. This paper proposes a new interpretation of the linear McCarthy’s relation. This relation can be interpreted only together with the storage equation, whereas the Muskingum equation can be derived directly from the system...
-
Numerical modelling and experimental verification of compressible squeeze film pressure
PublicationThe validity of using the Reynolds equation for compressible squeeze film pressure was tested with computational fluid dynamics (CFD). A squeeze film air bearing was instrumented with pressure sensors and non-contacting displacement probes to provide transient measurements of film thickness and pressure. The film thickness measurements also provided input parameters to the numerical prediction. However, numerical results showed...
-
Computational issues of solving the 1D steady gradually varied flow equation
PublicationIn this paper a problem of multiple solutions of steady gradually varied flow equation in the form of the ordinary differential energy equation is discussed from the viewpoint of its numerical solution. Using the Lipschitz theorem dealing with the uniqueness of solution of an initial value problem for the ordinary differential equation it was shown that the steady gradually varied flow equation can have more than one solution....
-
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...