В рассматривается задача Коши для линейного параболического уравнения второго порядка: . Матрица коэффициентов a(x) измерима, симметрична и 1-периодична по каждой переменной . Задача моделирует диффузию в неоднородной среде, имеющей периодическую структуру. Решение u(x, t) можно интерпретировать, например, как плотность распределения некоторой случайной величины в момент времени t и тогда подразумевается, что начальное распределение f(x) - неотрицательная функция, интеграл от которой по равен 1. Известно, что при большом значении времени t решение u(x,t) близко к решению задачи Коши: с постоянной матрицей диффузии . Иными словами, при больших t диффузия в периодической среде описывается эффективно через диффузию в однородной среде, которой соответствует постоянная матрица диффузии , называемая матрицей эффективной диффузии. Недавно была доказана оценка погрешности приближения функции u(x,t) решением в норме лебегова -пространства по сечению t=const для любого . Эта оценка порядка при и имеет операторный тип. В настоящей работе найдено приближение для решения исходной задачи с оценкой погрешности того же типа, но порядка при . Это приближение оказывается суммой нулевого приближения и некоторого корректора. Обоснование проведено при дополнительном условии липшицевости матрицы диффузии a(x). Описанный выше результат используется для построения аппроксимации решения уравнения диффузии с быстро осциллирующими ε-периодическими коэффициентами в норме - пространства на сечении t = 1. Погрешность аппроксимации имеет оценку порядка при ε → 0. Результаты работы можно использовать в разных областях, например, при расчете теплового потока в композиционной среде с мелкоячеистой периодической структурой или расчете плотности популяции бактерий в периодической питательной среде.
The Cauchy problem for a linear second order parabolic equation with 1-periodic measurable coefficients is considered Rd, d ≥ 2. The problem models diffusion in a nonhomogeneous periodic medium. The appropriate diffusion operator A is self-adjoint in L2 (Rd). The large-time behaviour of the solution of the Cauchy problem is of main interest for us. To this end, we study first of all the fundamental solution, in other words, the kernel of the exponential exp(-tA), more exactly, its large-time asymptotics. Its approximation is found with integral error estimate of order O(t-1), as time t tends to +∞. To construct this approximation and carry out its justification we use (i) the known fundamental solution to the homogenized diffusion equation (having constant coefficients); (ii) solutions of so-called auxiliary problems on a periodicity cell, which are formulated in a recurrent way. We substantiate this approximation under additional regularity condition on the diffusion matrix a(x): it should be Lipschitz continuous. The results of asymptotic behaviour of the fundamental solution are applied to obtain an approximation of order O(t-1) for the exponential exp(-tA) in operator Lp-norms, on the section t=const as t tends to +∞. There are also some corollaries of these results to operator estimates for a similar exponential exp(-tAε), Aε being a diffusion operator with quickly oscillating ε-periodic coefficients, as tends to zero. This exponential corresponds to the Cauchy problem: modelling diffusion in a strongly nonhomogeneous ε-periodic medium. Here ε is a small parameter. We construct approximations of order O(ε2) for the exponential exp(-tAε) in operator Lp - norm on the section t=const for arbitrary finite fixed t (say, t=1). This approximation is a sum of the exponential exp(-tA0) with a homogenized operator A0 and some correcting operator. The results have a broad range of applications, e.g., for computing heat flow in a periodic composite medium with a small periodicity cell or for bacterial density estimation in a periodic culture medium.
The authors declare that there are no conflicts of interest present.