Search results

The main purpose of this paper is to introduce the gradient discretisation method (GDM) to a system of reaction diffusion equations subject to non-homogeneous Dirichlet boundary conditions. Then, the authors show that the GDM provides a comprehensive convergence analysis of several numerical methods for the considered model. The convergence is established without non-physical regularity assumptions on the solutions.

In this paper, the authors use the GDM to discretise a system of reaction diffusion equations with non-homogeneous Dirichlet boundary conditions.

The authors provide a generic convergence analysis of a system of reaction diffusion equations. The authors introduce a specific example of numerical scheme that fits in the gradient discretisation method. The authors conduct a numerical test to measure the efficiency of the proposed method.

This work provides a unified convergence analysis of several numerical methods for a system of reaction diffusion equations. The generic convergence is proved under the classical assumptions on the solutions.

This study aims to develop an experimentally validated three-dimensional numerical model for predicting different flow patterns produced with a gas dynamic virtual nozzle (GDVN).

The physical model is posed in the mixture formulation and copes with the unsteady, incompressible, isothermal, Newtonian, low turbulent two-phase flow. The computational fluid dynamics numerical solution is based on the half-space finite volume discretisation. The geo-reconstruct volume-of-fluid scheme tracks the interphase boundary between the gas and the liquid. To ensure numerical stability in the transition regime and adequately account for turbulent behaviour, the k-ω shear stress transport turbulence model is used. The model is validated by comparison with the experimental measurements on a vertical, downward-positioned GDVN configuration. Three different combinations of air and water volumetric flow rates have been solved numerically in the range of Reynolds numbers for airflow 1,009–2,596 and water 61–133, respectively, at Weber numbers 1.2–6.2.

The half-space symmetry allows the numerical reconstruction of the dripping, jetting and indication of the whipping mode. The kinetic energy transfer from the gas to the liquid is analysed, and locations with locally increased gas kinetic energy are observed. The calculated jet shapes reasonably well match the experimentally obtained high-speed camera videos.

The model is used for the virtual studies of new GDVN nozzle designs and optimisation of their operation.

To the best of the authors’ knowledge, the developed model numerically reconstructs all three GDVN flow regimes for the first time.

Performing accurate numerical simulations of electrical drives, the precise knowledge of the local magnetic material properties is of utmost importance. Due to the various manufacturing steps, e.g. heat treatment or cutting techniques, the magnetic material properties can strongly vary locally, and the assumption of homogenized global material parameters is no longer feasible. This paper aims to present the general methodology and two different solution strategies for determining the local magnetic material properties using reference and simulation data.

The general methodology combines methods based on measurement, numerical simulation and solving an inverse problem. Therefore, a sensor-actuator system is used to characterize electrical steel sheets locally. Based on the measurement data and results from the finite element simulation, the inverse problem is solved with two different solution strategies. The first one is a quasi Newton method (QNM) using Broyden's update formula to approximate the Jacobian and the second is an adjoint method. For comparison of both methods regarding convergence and efficiency, an artificial example with a linear material model is considered.

The QNM and the adjoint method show similar convergence behavior for two different cutting-edge effects. Furthermore, considering a priori information improved the convergence rate. However, no impact on the stability and the remaining error is observed.

The presented methodology enables a fast and simple determination of the local magnetic material properties of electrical steel sheets without the need for a large number of samples or special preparation procedures.

This paper focuses on the unconditionally optimal error estimates of a linearized second-order scheme for a nonlocal nonlinear parabolic problem. The first step of the scheme is based on Crank–Nicholson method while the second step is the second-order BDF method.

A rigorous error analysis is done, and optimal L² error estimates are derived using the error splitting technique. Some numerical simulations are presented to confirm the study’s theoretical analysis.

Optimal L² error estimates and energy norm.

The goal of this research article is to present and establish the unconditionally optimal error estimates of a linearized second-order BDF finite element scheme for the reaction-diffusion problem. An optimal error estimate for the proposed methods is derived by using the temporal-spatial error splitting techniques, which split the error between the exact solution and the numerical solution into two parts, that is, the temporal error and the spatial error. Since the spatial error is not dependent on the time step, the boundedness of the numerical solution in L∞-norm follows an inverse inequality immediately without any restriction on the grid mesh.

In fire resistance tests (FRTs) of building materials, a crucial criterion to pass the test procedure is to avoid the leakage of the hot flue gases caused by gaps and cracks occurring due to the thermal exposure. The present study's aim is to calculate the deformation of a steel door, which is embedded within a wall made of bricks, and qualitatively determine the flue gas leakage.

A computational fluid dynamics/finite element method (CFD/FEM) coupling was introduced representing an intermediate approach between a one-way and a full two-way coupling methodology, leading to a simplified two-way coupling (STWC). In contrast to a full two way-coupling, the heat transfer through the steel door was simulated based on a one-way approach. Subsequently, the predicted temperatures at the door from the one-way simulation were used in the following CFD/FEM simulation, where the fluid flow inside and outside the furnace as well as the deformation of the door were calculated simultaneously.

The simulation showed large gaps and flue gas leakage above the door lock and at the upper edge of the door, which was in close accordance to the experiment. Furthermore, it was found that STWC predicted similar deformations compared to the one-way coupling.

Since two-way coupling approaches for fluid/structure interaction in fire research are computationally demanding, the number of studies is low. Only a few are dealing with the flue gas exit from rooms due to destruction of solid components. Thus, the present study is the first two-way approach dealing with flue gas leakage due to gap formation.

This study aims to reduce the redundant weight of the anti-roll torsion bar brought by the traditional empirical design and improving its strength and stiffness.

Based on the finite element approach coupled with the improved beluga whale optimization (IBWO) algorithm, a collaborative optimization method is suggested to optimize the design of the anti-roll torsion bar structure and weight. The dimensions and material properties of the torsion bar were defined as random variables, and the torsion bar's mass and strength were investigated using finite elements. Then, chaotic mapping and differential evolution (DE) operators are introduced to improve the beluga whale optimization (BWO) algorithm and run case studies.

The findings demonstrate that the IBWO has superior solution set distribution uniformity, convergence speed, solution correctness and stability than the BWO. The IBWO algorithm is used to optimize the anti-roll torsion bar design. The error between the optimization and finite element simulation results was less than 1%. The weight of the optimized anti-roll torsion bar was lessened by 4%, the maximum stress was reduced by 35% and the stiffness was increased by 1.9%.

The study provides a methodological reference for the simulation optimization process of the lateral anti-roll torsion bar.

The “chiralisation” of Euclidean polygonal tessellations is a novel, recent method which has been used to design new auxetic metamaterials with complex topologies and improved geometric versatility over traditional chiral honeycombs. This paper aims to design and manufacture chiral honeycombs representative of four distinct classes of 2D Euclidean tessellations with hexagonal rotational symmetry using fused-deposition additive manufacturing and experimentally analysed the mechanical properties and failure modes of these metamaterials.

Finite Element simulations were also used to study the high-strain compressive performance of these systems under both periodic boundary conditions and realistic, finite conditions. Experimental uniaxial compressive loading tests were applied to additively manufactured prototypes and digital image correlation was used to measure the Poisson’s ratio and analyse the deformation behaviour of these systems.

The results obtained demonstrate that these systems have the ability to exhibit a wide range of Poisson’s ratios (positive, quasi-zero and negative values) and stiffnesses as well as unusual failure modes characterised by a sequential layer-by-layer collapse of specific, non-adjacent ligaments. These findings provide useful insights on the mechanical properties and deformation behaviours of this new class of metamaterials and indicate that these chiral honeycombs could potentially possess anomalous characteristics which are not commonly found in traditional chiral metamaterials based on regular monohedral tilings.

To the best of the authors’ knowledge, the authors have analysed for the first time the high strain behaviour and failure modes of chiral metamaterials based on Euclidean multi-polygonal tessellations.

This paper comprehensively addresses the influence of chopped strand mat glass fiber-reinforced polymer (GFRP) patch configurations such as geometry, dimensions, position and the number of layers of patches, whether a single or double patch is used and how well debonding the area under the patch improves the strength of the cracked aluminum plates with different crack lengths.

Single-edge cracked aluminum specimens of 150 mm in length and 50 mm in width were tested using the tensile test. The cracked aluminum specimens were then repaired using GFRP patches with various configurations. A three-dimensional (3D) finite element method (FEM) was adopted to simulate the repaired cracked aluminum plates using composite patches to obtain the stress intensity factor (SIF). The numerical modeling and validation of ABAQUS software and the contour integral method for SIF calculations provide a valuable tool for further investigation and design optimization.

The width of the GFRP patches affected the efficiency of the rehabilitated cracked aluminum plate. Increasing patch width WP from 5 mm to 15 mm increases the peak load by 9.7 and 17.5%, respectively, if compared with the specimen without the patch. The efficiency of the GFRP patch in reducing the SIF increased as the number of layers increased, i.e. the maximum load was enhanced by 5%.

This study assessed repairing metallic structures using the chopped strand mat GFRP. Furthermore, it demonstrated the superiority of rectangular patches over semicircular ones, along with the benefit of using double patches for out-of-plane bending prevention and it emphasizes the detrimental effect of defects in the bonding area between the patch and the cracked component. This underlines the importance of proper surface preparation and bonding techniques for successful repair.

This paper empirically investigates whether trade openness (TO) in Brazil, Russia, India, China and South Africa (BRICS) countries affects how banks might employ loan loss provisions (LLPs) to smooth out their earnings and how adopting the International Financial Reporting Standards (IFRS) can mitigate it.

The analysis includes 78 commercial banks from five BRICS nations and spans 2014 through 2020. To test these hypotheses, the authors utilized a fixed-effect and two-step system panel generalized methods of moments (GMM) estimator.

TO positively affects income smoothing (earnings management) across BRICS commercial banks. The effect is clearer in banks that make financial reports under the IFRS. Path analysis reveals that the effect of TO is driven by nonperforming loans (NPLs). Additionally, the IFRS restricts earnings management in the BRICS banking sector when a better institutional environment is present. The authors found that accounting rules (IFRS) and enforcement (better institutional settings) interact to enhance earnings’ quality.

The relationship between TO and bank earnings management practices is important for understanding the complex interplay between trade and finance and ensuring financial stability, investor confidence and regulatory compliance. This study recommends better regulations and governance mechanisms for financial reports in emerging nations like BRICS. Additionally, macro-prudential regulators and banking supervisors should work closely to ensure transparent TO decisions with improved discipline, institutional quality and regulatory support to enhance bank stability.

The study finds evidence of bank income smoothing in the BRICS and introduces TO as a determinant. It also identifies the evolving role of IFRS in the presence of higher institutional quality and TO, thereby expanding the financial reporting literature.

The vibration of the rails is a significant source of railway rolling noise, often forming the dominant component of noise in the important frequency region between 400 and 2000 Hz. The purpose of the paper is to investigate the influence of the ground profile and the presence of the train body on the sound radiation from the rail.

Two-dimensional boundary element calculations are used, in which the rail vibration is the source. The ground profile and various different shapes of train body are introduced in the model, and results are observed in terms of sound power and sound pressure. Comparisons are also made with vibro-acoustic measurements performed with and without a train present.

The sound radiated by the rail in the absence of the train body is strongly attenuated by shielding due to the ballast shoulder. When the train body is present, the sound from the vertical rail motion is reflected back down toward the track where it is partly absorbed by the ballast. Nevertheless, the sound pressure at the trackside is increased by typically 0–5 dB. For the lateral vibration of the rail, the effects are much smaller. Once the sound power is known, the sound pressure with the train present can be approximated reasonably well with simple line source directivities.

Numerical models used to predict the sound radiation from railway rails have generally neglected the influence of the ground profile and reflections from the underside of the train body on the sound power and directivity of the rail. These effects are studied in a systematic way including comparisons with measurements.