Search results

1 – 10 of 160
Article
Publication date: 1 February 2016

S. Ali Faghidian

The linear regression technique is widely used to determine empirical parameters of fatigue life profile while the results may not continuously depend on experimental data. Thus…

Abstract

Purpose

The linear regression technique is widely used to determine empirical parameters of fatigue life profile while the results may not continuously depend on experimental data. Thus Tikhonov-Morozov method is utilized here to regularize the linear regression results and consequently reduces the influence of measurement noise without notably distorting the fatigue life distribution. The paper aims to discuss these issues.

Design/methodology/approach

Tikhonov-Morozov regularization method would be shown to effectively reduce the influences of measurement noise without distorting the fatigue life distribution. Moreover since iterative regularization methods are known to be an attractive alternative to Tikhonov regularization, four gradient iterative methods called as simple iteration, minimum error, steepest descent and conjugate gradient methods are examined with an appropriate initial guess of regularized coefficients.

Findings

It has been shown that in case of sparse fatigue life measurements, linear regression results may not have continuous dependence on experimental data and measurement error could lead to misinterpretations of the solution. Therefore from engineering safety point of view, utilizing regularization method could successfully reduce the influence of measurement noise without significantly distorting the fatigue life distribution.

Originality/value

An excellent initial guess for mixed iterative-direct algorithm is introduced and it has been shown that the combination of Newton iterative approach and Morozov discrepancy principle is one of the interesting strategies for determination of regularization parameter having an excellent rate of convergence. Moreover since iterative methods are known to be an attractive alternative to Tikhonov regularization, four gradient descend methods are examined here for regularization of the linear regression problem. It has been found that all of gradient decent methods with an appropriate initial guess of regularized coefficients have an excellent convergence to Tikhonov-Morozov regularization results.

Details

International Journal of Structural Integrity, vol. 7 no. 1
Type: Research Article
ISSN: 1757-9864

Keywords

Open Access
Article
Publication date: 16 January 2020

Magda Joachimiak

In this paper, the Cauchy-type problem for the Laplace equation was solved in the rectangular domain with the use of the Chebyshev polynomials. The purpose of this paper is to…

4380

Abstract

Purpose

In this paper, the Cauchy-type problem for the Laplace equation was solved in the rectangular domain with the use of the Chebyshev polynomials. The purpose of this paper is to present an optimal choice of the regularization parameter for the inverse problem, which allows determining the stable distribution of temperature on one of the boundaries of the rectangle domain with the required accuracy.

Design/methodology/approach

The Cauchy-type problem is ill-posed numerically, therefore, it has been regularized with the use of the modified Tikhonov and Tikhonov–Philips regularization. The influence of the regularization parameter choice on the solution was investigated. To choose the regularization parameter, the Morozov principle, the minimum of energy integral criterion and the L-curve method were applied.

Findings

Numerical examples for the function with singularities outside the domain were solved in this paper. The values of results change significantly within the calculation domain. Next, results of the sought temperature distributions, obtained with the use of different methods of choosing the regularization parameter, were compared. Methods of choosing the regularization parameter were evaluated by the norm Nmax.

Practical implications

Calculation model described in this paper can be applied to determine temperature distribution on the boundary of the heated wall of, for instance, a boiler or a body of the turbine, that is, everywhere the temperature measurement is impossible to be performed on a part of the boundary.

Originality/value

The paper presents a new method for solving the inverse Cauchy problem with the use of the Chebyshev polynomials. The choice of the regularization parameter was analyzed to obtain a solution with the lowest possible sensitivity to input data disturbances.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 30 no. 10
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 3 October 2019

Magda Joachimiak, Michał Ciałkowski and Andrzej Frąckowiak

The purpose of this paper is to present the method for solving the inverse Cauchy-type problem for the Laplace’s equation. Calculations were made for the rectangular domain with…

Abstract

Purpose

The purpose of this paper is to present the method for solving the inverse Cauchy-type problem for the Laplace’s equation. Calculations were made for the rectangular domain with the target temperature on three boundaries and, additionally, on one of the boundaries, the heat flux distribution was selected. The purpose of consideration was to determine the distribution of temperature on a section of the boundary of the investigated domain (boundary Γ1) and find proper method for the problem regularization.

Design/methodology/approach

The solution of the direct and the inverse (Cauchy-type) problems for the Laplace’s equation is presented in the paper. The form of the solution is noted as the linear combination of the Chebyshev polynomials. The collocation method in which collocation points had been determined based on the Chebyshev nodes was applied. To solve the Cauchy problem, the minimum of functional describing differences between the target and the calculated values of temperature and the heat flux on a section of the domain’s boundary was sought. Various types of the inverse problem regularization, based on Tikhonov and Tikhonov–Philips regularizations, were analysed. Regularization parameter α was chosen with the use of the Morozov discrepancy principle.

Findings

Calculations were performed for random disturbances to the heat flux density of up to 0.01, 0.02 and 0.05 of the target value. The quality of obtained results was next estimated by means of the norm. Effect of the disturbance to the heat flux density and the type of regularization on the sought temperature distribution on the boundary Γ1 was investigated. Presented methods of regularization are considerably less sensitive to disturbances to measurement data than to Tikhonov regularization.

Practical implications

Discussed in this paper is an example of solution of the Cauchy problem for the Laplace’s equation in the rectangular domain that can be applied for determination of the temperature distribution on the boundary of the heated element where it is impossible to measure temperature or the measurement is subject to a great error, for instance on the inner wall of the boiler. Authors investigated numerical examples for functions with singularities outside the domain, for which values of gradients change significantly within the calculation domain what corresponds to significant changes in temperature on the wall of the boiler during the fuel combustion.

Originality/value

In this paper, a new method for solving the Cauchy problem for the Laplace’s equation is described. To solve this problem, the Chebyshev polynomials and nodes were used. Various types of regularization of this problem were considered.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 30 no. 3
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 6 December 2022

M.J. Huntul and Mohammad Tamsir

The purpose of this paper is to reconstruct the potential numerically in the fourth-order Rayleigh–Love equation with boundary and nonclassical boundary conditions, from…

Abstract

Purpose

The purpose of this paper is to reconstruct the potential numerically in the fourth-order Rayleigh–Love equation with boundary and nonclassical boundary conditions, from additional measurement.

Design/methodology/approach

Although, the aforesaid inverse identification problem is ill-posed but has a unique solution. The authors use the cubic B-spline (CBS) collocation and Tikhonov regularization techniques to discretize the direct problem and to obtain stable as well as accurate solutions, respectively. The stability, for the discretized system of the direct problem, is also carried out by means of the von Neumann method.

Findings

The acquired results demonstrate that accurate as well as stable solutions for the a(t) are accessed for λ {10–8, 10–7, 10–6, 10–5}, when p {0.01%, 0.1%} for both linear and nonlinear potential coefficient a(t). The stability analysis shows that the discretized system of the direct problem is unconditionally stable.

Research limitations/implications

Since the noisy data are introduced, the investigation and analysis model real circumstances where the practical quantities are naturally infested with noise.

Practical implications

The acquired results demonstrate that accurate as well as stable solutions for the a(t) are accessed for λ {10–8, 10–7, 10–6, 10–5}, when p {0.01%, 0.1%} for both linear and nonlinear potential coefficient a(t). The stability analysis shows that the discretized system of the direct problem is unconditionally stable.

Originality/value

The potential term in the fourth-order Rayleigh–Love equation from additional measurement is reconstructed numerically, for the first time. The technique establishes that accurate, as well as stable solutions are obtained.

Article
Publication date: 21 May 2021

M.J. Huntul, Mohammad Tamsir and Abdullah Ahmadini

The paper aims to numerically solve the inverse problem of determining the time-dependent potential coefficient along with the temperature in a higher-order Boussinesq-Love…

Abstract

Purpose

The paper aims to numerically solve the inverse problem of determining the time-dependent potential coefficient along with the temperature in a higher-order Boussinesq-Love equation (BLE) with initial and Neumann boundary conditions supplemented by boundary data, for the first time.

Design/methodology/approach

From the literature, the authors already know that this inverse problem has a unique solution. However, the problem is still ill-posed by being unstable to noise in the input data. For the numerical realization, the authors apply the generalized finite difference method (GFDM) for solving the BLE along with the Tikhonov regularization to find stable and accurate numerical solutions. The regularized nonlinear minimization is performed using the MATLAB subroutine lsqnonlin. The stability analysis of solution of the BLE is proved using the von Neumann method.

Findings

The present numerical results demonstrate that obtained solutions are stable and accurate.

Practical implications

Since noisy data are inverted, the study models real situations in which practical measurements are inherently contaminated with noise.

Originality/value

The knowledge of this physical property coefficient is very important in various areas of human activity such as seismology, mineral exploration, biology, medicine, quality control of industrial products, etc. The originality lies in the insight gained by performing the numerical simulations of inversion to find the potential co-efficient on time in the BLE from noisy measurement.

Article
Publication date: 24 May 2021

M.J. Huntul and Mohammad Tamsir

The inverse problem of identifying the time-dependent potential coefficient along with the temperature in the fourth-order Boussinesq–Love equation (BLE) with initial and boundary…

Abstract

Purpose

The inverse problem of identifying the time-dependent potential coefficient along with the temperature in the fourth-order Boussinesq–Love equation (BLE) with initial and boundary conditions supplemented by mass measurement is, for the first time, numerically investigated. From the literature, the authors already know that this inverse problem has a unique solution. However, the problem is still ill-posed by being unstable to noise in the input data.

Design/methodology/approach

For the numerical discretization, the authors apply the Crank–Nicolson finite difference method along with the Tikhonov regularization for finding a stable and accurate approximate solution. The resulting nonlinear minimization problem is solved using the MATLAB routine lsqnonlin. Both exact and numerically simulated noisy input data are inverted.

Findings

The present computational results demonstrate that obtained solutions are stable and accurate.

Originality/value

The inverse problem presented in this paper was already showed to be locally uniquely solvable, but no numerical identification has been studied yet. Therefore, the main aim of the present work is to undertake the numerical realization. The von Neumann stability analysis is also discussed.

Article
Publication date: 6 March 2017

Fabio Bozzoli, Luca Cattani, Sara Rainieri, Fermín S.V. Bazán and Leonardo S. Borges

Most of the passive techniques for enhancing heat transfer inside pipes (e.g. rough surfaces, swirl-flow devices and coiled tubes) give origin to an irregular distribution of the…

Abstract

Purpose

Most of the passive techniques for enhancing heat transfer inside pipes (e.g. rough surfaces, swirl-flow devices and coiled tubes) give origin to an irregular distribution of the heat transfer coefficient at the fluid–wall interface along the wall perimeter. This irregular distribution could be critical in some industrial applications, but most of the available research papers, mainly due to the practical difficulty of local measuring heat flux on the internal wall surface of a pipe, present the results only in terms of Nusselt number averaged along the wall circumference. This paper aims to study the application of inverse problem solution techniques, which could overcome this limitation.

Design/methodology/approach

With regard to the estimation of the local convective heat transfer coefficient in coiled tubes, two different inverse heat conduction problem solution techniques were considered and compared both by synthetic and experimental data.

Findings

The paper shows the success of two inverse problem solution techniques in the estimation of the local convective heat transfer coefficient in coiled tubes.

Originality/value

This paper fulfills an identified need because most of the available research papers present the results only in terms of average thermal performance, neglecting local behavior.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 27 no. 3
Type: Research Article
ISSN: 0961-5539

Keywords

Open Access
Article
Publication date: 13 December 2022

Marcelo Colaço, Fabio Bozzoli, Luca Cattani and Luca Pagliarini

The purpose of this paper is to apply the conjugate gradient (CG) method, together with the adjoint operator (AO) to the pulsating heat pipe problem, including some quite…

379

Abstract

Purpose

The purpose of this paper is to apply the conjugate gradient (CG) method, together with the adjoint operator (AO) to the pulsating heat pipe problem, including some quite interesting experimental results. The CG method, together with the AO, was able to estimate the unknown functions more efficiently than the other techniques presented in this paper. The estimation of local heat transfer coefficients, rather than the global ones, in pulsating heat pipes is a relatively new subject and presenting a robust, efficient and self-regularized inverse tool to estimate it, supported also by some experimental results, is the main purpose of this paper. To also increase the visibility and the general use of the paper to the heat transfer community, the authors include, as supplemental material, all numerical and experimental data used in this paper.

Design/methodology/approach

The approach was established on the solution of the inverse heat conduction problem in the wall by using as starting data the temperature measurements on the outer surface. The procedure is based on the CG method with AO. The here proposed approach was first verified adopting synthetic data and then it was validated with real cases regarding pulsating heat pipes.

Findings

An original fast methodology to estimate local convective heat flux is proposed. The procedure has been validated both numerically and experimentally. The procedure has been compared to other classical methods presenting some peculiar benefits.

Practical implications

The approach is suitable for pulsating heat pipes performance evaluation because these devices present a local heat flux distribution characterized by an important variation both in time and in space as a result of the complex flow patterns that are generated in this type of devices.

Originality/value

The procedure here proposed shows these benefits: it affords a general model of the heat conduction problem that is effortlessly customized for the particular case, it can be applied also to large datasets and it presents reduced computational expense.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 33 no. 5
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 30 May 2023

Everton Boos, Fermín S.V. Bazán and Vanda M. Luchesi

This paper aims to reconstruct the spatially varying orthotropic conductivity based on a two-dimensional inverse heat conduction problem described by a partial differential…

23

Abstract

Purpose

This paper aims to reconstruct the spatially varying orthotropic conductivity based on a two-dimensional inverse heat conduction problem described by a partial differential equation (PDE) model with mixed boundary conditions. The proposed discretization uses a highly accurate technique and allows simple implementations. Also, the authors solve the related inverse problem in such a way that smoothness is enforced on the iterations, showing promising results in synthetic examples and real problems with moving heat source.

Design/methodology/approach

The discretization procedure applied to the model for the direct problem uses a pseudospectral collocation strategy in the spatial variables and Crank–Nicolson method for the time-dependent variable. Then, the related inverse problem of recovering the conductivity from temperature measurements is solved by a modified version of Levenberg–Marquardt method (LMM) which uses singular scaling matrices. Problems where data availability is limited are also considered, motivated by a face milling operation problem. Numerical examples are presented to indicate the accuracy and efficiency of the proposed method.

Findings

The paper presents a discretization for the PDEs model aiming on simple implementations and numerical performance. The modified version of LMM introduced using singular scaling matrices shows the capabilities on recovering quantities with precision at a low number of iterations. Numerical results showed good fit between exact and approximate solutions for synthetic noisy data and quite acceptable inverse solutions when experimental data are inverted.

Originality/value

The paper is significant because of the pseudospectral approach, known for its high precision and easy implementation, and usage of singular regularization matrices on LMM iterations, unlike classic implementations of the method, impacting positively on the reconstruction process.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 33 no. 8
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 2 January 2018

Andrzej Frąckowiak and Michał Ciałkowski

This paper aims to present the Cauchy problem for the Laplace’s equation for profiles of gas turbine blades with one and three cooling channels. The distribution of heat transfer…

Abstract

Purpose

This paper aims to present the Cauchy problem for the Laplace’s equation for profiles of gas turbine blades with one and three cooling channels. The distribution of heat transfer coefficient and temperature on the outer boundary of the blade are known. On this basis, the temperature on inner surfaces of the blade (the walls of cooling channels) is determined.

Design/methodology/approach

Such posed inverse problem was solved using the finite element method in the domain of the discrete Fourier transform (DFT).

Findings

Calculations indicate that the regularization in the domain of the DFT enables obtaining a stable solution to the inverse problem. In the example under consideration, problems with reconstruction constant temperature, assumed on the outer boundary of the blade, in the vicinity of the trailing and leading edges occurred.

Originality/value

The application of DFT in connection with regularization is an original achievement presented in this study.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 28 no. 1
Type: Research Article
ISSN: 0961-5539

Keywords

1 – 10 of 160