Collaborative Research: Effective Numerical Schemes for Fundamental Problems Related to Incompressible Fluids

该项目旨在开发一种新的计算方法来研究浮力驱动的流体和由于导电而产生的湍流的稳定性，这种湍流被称为磁流体动力学(MHD)湍流。通过准确模拟这些现象，这项研究将为改进龙卷风、天文事件、北极光和太阳耀斑等现象以及等离子体和液态金属导电流体的建模和预测提供见解。新的计算方法将成为科学计算界的宝贵工具。研究生，包括那些来自代表性不足的群体的学生，将接受理论和计算机领域的培训。这项研究还将吸引本科生和K-12学生参与，使当地学校和社区受益。该项目旨在研究求解不可压缩的Navier-Stokes方程的数值方法，提供散度自由速度、可证明的稳定性、高雷诺数下的稳健性和长时间计算的高效率。本研究将致力于改进经典投影法的不足，目标如下。首先，该方法将达到无发散的条件。该功能对于精确模拟浮力驱动的流体至关重要。这些流体中的温度将被速度场传输、重新排列和分层，而无散度条件的实施将更准确地描述演化和最终状态。其次，该算法将能够以高精度处理任意大的雷诺数。对高雷诺数的三维N-S流动进行数值模拟，可以得到最优的收敛，避免非物理振荡。第三，将建立该格式的稳定性和误差估计，其中结果将独立于雷诺数。研究人员将应用这种数值方法来模拟流体静力平衡附近的浮力驱动流体和背景磁场附近的导电流体。在第一个问题中，将求解一个控制流体静力平衡附近摄动的Boussinesq-Navier-Stokes系统。Boussinesq系统将浮力作用下的Navier-Stokes方程与温度输运方程相耦合。数值方法还将扩展到模拟垂直粘性远小于水平粘性的各向异性流动。在对大气和海洋中发生的Ekman层中的湍流进行建模时，出现了这个问题。在第二个问题中，将计算和分析导电流体在引导磁场下的稳定性和长期行为。当与严格的分析相辅相成时，这些模拟将有助于加速解决Boussinesq和MHD方程上的几个公开的稳定性问题。该项目由计算数学计划和既定的激励竞争研究计划(EPSCoR)联合资助。该奖项反映了NSF的法定使命，并通过使用基金会的智力优势和更广泛的影响审查标准进行评估，被认为值得支持。

This project aims to develop a novel computational method to investigate the stability of buoyancy-driven fluids and turbulent flows due to electrical conduction, known as magnetohydrodynamic (MHD) turbulence. By accurately simulating these phenomena, the research will provide insights into improving modeling and prediction of extreme weather events such as tornados, astronomical occurrences, phenomena like Northern lights and solar flares, and electrically conducting fluid of plasma and liquid metals. The new computational method will be a valuable tool for the scientific computing community. Graduate students, including those from underrepresented groups, will be trained in both theoretical and computer fields. The research will also engage undergraduates and K-12 students, benefiting local schools and communities.The project aims to investigate numerical methods to solve the incompressible Navier-Stokes equations, delivering the divergence free velocity, provable stability, robustness in high Reynolds numbers, and high efficiency in long-time computations. The research will focus on improving the shortcomings of the classical projection method with the following goals. First, the method will achieve the divergence-free condition. This feature is critical in the accurate simulation of buoyancy-driven fluids. The temperature in these fluids will be transported, rearranged, and stratified by the velocity field, and the enforcement of the divergence-free condition will give a more accurate account of the evolution and the eventual states. Second, this algorithm will be able to handle arbitrarily large Reynolds numbers with high accuracy. The simulations of the 3D Navier-Stokes flows with a high Reynolds number will result in optimal convergence and avoid nonphysical oscillations. Third, stability and error estimates will be established for this scheme, where the results will be independent of the Reynolds numbers. The investigators will apply this numerical method to simulate the buoyancy-driven fluids near the hydrostatic equilibrium and the electrically conducting fluids near a background magnetic field. In the first problem, a Boussinesq-Navier-Stokes system governing the perturbations near the hydrostatic equilibrium will be solved. The Boussinesq system couples the Navier-Stokes equations forced by buoyancy with the temperature transport equation. The numerical method will also be extended to simulate anisotropic flows for which the vertical viscosity is much smaller than the horizontal one. This problem arises in modeling turbulent flows in Ekman layers occurring in the atmosphere and the ocean. In the second problem, the stability and long-time behavior of electrically conducting fluids under a guiding magnetic field will be computed and analyzed. When complemented with rigorous analysis, these simulations will help accelerate the resolution of several open stability problems on the Boussinesq and the MHD equations.This project is jointly funded by the Computational Mathematics program and the Established Program to Stimulate Competitive Research (EPSCoR).This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.