Search

Reduced-order models (ROMs) have achieved a lot of success in reducing the computational cost of traditional numerical methods across many disciplines. In fluid dynamics, ROMs have been successful in providing efficient and relatively accurate solutions for the numerical simulation of laminar flows. For convection-dominated (e.g., turbulent) flows, however, standard ROMs generally yield inaccurate results, usually affected by spurious oscillations. Thus, ROMs are usually equipped with numerical stabilization or closure models in order to account for the effect of the discarded modes. The literature on ROM closures and stabilizations is large and growing fast. In this paper, instead of reviewing all the ROM closures and stabilizations, we took a more modest step and focused on one particular type of ROM closure and stabilization that is inspired by large eddy simulation (LES), a classical strategy in computational fluid dynamics (CFD). These ROMs, which we call LES-ROMs, are extremely easy to implement, very efficient, and accurate. Indeed, LES-ROMs are modular and generally require minimal modifications to standard (“legacy”) ROM formulations. Furthermore, the computational overhead of these modifications is minimal. Finally, carefully tuned LES-ROMs can accurately capture the average physical quantities of interest in challenging convection-dominated flows in science and engineering applications. LES-ROMs are constructed by leveraging spatial filtering, which is the same principle used to build classical LES models. This ensures a modeling consistency between LES-ROMs and the approaches that generated the data used to train them. It also “bridges” two distinct research fields (LES and ROMs) that have been disconnected until now. This paper is a review of LES-ROMs, with a particular focus on the LES concepts and models that enable the construction of LES-inspired ROMs and the bridging of LES and reduced-order modeling. This paper starts with a description of a versatile LES strategy called evolve–filter–relax (EFR) that has been successfully used as a full-order method for both incompressible and compressible convection-dominated flows. We present evidence of this success. We then show how the EFR strategy, and spatial filtering in general, can be leveraged to construct LES-ROMs (e.g., EFR-ROM). Several applications of LES-ROMs to the numerical simulation of incompressible and compressible convection-dominated flows are presented. Finally, we draw conclusions and outline several research directions and open questions in LES-ROM development. While we do not claim this review to be comprehensive, we certainly hope it serves as a brief and friendly introduction to this exciting research area, which we believe has a lot of potential in the practical numerical simulation of convection-dominated flows in science, engineering, and medicine.

The COVID-19 pandemic is still ongoing, even if the emergency is over, and we now have enough data to analyze the outbreak over a long timeline. There is evidence that the outbreak alternates periods of high and low infections. Retrospectively, this can help in understanding the nature of an appropriate mathematical model for this dramatic infection. The periodic behavior may be the consequence of time-dependent coefficients related to seasonal effects and specific political actions, or an intrinsic feature of the model. The present paper relies on the assumption that the periodic spikes are an intrinsic feature of the disease, and, as such, it should be properly reflected in the mathematical model. Based on the concept of waning immunity proposed for other pathologies, we introduce a new model with (i) a compartment for weakly immune people subject to immunity booster, represented by a non-linear term; (ii) discrimination between individuals infected/vaccinated for the first time, and individuals already infected/vaccinated, undergoing to new infections/doses. We analyze some preliminary properties of our model, called SIRW2, and provide a proof-of-concept that it is capable of reproducing qualitatively the long-term oscillatory behavior of COVID-19 infection.

Kidney dialysis is the most widespread treatment method for end-stage renal disease, a debilitating health condition common in industrialized societies. While ubiquitous, kidney dialysis suffers from an inability to remove larger toxins, resulting in a gradual buildup of these toxins in dialysis patients, ultimately leading to further health complications. To improve dialysis, hollow fibers incorporating a cell-monolayer with cultured kidney cells have been proposed; however, the design of such a fiber is nontrivial. In particular, the effects of fluid wall-shear stress have an important influence on the ability of the cell layer to transport toxins. In the present work, we introduce a model for cell-transport aided dialysis, incorporating the effects of the shear stress. We analyze the model mathematically and establish its well-posedness. We then present a series of numerical results, which suggest that a hollow-fiber design with a wavy profile may increase the efficiency of the dialysis treatment. We investigate numerically the shape of the wavy channel to maximize the toxin clearance. These results demonstrate the potential for the use of computational models in the study and advancement of renal therapies.

Transcatheter aortic valve replacement (TAVR) is a rapidly growing field enabling replacement of diseased aortic valves without the need for open heart surgery. However, due to the nature of the procedure and nonremoval of the diseased tissue, there are rates of complications ranging from tissue rupture and coronary obstruction to paravalvular leak, valve thrombosis, and permanent pacemaker implantation. In recent years, computational modeling has shown a great deal of promise in its capabilities to understand the biomechanical implications of TAVR as well as help preoperatively predict risks inherent to device–patient-specific anatomy biomechanical interaction. This includes intricate replication of stent and leaflet designs and tested and validated simulated deployments with structural and fluid mechanical simulations. This review outlines current biomechanical understanding of device-related complications from TAVR and related predictive strategies using computational modeling. An outlook on future modeling strategies highlighting reduced order modeling which could significantly reduce the high time and cost that are required for computational prediction of TAVR outcomes is presented in this review paper. A summary of current commercial/in-development software is presented in the final section.

Total cavopulmonary connection (TCPC) hemodynamics has been hypothesized to be associated with long-term complications in single ventricle heart defect patients. Rigid wall assumption has been commonly used when evaluating TCPC hemodynamics using computational fluid dynamics (CFD) simulation. Previous study has evaluated impact of wall compliance on extra-cardiac TCPC hemodynamics using fluid-structure interaction (FSI) simulation. However, the impact of ignoring wall compliance on the presumably more compliant intra-atrial TCPC hemodynamics is not fully understood. To narrow this knowledge gap, this study aims to investigate impact of wall compliance on an intra-atrial TCPC hemodynamics. A patient-specific model of an intra-atrial TCPC is simulated with an FSI model. Patient-specific 3D TCPC anatomies were reconstructed from transverse cardiovascular magnetic resonance images. Patient-specific vessel flow rate from phase-contrast magnetic resonance imaging (MRI) at the Fontan pathway and the superior vena cava under resting condition were prescribed at the inlets. From the FSI simulation, the degree of wall deformation was compared with in vivo wall deformation from phase-contrast MRI data as validation of the FSI model. Then, TCPC flow structure, power loss and hepatic flow distribution (HFD) were compared between rigid wall and FSI simulation. There were differences in instantaneous pressure drop, power loss and HFD between rigid wall and FSI simulations, but no difference in the time-averaged quantities. The findings of this study support the use of a rigid wall assumption on evaluation of time-averaged intra-atrial TCPC hemodynamic metric under resting breath-held condition.

Although plaque erosion causes approximately 40% of all coronary thrombi and disproportionally affects women more than men, its mechanism is not well understood. The role of tissue mechanics in plaque rupture and regulation of mechanosensitive inflammatory proteins is well established, but their role in plaque erosion is unknown. Given obvious differences in morphology between plaque erosion and rupture, we hypothesized that inflammation in general as well as the association between local mechanical strain and inflammation known to exist in plaque rupture may not occur in plaque erosion. Therefore, our objective was to determine if similar mechanisms underlie plaque rupture and plaque erosion.We studied a total of 74 human coronary plaque specimens obtained at autopsy. Using lesion-specific computer modeling of solid mechanics, we calculated the stress and strain distribution for each plaque and determined if there were any relationships with markers of inflammation. Consistent with previous studies, inflammatory markers were positively associated with increasing strain in specimens with rupture and thin-cap fibroatheromas. Conversely, overall staining for inflammatory markers and apoptosis were significantly lower in erosion, and there was no relationship with mechanical strain. Samples with plaque erosion most closely resembled those with the stable phenotype of thick-cap fibroatheromas.In contrast to classic plaque rupture, plaque erosion was not associated with markers of inflammation and mechanical strain. These data suggest that plaque erosion is a distinct pathophysiological process with a different etiology and therefore raises the possibility that a different therapeutic approach may be required to prevent plaque erosion.