Все выпуски

В результате численного моделирования впервые определен механизм инициированного кавитацией подъема уровня жидкости в трубках и капиллярах, известного как лазерно-индуцированный светокапиллярный эффект, а также его аналоги — звукокапиллярный и плазмокапиллярный эффекты. Показано, что ключевым условием возникновения подъема жидкости является асимметричное схлопывание относительно крупного одиночного кавитационного пузырька внутри вертикально ориентированной трубки или капилляра. Близость границ (стенки трубки, торца оптоволокона и другие) нарушает сферическую симметрию пузырька при его коллапсе, что приводит к появлению потока жидкости, который замыкается в долгоживущее тороидальное вихревое кольцо (тороидальный вихрь). Вихрь за счет вязкого увлечения окружающей среды генерирует направленное течение жидкости вверх, а также обеспечивает всасывание новой порции жидкости через открытый нижний конец трубки. Результаты моделирования показывают, что характерное время жизни тороидального вихря значительно превышает время стадий роста и схлопывания кавитационного пузырька, который его породил. Показано, что подъем уровня жидкости в трубке начинается не в момент расширения пузырька, а после его полного исчезновения и за счет инерции вихревого движения продолжается на протяжении относительно длительного промежутка времени. Этот результат полностью согласуется с экспериментальными данными, что подтверждает достоверность предложенного механизма. В работе исследована практически значимая конфигурация лазерно-индуцированного светокапиллярного эффекта с использованием оптоволокна. Такая конфигурация открывает широкие перспективы для технических и медицинских приложений, в частности в лазерной хирургии. Исследованный механизм может быть использован для создания кавитационных насосов — эффективных инструментов очистки технических поверхностей и поверхностей ран, где процесс удаления поврежденных тканей и инородных тел в результате температурного воздействия будет сопровождаться удалением продуктов очистки через трубку, что существенно повышает эффективность и безопасность процедуры. Полученные результаты представляют собой первое непротиворечивое объяснение класса кавитационно-индуцированных капиллярных явлений и создают основу для их контролируемого применения в биомедицинских и микрофлюидных технологиях.

For the first time, numerical modeling has determined the mechanism of the initiated-by-cavitation rise of the liquid level in tubes and capillaries, known as the laser-induced optocapillary effect, as well as its analogues — the acoustocapillary and plasmocapillary effects. It is shown that the key condition for the occurrence of the liquid rise is the asymmetric collapse of a single relatively large cavitation bubble inside a vertically oriented tube or capillary. The proximity of boundaries (the tube wall, the fiber optic tip, and others) disrupts the spherical symmetry of the bubble during its collapse, leading to the appearance of a liquid flow that rolls up into a long-lived toroidal vortex ring. Due to viscous entrainment of the surrounding medium, the vortex generates a directed liquid flow upward and also ensures the suction of a new portion of liquid through the open lower end of the tube. The simulation results show that the characteristic lifetime of the toroidal vortex significantly exceeds the duration of the growth and collapse stages of the cavitation bubble that generated it. It is demonstrated that the rise of the liquid level in the tube does not begin at the moment of bubble expansion, but after its complete disappearance, and continues over a relatively long period due to the inertia of the vortex motion. This result is in complete agreement with experimental data, confirming the validity of the proposed mechanism.

The study investigated the practically significant configuration of the laser-induced optocapillary effect using an optical fiber. This configuration opens broad prospects for technical and medical applications, particularly in laser surgery. The investigated mechanism can be used to create cavitation pumps — effective tools for cleaning technical surfaces and wound surfaces, where the process of removing damaged tissue and foreign bodies due to thermal exposure will be accompanied by the removal of debris through the tube, significantly increasing the efficiency and safety of the procedure.

The obtained results represent the first consistent explanation for a class of cavitation-induced capillary phenomena and create a foundation for their controlled application in biomedical and microfluidic technologies.

Для решения уравнений Навье – Стокса в случае несжимаемых течений разработано большое количество методов, наиболее популярными из которых являются методы с коррекцией скорости по алгоритму SIMPLE, аналогом которого является метод расщепления по физическим переменным. Данные методы, разработанные еще в прошлом веке, использовались для решения достаточно простых задач — расчета как стационарных течений, так и нестационарных, в которых границы расчетной области были неподвижны. В настоящее время задачи вычислительной гидродинамики существенно усложнились. Интерес представляют задачи с движением тел в расчетной области, движением контактных границ, кавитацией и задачи с динамической локальной адаптацией расчетной сетки. При этом расчетная сетка меняется, что приводит к нарушению условия дивергентности скорости на ней. Поскольку дивергентные скорости используются не только для уравнений Навье – Стокса, но и для всех остальных уравнений математической модели движения жидкости — моделей турбулентности, массопереноса и сохранения энергии, нарушение этого условия ведет к численным ошибкам и, зачастую, к расхождению вычислительного алгоритма.

В статье представлен неявный метод расщепления по физическим переменным, который использует дивергентные скорости с данного шага по времени для решения несжимаемых уравнений Навье – Стокса. Метод разработан для расчета течений при наличии подвижных и контактных границ, моделируемых в постановке Эйлера. Метод позволяет проводить расчеты с шагом интегрирования, на порядки превышающем явный шаг по времени (число Куранта – Фридрихcа – Леви $CFL\gg1$). В данной статье представлен вариант метода для несжимаемых течений. Вариант метода, позволяющий рассчитывать движение жидкости и газа при любых числах Маха, будет опубликован в ближайшее время. Метод для полностью сжимаемых течений реализован в программном комплексе FlowVision.

В статье приводятся результаты численного решения классической задачи обтекания кругового цилиндра при малых числах Рейнольдса ($50<Re<140$), при которых ламинарное обтекание цилиндра становиться нестационарным и образуется дорожка Кармана. Показано хорошее совпадение расчетов с экспериментальными данными, опубликованными в классических работах Ван-Дайка и Танеды.

A large number of methods have been developed to solve the Navier – Stokes equations in the case of incompressible flows, the most popular of which are methods with velocity correction by the SIMPLE algorithm and its analogue — the method of splitting by physical variables. These methods, developed more than 40 years ago, were used to solve rather simple problems — simulating both stationary flows and non-stationary flows, in which the boundaries of the calculation domain were stationary. At present, the problems of computational fluid dynamics have become significantly more complicated. CFD problems are involving the motion of bodies in the computational domain, the motion of contact boundaries, cavitation and tasks with dynamic local adaptation of the computational mesh. In this case the computational mesh changes resulting in violation of the velocity divergence condition on it. Since divergent velocities are used not only for Navier – Stokes equations, but also for all other equations of the mathematical model of fluid motion — turbulence, mass transfer and energy conservation models, violation of this condition leads to numerical errors and, often, to undivergence of the computational algorithm.

This article presents an implicit method of splitting by physical variables that uses divergent velocities from a given time step to solve the incompressible Navier – Stokes equations. The method is developed to simulate flows in the case of movable and contact boundaries treated in the Euler paradigm. The method allows to perform computations with the integration step exceeding the explicit time step by orders of magnitude (Courant – Friedrichs – Levy number $CFL\gg1$). This article presents a variant of the method for incompressible flows. A variant of the method that allows to calculate the motion of liquid and gas at any Mach numbers will be published shortly. The method for fully compressible flows is implemented in the software package FlowVision.

Numerical simulating classical fluid flow around circular cylinder at low Reynolds numbers ($50 < Re < 140$), when laminar flow is unsteady and the Karman vortex street is formed, are presented in the article. Good agreement of calculations with the experimental data published in the classical works of Van Dyke and Taneda is demonstrated.