The results of mathematical modeling of stationary physical processes in the electron-hole plasma of the active region (i-region) of integrated p - i - n-structures are presented. The mathematical model is written in the framework of the hydrodynamic thermal approximation, taking into account the phenomenological data on the effect on the dynamic characteristics of charge carriers of heating of the electron-hole plasma as a result of the release of Joule heat in the volume of the i-th region and the release of recombination energy. The model is based on a nonlinear boundary value problem on a given spatial domain with curvilinear sections of the boundary for the system of equations for the continuity of the current of charge carriers, Poisson, and thermal conductivity. The statement of the problem contains a naturally formed small parameter, which made it possible to use asymptotic methods for its analytical-numerical solution. A model nonlinear boundary value problem with a small parameter is reduced to a sequence of linear boundary value problems by the methods of perturbation theory, and the physical domain of the problem with curvilinear sections of the boundary is reduced to the canonical form by the method of conformal mappings. Stationary distributions of charge carrier concentrations and the corresponding temperature field in the active region of p - i - n-structures are obtained in the form of asymptotic series in powers of a small parameter. The process of refining solutions is iterative, with the alternate fixation of unknown tasks at different stages of the iterative process. The asymptotic series describing the behavior of the plasma concentration and potential in the region under study, in contrast to the classical ones, contain boundary layer corrections. It was found that boundary functions play a key role in describing the electrostatic plasma field. The proposed approach to solving the corresponding nonlinear problem can significantly save computing resources.