Search
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...
-
Database of the thermal ablation model
Open Research DataThermal ablation is a low invasive technique which eliminates cancerous tissue using high temperature. The presented database was used to show the temperature distribution for t=600[s] in two cases: when the value of the thermal conductivity of tissue k(x;T) is constant and for the variable k(x;T). In addition, using these data we showed the difference...
-
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...
-
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...
-
Wioletta Gorczewska-Langner dr inż.
People -
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...
-
Karolina Lademann mgr
PeopleCurriculum vitae
-
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...
-
Michał Michna dr hab. inż.
PeopleMichal Michna received the M.Sc. and Ph.D. degrees in electrical engineering from the Gdansk University of Technology (GUT), Gdansk, Poland, in 1998 and 2005, respectively. Since 2004, he was employed at the Department of Power Electronics and Electrical Machines of the Gdańsk University of Technology (assistant, assistant professor, senior lecturer). In 2010-2015 he was a deputy of head of the Department of Power Electronics and...
-
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...
-
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...
-
Karolina Lademann Mgr
People -
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...
-
Dataset of phase portraits of the fractional prey-predator model with Holling type-II interaction (without predator harvesting)
Open Research DataThe need for a fractional generalization of a given classical model is often due to new behaviors which cannot be taken into account by the model. In this situation, it can be useful to look for a fractional deformation of the initial system, trying to fit the fractional exponent of differentiation in order to catch properly the data.
-
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)...
-
Marek Czachor prof. dr hab.
People -
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...