scholarly journals Numerical Algorithms for Diffusion Coefficient Identification in Problems of Tissue Engineering

2016 ◽  
Vol 11 (2) ◽  
pp. 426-444 ◽  
Author(s):  
А.В. Пененко ◽  
A.V. Penenko
Keyword(s):  
Inverse Problem ◽  
Process Model ◽  
Numerical Algorithms ◽  
Adjoint Equations ◽  
Problem Solution ◽  
Type Method ◽  
Tomographic Images ◽  
Gradient Based ◽  
Pseudo Inverse ◽  
Sensitivity Operator

Identification algorithms of diffusion coefficients in a specimen with tomographic images of the solution penetration dynamics are considered. With the sensitivity operator, built on the basis of adjoint equations for diffusion process model, the corresponding coefficient inverse problem is reduced to the quasilinear operator equation which is then solved by the Newton-type method with successive evaluation of r-pseudo inverse operators of increasing dimensionality. The efficiency of the constructed algorithm is tested in numerical experiments. For comparison, a gradient-based algorithm for the inverse problem solution is considered.

scholarly journals The inverse coefficient problem of heat transfer in layered nanostructures

2017 ◽  
Vol 20 (3) ◽  
pp. 213-219
Author(s):  
K. K. Abgarian ◽  
R. G. Noskov ◽  
D. L. Reviznikov
Keyword(s):  
Heat Transfer ◽  
Inverse Problem ◽  
Thermal Resistance ◽  
Transfer Problem ◽  
Dominant Role ◽  
Rapid Development ◽  
Problem Solution ◽  
Microwave Range ◽  
Contact Thermal Resistance ◽  
Resistance Coefficients

The rapid development of electronics leads to the creation and use of electronic components of small dimensions, including nanoelements of complex, layered structure. The search for effective methods for cooling electronic systems dictates the need for the development of methods for the numerical analysis of heat transfer in nanostructures. A characteristic feature of energy transfer in such systems is the dominant role of contact thermal resistance at interlayer interfaces. Since the contact resistance depends on a number of factors associated with the technology of heterostructures manufacturing, it is of great importance to determine the corresponding coefficients from the results of temperature measurements.The purpose of this paper is to evaluate the possibility of reconstructing the thermal resistance coefficients at the interfaces between layers by solving the inverse problem of heat transfer.The complex of algorithms includes two major blocks — a block for solving the direct heat transfer problem in a layered nanostructure and an optimization block for solving the inverse problem. The direct problem was formulated in an algebraic (finite difference) form under the assumption of a constant temperature within each layer, which is due to the small thickness of the layers. The inverse problem was solved in the extreme formulation, the optimization was carried out using zero-order methods that do not require the calculation of the derivatives of the optimized function. As a basic optimization algorithm, the Nelder—Mead method was used in combination with random restarts to search for a global minimum.The results of the identification of the contact thermal resistance coefficients obtained in the framework of a quasi-real experiment are presented. The accuracy of the identification problem solution is estimated as a function of the number of layers in the heterostructure and the «measurements» error.The obtained results are planned to be used in the new technique of multiscale modeling of thermal regimes of the electronic component base of the microwave range, when identifying the coefficients of thermal conductivity of heterostructure.

scholarly journals Adjoint Numerical Method for a Multiphysical Inverse Problem of Two-Phase Well Testing in Petroleum Reservoirs

2021 ◽  
Vol 2090 (1) ◽  
pp. 012139
Author(s):  
OA Shishkina ◽  
I M Indrupskiy
Keyword(s):  
Inverse Problem ◽  
Inverse Problems ◽  
Objective Function ◽  
Petroleum Reservoirs ◽  
Forward Problem ◽  
Adjoint Equations ◽  
Well Testing ◽  
Adjoint Methods ◽  
Two Phase ◽  
Gradient Calculation

Abstract Inverse problem solution is an integral part of data interpretation for well testing in petroleum reservoirs. In case of two-phase well tests with water injection, forward problem is based on the multiphase flow model in porous media and solved numerically. The inverse problem is based on a misfit or likelihood objective function. Adjoint methods have proved robust and efficient for gradient calculation of the objective function in this type of problems. However, if time-lapse electrical resistivity measurements during the well test are included in the objective function, both the forward and inverse problems become multiphysical, and straightforward application of the adjoint method is problematic. In this paper we present a novel adjoint algorithm for the inverse problems considered. It takes into account the structure of cross dependencies between flow and electrical equations and variables, as well as specifics of the equations (mixed parabolic-hyperbolic for flow and elliptic for electricity), numerical discretizations and grids, and measurements in the inverse problem. Decomposition is proposed for the adjoint problem which makes possible step-wise solution of the electric adjoint equations, like in the forward problem, after which a cross-term is computed and added to the right-hand side of the flow adjoint equations at this timestep. The overall procedure provides accurate gradient calculation for the multiphysical objective function while preserving robustness and efficiency of the adjoint methods. Example cases of the adjoint gradient calculation are presented and compared to straightforward difference-based gradient calculation in terms of accuracy and efficiency.

scholarly journals Estimation of the Directions of Alpha Rhythm Elementary Sources Using the Method of Human Brain Functional Tomography Based On the Magnetic Encephalography Data

2018 ◽  
Vol 13 (2) ◽  
pp. 426-436 ◽  
Author(s):  
M.N. Ustinin ◽  
S.D. Rykunov ◽  
A.I. Boyko ◽  
O.A. Maslova ◽  
K.D. Walton ◽  
...  
Keyword(s):  
Inverse Problem ◽  
Time Series ◽  
Human Brain ◽  
Alpha Rhythm ◽  
Brain Activity ◽  
Problem Solution ◽  
New York University ◽  
Magnetic Encephalography ◽  
The Brain ◽  
Calculation Grid

New method for the magnetic encephalography data analysis was proposed. The method transforms multichannel time series into the spatial structure of the human brain activity. In this paper we further develop this method to determine the dominant direction of the electrical sources of brain activity at each node of the calculation grid. We have considered the experimental data, obtained with three 275-channel magnetic encephalographs in New York University, McGill University and Montreal University. The human alpha rhythm phenomenon was selected as a model object. Magnetic encephalograms of the brain spontaneous activity were registered for 5-7 minutes in magnetically shielded room. Detailed multichannel spectra were obtained by the Fourier transform of the whole time series. For all spectral components, the inverse problem was solved in elementary current dipole model and the functional structure of the brain activity was calculated in the frequency band 8-12 Hz. In order to estimate the local activity direction, at the each node of calculation grid the vector of the inverse problem solution was selected, having the maximal spectral power. So, the 3D-map of the brain activity vector field was produced – the directional functional tomogram. Such maps were generated for 15 subjects and some common patterns were revealed in the directions of the alpha rhythm elementary sources. The proposed method can be used to study the local properties of the brain activity in any spectral band and in any brain compartment.

scholarly journals On investigation of the inverse problem for a parabolic integro-differential equation with a variable coefficient of thermal conductivity

2020 ◽  
Vol 30 (4) ◽  
pp. 572-584
Author(s):  
D.K. Durdiev ◽  
J.Z. Nuriddinov
Keyword(s):  
Thermal Conductivity ◽  
Inverse Problem ◽  
Variable Coefficient ◽  
Direct Problem ◽  
Local Existence ◽  
Problem Solution ◽  
Volterra Type ◽  
Integral Term ◽  
Local Existence And Uniqueness ◽  
Coefficient Of Thermal Conductivity

The inverse problem of determining a multidimensional kernel of an integral term depending on a time variable $t$ and $ (n-1)$-dimensional spatial variable $x'=\left(x_1,\ldots, x_ {n-1}\right)$ in the $n$-dimensional heat equation with a variable coefficient of thermal conductivity is investigated. The direct problem is the Cauchy problem for this equation. The integral term has the time convolution form of kernel and direct problem solution. As additional information for solving the inverse problem, the solution of the direct problem on the hyperplane $x_n = 0$ is given. At the beginning, the properties of the solution to the direct problem are studied. For this, the problem is reduced to solving an integral equation of the second kind of Volterra-type and the method of successive approximations is applied to it. Further the stated inverse problem is reduced to two auxiliary problems, in the second one of them an unknown kernel is included in an additional condition outside integral. Then the auxiliary problems are replaced by an equivalent closed system of Volterra-type integral equations with respect to unknown functions. Applying the method of contraction mappings to this system in the Hölder class of functions, we prove the main result of the article, which is a local existence and uniqueness theorem of the inverse problem solution.

scholarly journals Numerical Algorithms for Simulation of a Fluid-Filed Fracture Evolution in a Poroelastic Medium

2021 ◽  
pp. 24-35
Author(s):  
V. E Borisov ◽  
A. V Ivanov ◽  
B. V Kritsky ◽  
E. B Savenkov
Keyword(s):  
Finite Element ◽  
Fluid Flow ◽  
Geometric Model ◽  
Numerical Algorithms ◽  
Problem Solution ◽  
Poroelastic Medium ◽  
Heterogeneous Distribution ◽  
Coupled Problem ◽  
Fracture Evolution ◽  
Point Projection

The paper deals with the computational framework for the numerical simulation of the three dimensional fluid-filled fracture evolution in a poroelastic medium. The model consists of several groups of equations including the Biot poroelastic model to describe a bulk medium behavior, Reynold’s lubrication equations to describe a flow inside fracture and corresponding bulk/fracture interface conditions. The geometric model of the fracture assumes that it is described as an arbitrary sufficiently smooth surface with a boundary. Main attention is paid to describing numerical algorithms for particular problems (poroelasticity, fracture fluid flow, fracture evolution) as well as an algorithm for the coupled problem solution. An implicit fracture mid-surface representation approach based on the closest point projection operator is a particular feature of the proposed algorithms. Such a representation is used to describe the fracture mid-surface in the poroelastic solver, Reynold’s lubrication equation solver and for simulation of fracture evolutions. The poroelastic solver is based on a special variant of X-FEM algorithms, which uses the closest point representation of the fracture. To solve Reynold’s lubrication equations, which model the fluid flow in fracture, a finite element version of the closet point projection method for PDEs surface is used. As a result, the algorithm for the coupled problem is purely Eulerian and uses the same finite element mesh to solve equations defined in the bulk and on the fracture mid-surface. Finally, we present results of the numerical simulations which demonstrate possibilities of the proposed numerical techniques, in particular, a problem in a media with a heterogeneous distribution of transport, elastic and toughness properties.

Sign in / Sign up

Export Citation Format

Share Document