More Like this

1. Modeling epidemic flow with fluid dynamics

   https://doi.org/10.3934/mbe.2022388  

   Cheng, Ziqiang; Wang, Jin (January 2022, Mathematical Biosciences and Engineering)

   In this paper, a new mathematical model based on partial differential equations is proposed to study the spatial spread of infectious diseases. The model incorporates fluid dynamics theory and represents the epidemic spread as a fluid motion generated through the interaction between the susceptible and infected hosts. At the macroscopic level, the spread of the infection is modeled as an inviscid flow described by the Euler equation. Nontrivial numerical methods from computational fluid dynamics (CFD) are applied to investigate the model. In particular, a fifth-order weighted essentially non-oscillatory (WENO) scheme is employed for the spatial discretization. As an application, this mathematical and computational framework is used in a simulation study for the COVID-19 outbreak in Wuhan, China. The simulation results match the reported data for the cumulative cases with high accuracy and generate new insight into the complex spatial dynamics of COVID-19. 

   more » « less
2. An Improved Level-Set-Based Immersed Boundary Reconstruction Method for Computing Bio-Inspired Underwater Propulsion

   Pan, Y. (August 2021, ASME Fluids Engineering Division Summer Meeting)

   The immersed boundary method (IBM) has been widely employed to study bio-inspired underwater propulsion which often involves the high Reynolds number, complex body morphologies and large computational domain. Due to these problems, the immersed boundary (IB) reconstruction can be very costly in a simulation. Based on our previous work, an improved level-set-based immersed boundary method (LS-IBM) has been developed in this paper by introducing the narrow-band technique. Comparing with the previous LS-IBM, the narrowband level-set-based immersed boundary method (NBLS-IBM) is only required to propagate the level set values from the points near the boundaries to all the points in the narrow band. By simulating a steadyswimming Jackfish-like body, the consistency and stability of the new reconstruction method in the flow solver have been verified. Applications to a dolphin-like body swimming and a shark-like body swimming are used to demonstrate the efficiency and accuracy of the NBLS-IBM. The time for reconstructions shows that the reconstruction efficiency can increase up to 64.6% by using the NBLS-IBM while keeping the accuracy and robustness of the original LS-IBM. The vortex wake of the shark-like body in steady swimming shows the robustness, fastness and compatibility of the NBLS-IBM to our current flow solver 

   more » « less
3. Von Neumann Stable, Implicit, High Order, Finite Volume WENO Schemes

   https://doi.org/10.2118/193817-MS  

   Arbogast, Todd; Huang, Chieh-Sen; Zhao, Xikai (April 2019, SPE Reservoir Simulation Conference 2019)

   Simulation of flow and transport in petroleum reservoirs involves solving coupled systems of advection-diffusion-reaction equations with nonlinear flux functions, diffusion coefficients, and reactions/wells. It is important to develop numerical schemes that can approximate all three processes at once, and to high order, so that the physics can be well resolved. In this paper, we propose an approach based on high order, finite volume, implicit, Weighted Essentially NonOscillatory (iWENO) schemes. The resulting schemes are locally mass conservative and, being implicit, suited to systems of advection-diffusion-reaction equations. Moreover, our approach gives unconditionally L-stable schemes for smooth solutions to the linear advection-diffusion-reaction equation in the sense of a von Neumann stability analysis. To illustrate our approach, we develop a third order iWENO scheme for the saturation equation of two-phase flow in porous media in two space dimensions. The keys to high order accuracy are to use WENO reconstruction in space (which handles shocks and steep fronts) combined with a two-stage Radau-IIA Runge-Kutta time integrator. The saturation is approximated by its averages over the mesh elements at the current time level and at two future time levels; therefore, the scheme uses two unknowns per grid block per variable, independent of the spatial dimension. This makes the scheme fairly computationally efficient, both because reconstructions make use of local information that can fit in cache memory, and because the global system has about as small a number of degrees of freedom as possible. The scheme is relatively simple to implement, high order accurate, maintains local mass conservation, applies to general computational meshes, and appears to be robust. Preliminary computational tests show the potential of the scheme to handle advection-diffusion-reaction processes on meshes of quadrilateral gridblocks, and to do so to high order accuracy using relatively long time steps. The new scheme can be viewed as a generalization of standard cell-centered finite volume (or finite difference) methods. It achieves high order in both space and time, and it incorporates WENO slope limiting. 

   more » « less
4. Essentially non-oscillatory and weighted essentially non-oscillatory schemes

   https://doi.org/10.1017/S0962492920000057  

   Shu, Chi-Wang (May 2020, Acta Numerica)

   null (Ed.)

   Essentially non-oscillatory (ENO) and weighted ENO (WENO) schemes were designed for solving hyperbolic and convection–diffusion equations with possibly discontinuous solutions or solutions with sharp gradient regions. The main idea of ENO and WENO schemes is actually an approximation procedure, aimed at achieving arbitrarily high-order accuracy in smooth regions and resolving shocks or other discontinuities sharply and in an essentially non-oscillatory fashion. Both finite volume and finite difference schemes have been designed using the ENO or WENO procedure, and these schemes are very popular in applications, most noticeably in computational fluid dynamics but also in other areas of computational physics and engineering. Since the main idea of the ENO and WENO schemes is an approximation procedure not directly related to partial differential equations (PDEs), ENO and WENO schemes also have non-PDE applications. In this paper we will survey the basic ideas behind ENO and WENO schemes, discuss their properties, and present examples of their applications to different types of PDEs as well as to non-PDE problems. 

   more » « less
5. Multidimensional WENO-AO Reconstructions Using a Simplified Smoothness Indicator and Applications to Conservation Laws

   https://doi.org/10.1007/s10915-023-02319-x  

   Huang, Chieh-Sen; Arbogast, Todd; Tian, Chenyu (August 2023, Journal of Scientific Computing)

   Abstract Finite volume, weighted essentially non-oscillatory (WENO) schemes require the computation of a smoothness indicator. This can be expensive, especially in multiple space dimensions. We consider the use of the simple smoothness indicator$$\sigma ^{\textrm{S}}= \frac{1}{N_{\textrm{S}}-1}\sum _{j} ({\bar{u}}_{j} - {\bar{u}}_{m})^2$$ ${\sigma }^{\text{S}}=\frac{1}{N_{\text{S}}-1}{\sum }_j{({\overline{u}}_j-{\overline{u}}_m)}^2$, where$$N_{\textrm{S}}$$ $N_{\text{S}}$is the number of mesh elements in the stencil,$${\bar{u}}_j$$ ${\overline{u}}_j$is the local function average over mesh elementj, and indexmgives the target element. Reconstructions utilizing standard WENO weighting fail with this smoothness indicator. We develop a modification of WENO-Z weighting that gives a reliable and accurate reconstruction of adaptive order, which we denote as SWENOZ-AO. We prove that it attains the order of accuracy of the large stencil polynomial approximation when the solution is smooth, and drops to the order of the small stencil polynomial approximations when there is a jump discontinuity in the solution. Numerical examples in one and two space dimensions on general meshes verify the approximation properties of the reconstruction. They also show it to be about 10 times faster in two space dimensions than reconstructions using the classic smoothness indicator. The new reconstruction is applied to define finite volume schemes to approximate the solution of hyperbolic conservation laws. Numerical tests show results of the same quality as standard WENO schemes using the classic smoothness indicator, but with an overall speedup in the computation time of about 3.5–5 times in 2D tests. Moreover, the computational efficiency (CPU time versus error) is noticeably improved. 

   more » « less