Search
Search results for: numerical solution
-
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...
-
Comparison of strain results at a laser weld notch obtained by numerical calculations and experimental measurements
PublicationIn the development of ship structures applying new materials and it’s purposeful placement play an important role. During the last years, especially in a construction of ro-ro type vessels, the usage of novel sandwich structures in cargo decks is profitable. Steel sandwich panel is an innovative solution which at a todays state of development can be used for the construction of any members not taking part in a global bending of...
-
Examples of numerical simulations of two-dimensional unsaturated flow with VS2DI code using different interblock conductivity averaging schemes
PublicationFlow in unsaturated porous media is commonly described by the Richards equation. This equation is strongly nonlinear due to interrelationships between water pressure head (negative in unsaturated conditions), water content and hydraulic conductivity. The accuracy of numerical solution of the Richards equation often depends on the method used to estimate average hydraulic conductivity between neighboring nodes or cells of the numerical...
-
Advanced numerical modelling for predicting residual compressive strength of corroded stiffened plates
PublicationAn advanced methodology for predicting the residual compressive strength of corroded stiffened plates is developed here using the non-linear finite element method. The non-uniform loss of a plate thickness is accounted for on a macro-scale. In contrast, mechanical properties are changed using the constitutive model to reflect the corrosion degradation impact on a micro-scale. Three different stiffened plate thicknesses are considered,...
-
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...
-
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...
-
Numerical analysis of the flow rig for UWS spray examination in exhaust system-relevant conditions
PublicationIn the present study, a flow rig with optical access intended for spray investigations in exhaust system-relevant conditions was analysed in terms of flow and temperature in the spray area using numerical simulations. The operation of the rig was examined for a wide range of exhaust mass flow rates, temperatures and various forms of UWS (urea-water solution) spray plumes. The locations of the injector and thermocouple were verified....
-
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....
-
Cavity-expansion approximation for projectile impact and penetration into sand
PublicationA one-dimensional problem of a spherical cavity expanding at a constant velocity from zero initial radius in an infinite granular medium, which has the first-kind self-similar solution, is considered. We are solving this dynamic spherical cavity-expansion problem to model rigid spheres penetrating into a granular media. Elastic–plastic deformation of the granular media is described in a barotropic approximation, using the high-pressure...
-
The calculation long-term physical field created by water flow around the ship
PublicationThe paper presents a description of the problems with fluid flow around a ship. Using the described solution of the problem were performed numerical calculations using the boundary element method. Were also presented preliminary results of the calculation of pressure fields at the bottom of the ship, taking into account the impact of the bottom.
-
On neutral differential equations and the monotone iterative method
PublicationThe application of the monotone iterative method to neutral differential equations with deviating arguments is considered in this paper. We formulate existence results giving sufficient conditions which guarantee that such problems have solutions. This approach is new and to the Authors' knowledge, this is the first paper when the monotone iterative method is applied to neutral first-order differential equations with deviating...
-
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...
-
Application of the Boundary Element Method for the Simulation of Two-dimensional Viscous Incompressible Flow
PublicationThe paper presents the application of an indirect variant of the boundary element method (BEM) to solve the two-dimensional steady flow of a Stokes liquid. In the BEM, a system of differential equations is transformed into integral equations. Thi smakes it possible to limit discretization to the border of the solution. Numerical discretization of the computational domain was performed with linear boundary elements, for which a...
-
Theoretical analysis of the cable-stayed bridge over brda river in Bydgoszcz
PublicationThe topic of this article is the numerical analysis carried out for the cable-stayed tram (road) bridge over Brda river in Bydgoszcz. The goal of numerical studies was the verification of project assumptions and construction solutions elaborated by the autonomous team. The bridge calculations were conducted on three independent FEM models The internal forces and the elements effort under different loading schemes were examined....
-
Solving the Problem of Dynamic Adaptability of Artificial Intelligence Systems that Control Dynamic Technical Objects
PublicationThis paper investigates the increase in the response speed and stability of artificial intelligence systems that control dynamic technical objects. The problem of calculating the optimal time of switching an artificial intelligence system between software classes by the criterion of the rigidity degree of the model of a control object is considered. The solution of this problem is proposed for the general case of the control object...
-
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...
-
Selected local stability problems of channel section flanges made of aluminium alloys
PublicationThe paper addresses the issue of local buckling of compressed flanges of cold-formed thin-walled channel columns and beams with nonstandard flanges composed of aluminium alloys. The material behaviour follows the Ramberg–Osgood law. It should be noted that the proposed solution may be also applied for other materials, for example: stainless steel, carbon steel. The paper is motivated by an increasing interest in nonstandard cold-formed...
-
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...
-
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...
-
Application and mechanical properties of aluminium alloys
PublicationThe paper presents a short survey of the application and mechanical properties of aluminium alloys. This area of scientific research is still developing, the authors restricted it to about 30 papers. Numerical simulations of the aluminium bracket are also investigated. The authors propose a solution to increase the load capacity of aluminium brackets with application of the additional cover plates. This study is proposed as an...
-
Torsional buckling and post-buckling of columns made of aluminium alloy
PublicationThe paper concerns torsional buckling and the initial post-buckling of axially compressed thin-walled aluminium alloy columns with bisymmetrical cross-section. It is assumed that the column material behaviour is described by the Ramberg–Osgood constitutive equation in non-linear elastic range. The stationary total energy principle is used to derive the governing non-linear differential equation. An approximate solution of the equation...
-
On Wrinkling in Sandwich Panels with an Orthotropic Core
PublicationThis paper deals with the local loss of stability (wrinkling) problem of a thin facing of a sandwich panel. Classical solutions to the problem of facing instability resting on a homogeneous and isotropic substructure (a core) are compared. The relations between strain energy components associated with different forms of core deformations are discussed. Next, a new solution for the orthotropic core is presented in detail, which...
-
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...
-
Investigation of the Aerodynamics of an Innovative Vertical-Axis Wind Turbine
PublicationThe paper presents the results of an experimental investigation and a three dimensional numerical analysis of the transient aerodynamic phenomena occurring in the innovative modification of classic Savonius wind turbine. An attempt to explain the increased efficiency of the innovative design a comparison with the traditional solution is undertaken. A vorticity measure based on the integral of the velocity gradient tensor second...
-
Reduced-order models in the finite element analysis
PublicationA novel technique of incorporating macromodels into finite element electromagnetic analysis of waveguide components is presented. Macromodels are generated by using a model order reduction algorithm (ENOR), which results in significant decrease of the number of variables, that describe the computational region. Proposed technique allows for using a few independent macromodels as well as to duplicating one macromodel in many subregions...
-
Results of the application of tropospheric corrections from different troposphere models for precise GPS rapid static positioning
PublicationIn many surveying applications, determination of accurate heights is of significant interest. The delay caused by the neutral atmosphere is one of the main factors limiting the accuracy of GPS positioning and affecting mainly the height coordinate component rather than horizontal ones. Estimation of the zenith total delay is a commonly used technique for accounting for the tropospheric delay in static positioning. However, in the...
-
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...
-
CONSTRUCTION TECHNOLOGY OF TIMBER-FRAME HOUSES RESISTANT TO DYNAMIC LOADS – STUDY ON MODELS OF EXTERIOR WALLS
PublicationThe aim of this paper is to show the numeric representation of experimental studies concerning the behaviour of exterior wall models of a timber-frame house under harmonic loading. A single wall model according to traditional technology of timberframe house walls (filling with mineral wool) was tested. The analysis was conducted for the following frequencies: 0.5 Hz, 1.0 Hz, 2.0 Hz and 5.0 Hz for various values of the specified...
-
The smoothness test for a density function
PublicationThe problem of testing hypothesis that a density function has no more than μ derivatives versus it has more than μ derivatives is considered. For a solution, the L2 norms of wavelet orthogonal projections on some orthogonal ‘‘differences’’ of spaces from a multiresolution analysis is used. For the construction of the smoothness test an asymptotic distribution of a smoothness estimator is used. To analyze that asymptotic distribution,...
-
The Optimal Location of Ground-Based GNSS Augmentation Transceivers
PublicationModern Global Navigation Satellite Systems (GNSS) allow for positioning with accuracies ranging from tens of meters to single millimeters depending on user requirements and available equipment. A major disadvantage of these systems is their unavailability or limited availability when the sky is obstructed. One solution is to use additional range measurements from ground-based nodes located in the vicinity of the receiver. The highest...
-
Two-criteria optimisation problem for ventral hernia repair
PublicationTwo-criteria optimisation problem related to laparoscopic ventral hernia repair is formulated in this paper. An optimal implant from a given set and its orientation is sought. The implant is subjected to kinematic extortions due to a patient’s body movement and intra-abdominal pressure. The first criterion of the optimisation problem deals with the reaction force in the implant fastener, while the deflection of the implant constitutes...
-
Novel Analytic-Numerical Model of Free Convection: with Leading Edge Considered
PublicationA novel solution of the free convection boundary problem is represented in analytical form for velocity and temperature for an isothermal vertical plate, as an example. These fields are built as a Taylor Series in the x coordinate with coefficients as functions of the vertical coordinate (y). We restrict ourselves by cubic approximation for both functions. The basic Navier-Stokes and Fourier-Kirchhoff equations and boundary conditions...
-
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...
-
On rotational instability within the nonlinear six-parameter shell theory
PublicationWithin the six-parameter nonlinear shell theory we analyzed the in-plane rotational instability which oc- curs under in-plane tensile loading. For plane deformations the considered shell model coincides up to notations with the geometrically nonlinear Cosserat continuum under plane stress conditions. So we con- sidered here both large translations and rotations. The constitutive relations contain some additional mi- cropolar parameters...
-
Non-Least Square GNSS Positioning Algorithm for Densely Urbanized Areas
PublicationThe paper introduces an essentially new algorithm for calculating the GNSS position as an alternative to the least-square method. The proposed approach can be widely applied to any positioning method that uses multiple position lines for position calculation and is an example ofhow using a numerical solution can improve position accuracy without access to historical data. In essence, the method is based on the adaptation of the...
-
Numerical Methods for Partial Differential Equations
e-Learning CoursesCourse description: This course focuses on modern numerical techniques for linear and nonlinear elliptic, parabolic and hyperbolic partial differential equations (PDEs), and integral equations fundamental to a large variety of applications in science and engineering. Topics include: formulations of problems in terms of initial and boundary value problems; finite difference and finite element discretizations; boundary element approach;...
-
GIS Solution for Weather Forecast Data Analysis
PublicationIn this paper authors present the GIS system for the analysis of the numerical weather prediction data. This kind of data has multidimensional character (three dimensions and time) and its analysis should consider all the available factors. Proposed GIS system consists of RASDAMAN application with implemented OLAP cube mechanism, which enables the user to process data in the spatial-time domain. It also simplifies the meteorological...
-
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...
-
Diagonalized Macromodels in Finite Element Method for Fast Electromagnetic Analysis of Waveguide Components
PublicationA new technique of local model-order reduction (MOR) in 3-D finite element method (FEM) for frequency-domain electromagnetic analysis of waveguide components is proposed in this paper. It resolves the problem of increasing solution time of the reduced-order system assembled from macromodels created in the subdomains, into which an analyzed structure is partitioned. This problem becomes particularly relevant for growing size and...
-
Modal Adjustment of Rayleigh Based Structural Damping and Coordinate-Partitioning Algorithm Dedicated to Frictionless Contact Constraints between Multibody System and Structure Modelled with Finite Elements
PublicationThe paper presents a dedicated numerical algorithm. The algorithm is advantageous during investigations of the dynamics of a hybrid multibody / finite-elements system. We focus our attention on interactions resulting from mechanical contact. Pointwise contact connects a vertex of the multibody structure and a surface of the elastic reference body. Instead of a positive value of the relative penetration factor, constraint equations...
-
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...
-
Control of the bridge span vibration with high coefficient passive damper. Theoretical consideration and application
PublicationThe research was carried out due to the problem of vibration on the lively pedestrian drawbridge across the Motlawa River in the city of Gdansk. In the design stage, the main span of the footbridge showed unfavorable dynamic properties, which may create a comfort problem for pedestrians. The first vertical bending eigenfrequency was recognized as 1.64 Hz. The original design of the footbridge was equipped with a driving cylinder...
-
Damaged plasticity model - seismic analysis of structure of the Petronas Twin Towers (Kuala Lumpur) in Malaysia
PublicationThe article presents numerical simulations for the modelling of seismic impact on the structure of Petronas Twin Towers (Kuala Lumpur) in Malaysia with application of two methods. The Response Spectrum method, in which a spectrum of the structure's responses to an earthquake's impact is generated and the Acceletogram (BETA) method, in which we generate dynamic load in the form of a diagram of the connection between acceleration...
-
MULTI-OBJECTIVE OPTIMIZATION PROBLEM IN THE OptD-MULTI METHOD
PublicationNew measurement technologies, e.g. Light Detection And Ranging (LiDAR), generate very large datasets. In many cases, it is reasonable to reduce the number of measuring points, but in such a way that the datasets after reduction satisfy specific optimization criteria. For this purpose the Optimum Dataset (OptD) method proposed in [1] and [2] can be applied. The OptD method with the use of several optimization criteria is called...
-
Multiscale model for blood flow after a bileaflet artificial aortic valve implantation
PublicationCardiovascular diseases are the leading cause of mortality in the world, mainly due to atherosclerosis and its consequences. The article presents the numerical model of the blood flow through artificial aortic valve. The overset mesh approach was applied to simulate the valve leaflets motion and to realize the moving mesh, in the aortic arch and the main branches of cardiovascular system. To capture the cardiac system’s response...
-
Four Degree-of-Freedom Hydrodynamic Maneuvering Model of a Small Azipod-Actuated Ship With Application to Onboard Decision Support Systems
PublicationThe main contribution of this paper is a numerical ship motion model of NTNU’s research vessel Gunnerus, capturing the surge, sway, roll, and yaw dynamics when sailing in uniform and steady currents. The model utilizes a crossflow drag formulation for the transverse viscous loads, and it includes a nonlinear formulation for the propulsion and steering loads provided by two azipod thrusters. A wide range of experimental data obtained...
-
EVALUATION OF EFFECTIVENESS OF WATERJET PROPULSOR FOR A SMALL UNDERWATER VEHICLE
PublicationT he goal of the project described is to replace the existing propulsion system of a small underwater vehicle with a solution less prone to mechanical damage and ensuring a lower risk of the entanglement of fibrous objects suspended in the body of water. Four typical marine screws are utilised in the current design of the vehicle. One possible solution of the problem is the application of waterjet propulsors located...
-
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...
-
Study on 19th-century cast iron columns from the former financial office building in Kwidzyn
PublicationThe article presents the results for cast iron columns used as floor supports in the former financial office building in Kwidzyn (Poland).The stocktaking drawings provided the general and detailed data about the 19th century designer’s decisions. During the research, architectural and structural analysis was carried out. The graphics were used as a starting material for preparing 3D calculation models. The FEM-based numerical simulations...
-
GPU-Accelerated 3D Mesh Deformation for Optimization Based on the Finite Element Method
PublicationThis paper discusses a strategy for speeding up the mesh deformation process in the design-byoptimization of high-frequency components involving electromagnetic field simulations using the 3D finite element method (FEM). The mesh deformation is assumed to be described by a linear elasticity model of a rigid body; therefore, each time the shape of the device is changed, an auxiliary elasticity finite-element problem must be solved....