An efficient implicit spectral element method for time-dependent nonlinear diffusion equations by evaluating integrals at one quadrature point

The influence of the temperature dependence of the thermal conductivity and heat transfer coefficients, as well as of the emissivity of a dielectric film and of an insulating base on formation of a stationary temperature field of the dielectric–base system exposed to irradiation by a continuous ion beam in air, has been studied theoretically. Analytical expressions have been obtained for the temperature fields of a sample and base, as well as a system of interrelated nonlinear algebraic equations for characteristic temperatures. With the aid of numerical solution of a system of nonlinear algebraic equations, the dependences of the temperature of the irradiated surface of the sample–base system and of the rear side of this system on the incident ion beam parameters have been found. It has been established that these dependences are nonlinear. It has been found that heating of a sample, including its irradiated surface, decreases substantially with increase in the thermal conductivity of the base.

Modified bounded homotopies to enable a narrow bounding zone

CANM, a program for numerical solution of a system of nonlinear equations using the continuous analog of Newton's method

Systems of nonlinear algebraic equations have a special place in solving modeling and calculation problems in various areas of construction Sciences. These are tasks for calculating building structures, calculating distributed loads, mass and heat transfer tasks, calculating engineering networks, and many other computational and optimization tasks related to the construction industry. As is known, there are a sufficient number of fundamentally different methods for solving systems of nonlinear algebraic equations, which is related to the specifics of the systems, their degree of nonlinearity and conditionality. Non-linear ill-conditioned systems are particularly difficult; an erroneous choice of the method for solving such a system can lead to significantly distorted or simply incorrect results, which is absolutely unacceptable when performing the calculated stages of solving construction problems. In this paper, we propose to use a modification of the well – known Levenberg – Marquardt method developed by the authors, based on the regularization of the Jacobi matrix used in the classical Newton method, to solve complex, ill-conditioned systems of nonlinear equations. The modified method allows solving poorly defined systems and can significantly reduce the amount of calculations needed to ensure the required accuracy by reducing the number of iterative procedures. The paper presents a detailed description of the algorithm, given the solution of the model problem and, as an example the task of modeling and optimization of the port berthing facilities – mooring wall of massive masonry with various add-ons. When composing mathematical equations, the method of limit States was used, which is generally accepted for calculating such port hydraulic structures. The choice of the optimum scheme was carried out by minimizing the cost of structure during limit conditions the reliability of structures when using existing construction materials. The most significant geometric dimensions of the structure were selected using independent parameters to be optimized.

Aim. To develop an effi cient method of building a numerical model for the process of moisture transfer under trickle irrigation, with the mathematical modeling of the method involving the system of differential equations in partial derivatives of Klute-Richards, and to perform computing experiments regarding fl at-vertical profi le moisture transfer with point sources. Methods. The mathematical apparatus of the theory of differential schemes of solving differential equations in partial derivatives, and Newton’s method of iterative approximate solving of non-linear equations. Results. A stable differential two-step symmetrized algorithm (TS-algorithm) along with the corresponding scheme of the method of numerical solution for initially-boundary task for Richards’ equation was created. The method was realized in the form of a computer program in C++ language, the computing experiments were performed with three deeper points, the humidity zones for volume moisture and potential were obtained. Conclusions. The numerical method was suggested, ensuring the effi cient solution to Richards’ non-linear equation in conditions of several deep point sources. The algorithm structure allows reducing the system of non-linear algebraic equations with many unknowns to solving independent non-linear equations with one unknown. The presented method may easily be expanded for three-dimensional cases. The results of computing experiments are in agreement with natural observations.

: Complex physical phenomena involving chemically reacting systems or the transport of heat or fluids are often modeled by coupled systems of time-dependent, nonlinear partial differential equations. The difficulties in understanding the stability of the differential equation systems and in designing efficient, accurate numerical methods for their solution are widely recognized and were the focus of this research. We have worked on four general aspects of the analysis and numerical approximation of systems of partial differential equations. These areas of research are: (1) modeling aspects and stability analysis for nonlinear time-dependent partial differential equations; (2) use and analysis of finite element or finite difference methods to discretize coupled systems of nonlinear differential equations; (3) development of adaptive or local grid refinement capabilities to resolve local phenomena in large-scale applications; and (4) development of data structures, preconditioners, and efficient solution algorithms for large-scale problems on new computer architectures. Emphasis has been placed upon multiphase or multicomponent, transport-dominated flow processes with dynamic local phenomena. The research also involved a mix of analysis, algorithm development, and large-scale computation using newer computer architectures.

A nonlinear loaded differential equation with a parameter on a finite interval is studied. The interval is partitioned by the load points, at which the values of the solution to the equation are set as additional parameters. A nonlinear boundary value problem for the considered equation is reduced to a nonlinear multipoint boundary value problem for the system of nonlinear ordinary differential equations with parameters. For fixed parameters, we obtain the Cauchy problems for ordinary differential equations on the subintervals. Substituting the values of the solutions to these problems into the boundary condition and continuity conditions at the partition points, we compose a system of nonlinear algebraic equations in parameters. A method of solving the boundary value problem with a parameter is proposed. The method is based on finding the solution to the system of nonlinear algebraic equations composed.

A hybrid approximation method for solving Hutchinson’s equation

97/04852 Coal pulsating combustion model: Xu, W. Ranshao Kexue Yu Jishu., 1996, 2, (3), 260–269. (In Chinese)

98/03268 Reopening of pulverized coal injection in blast furnace after a pause

99/00738 On the combustion of heat pump with solar steam turbine

99/00035 ESR study of the macromolecular structure of copper-porphyrin-labeled coal

The global asymptotic nonlinear behavior of some standard iterative procedures in solving nonlinear systems of algebraic equations arising from four implicit linear multistep methods (LMMs) in discretizing 2 x 2 systems of first-order autonomous nonlinear ordinary differential equations is analyzed using the theory of dynamical systems. With the aid of parallel Connection Machines (CM-2 and CM-5), the associated bifurcation diagrams as a function of the time step, and the complex behavior of the associated 'numerical basins of attraction' of these iterative implicit schemes are revealed and compared. Studies showed that all of the four implicit LMMs exhibit a drastic distortion and segmentation but less shrinkage of the basin of attraction of the true solution than standard explicit methods. The numerical basins of attraction of a noniterative implicit procedure mimic more closely the basins of attraction of the differential equations than the iterative implicit procedures for the four implicit LMMs.

The solution of a system of nonlinear algebraic or transcendental equations is probably one of the most difficult problems in numerical analysis. In general, it is almost impossible to find the solution, unless an approximate value is already known. Linear equations occur more frequently in numerical calculations, simply because it is difficult and sometimes impossible to solve the corresponding nonlinear problem. Most of the equations governing the physical world are nonlinear, but very often for simplicity and tractability linear models are considered. Solution of boundary value problems for differential equations also lead to a nonlinear equation or a system of nonlinear equations. Solution of a system of nonlinear equations also arises in numerical analysis, while deriving certain formulae. As we have seen in Chapter 6, the Gaussian quadrature formulae can be obtained by solving a system of nonlinear equations.

The problem of numerical solution of systems of nonlinear algebraic equations with even powers is studied in this paper. Nondifferentiable (nonsmooth) unconstrained minimization problems are associated with the system under consideration, with objective functions, based on discrete l 1 - and l ∞ - norm, respectively, by using these two norms as proximity criteria, and the considered system is solved by minimizing the sum of absolute values of residuals, in the case of l 1 - norm, or by minimizing the maximal residual, in the case of l ∞ - norm, respectively. Also, a “differentiable” unconstrained minimization problem (least squares problem) is associated with the considered system of nonlinear equations, with objective function, based on the discrete l 2 -norm, and the considered system is solved by minimizing the sum of squares of the residuals. Subgradients of the objective functions of unconstrained minimization problems are calculated, and a subgradient method for solving the associated problems is presented. The iterative gradient method for solving the corresponding “differentiable” discrete least squares problem, based on l 2 - norm, is also outlined, and gradients of the objective function of the corresponding “differentiable” unconstrained minimization problems are calculated.