Search results for: BIOHEAT EQUATION, IMPLICIT NUMERICAL SCHEME
-
A graphical approach to yield and boundary surfaces of selected hypoplastic constitutive equations
PublicationThe article describes how to identify the boundary and yield surface for hypoplastic constitutive equations proposed by Wu, Gudehus and Bauer. It is shown how to identify and plot the surfaces for any equation in this class. Calculation errors are analyzed characteristic for appleid set of numerical formulas. In the paper there are computer links to the source code prepared in the MATLAB system, based on istructions in the article....
-
A New Approach to Stability Evaluation of Digital Filters
PublicationIn this paper, a new numerical method of evaluating digital filter stability is presented. This approach is based on novel root-finding algorithms at the complex plane using the Delaunay triangulation and Cauchy's Argument Principle. The presented algorithm locates unstable zeros of the characteristic equation with their multiplicities. The proposed method is generic and can be applied to a vast range of systems. Verification of...
-
On local buckling of cold-formed channel members
PublicationThe paper deals with local buckling of the compressed flanges of cold-formed thin-walled channel beams subjected to pure bending or axially compressed columns. Arbitrarily shaped flanges of open cross-sections and the web-flange interactions are taken into account. Buckling deformation of a beam flange is described by displacement related to torsion of the flange about the line of its connection with the web. Total potential energy...
-
Long-term hindcast simulation of sea level in the Baltic Sea
Open Research DataThe dataset contains the results of numerical modelling of sea level fluctuations over a period of 50 years (1958-2007) in the Baltic Sea. A long-term hindcast simulation was performed using a three-dimensional hydrodynamic model PM3D (Kowalewski and Kowalewska-Kalkowska, 2017), a new version of the M3D model (Kowalewski, 1997). The hydrodynamic model...
-
An extended basis inexact shift–invert Lanczos for the efficient solution of large-scale generalized eigenproblems
PublicationThis paper proposes a technique, based on the Inexact Shift–Invert Lanczos (ISIL) method with Inexact Jacobi Orthogonal Component Correction (IJOCC) refinement, and a preconditioned conjugate-gradient (PCG) linear solver with multilevel preconditioner, for finding several eigenvalues for generalized symmetric eigenproblems. Several eigenvalues are found by constructing (with the ISIL process) an extended projection basis. Presented...
-
Modelling of FloodWave Propagation with Wet-dry Front by One-dimensional 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 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ünwald-Letnikov definition of the fractional derivative in the model. The second involves building an analog scheme in the FOMCON Toolbox...
-
Model Order Reduction for Problems With Dispersive Surface Boundary Conditions
PublicationThis letter proposes a new scheme for reduced-order finite-element modeling of electromagnetic structures with nonlinear, dispersive surface boundary conditions, which optimally exploits the numerically stable and efficient MOR framework for second-order systems provided by SAPOR method. The presented results of numerical experiments for an example of a waveguide filter demonstrate the superior accuracy of the resulting reduced models...
-
Simulating propagation of coherent light in random media using the Fredholm type integral equation
PublicationStudying propagation of light in random scattering materials is important for both basic and applied research. Such studies often require usage of numerical method for simulating behavior of light beams in random media. However, if such simulations require consideration of coherence properties of light, they may become a complex numerical problems. There are well established methods for simulating multiple scattering of light (e.g....
-
Parametric method applicable in calculating breakout force and time for lifting axisymmetric objects from seabed.
PublicationThe issue of calculating parameters for lifting objects of axisymmetric shapes from the seabed is studied. The article presents the results of numerical simulation examination of the equation formulated by Foda for the force needed to extricate the object from the seabed depending on soil and water parameters, extrication speed, and object diameter. The simulations were performed within the range of parameters characteristic for...
-
Local buckling and initial post-buckling behaviour of channel member flange - analytical approach
PublicationThe local buckling and initial post-buckling behaviour of the cold-formed channel member flange is investigated. The governing nonlinear differential equation for axially compressed columns and beams undergoing pure bending is derived using the stationary total potential energy principle. The critical stress and initial post-buckling equilibrium path is determined by means of a perturbation approach. The results obtained allow...
-
Prediction of coking dynamics for wet coal charge
PublicationA one-dimensional transient mathematical model describing thermal and flow phenomena during coal coking in an oven chamber was studied in the paper. It also accounts for heat conduction in the ceramic oven wall when assuming a constant temperature at the heating channel side. The model was solved numerically using partly implicit methods for gas flow and heat transfer problems. The histories of temperature, gas evolution and internal...
-
Optical Magnetometry Based on Nanodiamonds with Nitrogen-Vacancy Color Centers
PublicationNitrogen-vacancy color centers in diamond are a very promising medium for many sensing applications such as magnetometry and thermometry. In this work, we study nanodiamonds deposited from a suspension onto glass substrates. Fluorescence and optically detected magnetic resonance spectra recorded with the dried-out nanodiamond ensembles are presented and a suitable scheme for tracking the magnetic-field value using a continuous...
-
Parametrized Local Reduced-Order Models With Compressed Projection Basis for Fast Parameter-Dependent Finite-Element Analysis
PublicationThis paper proposes an automated parametric local model-order reduction scheme for the expedited design of microwave devices using the full-wave finite-element method (FEM). The approach proposed here results in parameterized reduced-order models (ROMs) that account for the geometry and material variation in the selected subregion of the structure. In each subregion, a parameter-dependent projection basis is generated by concatenating...
-
Flow structure, heat transfer and scaling analysis in the case of thermo-magnetic convection in a differentially heated cylindrical enclosure
PublicationThe experimental, numerical and scaling analysis in the case of thermo-magnetic convection in a thermosyphon-like enclosure filled with a paramagnetic fluid is presented. Visualization of temperature field together with the numerical simulation gave an information about the flow structure, which indicated “finger-like” structures of hot and cold streams advecting each other. Their number depended on the Rayleigh number and also...
-
Greedy Multipoint Model-Order Reduction Technique for Fast Computation of Scattering Parameters of Electromagnetic Systems
PublicationThis paper attempts to develop a new automated multipoint model-order reduction (MOR) technique, based on matching moments of the system input–output function, which would be suited for fast and accurate computation of scattering parameters for electromagnetic (EM) systems over a wide frequency band. To this end, two questions are addressed. Firstly, the cost of the wideband reduced model generation is optimized by automating a...
-
Flow through a prosthetic mechanical aortic valve: Numerical model and experimental study
PublicationThis research presents a numerical model dedicated for virtual patient diagnostics in the field of synthetic valve implantation. The model operates based on computational fluid dynamics solver with implemented rigid body motion solver. Characteristic indicators related to the prosthetic valve were determined to assess the correctness of cardiac system operation after implantation. A novel approach for dynamic time discretization...
-
DL_MG: A Parallel Multigrid Poisson and Poisson–Boltzmann Solver for Electronic Structure Calculations in Vacuum and Solution
PublicationThe solution of the Poisson equation is a crucial step in electronic structure calculations, yielding the electrostatic potential -- a key component of the quantum mechanical Hamiltonian. In recent decades, theoretical advances and increases in computer performance have made it possible to simulate the electronic structure of extended systems in complex environments. This requires the solution of more complicated variants of the...
-
Numerical Analysis of Steady Gradually Varied Flow in Open Channel Networks with Hydraulic Structures
PublicationIn this paper, a method for numerical analysis of steady gradually varied fl ow in channel networks with hydraulic structures is considered. For this purpose, a boundary problem for the system of ordinary differential equations consisting of energy equation and mass conservation equations is formulated. The boundary problem is solved using fi nite difference technique which leads to the system of non-linear algebraic equations....
-
Nonlinear planar modeling of massive taut strings travelled by a force-driven point-mass
PublicationThe planar response of horizontal massive taut strings, travelled by a heavy point-mass, either driven by an assigned force, or moving with an assigned law, is studied. A kinematically exact model is derived for the free boundary problem via a variational approach, accounting for the singularity in the slope of the deflected string. Reactive forces exchanged between the point-mass and the string are taken into account via Lagrange...
-
Long-term hindcast simulation of water temperature and salinity in the Baltic Sea
Open Research DataThe dataset contains the results of numerical modelling of water temperature and salinity over a period of 50 years (1958-2007) in the Baltic Sea. A long-term hindcast simulation was performed using a three-dimensional hydrodynamic model PM3D (Kowalewski and Kowalewska-Kalkowska, 2017), a new version of the M3D model (Kowalewski, 1997). The hydrodynamic...
-
Long-term hindcast simulation of currents in the Baltic Sea
Open Research DataThe dataset contains the results of numerical modelling of currents over a period of 50 years (1958-2007) in the Baltic Sea . A long-term hindcast simulation was performed using a three-dimensional hydrodynamic model PM3D (Kowalewski and Kowalewska-Kalkowska, 2017), a new version of the M3D model (Kowalewski, 1997). The hydrodynamic model was coupled...
-
Method of reconstructing two-dimensional velocity fields on the basis of temperature field values measured with a thermal imaging camera
PublicationThis paper describes a novel numerical reconstruction procedure (NRP) of the velocity field during natural convective heat transfer from a two-sided, isothermal, heated vertical plate based only on the known temperature field obtained, e.g. with a thermal imaging camera. It has been demonstrated that with a knowledge of temperature distributions, the NRP enables the reconstruction of velocity fields by solving the Navier-Stokes...
-
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....
-
Hydrodynamic reanalysis of sea level in the Baltic Sea using the PM3D model
Open Research DataThe data set contains the results of numerical modelling of sea level fluctuations in the Baltic Sea in the Baltic Sea since 1998. A long-term reanalysis was performed using a three-dimensional hydrodynamic model PM3D (Kowalewski and Kowalewska-Kalkowska, 2017), a new version of the M3D model (Kowalewski, 1997).
-
Numerical Test for Stability Evaluation of Analog Circuits
PublicationIn this contribution, a new numerical test for the stability evaluation of analog circuits is presented. Usually, if an analog circuit is unstable then the roots of its characteristic equation are localized on the right half-plane of the Laplace s- plane. Because this region is unbounded, we employ the bilinear transformation to map it into the unit disc on the complex plane. Hence, the existence of any root inside the unit disc...
-
Numerical Simulations of Seepage in Dikes Using unsaturated and Two-Phase Flow Models
PublicationModeling of water flow in variably saturated porous media, including flood dikes, is often based on the Richards equation, which neglects the flow of pore air, assuming that it remains at constant atmospheric pressure. However, there is also evidence that the air flow can be important, especially when the connectivity between the pore air and atmospheric air is lost. In such cases a full two-phase air-water flow model should be...
-
Simulations of air and water flow in a model dike during overflow experiments
PublicationFlow in flood dikes, earth dams, and embankments occurs in variably saturated conditions, with pores of the earth material filled partly with water and partly with air. In routine engineering analysis, the influence of pore air is neglected and the air pressure is assumed equal to atmospheric. In some circumstances, for example, during overtopping of the dike by water, the effect of pore air on water flow and stability of the structure...
-
Local buckling of composite channel columns
PublicationThe investigation concerns local buckling of compressed flanges of axially compressed composite channel columns. Cooperation of the member flange and web is taken into account here. The buckling mode of the member flange is defined by rotation angle a flange about the line of its connection with the web. The channel column under investigation is made of unidirectional fibre-reinforced laminate. Two approaches to member orthotropic...
-
Air trapping problem during infiltration on the large areas
PublicationThe process of flow modeling in unsaturated porous medium is often found in many fields of sciences: geology, fluid mechanics, thermodynamics, microbiology or chemistry. Problem is relatively complicated due to complexity of the system which contains three phases: water, air and soil skeleton. The flow of water in such a medium can be described using two-phase (2PH) flow formulation, which accounts the inflow of air and water phases,...
-
Guessing Intrinsic Forwarding Trustworthiness of Wireless Ad Hoc Network Nodes
PublicationA novel node misbehavior detection system called GIFTED is proposed for a multihop wireless ad hoc network (WAHN) whose nodes may selfishly refuse to forward transit packets. The system guesses the nodes’ intrinsic forwarding trustworthiness (IFT) by analyzing end-to-end path performance rather than utilizing unreliable and incentive incompatible low-layer mechanisms. It can work with occasional IFT jumps, directional antennae,...
-
Absorbing Boundary Conditions Derived Based on Pauli Matrices Algebra
PublicationIn this letter, we demonstrate that a set of absorbing boundary conditions (ABCs) for numerical simulations of waves, proposed originally by Engquist and Majda and later generalized by Trefethen and Halpern, can alternatively be derived with the use of Pauli matrices algebra. Hence a novel approach to the derivation of one-way wave equations in electromagnetics is proposed. That is, the classical wave equation can be factorized...
-
Optimal shape and stress control of geometrically nonlinear structures with exact gradient with respect to the actuation inputs
PublicationThis paper presents an efficient and robust optimization methodology for stress and shape control of actuated geometrically nonlinear elastic structures, applied to 3D trusses. The actuation inputs, modeled as prescribed strains, serve as the optimization variables. The objective is to minimize total actuation while satisfying several constraints: (i) actuation bounds in each actuated element and (ii) target ranges for nodal displacements...
-
NUMERYCZNE MODELOWANIE HYDRODYNAMIKI JAKO NARZĘDZIE WSPOMAGANIA PROJEKTOWANIA REKREACYJNYCH I SPORTOWYCH OBIEKTÓW WODNYCH
PublicationW artykule opisano zastosowanie obliczeń numerycznych w procesie projektowania, analizy działania oraz oceny bezpieczeństwa użytkowania wodnych obiektów rekreacyjnych i sportowych. Użytkownicy aquaparków, sportowcy, a szczególnie kibice, oczekują wyjątkowych doznań na obiektach wodnych z jednoczesną gwarancją bezpieczeństwa. Na te kwestie duży wpływ ma hydrodynamika przepływu. Aby właściwie dobrać parametry geometryczne i hydrauliczne...
-
On the Structure of Time in Computational Semantics of a Variable-Step Solver for Hybrid Behavior Analysis
PublicationHybrid dynamic systems combine continuous and discrete behavior. Often, computational approaches are employed to derive behaviors that approximate the analytic solution. An important part of this is the approximation of differential equation behavior by numerical integration. The accuracy and computational efficiency of the integration usually depend on the complexity of the method and its implicated approximation errors, especially...
-
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 two-dimensional 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...
-
Assessing groundwater vulnerability to pollution in the Puck region (denudation moraine upland) using vertical seepage method
PublicationDegradation of groundwater quality can cause a serious water supply and environmental problems. The identify of potential groundwater pollution can be determined by assessment of groundwater vulnerability method. The assessment of groundwater vulnerability to pollution was based on estimation of migration time of potential conservative contamination through the vadose zone. Area of investigation is a type of denudation moraine...
-
Finite-difference time-domain analyses of active cloaking for electrically-large objects
PublicationInvisibility cloaking devices constitute a unique and potentially disruptive technology, but only if they can work over broad bandwidths for electrically-large objects. So far, the only known scheme that allows for broadband scattering cancellation from an electrically-large object is based on an active implementation where electric and magnetic sources are deployed over a surface surrounding the object, but whose ‘switching on’...
-
Modeling of the Two-Dimensional Flow Caused by Sea Conditions and Wind Stresses on the Example of Dead Vistula
PublicationThe article presents the results of two-dimensional modeling of flows caused by the sea conditions and wind stresses on the example of Dead Vistula. Based on the available bathymetric data, a numerical model of the river section was created, which was supplemented with data on the position of the water table depending on hydrometeorological conditions. To describe the flow field in steady conditions, a simplified model of two-dimensional...
-
METHOD FOR SHIP'S ROLLING PERIOD PREDICTION WITH REGARD TO NON-LINEARITY OF GZ CURVE
PublicationThe paper deals with the problem of prediction of the rolling period. A special emphasis is put on the practical application of the new method for rolling period prediction with regard to non-linearity of the GZ curve. The one degree-of-freedom rolling equation is applied with using the non-linear stiffness moment and linear damping moment formulas. A number of ships are considered to research the discrepancies between the pending...
-
Torsion of restrained thin-walled bars of open constant bisymmetric cross-section
PublicationElastic and geometric stiffness matrices were derived using Castigliano's first theorem, for the case of torsion of restrained thin-walled bars of open constant bisymmetric cross-section. Functions which describe the angles of torsion were adopted from the solutions of thedifferential equation for restrained torsion. The exact solutions were simplified by expanding them in a power series. Numerical examples were taken from Kujawa...
-
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...
-
Distortional buckling of composite thin-walled columns of a box-type cross section with diaphragms
PublicationDistortional buckling of axially compressed columns of box-like composite cross sections with andwithout internal diaphragms is investigated in the framework of one-dimensional theory. The channel membersare composed of unidirectional fibre-reinforced laminate. Two approaches to the member orthotropic materialare applied: homogenization based on the theory of mixture and periodicity cells, and homogenization basedon the Voigt–Reuss...
-
Nonlinear phenomena of small-scale sound in a gas with exponential stratification
PublicationThe nonlinear dynamics of perturbations, quickly varying in space, with comparatively large characteristic wavenumbers k: k>1/H, is considered. H is the scale of density and pressure reduction in unperturbed gas, as the coordinate (H is the so-called height of the uniform equilibrium gas). Coupling nonlinear equations which govern the sound and the entropy mode in a weakly nonlinear flow are derived. They describe the dynamics...
-
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...
-
Coherent-wave Monte Carlo method for simulating light propagation in tissue
PublicationSimulating propagation and scattering of coherent light in turbid media, such as biological tissues, is a complex problem. Numerical methods for solving Helmholtz or wave equation (e.g. finite-difference or finite-element methods) require large amount of computer memory and long computation time. This makes them impractical for simulating laser beam propagation into deep layers of tissue. Other group of methods, based on radiative...
-
A Novel Trust-Region-Based Algorithm with Flexible Jacobian Updates for Expedited Optimization of High-Frequency Structures
PublicationSimulation-driven design closure is mandatory in the design of contemporary high-frequency components. It aims at improving the selected performance figures through adjustment of the structure’s geometry (and/or material) parameters. The computational cost of this process when employing numerical optimization is often prohibitively high, which is a strong motivation for the development of more efficient methods. This is especially...
-
Acceleration of Electromagnetic Simulations on Reconfigurable FPGA Card
PublicationIn this contribution, the hardware acceleration of electromagnetic simulations on the reconfigurable field-programmable-gate-array (FPGA) card is presented. In the developed implementation of scientific computations, the matrix-assembly phase of the method of moments (MoM) is accelerated on the Xilinx Alveo U200 card. The computational method involves discretization of the frequency-domain mixed potential integral equation using...
-
Determination of Thermophysical Parameters Involved in The Numerical Model to Predict the Temperature Field of Cast-In-Place Concrete Bridge Deck
PublicationThe paper dealswith a concept of a practical computationmethod to simulate the temperature distribution in an extradosed bridge deck. The main goal of the study is to develop a feasible model of hardening of concrete consistent with in-situ measurement capabilities. The presented investigations include laboratory tests of high performance concrete, measurements of temperature evolution in the bridge deck and above all, numerical...
-
Rapid Simulation-Driven Multiobjective Design Optimization of Decomposable Compact Microwave Passives
PublicationIn this paper, a methodology for fast multiobjective optimization of the miniaturized microwave passives has been presented. Our approach is applicable to circuits that can be decomposed into individual cells [e.g., compact microstrip resonant cells (CMRCs)]. The structures are individually modeled using their corresponding equivalent circuits and aligned with their accurate, EM simulated...