We derive macroscopic dynamics for self-propelled particles in a fluid. The starting point is a coupled Vicsek–Stokes system. The Vicsek model describes self-propelled agents interacting through alignment. It provides a phenomenological description of hydrodynamic interactions between agents at high density. Stokes equations describe a low Reynolds number fluid. These two dynamics are coupled by the interaction between the agents and the fluid. The fluid contributes to rotating the particles through Jeffery’s equation. Particle self-propulsion induces a force dipole on the fluid. After coarse-graining we obtain a coupled Self-Organised Hydrodynamics–Stokes system. We perform a linear stability analysis for this system which shows that both pullers and pushers have unstable modes. We conclude by providing extensions of the Vicsek–Stokes model including short-distance repulsion, finite particle inertia and finite Reynolds number fluid regime.

The onset of thermal convection in a rapidly rotating spherical shell is studied by linear stability analysis based on the fully compressible Navier-Stokes equations. Compressibility is quantified by the number of density scale heights N-rho, which measures the intensity of density stratification of the motionless, polytropic base state. The nearly adiabatic flow with polytropic index n = 1.499 < na = 1.5 is considered, where na is the adiabatic polytropic index. By investigating the stability of the base state with respect to the disturbance of specified wavenumber, the instability process is found to be sensitive to the Prandtl number Pr and to N-rho. For large Pr and small N-rho, the quasi-geostrophic columnar mode loses stability first; while for relatively small Pr a new quasi-geostrophic compressible mode is identified, which becomes unstable first under strong density stratification. The inertial mode can also occur first for relatively small Pr and a certain intensity of density stratification in the parameter range considered. Although the Rayleigh numbers Ra for the onsets of the quasi-geostrophic compressible mode and columnar mode are different by several orders of magnitude, we find that they follow very similar scaling laws with the Taylor number. The critical Ra for convection onset is found to be always positive, in contrast with previous results based on the widely used anelastic model that convection can occur at negative Ra. By evaluating the relative magnitude of the time derivative of density perturbation in the continuity equation, we show that the anelastic approximation in the present system cannot be applied in the small-Ra and large-N-rho regime.

Gaseous fluids may move slowly, as smoke does, or at high speed, such as occurs with explosions. High-speed gas flow is always accompanied by low-speed gas flow, which produces rich visual details in the fluid motion. Realistic visualization involves a complex dynamic flow field with both low and high speed fluid behavior. In computer graphics, algorithms to simulate gaseous fluids address either the low speed case or the high speed case, but no algorithm handles both efficiently. With the aim of providing visually pleasing results, we present a hybrid algorithm that efficiently captures the essential physics of both low- and high-speed gaseous fluids. We model the low speed gaseous fluids by a grid approach and use a particle approach for the high speed gaseous fluids. In addition, we propose a physically sound method to connect the particle model to the grid model. By exploiting complementary strengths and avoiding weaknesses of the grid and particle approaches, we produce some animation examples and analyze their computational performance to demonstrate the effectiveness of the new hybrid method.

Surface flow phenomena, such as rain water flowing down a tree trunk and progressive water front in a shower room, are common in real life. However, compared with the 3D spatial fluid flow, these surface flow problems have been much less studied in the graphics community. To tackle this research gap, we present an efficient, robust and high-fidelity simulation approach based on the shallow-water equations. Specifically, the standard shallow-water flow model is extended to general triangle meshes with a feature-based bottom friction model, and a series of coherent mathematical formulations are derived to represent the full range of physical effects that are important for real-world surface flow phenomena. In addition, by achieving compatibility with existing 3D fluid simulators and by supporting physically realistic interactions with multiple fluids and solid surfaces, the new model is flexible and readily extensible for coupled phenomena. A wide range of simulation examples are presented to demonstrate the performance of the new approach.

In this paper, our aim is to present (1) an embedded fracture model (EFM) for coupled flow and mechanics problem based on the dual continuum approach on the fine grid and (2) an upscaled model for the resulting fine grid equations. The mathematical model is described by the coupled system of equation for displacement, fracture and matrix pressures. For a fine grid approximation, we use the finite volume method for flow problem and finite element method for mechanics. Due to the complexity of fractures, solutions have a variety of scales, and fine grid approximation results in a large discrete system. Our second focus in the construction of the upscaled coarse grid poroelasticity model for fractured media. Our upscaled approach is based on the nonlocal multicontinuum (NLMC) upscaling for coupled flow and mechanics problem, which involves computations of local basis functions via an energy minimization