The problem of determining the thermal conductivity coefficient that depends on temperature is studied. The consideration is based on the initial-boundary value problem for the one-dimensional unsteady heat equation. ...
详细信息
The problem of determining the thermal conductivity coefficient that depends on temperature is studied. The consideration is based on the initial-boundary value problem for the one-dimensional unsteady heat equation. The mean-root-square deviation of the temperature distribution field and the heat flux from the experimental data on the left boundary of the domain is used as the objective functional. An analytical expression for the gradient of the objective functional is obtained. An algorithm for the numerical solution of the problem based on the modern fast automatic differentiation technique is proposed. Examples of solving the problem are discussed.
An efficient numerical method for solution of boundary value problems with additional condition is presented. The approach is based on the shooting method but the procedure of seeking "the proper shot" allow...
详细信息
An efficient numerical method for solution of boundary value problems with additional condition is presented. The approach is based on the shooting method but the procedure of seeking "the proper shot" allows one to satisfy "additional" boundary conditions. General considerations are illustrated by a real example. The computational example concerns the "dead zone" regime for the non-linear diffusion-reaction equation in heterogeneous catalysis. Accuracy and efficiency of the presented method confirm results obtained for a wide range of changes of process parameters, including the vicinity of a critical point. Calculations were performed with the use of the Maple (R) program.
The aim of this paper is to develop an effective finite volume method for numerical simulation of the adiabatic shear bands (ASB) formation processes. A formation of ASB happens at high-speed shear strains of ductile ...
详细信息
The aim of this paper is to develop an effective finite volume method for numerical simulation of the adiabatic shear bands (ASB) formation processes. A formation of ASB happens at high-speed shear strains of ductile materials. A numerical simulation of such problems using Lagrangian approach is associated with some problems, the main one of which is a mesh distortion at large deformations. We use Eulerian approach to describe a motion of the non-linear elasto-plastic material. More specifically, we consider a modification of a well-known hypoelastic Wilkins model. In this paper we suggest a numerical method for modeling of high-speed shear deformations on two-dimensional meshes. The method is verified on the three test problems suggested by other authors. (c) 2021 Elsevier B.V. All rights reserved.
The paper presents a numerical method for determining the contact area in three-dimensional elastostatic normal contact without friction. The method makes use of the theorem developed by Barber, the contact area is th...
详细信息
The paper presents a numerical method for determining the contact area in three-dimensional elastostatic normal contact without friction. The method makes use of the theorem developed by Barber, the contact area is that over which the total indentation force achieves its maximum value. By approximating the punch by linear interpolation, the analytical expression for the indentation force is derived by virtue of the reciprocal theorem. The physical meaning of the parameter which determines the contact boundary is discussed, and its feasible range corresponding to the contact area is found. Then, the numerical algorithm for determining the parameter is developed and applied to solve several normal contact problems. The results show that the proposed numerical method possesses a good property on accuracy and convergency.
Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust num...
详细信息
Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of dis-cretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).
We show that a constant-potential time-independent Schrodinger equation with Dirichlet boundary data can be reformulated as a Laplace equation with Dirichlet boundary data. With this reformulation, which we call the D...
详细信息
We show that a constant-potential time-independent Schrodinger equation with Dirichlet boundary data can be reformulated as a Laplace equation with Dirichlet boundary data. With this reformulation, which we call the Duffin correspondence, we provide a classical Walk On Spheres (WOS) algorithm for Monte Carlo simulation of the solutions of the boundary value problem. We compare the obtained Duffin WOS algorithm with existing modified WOS algorithms.
Results obtained by the authors in solving inverse coefficient problems are overviewed. The inverse problem under consideration is to determine a temperature-dependent thermal conductivity coefficient from experimenta...
详细信息
Results obtained by the authors in solving inverse coefficient problems are overviewed. The inverse problem under consideration is to determine a temperature-dependent thermal conductivity coefficient from experimental observations of the temperature field in the studied substance and (or) the heat flux on the surface of the object. The study is based on the Dirichlet boundary value problem for the nonstationary heat equation stated in the general -dimensional formulation. For this general case, an analytical expression for the cost functional gradient is obtained. The features of solving the inverse problem and the difficulties encountered in the solution process are discussed.
In this article we investigate the hitting time of some given boundaries for Bessel processes. The main motivation comes from mathematical finance when dealing with volatility models, but the results can also be used ...
详细信息
In this article we investigate the hitting time of some given boundaries for Bessel processes. The main motivation comes from mathematical finance when dealing with volatility models, but the results can also be used in optimal control problems. The aim here is to construct a new and efficient algorithm in order to approach this hitting time. As an application we will consider the hitting time of a given level for the Cox-Ingersoll-Ross process. The main tools we use are on one side, an adaptation of the method of images to this particular situation and on the other side, the connection that exists between Cox-Ingersoll-Ross processes and Bessel processes.
The accuracy of stability assessment provided by the transient energy function (TEF) method depends on the determination of the controlling unstable equilibrium point (UEP), The technique that currently determines the...
详细信息
The accuracy of stability assessment provided by the transient energy function (TEF) method depends on the determination of the controlling unstable equilibrium point (UEP), The technique that currently determines the controlling UEP in the TEF method is based on the so-called exit point method and has also been recently labeled the BCU method, The exit point method consists of two basic steps, First, the exit point is approximated by the point theta(egsa), where the first maximum of the potential energy along the fault-on trajectory is encountered, Second, the minimum gradient point theta(mgp) along the trajectory from theta(egsa) is computed, The controlling UEP is then obtained by solving a system of nonlinear algebraic equations with theta(mgp) as an initial guess, It has been observed that this method lacks robustness in the sense that the following two problems may occur, 1) There may be no detection of the minimum gradient point theta(mgp) and hence, no determination of the controlling UEP, 2) if theta(mgp) is found, then based on the definition of the controlling UEP, it may not be in the domain of convergence of the controlling UEP for the particular solving algorithm used, Hence, another equilibrium point, possibly a stable equilibrium paint, not the controlling UEP will be located, This results in a flawed transient stability assessment, The result of this research has been the development of a new numerical technique for determining the controlling UEP, With an initial starting point that is close to the exit point this technique efficiently produces a sequence of points, An analytical foundation for this method is given which shows that under certain assumptions this sequence will converge to the controlling UEP, Hence this new technique exhibits a substantial improvement over the exit point method because of the following reasons: (1) the technique does not attempt to detect the point theta(mgp), (2) the technique can produce a point that is close to t
A previously formulated new approach to the consideration of systems of quasilinear hyperbolic equations on the basis of variational principles is described in more detail in the case of special systems of three equat...
详细信息
A previously formulated new approach to the consideration of systems of quasilinear hyperbolic equations on the basis of variational principles is described in more detail in the case of special systems of three equations. It is shown that each field of characteristics can be represented as a solution of a variational problem. Moreover, the Rankine-Hugoniot relations at the corner points of the characteristics or at the intersections of the characteristics of a single family hold automatically. In the simplest case of the Hopf equation, a numerical algorithm is constructed on the basis of a variational principle.
暂无评论