Showing 1–50 of 58 results for author: Bagheri, S

BraTS-PEDs: Results of the Multi-Consortium International Pediatric Brain Tumor Segmentation Challenge 2023

Authors: Anahita Fathi Kazerooni, Nastaran Khalili, Xinyang Liu, Debanjan Haldar, Zhifan Jiang, Anna Zapaishchykova, Julija Pavaine, Lubdha M. Shah, Blaise V. Jones, Nakul Sheth, Sanjay P. Prabhu, Aaron S. McAllister, Wenxin Tu, Khanak K. Nandolia, Andres F. Rodriguez, Ibraheem Salman Shaikh, Mariana Sanchez Montano, Hollie Anne Lai, Maruf Adewole, Jake Albrecht, Udunna Anazodo, Hannah Anderson, Syed Muhammed Anwar, Alejandro Aristizabal, Sina Bagheri , et al. (54 additional authors not shown)

Abstract: Pediatric central nervous system tumors are the leading cause of cancer-related deaths in children. The five-year survival rate for high-grade glioma in children is less than 20%. The development of new treatments is dependent upon multi-institutional collaborative clinical trials requiring reproducible and accurate centralized response assessment. We present the results of the BraTS-PEDs 2023 cha… ▽ More Pediatric central nervous system tumors are the leading cause of cancer-related deaths in children. The five-year survival rate for high-grade glioma in children is less than 20%. The development of new treatments is dependent upon multi-institutional collaborative clinical trials requiring reproducible and accurate centralized response assessment. We present the results of the BraTS-PEDs 2023 challenge, the first Brain Tumor Segmentation (BraTS) challenge focused on pediatric brain tumors. This challenge utilized data acquired from multiple international consortia dedicated to pediatric neuro-oncology and clinical trials. BraTS-PEDs 2023 aimed to evaluate volumetric segmentation algorithms for pediatric brain gliomas from magnetic resonance imaging using standardized quantitative performance evaluation metrics employed across the BraTS 2023 challenges. The top-performing AI approaches for pediatric tumor analysis included ensembles of nnU-Net and Swin UNETR, Auto3DSeg, or nnU-Net with a self-supervised framework. The BraTSPEDs 2023 challenge fostered collaboration between clinicians (neuro-oncologists, neuroradiologists) and AI/imaging scientists, promoting faster data sharing and the development of automated volumetric analysis techniques. These advancements could significantly benefit clinical trials and improve the care of children with brain tumors. △ Less

Submitted 11 July, 2024; originally announced July 2024.

Drag prediction of rough-wall turbulent flow using data-driven regression

Authors: Zhaoyu Shi, Seyed Morteza Habibi Khorasani, Heesoo Shin, Jiasheng Yang, Sangseung Lee, Shervin Bagheri

Abstract: Efficient tools for predicting the drag of rough walls in turbulent flows would have a tremendous impact. However, methods for drag prediction rely on experiments or numerical simulations which are costly and time-consuming. Data-driven regression methods have the potential to provide a prediction that is accurate and fast. We assess the performance and limitations of linear regression, kernel met… ▽ More Efficient tools for predicting the drag of rough walls in turbulent flows would have a tremendous impact. However, methods for drag prediction rely on experiments or numerical simulations which are costly and time-consuming. Data-driven regression methods have the potential to provide a prediction that is accurate and fast. We assess the performance and limitations of linear regression, kernel methods and neural networks for drag prediction using a database of 1000 homogeneous rough surfaces. Model performance is evaluated using the roughness function obtained at friction-scaled Reynolds number 500. With two trainable parameters, the kernel method can fully account for nonlinear relations between $ΔU^+$ and surface statistics (roughness height, effective slope, skewness, etc). In contrast, linear regression cannot account for nonlinear correlations and display large errors and high uncertainty. Multilayer perceptron and convolutional neural networks demonstrate performance on par with the kernel method but have orders of magnitude more trainable parameters. For the current database size, the networks' capacity cannot be fully exploited, resulting in reduced generalizability and reliability. Our study provides insight into the appropriateness of different regression models for drag prediction. We also discuss the remaining steps before data-driven methods emerge as useful tools in applications. △ Less

Submitted 15 May, 2024; originally announced May 2024.

Comments: This manuscript consists of 18 pages, where 2 appendices, 11 figures and 3 tables are included. It is currently under review at FLOW journal. Dr. Zhaoyu Shi developed the machine learning models and conducted the data analysis as well as the draft writing. The direct numerical simulations (see Appendix A) were conducted by Dr S.M.H. Khorasani

Data-driven discovery of drag-inducing elements on a rough surface through convolutional neural networks

Authors: Heesoo Shin, Seyed Morteza Habibi Khorasani, Zhaoyu Shi, Jiasheng Yang, Sangseung Lee, Shervin Bagheri

Abstract: Understanding the influence of surface roughness on drag forces remains a significant challenge in fluid dynamics. This paper presents a convolutional neural network (CNN) that predicts drag solely by the topography of rough surfaces and is capable of discovering spatial patterns linked to drag-inducing structures. A CNN model was developed to analyze spatial information from the topography of a r… ▽ More Understanding the influence of surface roughness on drag forces remains a significant challenge in fluid dynamics. This paper presents a convolutional neural network (CNN) that predicts drag solely by the topography of rough surfaces and is capable of discovering spatial patterns linked to drag-inducing structures. A CNN model was developed to analyze spatial information from the topography of a rough surface and predict the roughness function, $ΔU^+$, obtained from direct numerical simulation. This model enables the prediction of drag from rough surface data alone, which was not possible with previous methods owing to the large number of surface-derived parameters. Additionally, the retention of spatial information by the model enables the creation of a feature map that accentuates critical areas for drag prediction on rough surfaces. By interpreting the feature maps, we show that the developed CNN model is able to discover spatial patterns associated with drag distributions across rough surfaces, even without a direct training on drag distribution data. The analysis of the feature map indicates that, even without flow field information, the CNN model extracts the importance of the flow-directional slope and height of roughness elements as key factors in inducing pressure drag. This study demonstrates that CNN-based drag prediction is grounded in physical principles of fluid dynamics, underscoring the utility of CNNs in both predicting and understanding drag on rough surfaces. △ Less

Submitted 14 May, 2024; originally announced May 2024.

The role of fluid friction in streamer formation and biofilm growth

Authors: Cornelius Wittig, Michael Wagner, Romain Vallon, Thomas Crouzier, Wouter van der Wijngaart, Harald Horn, Shervin Bagheri

Abstract: Bacillus subtilis biofilms were grown in laminar channel flow at wall shear stress spanning one order of magnitude ($τ_w = 0.068$ Pa to $τ_w = 0.67$ Pa). We monitor, non-invasively, the evolution of the three-dimensional distribution of biofilm over seven days using optical coherence tomography (OCT). The obtained biofilms consist of many microcolonies where the characteristic colony has a base st… ▽ More Bacillus subtilis biofilms were grown in laminar channel flow at wall shear stress spanning one order of magnitude ($τ_w = 0.068$ Pa to $τ_w = 0.67$ Pa). We monitor, non-invasively, the evolution of the three-dimensional distribution of biofilm over seven days using optical coherence tomography (OCT). The obtained biofilms consist of many microcolonies where the characteristic colony has a base structure in the form of a leaning pillar and a streamer in the form of a thin filament that originates near the tip of the pillar. While the shape, size and distribution of these microcolonies depend on the imposed shear stress, the same structural features appear consistently for all shear stress values. The formation of streamers seems to occur after the development of a base structure, suggesting that the latter induces a curved secondary flow that triggers the formation of the streamers. Moreover, we observe that the biofilm volume grows approximately linearly over seven days for all the shear stress values, with a growth rate that is inversely proportional to the wall shear stress. We develop a simple model of friction-limited growth, which agrees with the experimental observations. The model provides physical insight into growth mechanisms and can be used to develop accurate continuum models of bacterial biofilm growth. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Comments: 22 pages, 20 figures

Training and Comparison of nnU-Net and DeepMedic Methods for Autosegmentation of Pediatric Brain Tumors

Authors: Arastoo Vossough, Nastaran Khalili, Ariana M. Familiar, Deep Gandhi, Karthik Viswanathan, Wenxin Tu, Debanjan Haldar, Sina Bagheri, Hannah Anderson, Shuvanjan Haldar, Phillip B. Storm, Adam Resnick, Jeffrey B. Ware, Ali Nabavizadeh, Anahita Fathi Kazerooni

Abstract: Brain tumors are the most common solid tumors and the leading cause of cancer-related death among children. Tumor segmentation is essential in surgical and treatment planning, and response assessment and monitoring. However, manual segmentation is time-consuming and has high inter-operator variability, underscoring the need for more efficient methods. We compared two deep learning-based 3D segment… ▽ More Brain tumors are the most common solid tumors and the leading cause of cancer-related death among children. Tumor segmentation is essential in surgical and treatment planning, and response assessment and monitoring. However, manual segmentation is time-consuming and has high inter-operator variability, underscoring the need for more efficient methods. We compared two deep learning-based 3D segmentation models, DeepMedic and nnU-Net, after training with pediatric-specific multi-institutional brain tumor data using based on multi-parametric MRI scans.Multi-parametric preoperative MRI scans of 339 pediatric patients (n=293 internal and n=46 external cohorts) with a variety of tumor subtypes, were preprocessed and manually segmented into four tumor subregions, i.e., enhancing tumor (ET), non-enhancing tumor (NET), cystic components (CC), and peritumoral edema (ED). After training, performance of the two models on internal and external test sets was evaluated using Dice scores, sensitivity, and Hausdorff distance with reference to ground truth manual segmentations. Dice score for nnU-Net internal test sets was (mean +/- SD (median)) 0.9+/-0.07 (0.94) for WT, 0.77+/-0.29 for ET, 0.66+/-0.32 for NET, 0.71+/-0.33 for CC, and 0.71+/-0.40 for ED, respectively. For DeepMedic the Dice scores were 0.82+/-0.16 for WT, 0.66+/-0.32 for ET, 0.48+/-0.27, for NET, 0.48+/-0.36 for CC, and 0.19+/-0.33 for ED, respectively. Dice scores were significantly higher for nnU-Net (p<=0.01). External validation of the trained nnU-Net model on the multi-institutional BraTS-PEDs 2023 dataset revealed high generalization capability in segmentation of whole tumor and tumor core with Dice scores of 0.87+/-0.13 (0.91) and 0.83+/-0.18 (0.89), respectively. Pediatric-specific data trained nnU-Net model is superior to DeepMedic for whole tumor and subregion segmentation of pediatric brain tumors. △ Less

Submitted 30 January, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Extreme types and extremal models

Authors: Seyed-Mohammad Bagheri

Abstract: In the affine fragment of continuous logic, type spaces are compact convex sets. I study some model theoretic properties of extreme types. It is proved that every complete theory $T$ has an extremal model, i.e. a model which realizes only extreme types. Extremal models form an elementary class in the full continuous logic sense if and only if the set of extreme $n$-types is closed in $S_n(T)$ for… ▽ More In the affine fragment of continuous logic, type spaces are compact convex sets. I study some model theoretic properties of extreme types. It is proved that every complete theory $T$ has an extremal model, i.e. a model which realizes only extreme types. Extremal models form an elementary class in the full continuous logic sense if and only if the set of extreme $n$-types is closed in $S_n(T)$ for each $n$. Also, some applications are given in the special cases where the theory has a compact or first order model. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Definability in affine continuous logic

Authors: Seyed-Mohammad Bagheri

Abstract: I study definable sets in affine continuous logic. Let $T$ be an affine theory. After giving some general results, it is proved that if $T$ has a first order model, its extremal theory is a complete first order theory and first order definable sets are affinely definable. In this case, the type spaces of $T$ are Bauer simplices and they coincide with the sets of Keisler measures of the extremal th… ▽ More I study definable sets in affine continuous logic. Let $T$ be an affine theory. After giving some general results, it is proved that if $T$ has a first order model, its extremal theory is a complete first order theory and first order definable sets are affinely definable. In this case, the type spaces of $T$ are Bauer simplices and they coincide with the sets of Keisler measures of the extremal theory. In contrast, if $T$ has a compact model, definable sets are exactly the end-sets of definable predicates. As an example, it is proved in the theory of probability algebras that one dimensional definable sets are exactly the intervals $[a,b]$. △ Less

Submitted 12 March, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Contact-angle hysteresis provides resistance to drainage of liquid-infused surfaces in turbulent flows

Authors: Sofia Saoncella, Si Suo, Johan Sundin, Agastya Parikh, Marcus Hultmark, Wouter Metsola van der Wijngaart, Fredrik Lundell, Shervin Bagheri

Abstract: Lubricated textured surfaces immersed in liquid flows offer tremendous potential for reducing fluid drag, enhancing heat and mass transfer, and preventing fouling. According to current design rules, the lubricant must chemically match the surface to remain robustly trapped within the texture. However, achieving such chemical compatibility poses a significant challenge for large-scale flow systems,… ▽ More Lubricated textured surfaces immersed in liquid flows offer tremendous potential for reducing fluid drag, enhancing heat and mass transfer, and preventing fouling. According to current design rules, the lubricant must chemically match the surface to remain robustly trapped within the texture. However, achieving such chemical compatibility poses a significant challenge for large-scale flow systems, as it demands advanced surface treatments or severely limits the range of viable lubricants. In addition, chemically tuned surfaces often degrade over time in harsh environments. Here, we demonstrate that a lubricant-infused surface (LIS) can resist drainage in the presence of external shear flow without requiring chemical compatibility. Surfaces featuring longitudinal grooves can retain up to 50% of partially wetting lubricants in fully developed turbulent flows. The retention relies on contact-angle hysteresis, where triple-phase contact lines are pinned to substrate heterogeneities, creating capillary resistance that prevents lubricant depletion. We develop an analytical model to predict the maximum length of pinned lubricant droplets in microgrooves. This model, validated through a combination of experiments and numerical simulations, can be used to design chemistry-free LISs for applications where the external environment is continuously flowing. Our findings open up new possibilities for using functional surfaces to control transport processes in large systems. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 20 pages, 12 figures

Dewetting of a corner film wrapping a wall-mounted cylinder

Authors: Si Suo, Seyed Morteza Habibi Khorasani, Shervin Bagheri

Abstract: In this study, we investigate the stability of a liquid film that partially wets a corner between a cylinder and a substrate, using a combination of theoretical and numerical approaches. The film stability, which depends on the film size and the wall wettability, is firstly predicted by a standard linear stability analysis (LSA) within the long-wave theoretical framework. We find that the film siz… ▽ More In this study, we investigate the stability of a liquid film that partially wets a corner between a cylinder and a substrate, using a combination of theoretical and numerical approaches. The film stability, which depends on the film size and the wall wettability, is firstly predicted by a standard linear stability analysis (LSA) within the long-wave theoretical framework. We find that the film size plays the most important role in controlling the film stability. Specifically, the thinner the film is, the more sensitive it becomes to the large-wavenumber perturbation. The wall wettability mainly impacts the growth rates of perturbations and slightly influences the marginal stability and post-instability patterns of wrapping films. We compare the LSA predictions with numerical results obtained from a disjoining pressure model (DPM) and Volume-of-Fluid (VOF) simulations, which provide more insights into the film breakup process. At the early stage there is a strong agreement between the LSA predictions and the DPM results. Notably, as the perturbation grows, thin film regions connecting two neighboring satellite droplets form which may eventually lead to a stable or temporary secondary droplet, an aspect which the LSA is incapable of capturing. In addition, the VOF simulations suggest that beyond a critical film size the crest coalescence mechanism becomes involved during the breakup stage. Therefore, the LSA predictions are able to provide only an upper limit on the final number of satellite droplets. △ Less

Submitted 9 October, 2023; originally announced October 2023.

A multi-institutional pediatric dataset of clinical radiology MRIs by the Children's Brain Tumor Network

Authors: Ariana M. Familiar, Anahita Fathi Kazerooni, Hannah Anderson, Aliaksandr Lubneuski, Karthik Viswanathan, Rocky Breslow, Nastaran Khalili, Sina Bagheri, Debanjan Haldar, Meen Chul Kim, Sherjeel Arif, Rachel Madhogarhia, Thinh Q. Nguyen, Elizabeth A. Frenkel, Zeinab Helili, Jessica Harrison, Keyvan Farahani, Marius George Linguraru, Ulas Bagci, Yury Velichko, Jeffrey Stevens, Sarah Leary, Robert M. Lober, Stephani Campion, Amy A. Smith , et al. (15 additional authors not shown)

Abstract: Pediatric brain and spinal cancers remain the leading cause of cancer-related death in children. Advancements in clinical decision-support in pediatric neuro-oncology utilizing the wealth of radiology imaging data collected through standard care, however, has significantly lagged other domains. Such data is ripe for use with predictive analytics such as artificial intelligence (AI) methods, which… ▽ More Pediatric brain and spinal cancers remain the leading cause of cancer-related death in children. Advancements in clinical decision-support in pediatric neuro-oncology utilizing the wealth of radiology imaging data collected through standard care, however, has significantly lagged other domains. Such data is ripe for use with predictive analytics such as artificial intelligence (AI) methods, which require large datasets. To address this unmet need, we provide a multi-institutional, large-scale pediatric dataset of 23,101 multi-parametric MRI exams acquired through routine care for 1,526 brain tumor patients, as part of the Children's Brain Tumor Network. This includes longitudinal MRIs across various cancer diagnoses, with associated patient-level clinical information, digital pathology slides, as well as tissue genotype and omics data. To facilitate downstream analysis, treatment-naïve images for 370 subjects were processed and released through the NCI Childhood Cancer Data Initiative via the Cancer Data Service. Through ongoing efforts to continuously build these imaging repositories, our aim is to accelerate discovery and translational AI models with real-world data, to ultimately empower precision medicine for children. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Graph Sparsification for GCN Towards Optimal Crop Yield Predictions

Authors: Saghar Bagheri, Gene Cheung, Tim Eadie

Abstract: In agronomics, predicting crop yield at a per field/county granularity is important for farmers to minimize uncertainty and plan seeding for the next crop cycle. While state-of-the-art prediction techniques employ graph convolutional nets (GCN) to predict future crop yields given relevant features and crop yields of previous years, a dense underlying graph kernel requires long training and executi… ▽ More In agronomics, predicting crop yield at a per field/county granularity is important for farmers to minimize uncertainty and plan seeding for the next crop cycle. While state-of-the-art prediction techniques employ graph convolutional nets (GCN) to predict future crop yields given relevant features and crop yields of previous years, a dense underlying graph kernel requires long training and execution time. In this paper, we propose a graph sparsification method based on the Fiedler number to remove edges from a complete graph kernel, in order to lower the complexity of GCN training/execution. Specifically, we first show that greedily removing an edge at a time that induces the minimal change in the second eigenvalue leads to a sparse graph with good GCN performance. We then propose a fast method to choose an edge for removal per iteration based on an eigenvalue perturbation theorem. Experiments show that our Fiedler-based method produces a sparse graph with good GCN performance compared to other graph sparsification schemes in crop yield prediction. △ Less

Submitted 2 June, 2023; originally announced June 2023.

The Brain Tumor Segmentation (BraTS) Challenge 2023: Focus on Pediatrics (CBTN-CONNECT-DIPGR-ASNR-MICCAI BraTS-PEDs)

Authors: Anahita Fathi Kazerooni, Nastaran Khalili, Xinyang Liu, Debanjan Haldar, Zhifan Jiang, Syed Muhammed Anwar, Jake Albrecht, Maruf Adewole, Udunna Anazodo, Hannah Anderson, Sina Bagheri, Ujjwal Baid, Timothy Bergquist, Austin J. Borja, Evan Calabrese, Verena Chung, Gian-Marco Conte, Farouk Dako, James Eddy, Ivan Ezhov, Ariana Familiar, Keyvan Farahani, Shuvanjan Haldar, Juan Eugenio Iglesias, Anastasia Janas , et al. (48 additional authors not shown)

Abstract: Pediatric tumors of the central nervous system are the most common cause of cancer-related death in children. The five-year survival rate for high-grade gliomas in children is less than 20\%. Due to their rarity, the diagnosis of these entities is often delayed, their treatment is mainly based on historic treatment concepts, and clinical trials require multi-institutional collaborations. The MICCA… ▽ More Pediatric tumors of the central nervous system are the most common cause of cancer-related death in children. The five-year survival rate for high-grade gliomas in children is less than 20\%. Due to their rarity, the diagnosis of these entities is often delayed, their treatment is mainly based on historic treatment concepts, and clinical trials require multi-institutional collaborations. The MICCAI Brain Tumor Segmentation (BraTS) Challenge is a landmark community benchmark event with a successful history of 12 years of resource creation for the segmentation and analysis of adult glioma. Here we present the CBTN-CONNECT-DIPGR-ASNR-MICCAI BraTS-PEDs 2023 challenge, which represents the first BraTS challenge focused on pediatric brain tumors with data acquired across multiple international consortia dedicated to pediatric neuro-oncology and clinical trials. The BraTS-PEDs 2023 challenge focuses on benchmarking the development of volumentric segmentation algorithms for pediatric brain glioma through standardized quantitative performance evaluation metrics utilized across the BraTS 2023 cluster of challenges. Models gaining knowledge from the BraTS-PEDs multi-parametric structural MRI (mpMRI) training data will be evaluated on separate validation and unseen test mpMRI dataof high-grade pediatric glioma. The CBTN-CONNECT-DIPGR-ASNR-MICCAI BraTS-PEDs 2023 challenge brings together clinicians and AI/imaging scientists to lead to faster development of automated segmentation techniques that could benefit clinical trials, and ultimately the care of children with brain tumors. △ Less

Submitted 23 May, 2024; v1 submitted 26 May, 2023; originally announced May 2023.

Prediction of equivalent sand-grain size and identification of drag-relevant scales of roughness -- a data driven approach

Authors: Jiasheng Yang, Alexander Stroh, Sangseung Lee, Shervin Bagheri, Bettina Frohnapfel, Pourya Forooghi

Abstract: The purpose of the present work is to examine two possibilities; firstly, predicting equivalent sand-grain roughness size $k_s$ based on the roughness height probability density function and power spectrum leveraging machine learning as a regression tool, and secondly, extracting information about relevance of different roughness scales to skin-friction drag by interpreting the output of the train… ▽ More The purpose of the present work is to examine two possibilities; firstly, predicting equivalent sand-grain roughness size $k_s$ based on the roughness height probability density function and power spectrum leveraging machine learning as a regression tool, and secondly, extracting information about relevance of different roughness scales to skin-friction drag by interpreting the output of the trained data-driven model. The model is an ensemble neural network consisting of 50 deep neural networks. The data for the training of the model is obtained from direct numerical simulations (DNSs) of turbulent flow in plane channels over 85 irregular multi-scale roughness samples at friction Reynolds number Re$_τ=800$. The 85 roughness samples are selected from a repository of 4200 samples, covering a wide parameter space, through an active learning (AL) framework. The selection is made in several iterations, based on the informativeness of samples in the repository, quantified by the variance of ENN predictions. This AL framework aims to maximize the generalizability of the predictions with a certain amount of data. This is examined using three different testing data sets with different types of roughness, including 21 surfaces from the literature. The model yields an overall mean error of 5\% to 10\% on different testing data sets. Subsequently, a data interpretation technique, known as layer-wise relevance propagation, is applied to measure the contributions of different roughness wave-lengths to the predicted $k_s$. High-pass filtering is then applied to the roughness PS to exclude the wave-numbers identified as drag-irrelevant. The filtered rough surfaces are investigated using DNS, and it is demonstrated that, despite significant impact of filtering on the roughness topographical appearance and statistics, the skin-friction coefficient of the original roughness is successfully preserved. △ Less

Submitted 28 August, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 24 pages, 11 figures, submitted to JFM

Turbulent flows over porous lattices: alteration of near-wall turbulence and pore-flow amplitude modulation

Authors: Seyed Morteza Habibi Khorasani, Mitul Luhar, Shervin Bagheri

Abstract: Turbulent flows over porous lattices consisting of rectangular cuboid pores are investigated using scale-resolving direct numerical simulations. Beyond a certain threshold which is primarily determined by the wall-normal Darcy permeability, ${K_y}^+$, near-wall turbulence transitions from its canonical regime, marked by the presence of streak-like structures, to another marked by the presence of s… ▽ More Turbulent flows over porous lattices consisting of rectangular cuboid pores are investigated using scale-resolving direct numerical simulations. Beyond a certain threshold which is primarily determined by the wall-normal Darcy permeability, ${K_y}^+$, near-wall turbulence transitions from its canonical regime, marked by the presence of streak-like structures, to another marked by the presence of spanwise coherent structures reminiscent of the Kelvin-Helmholtz (K-H) type of instability. This permeability threshold agrees well with that previously established in studies where permeable-wall boundary conditions had been used as surrogates for a porous substrate. None of the substrates investigated demonstrate any drag reduction relative to smooth-wall turbulent flow. At the permeable surface, a significant component of the flow is that which adheres to the pore geometry and undergoes amplitude modulation (AM). This pore-coherent flow remains notable within the substrates, highlighting the importance of the porous substrate's microstructure when the overlying flow is turbulent, an aspect which cannot be accounted for when using continuum-based approaches to model porous media flows or effective representations such as wall boundary conditions. The severity of the AM is enhanced in the K-H-like regime, which has implications when designing porous substrates for transport processes. This suggests that the surface of the substrate can have a geometry which is different than the rest of it and tailored to influence the overlying flow in a particular way. △ Less

Submitted 14 August, 2023; v1 submitted 18 October, 2022; originally announced October 2022.

Rapid wetting of shear-thinning fluids

Authors: Susumu Yada, Kazem Bazesefidpar, Outi Tammisola, Gustav Amberg, Shervin Bagheri

Abstract: Using experiments and numerical simulations, we investigate the spontaneous spreading of droplets of aqueous glycerol (Newtonian) and aqueous polymer (shear-thinning) solutions on smooth surfaces. We find that in the first millisecond the spreading of the shear-thinning solutions is identical to the spreading of water, regardless of the polymer concentration. In contrast, aqueous glycerol solution… ▽ More Using experiments and numerical simulations, we investigate the spontaneous spreading of droplets of aqueous glycerol (Newtonian) and aqueous polymer (shear-thinning) solutions on smooth surfaces. We find that in the first millisecond the spreading of the shear-thinning solutions is identical to the spreading of water, regardless of the polymer concentration. In contrast, aqueous glycerol solutions show a different behavior, namely, significantly slower spreading rate than water. In the initial rapid spreading phase, the dominating forces that can resist the wetting are inertial forces and contact-line friction. For the glycerol solutions, an increase in glycerol concentration effectively increases the contact-line friction, resulting in increased resistance to wetting. For the polymeric solutions, however, an increase in polymer concentration does not modify contact-line friction. As a consequence, the energy dissipation at the contact line can not be controlled by varying the amount of additives for shear-thinning fluids. The reduction of the spreading rate of shear-thinning fluids on smooth surfaces in the rapid wetting regime can only be achieved by increasing solvent viscosity. Our results have implications for phase-change applications where the control of the rapid spreading rate is central, such as anti-icing and soldering △ Less

Submitted 24 August, 2022; originally announced August 2022.

Efficient Signed Graph Sampling via Balancing & Gershgorin Disc Perfect Alignment

Authors: Chinthaka Dinesh, Gene Cheung, Saghar Bagheri, Ivan V. Bajic

Abstract: A basic premise in graph signal processing (GSP) is that a graph encoding pairwise (anti-)correlations of the targeted signal as edge weights is exploited for graph filtering. However, existing fast graph sampling schemes are designed and tested only for positive graphs describing positive correlations. In this paper, we show that for datasets with strong inherent anti-correlations, a suitable gra… ▽ More A basic premise in graph signal processing (GSP) is that a graph encoding pairwise (anti-)correlations of the targeted signal as edge weights is exploited for graph filtering. However, existing fast graph sampling schemes are designed and tested only for positive graphs describing positive correlations. In this paper, we show that for datasets with strong inherent anti-correlations, a suitable graph contains both positive and negative edge weights. In response, we propose a linear-time signed graph sampling method centered on the concept of balanced signed graphs. Specifically, given an empirical covariance data matrix $\bar{\bf{C}}$, we first learn a sparse inverse matrix (graph Laplacian) $\mathcal{L}$ corresponding to a signed graph $\mathcal{G}$. We define the eigenvectors of Laplacian $\mathcal{L}_B$ for a balanced signed graph $\mathcal{G}_B$ -- approximating $\mathcal{G}$ via edge weight augmentation -- as graph frequency components. Next, we choose samples to minimize the low-pass filter reconstruction error in two steps. We first align all Gershgorin disc left-ends of Laplacian $\mathcal{L}_B$ at smallest eigenvalue $λ_{\min}(\mathcal{L}_B)$ via similarity transform $\mathcal{L}_p = §\mathcal{L}_B §^{-1}$, leveraging a recent linear algebra theorem called Gershgorin disc perfect alignment (GDPA). We then perform sampling on $\mathcal{L}_p$ using a previous fast Gershgorin disc alignment sampling (GDAS) scheme. Experimental results show that our signed graph sampling method outperformed existing fast sampling schemes noticeably on various datasets. △ Less

Submitted 15 January, 2023; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: arXiv admin note: text overlap with arXiv:2103.06153

Unsupervised Graph Spectral Feature Denoising for Crop Yield Prediction

Authors: Saghar Bagheri, Chinthaka Dinesh, Gene Cheung, Timothy Eadie

Abstract: Prediction of annual crop yields at a county granularity is important for national food production and price stability. In this paper, towards the goal of better crop yield prediction, leveraging recent graph signal processing (GSP) tools to exploit spatial correlation among neighboring counties, we denoise relevant features via graph spectral filtering that are inputs to a deep learning predictio… ▽ More Prediction of annual crop yields at a county granularity is important for national food production and price stability. In this paper, towards the goal of better crop yield prediction, leveraging recent graph signal processing (GSP) tools to exploit spatial correlation among neighboring counties, we denoise relevant features via graph spectral filtering that are inputs to a deep learning prediction model. Specifically, we first construct a combinatorial graph with edge weights that encode county-to-county similarities in soil and location features via metric learning. We then denoise features via a maximum a posteriori (MAP) formulation with a graph Laplacian regularizer (GLR). We focus on the challenge to estimate the crucial weight parameter $μ$, trading off the fidelity term and GLR, that is a function of noise variance in an unsupervised manner. We first estimate noise variance directly from noise-corrupted graph signals using a graph clique detection (GCD) procedure that discovers locally constant regions. We then compute an optimal $μ$ minimizing an approximate mean square error function via bias-variance analysis. Experimental results from collected USDA data show that using denoised features as input, performance of a crop yield prediction model can be improved noticeably. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Shear-induced droplet mobility within porous surfaces

Authors: Si Suo, Haibo Zhao, Shervin Bagheri, Peng Yu, Yixiang Gan

Abstract: Droplet mobility under shear flows is important in a wide range of engineering applications, e.g., fog collection, and self-cleaning surfaces. For structured surfaces to achieve superhydrophobicity, the removal of stains adhered within the microscale surface features strongly determines the functional performance and durability. In this study, we numerically investigate the shear-induced mobility… ▽ More Droplet mobility under shear flows is important in a wide range of engineering applications, e.g., fog collection, and self-cleaning surfaces. For structured surfaces to achieve superhydrophobicity, the removal of stains adhered within the microscale surface features strongly determines the functional performance and durability. In this study, we numerically investigate the shear-induced mobility of the droplet trapped within porous surfaces. Through simulations covering a wide range of flow conditions and porous geometries, three droplet mobility modes are identified, i.e., the stick-slip, crossover, and slugging modes. To quantitatively characterise the droplet dynamics, we propose a droplet-scale capillary number that considers the driving force and capillary resistance. By comparing against the simulation results, the proposed dimensionless number presents a strong correlation with the leftover volume. The dominating mechanisms revealed in this study provide a basis for further research on enhancing surface cleaning and optimising design of anti-fouling surfaces. △ Less

Submitted 29 April, 2022; originally announced April 2022.

Slip of submerged two-dimensional liquid-infused surfaces in the presence of surfactants

Authors: Johan Sundin, Shervin Bagheri

Abstract: Using numerical simulations, we investigate the effects of Marangoni stresses induced by surfactants on the effective slip length of liquid-infused surfaces (LIS) with transverse grooves. The surfactants are assumed soluble in the external liquid shearing the surface and can adsorb onto the interfaces. Two different adsorption models are used: a classical Frumkin model and a more advanced model th… ▽ More Using numerical simulations, we investigate the effects of Marangoni stresses induced by surfactants on the effective slip length of liquid-infused surfaces (LIS) with transverse grooves. The surfactants are assumed soluble in the external liquid shearing the surface and can adsorb onto the interfaces. Two different adsorption models are used: a classical Frumkin model and a more advanced model that better describes the decrease of surface tension for minuscule concentrations. The simulations show that LIS may face even more severe effects of surfactants than previously investigated superhydrophobic surfaces. Constructing an analytical model for the effective slip length, we can predict the critical surfactant concentration for which the slip length decreases significantly. This analytical model describes both adsorption models of LIS on a unified framework if properly adjusted. We also advance the understanding of when surfactant advection gives rise to highly skewed interfacial concentrations - the so-called partial stagnant cap regime. To a good approximation, this regime can only exist below a specific surfactant concentration given by the Marangoni number and the strength of the surfactants. △ Less

Submitted 1 October, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 34 pages, 13 figures

Droplet impact on asymmetric hydrophobic microstructures

Authors: Susumu Yada, Ugis Lacis, Wouter van der Wijngaart, Fredrik Lundell, Gustav Amberg, Shervin Bagheri

Abstract: Textured hydrophobic surfaces that repel liquid droplets unidirectionally are found in nature such as butterfly wings and ryegrass leaves and are also essential in technological processes such as self-cleaning and anti-icing. However, droplet impact on such surfaces is not fully understood. Here, we study, using a high-speed camera, droplet impact on surfaces with inclined micropillars. We observe… ▽ More Textured hydrophobic surfaces that repel liquid droplets unidirectionally are found in nature such as butterfly wings and ryegrass leaves and are also essential in technological processes such as self-cleaning and anti-icing. However, droplet impact on such surfaces is not fully understood. Here, we study, using a high-speed camera, droplet impact on surfaces with inclined micropillars. We observed directional rebound at high impact speeds on surfaces with dense arrays of pillars. We attribute this asymmetry to the difference in wetting behavior of the structure sidewalls, causing slower retraction of the contact line in the direction against the inclination compared to with the inclination. The experimental observations are complemented with numerical simulations to elucidate the detailed movement of the drops over the pillars. These insights improve our understanding of droplet impact on hydrophobic microstructures and may be a useful for designing structured surfaces for controlling droplet mobility. △ Less

Submitted 9 March, 2022; originally announced March 2022.

Comments: 24 pages, 9 figures

Hybrid Model-based / Data-driven Graph Transform for Image Coding

Authors: Saghar Bagheri, Tam Thuc Do, Gene Cheung, Antonio Ortega

Abstract: Transform coding to sparsify signal representations remains crucial in an image compression pipeline. While the Karhunen-Loève transform (KLT) computed from an empirical covariance matrix $\bar{C}$ is theoretically optimal for a stationary process, in practice, collecting sufficient statistics from a non-stationary image to reliably estimate $\bar{C}$ can be difficult. In this paper, to encode an… ▽ More Transform coding to sparsify signal representations remains crucial in an image compression pipeline. While the Karhunen-Loève transform (KLT) computed from an empirical covariance matrix $\bar{C}$ is theoretically optimal for a stationary process, in practice, collecting sufficient statistics from a non-stationary image to reliably estimate $\bar{C}$ can be difficult. In this paper, to encode an intra-prediction residual block, we pursue a hybrid model-based / data-driven approach: the first $K$ eigenvectors of a transform matrix are derived from a statistical model, e.g., the asymmetric discrete sine transform (ADST), for stability, while the remaining $N-K$ are computed from $\bar{C}$ for performance. The transform computation is posed as a graph learning problem, where we seek a graph Laplacian matrix minimizing a graphical lasso objective inside a convex cone sharing the first $K$ eigenvectors in a Hilbert space of real symmetric matrices. We efficiently solve the problem via augmented Lagrangian relaxation and proximal gradient (PG). Using WebP as a baseline image codec, experimental results show that our hybrid graph transform achieved better energy compaction than default discrete cosine transform (DCT) and better stability than KLT. △ Less

Submitted 2 March, 2022; originally announced March 2022.

Nanoscale sheared droplet: Volume-of-Fluid, phase-field and no-slip molecular dynamics

Authors: Uǧis Lācis, Michele Pellegrino, Johan Sundin, Gustav Amberg, Stephané Zaleski, Berk Hess, Shervin Bagheri

Abstract: The motion of the three-phase contact line between two immiscible fluids and a solid surface arises in a variety of wetting phenomena and technological applications. One challenge in continuum theory is the effective representation of molecular phenomena close to the contact line. Here, we characterize the molecular processes of the moving contact line to assess the accuracy of two different conti… ▽ More The motion of the three-phase contact line between two immiscible fluids and a solid surface arises in a variety of wetting phenomena and technological applications. One challenge in continuum theory is the effective representation of molecular phenomena close to the contact line. Here, we characterize the molecular processes of the moving contact line to assess the accuracy of two different continuum two-phase models. Specifically, molecular dynamics (MD) simulations of a two-dimensional droplet between two moving plates are used to create reference data for different capillary numbers and contact angles. We use a simple-point-charge/extended (SPC/E) water model with particle-mesh Ewald electrostatics treatment. This model provides a very small slip and a more realistic representation of the molecular physics than Lennards-Jones models. The Cahn-Hilliard phase-field model and the Volume-of-Fluid model are calibrated against the drop displacement from MD reference data. It is demonstrated that the calibrated continuum models can accurately capture droplet displacement and droplet breakup for different capillary numbers and contact angles. However, we also observe differences between continuum and atomistic simulations in describing the transient and unsteady droplet behavior, in particular, close to dynamical wetting transitions. The molecular dynamics of the sheared droplet provide insight of the line friction experienced by the advancing and receding contact lines and evidence of large-scale temporal "stick-slip" like oscillations. The presented results will serve as a stepping stone towards developing accurate continuum models for nanoscale hydrodynamics. △ Less

Submitted 12 April, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: 45 pages, 23 figures, 4 tables, 1 supplementary PDF, 2 supplementary movies, 1 supplementary archive, accepted for publication in "Journal of Fluid Mechanics"

Heat transfer increase by convection in liquid-infused surfaces for laminar and turbulent flows

Authors: Johan Sundin, Umberto Ciri, Stefano Leonardi, Marcus Hultmark, Shervin Bagheri

Abstract: Liquid-infused surfaces (LIS) can reduce friction drag in both laminar and turbulent flows. However, the heat transfer properties of such multi-phase surfaces have still not been investigated to a large extent. We use numerical simulations to study conjugate heat transfer of liquid-filled grooves. It is shown that heat transfer can increase for both laminar and turbulent liquid flows due to recirc… ▽ More Liquid-infused surfaces (LIS) can reduce friction drag in both laminar and turbulent flows. However, the heat transfer properties of such multi-phase surfaces have still not been investigated to a large extent. We use numerical simulations to study conjugate heat transfer of liquid-filled grooves. It is shown that heat transfer can increase for both laminar and turbulent liquid flows due to recirculation in the surface texture. For the increase to be substantial, the thermal conductivity of the solid must be similar to the thermal conductivity of the fluids, and the recirculation in the grooves must be sufficiently strong (Péclet number larger than 1). The ratio of the surface cavity to the system height is an upper limit of the direct contribution from the recirculation. While this ratio can be significant for laminar flows in microchannels, it is limited for turbulent flows, where the system scale (e.g. channel height) usually is much larger than the texture height. However, heat transfer enhancement on the order of $10\%$ is observed (with a net drag reduction) in a turbulent channel flow at a friction Reynolds number $Re_τ\approx 180$. It is shown that the turbulent convection in the bulk can be enhanced indirectly from the recirculation in the grooves. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: 23 pages, 14 figures

Final Adaptation Reinforcement Learning for N-Player Games

Authors: Wolfgang Konen, Samineh Bagheri

Abstract: This paper covers n-tuple-based reinforcement learning (RL) algorithms for games. We present new algorithms for TD-, SARSA- and Q-learning which work seamlessly on various games with arbitrary number of players. This is achieved by taking a player-centered view where each player propagates his/her rewards back to previous rounds. We add a new element called Final Adaptation RL (FARL) to all these… ▽ More This paper covers n-tuple-based reinforcement learning (RL) algorithms for games. We present new algorithms for TD-, SARSA- and Q-learning which work seamlessly on various games with arbitrary number of players. This is achieved by taking a player-centered view where each player propagates his/her rewards back to previous rounds. We add a new element called Final Adaptation RL (FARL) to all these algorithms. Our main contribution is that FARL is a vitally important ingredient to achieve success with the player-centered view in various games. We report results on seven board games with 1, 2 and 3 players, including Othello, ConnectFour and Hex. In most cases it is found that FARL is important to learn a near-perfect playing strategy. All algorithms are available in the GBG framework on GitHub. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 23 pages

Predicting drag on rough surfaces by transfer learning of empirical correlations

Authors: Sangseung Lee, Jiasheng Yang, Pourya Forooghi, Alexander Stroh, Shervin Bagheri

Abstract: Recent developments in neural networks have shown the potential of estimating drag on irregular rough surfaces. Nevertheless, the difficulty of obtaining a large high-fidelity dataset to train neural networks is deterring their use in practical applications. In this study, we propose a transfer learning framework to model the drag on irregular rough surfaces even with a limited amount of direct nu… ▽ More Recent developments in neural networks have shown the potential of estimating drag on irregular rough surfaces. Nevertheless, the difficulty of obtaining a large high-fidelity dataset to train neural networks is deterring their use in practical applications. In this study, we propose a transfer learning framework to model the drag on irregular rough surfaces even with a limited amount of direct numerical simulations. We show that transfer learning of empirical correlations, reported in the literature, can significantly improve the performance of neural networks for drag prediction. This is because empirical correlations include `approximate knowledge' of the drag dependency in high-fidelity physics. The `approximate knowledge' allows neural networks to learn the surface statistics known to affect drag more efficiently. The developed framework can be applied to applications where acquiring a large dataset is difficult, but empirical correlations have been reported. △ Less

Submitted 30 September, 2021; v1 submitted 10 June, 2021; originally announced June 2021.

Journal ref: Journal of Fluid Mechanics (2022)

The Brinkman viscosity for porous media exposed to a free flow

Authors: Aidan Rinehart, Uǧis Lācis, Shervin Bagheri

Abstract: The Brinkman equation has found great popularity in modelling the interfacial flow between free fluid and a porous medium. However, it is still unclear how to determine an appropriate effective Brinkman viscosity without resolving the flow at the pore scale. Here, we propose to determine the Brinkman viscosity for rough porous media from the interface slip length and the interior permeability. Bot… ▽ More The Brinkman equation has found great popularity in modelling the interfacial flow between free fluid and a porous medium. However, it is still unclear how to determine an appropriate effective Brinkman viscosity without resolving the flow at the pore scale. Here, we propose to determine the Brinkman viscosity for rough porous media from the interface slip length and the interior permeability. Both slip and permeability can be determined from unit-cell analysis, thus enabling an a priori estimate of the effective viscosity. By comparing the velocity distribution in the porous material predicted from the Brinkman equation with that obtained from pore-scale resolved simulations, we show that modelling errors are $\sim 10\%$ and not larger than $40\%$. We highlight the physical origins of the obtained errors and demonstrate that the Brinkman model can be much more accurate for irregular porous structures. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 21 page, 14 figures, 1 table. Under consideration for the Journal of Fluid Mechanics

Near-wall turbulence alteration with the Transpiration-Resistance Model

Authors: Seyed Morteza Habibi Khorasani, Uǧis Lācis, Simon Pasche, Marco Edoardo Rosti, Shervin Bagheri

Abstract: A set of boundary conditions called the Transpiration-Resistance Model (TRM) are investigated in altering near-wall turbulence. The TRM has been previously proposed by \citet{Lacis2020} as a means of representing the net effect of surface micro-textures on their overlying bulk flows. It encompasses conventional Navier-slip boundary conditions relating the streamwise and spanwise velocities to thei… ▽ More A set of boundary conditions called the Transpiration-Resistance Model (TRM) are investigated in altering near-wall turbulence. The TRM has been previously proposed by \citet{Lacis2020} as a means of representing the net effect of surface micro-textures on their overlying bulk flows. It encompasses conventional Navier-slip boundary conditions relating the streamwise and spanwise velocities to their respective shears through the slip lengths $\ell_x$ and $\ell_z$. In addition, it features a transpiration condition accounting for the changes induced in the wall-normal velocity by expressing it in terms of variations of the wall-parallel velocity shears through the transpiration lengths $m_x$ and $m_z$. Greater levels of drag increase occur when more transpiration takes place at the boundary plane, with turbulent transpiration being predominately coupled to the spanwise shear component for canonical near-wall turbulence. The TRM can reproduce the effect of a homogeneous and structured roughness up to ${k^+}\,{\approx18}$. In this transitionally rough flow regime, the transpiration lengths of the TRM must be empirically determined. The \emph{transpiration factor} is defined as the product between the slip and transpiration lengths, i.e. $(m\ell)_{x,z}$. This factor contains the compound effect of the wall-parallel velocity occurring at the boundary plane and increased permeability, both of which lead to the transport of momentum in the wall-normal direction. A linear relation between the transpiration factor and the roughness function is observed for regularly textured surfaces in the transitionally rough regime of turbulence. The relations obtained between the transpiration factor and the roughness function show that such effective flow quantities can be suitable measures for characterizing rough surfaces in this flow regime. △ Less

Submitted 23 October, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Journal ref: Journal of Fluid Mechanics , Volume 942 , 10 July 2022 , A45

Roughness on liquid-infused surfaces induced by capillary waves

Authors: Johan Sundin, Stéphane Zaleski, Shervin Bagheri

Abstract: Liquid-infused surfaces (LIS) are a promising technique for reducing friction, fouling and icing in both laminar and turbulent flows. Previous work has demonstrated that these surfaces are susceptible to shear-driven drainage. Here, we report a different failure mode using direct numerical simulations of a turbulent channel flow with liquid-infused longitudinal grooves. When the liquid-liquid surf… ▽ More Liquid-infused surfaces (LIS) are a promising technique for reducing friction, fouling and icing in both laminar and turbulent flows. Previous work has demonstrated that these surfaces are susceptible to shear-driven drainage. Here, we report a different failure mode using direct numerical simulations of a turbulent channel flow with liquid-infused longitudinal grooves. When the liquid-liquid surface tension is small and/or grooves are wide, we observe traveling-wave perturbations on the interface with amplitudes larger than the viscous sublayer of the turbulent flow. These capillary waves induce a roughness effect that increases drag. The generation mechanism of these waves is explained using the theory developed by Miles for gravity waves. Energy is transferred from the turbulent flow to the LIS provided that there is a negative curvature of the mean flow at the critical layer. Given the groove width, the Weber number and an estimate of the friction Reynolds number, we provide relations to determine whether a LIS behaves as a smooth or rough surface in a turbulent flow. △ Less

Submitted 12 April, 2021; v1 submitted 2 February, 2021; originally announced February 2021.

Journal ref: J. Fluid Mech. 915 (2021) R6

Learning Sparse Graph Laplacian with K Eigenvector Prior via Iterative GLASSO and Projection

Authors: Saghar Bagheri, Gene Cheung, Antonio Ortega, Fen Wang

Abstract: Learning a suitable graph is an important precursor to many graph signal processing (GSP) pipelines, such as graph spectral signal compression and denoising. Previous graph learning algorithms either i) make some assumptions on connectivity (e.g., graph sparsity), or ii) make simple graph edge assumptions such as positive edges only. In this paper, given an empirical covariance matrix $\bar{C}$ co… ▽ More Learning a suitable graph is an important precursor to many graph signal processing (GSP) pipelines, such as graph spectral signal compression and denoising. Previous graph learning algorithms either i) make some assumptions on connectivity (e.g., graph sparsity), or ii) make simple graph edge assumptions such as positive edges only. In this paper, given an empirical covariance matrix $\bar{C}$ computed from data as input, we consider a structural assumption on the graph Laplacian matrix $L$: the first $K$ eigenvectors of $L$ are pre-selected, e.g., based on domain-specific criteria, such as computation requirement, and the remaining eigenvectors are then learned from data. One example use case is image coding, where the first eigenvector is pre-chosen to be constant, regardless of available observed data. We first prove that the subspace of symmetric positive semi-definite (PSD) matrices $H_{u}^+$ with the first $K$ eigenvectors being $\{u_k\}$ in a defined Hilbert space is a convex cone. We then construct an operator to project a given positive definite (PD) matrix $L$ to $H_{u}^+$, inspired by the Gram-Schmidt procedure. Finally, we design an efficient hybrid graphical lasso/projection algorithm to compute the most suitable graph Laplacian matrix $L^* \in H_{u}^+$ given $\bar{C}$. Experimental results show that given the first $K$ eigenvectors as a prior, our algorithm outperforms competing graph learning schemes using a variety of graph comparison metrics. △ Less

Submitted 18 February, 2021; v1 submitted 25 October, 2020; originally announced October 2020.

Droplet impact on surfaces with asymmetric microscopic features

Authors: Susumu Yada, Blandine Allais, Wouter van der Wijngaart, Fredrik Lundell, Gustav Amberg, Shervin Bagheri

Abstract: The impact of liquid drops on a rigid surface is central in cleaning, cooling and coating processes in both nature and industrial applications. However, it is not clear how details of pores, roughness and texture on the solid surface influence the initial stages of the impact dynamics. Here, we experimentally study drop impacting at low velocities onto surfaces textured with asymmetric (tilted) ri… ▽ More The impact of liquid drops on a rigid surface is central in cleaning, cooling and coating processes in both nature and industrial applications. However, it is not clear how details of pores, roughness and texture on the solid surface influence the initial stages of the impact dynamics. Here, we experimentally study drop impacting at low velocities onto surfaces textured with asymmetric (tilted) ridges. We define the line-friction capillary number $Ca_f={μ_f V_0}/σ$ (where $μ_f$, $V_0$ and $σ$ are the line friction, impact velocity and surface tension, respectively) as a measure of the importance of the topology of surface textures for the dynamics of droplet impact. We show that when $Ca_f \ll 1$, the contact line speed in the direction against the inclination of the ridges is set by line-friction, whereas in the direction with inclination the contact line is pinned at acute corners of the ridge. When $Ca_f \sim 1$, the pinning is only temporary until the liquid-vapor interface reaches to the next ridge where a new contact line is formed. Finally, when $Ca_f\gg 1$, the geometric details of non-smooth surfaces play little role. △ Less

Submitted 4 February, 2021; v1 submitted 3 July, 2020; originally announced July 2020.

Comments: 12 pages, 7 figures, under consideration for Physical Review Fluids. Revised February 04 2021

Steady moving contact line of water over a no-slip substrate

Authors: Uǧis Lācis, Petter Johansson, Tomas Fullana, Berk Hess, Gustav Amberg, Shervin Bagheri, Stephané Zaleski

Abstract: The movement of the triple contact line plays a crucial role in many applications such as ink-jet printing, liquid coating and drainage (imbibition) in porous media. To design accurate computational tools for these applications, predictive models of the moving contact line are needed. However, the basic mechanisms responsible for movement of the triple contact line are not well understood but stil… ▽ More The movement of the triple contact line plays a crucial role in many applications such as ink-jet printing, liquid coating and drainage (imbibition) in porous media. To design accurate computational tools for these applications, predictive models of the moving contact line are needed. However, the basic mechanisms responsible for movement of the triple contact line are not well understood but still debated. We investigate the movement of the contact line between water, vapour and a silica-like solid surface under steady conditions in low capillary number regime. We use molecular dynamics (MD) with an atomistic water model to simulate a nanoscopic drop between two moving plates. We include hydrogen bonding between the water molecules and the solid substrate, which leads to a sub-molecular slip length. We benchmark two continuum methods, the Cahn-Hilliard phase-field (PF) model and a volume-of-fluid (VOF) model, against MD results. We show that both continuum models can reproduce the statistical measures obtained from MD reasonably well, with a trade-off in accuracy. We demonstrate the importance of the phase field mobility parameter and the local slip length in accurately modelling the moving contact line. △ Less

Submitted 5 April, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 26 pages, 13 figures, 4 tables, post-review version, under consideration for publication in "European Physical Journal Special Topics"

Lift induced by slip inhomogeneities in lubricated contacts

Authors: Aidan Rinehart, Uǧis L{ā}cis, Thomas Salez, Shervin Bagheri

Abstract: Lubrication forces depend to a high degree on elasticity, texture, charge, chemistry, and temperature of the interacting surfaces. Therefore, by appropriately designing surface properties, we may tailor lubrication forces to reduce friction, adhesion and wear between sliding surfaces or control repulsion, assembly, and collision of interacting particles. Here, we show that variations of slippage… ▽ More Lubrication forces depend to a high degree on elasticity, texture, charge, chemistry, and temperature of the interacting surfaces. Therefore, by appropriately designing surface properties, we may tailor lubrication forces to reduce friction, adhesion and wear between sliding surfaces or control repulsion, assembly, and collision of interacting particles. Here, we show that variations of slippage on one of the contacting surfaces induce a lift force. We demonstrate the consequences of this force on the mobility of a cylinder traveling near a wall and show the emergence of particle oscillation and migration that would not otherwise occur in the Stokes flow regime.Our study has implications for understanding how inhomogeneous biological interfaces interact with their environment; it also reveals a new method of patterning surfaces for controlling the motion of nearby particles. △ Less

Submitted 30 April, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Journal ref: Phys. Rev. Fluids 5, 082001 (2020)

The isomorphism theorem for linear fragments of continuous logic

Authors: Seyed-Mohammad Bagheri

Abstract: The ultraproduct construction is generalized to $p$-ultramean constructions ($1\leqslant p<\infty$) by replacing ultrafilters with finitely additive measures. These constructions correspond to the linear fragments $\mathscr L^p$ of continuous logic. A powermean variant of Keisler-Shelah isomorphism theorem is proved for $\mathscr L^p$. It is then proved that $\mathscr L^p$-sentences (and their app… ▽ More The ultraproduct construction is generalized to $p$-ultramean constructions ($1\leqslant p<\infty$) by replacing ultrafilters with finitely additive measures. These constructions correspond to the linear fragments $\mathscr L^p$ of continuous logic. A powermean variant of Keisler-Shelah isomorphism theorem is proved for $\mathscr L^p$. It is then proved that $\mathscr L^p$-sentences (and their approximations) are exactly those sentences of continuous logic which are preserved by such constructions. Some other applications are also given. △ Less

Submitted 2 October, 2019; originally announced October 2019.

Continuous integration logic

Authors: Seyed-Mohammad Bagheri, Massoud Pourmahdian

Abstract: We combine continuous and integral logics and found a logical framework for metric measure spaces equipped with a family of continuous relations and operations. We prove the ultraproduct theorem and deduce compactness and other usual results. We also give applications of the compactness theorem in metric measure theory. We combine continuous and integral logics and found a logical framework for metric measure spaces equipped with a family of continuous relations and operations. We prove the ultraproduct theorem and deduce compactness and other usual results. We also give applications of the compactness theorem in metric measure theory. △ Less

Submitted 30 September, 2019; originally announced October 2019.

Droplet leaping governs microstructured surface wetting

Authors: Susumu Yada, Shervin Bagheri, Jonas Hansson, Minh Do-Quang, Fredrik Lundell, Wouter van der Wijngaart, Gustav Amberg

Abstract: Microstructured surfaces that control the direction of liquid transport are not only ubiquitous in nature, but they are also central to technological processes such as fog/water harvesting, oil-water separation, and surface lubrication. However, a fundamental understanding of the initial wetting dynamics of liquids spreading on such surfaces is lacking. Here, we show that three regimes govern micr… ▽ More Microstructured surfaces that control the direction of liquid transport are not only ubiquitous in nature, but they are also central to technological processes such as fog/water harvesting, oil-water separation, and surface lubrication. However, a fundamental understanding of the initial wetting dynamics of liquids spreading on such surfaces is lacking. Here, we show that three regimes govern microstructured surface wetting on short time scales: spread, stick, and contact line leaping. The latter involves establishing a new contact line downstream of the wetting front as the liquid leaps over specific sections of the solid surface. Experimental and numerical investigations reveal how different regimes emerge in different flow directions during wetting of periodic asymmetrically microstructured surfaces. These insights improve our understanding of rapid wetting in droplet impact, splashing, and wetting of vibrating surfaces and may contribute to advances in designing structured surfaces for the mentioned applications. △ Less

Submitted 16 September, 2019; originally announced September 2019.

Comments: 9 pages, 5 figures, under consideration for Soft Matter

Journal ref: Soft Matter, 2019,15, 9528-9536

Higher order homogenized boundary conditions for flows over rough and porous surfaces

Authors: Y. Sudhakar, Ugis Lacis, Simon Pasche, Shervin Bagheri

Abstract: We derive a homogenized macroscopic model for fluid flows over ordered homogeneous porous surfaces. The unconfined free-flow is described by the Navier-Stokes equation, and the Darcy equation governs the seepage flow within the porous domain. Boundary conditions that accurately capture mass and momentum transport at the contact surface between these two domains are derived using the multiscale hom… ▽ More We derive a homogenized macroscopic model for fluid flows over ordered homogeneous porous surfaces. The unconfined free-flow is described by the Navier-Stokes equation, and the Darcy equation governs the seepage flow within the porous domain. Boundary conditions that accurately capture mass and momentum transport at the contact surface between these two domains are derived using the multiscale homogenization technique. In addition to obtaining the generalized version of the widely used Beavers-Joseph slip condition for tangential velocities, the present work provides an accurate formulation for the transpiration velocity and pressure jump at fluid-porous interfaces; these two conditions are essential for handling two- and three-dimensional flows over porous media. All the constitutive parameters appearing in the interface conditions are computed by solving a set of Stokes problems on a much smaller computational domain, making the formulations free of empirical parameters. The tensorial form of the proposed interface conditions makes it possible to handle flows over isotropic, orthotropic, and anisotropic media. A subset of interface conditions, derived for porous media, can be used to model flows over rough walls. The accuracy of the proposed macroscopic model is numerically quantified for flows over porous and rough walls by comparing the results from our homogenized model with those obtained from geometry-resolved microscopic simulations. △ Less

Submitted 14 December, 2020; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: 47 pages, 19 figures, 8 tables, Under consideration for publication in "Transport in porous media"

SACOBRA with Online Whitening for Solving Optimization Problems with High Conditioning

Authors: Samineh Bagheri, Wolfgang Konen, Thomas Bäck

Abstract: Real-world optimization problems often have expensive objective functions in terms of cost and time. It is desirable to find near-optimal solutions with very few function evaluations. Surrogate-assisted optimizers tend to reduce the required number of function evaluations by replacing the real function with an efficient mathematical model built on few evaluated points. Problems with a high conditi… ▽ More Real-world optimization problems often have expensive objective functions in terms of cost and time. It is desirable to find near-optimal solutions with very few function evaluations. Surrogate-assisted optimizers tend to reduce the required number of function evaluations by replacing the real function with an efficient mathematical model built on few evaluated points. Problems with a high condition number are a challenge for many surrogate-assisted optimizers including SACOBRA. To address such problems we propose a new online whitening operating in the black-box optimization paradigm. We show on a set of high-conditioning functions that online whitening tackles SACOBRA's early stagnation issue and reduces the optimization error by a factor between 10 to 1e12 as compared to the plain SACOBRA, though it imposes many extra function evaluations. Covariance matrix adaptation evolution strategy (CMA-ES) has for very high numbers of function evaluations even lower errors, whereas SACOBRA performs better in the expensive setting (less than 1e03 function evaluations). If we count all parallelizable function evaluations (population evaluation in CMA-ES, online whitening in our approach) as one iteration, then both algorithms have comparable strength even on the long run. This holds for problems with dimension D <= 20. △ Less

Submitted 17 April, 2019; originally announced April 2019.

Comments: 20 pages, 10 figures

Modal Analysis of Fluid Flows: Applications and Outlook

Authors: Kunihiko Taira, Maziar S. Hemati, Steven L. Brunton, Yiyang Sun, Karthik Duraisamy, Shervin Bagheri, Scott T. M. Dawson, Chi-An Yeh

Abstract: We present applications of modal analysis techniques to study, model, and control canonical aerodynamic flows. To illustrate how modal analysis techniques can provide physical insights in a complementary manner, we selected four fundamental examples of cylinder wakes, wall-bounded flows, airfoil wakes, and cavity flows. We also offer brief discussions on the outlook for modal analysis techniques,… ▽ More We present applications of modal analysis techniques to study, model, and control canonical aerodynamic flows. To illustrate how modal analysis techniques can provide physical insights in a complementary manner, we selected four fundamental examples of cylinder wakes, wall-bounded flows, airfoil wakes, and cavity flows. We also offer brief discussions on the outlook for modal analysis techniques, in light of rapid developments in data science. △ Less

Submitted 26 July, 2019; v1 submitted 13 March, 2019; originally announced March 2019.

Comments: 37 pages, 19 figures, 2 tables

Transfer of mass and momentum at rough and porous surfaces

Authors: Uǧis Lācis, Y. Sudhakar, Simon Pasche, Shervin Bagheri

Abstract: The surface texture of materials plays a critical role in wettability, turbulence and transport phenomena. In order to design surfaces for these applications, it is desirable to characterise non-smooth and porous materials by their ability to exchange mass and momentum with flowing fluids. While the underlying physics of the tangential (slip) velocity at a fluid-solid interface is well understood,… ▽ More The surface texture of materials plays a critical role in wettability, turbulence and transport phenomena. In order to design surfaces for these applications, it is desirable to characterise non-smooth and porous materials by their ability to exchange mass and momentum with flowing fluids. While the underlying physics of the tangential (slip) velocity at a fluid-solid interface is well understood, the importance and treatment of normal (transpiration) velocity and normal stress is unclear. We show that, when slip velocity varies at an interface above the texture, a non-zero transpiration velocity arises from mass conservation. The ability of a given surface texture to accommodate for a normal velocity of this kind is quantified by a transpiration length. We further demonstrate that normal momentum transfer gives rise to a pressure jump. For a porous material, the pressure jump can be characterised by so called resistance coefficients. By solving five Stokes problems, the introduced measures of slip, transpiration and resistance can be determined for any anisotropic non-smooth surface consisting of regularly repeating geometric patterns. The proposed conditions are a subset of effective boundary conditions derived from formal multi-scale expansion. We validate and demonstrate the physical significance of the effective conditions on two canonical problems -- a lid-driven cavity and a turbulent channel flow, both with non-smooth bottom surfaces. △ Less

Submitted 31 May, 2019; v1 submitted 21 December, 2018; originally announced December 2018.

Comments: 31 page, 12 figures, 5 tables, under consideration for publication in Journal of Fluid Mechanics

Interaction between hairy surfaces and turbulence for different surface time scales

Authors: Johan Sundin, Shervin Bagheri

Abstract: Surfaces with filamentous structures are ubiquitous in nature on many different scales, ranging from forests to micrometer-sized cilia in organs. Hairy surfaces are elastic and porous, and it is not fully understood how they modify turbulence near a wall. The interaction between hairy surfaces and turbulent flows is here investigated numerically in a turbulent channel flow configuration at… ▽ More Surfaces with filamentous structures are ubiquitous in nature on many different scales, ranging from forests to micrometer-sized cilia in organs. Hairy surfaces are elastic and porous, and it is not fully understood how they modify turbulence near a wall. The interaction between hairy surfaces and turbulent flows is here investigated numerically in a turbulent channel flow configuration at $Re_τ\approx 180$. We show that a filamentous bed of a given geometry can modify a turbulent flow very differently depending on the resonance frequency of the surface, which is determined by elasticity and mass of the filaments. Filaments having resonance frequencies lower than the main frequency content of the turbulent wall-shear stress conform to slowly traveling elongated streaky structures, since they are too slow to adapt to fluid forces of higher frequencies. On the other hand, a bed consisting of stiff and low-mass filaments has a high resonance frequency and shows local regions of increased permeability, which results in large entrainment and a vast increase in drag. △ Less

Submitted 16 November, 2018; v1 submitted 25 June, 2018; originally announced June 2018.

Shapley Value Methods for Attribution Modeling in Online Advertising

Authors: Kaifeng Zhao, Seyed Hanif Mahboobi, Saeed R. Bagheri

Abstract: This paper re-examines the Shapley value methods for attribution analysis in the area of online advertising. As a credit allocation solution in cooperative game theory, Shapley value method directly quantifies the contribution of online advertising inputs to the advertising key performance indicator (KPI) across multiple channels. We simplify its calculation by developing an alternative mathematic… ▽ More This paper re-examines the Shapley value methods for attribution analysis in the area of online advertising. As a credit allocation solution in cooperative game theory, Shapley value method directly quantifies the contribution of online advertising inputs to the advertising key performance indicator (KPI) across multiple channels. We simplify its calculation by developing an alternative mathematical formulation. The new formula significantly improves the computational efficiency and therefore extends the scope of applicability. Based on the simplified formula, we further develop the ordered Shapley value method. The proposed method is able to take into account the order of channels visited by users. We claim that it provides a more comprehensive insight by evaluating the attribution of channels at different stages of user conversion journeys. The proposed approaches are illustrated using a real-world online advertising campaign dataset. △ Less

Submitted 15 April, 2018; originally announced April 2018.

Edge state modulation by mean viscosity gradients

Authors: Enrico Rinaldi, Philipp Schlatter, Shervin Bagheri

Abstract: Motivated by the relevance of edge state solutions as mediators of transition, we use direct numerical simulations to study the effect of spatially non-uniform viscosity on their energy and stability in minimal channel flows. What we seek is a theoretical support rooted in a fully non-linear framework that explains the modified threshold for transition to turbulence in flows with temperature-depen… ▽ More Motivated by the relevance of edge state solutions as mediators of transition, we use direct numerical simulations to study the effect of spatially non-uniform viscosity on their energy and stability in minimal channel flows. What we seek is a theoretical support rooted in a fully non-linear framework that explains the modified threshold for transition to turbulence in flows with temperature-dependent viscosity. Consistently over a range of subcritical Reynolds numbers, we find that decreasing viscosity away from the walls weakens the streamwise streaks and the vortical structures responsible for their regeneration. The entire self-sustained cycle of the edge state is maintained on a lower kinetic energy level with a smaller driving force, compared to a flow with constant viscosity. Increasing viscosity away from the walls has the opposite effect. In both cases, the effect is proportional to the strength of the viscosity gradient. The results presented highlight a local shift in the state space of the position of the edge state relative to the laminar attractor with the consequent modulation of its basin of attraction in the proximity of the edge state and of the surrounding manifold. The implication is that the threshold for transition is reduced for perturbations evolving in the neighbourhood of the edge state in case viscosity decreases away from the walls, and vice versa. △ Less

Submitted 14 December, 2017; originally announced December 2017.

Comments: 25 pages, 21 figures. Accepted for publication on the Journal of Fluid Mechanics

Targeted Advertising Based on Browsing History

Authors: Yong Zhang, Hongming Zhou, Nganmeng Tan, Saeed Bagheri, Meng Joo Er

Abstract: Audience interest, demography, purchase behavior and other possible classifications are ex- tremely important factors to be carefully studied in a targeting campaign. This information can help advertisers and publishers deliver advertisements to the right audience group. How- ever, it is not easy to collect such information, especially for the online audience with whom we have limited interaction… ▽ More Audience interest, demography, purchase behavior and other possible classifications are ex- tremely important factors to be carefully studied in a targeting campaign. This information can help advertisers and publishers deliver advertisements to the right audience group. How- ever, it is not easy to collect such information, especially for the online audience with whom we have limited interaction and minimum deterministic knowledge. In this paper, we pro- pose a predictive framework that can estimate online audience demographic attributes based on their browsing histories. Under the proposed framework, first, we retrieve the content of the websites visited by audience, and represent the content as website feature vectors; second, we aggregate the vectors of websites that audience have visited and arrive at feature vectors representing the users; finally, the support vector machine is exploited to predict the audience demographic attributes. The key to achieving good prediction performance is preparing representative features of the audience. Word Embedding, a widely used tech- nique in natural language processing tasks, together with term frequency-inverse document frequency weighting scheme is used in the proposed method. This new representation ap- proach is unsupervised and very easy to implement. The experimental results demonstrate that the new audience feature representation method is more powerful than existing baseline methods, leading to a great improvement in prediction accuracy. △ Less

Submitted 13 November, 2017; originally announced November 2017.

Revenue-based Attribution Modeling for Online Advertising

Authors: Kaifeng Zhao, Seyed Hanif Mahboobi, Saeed Bagheri

Abstract: This paper examines and proposes several attribution modeling methods that quantify how revenue should be attributed to online advertising inputs. We adopt and further develop relative importance method, which is based on regression models that have been extensively studied and utilized to investigate the relationship between advertising efforts and market reaction (revenue). Relative importance m… ▽ More This paper examines and proposes several attribution modeling methods that quantify how revenue should be attributed to online advertising inputs. We adopt and further develop relative importance method, which is based on regression models that have been extensively studied and utilized to investigate the relationship between advertising efforts and market reaction (revenue). Relative importance method aims at decomposing and allocating marginal contributions to the coefficient of determination (R^2) of regression models as attribution values. In particular, we adopt two alternative submethods to perform this decomposition: dominance analysis and relative weight analysis. Moreover, we demonstrate an extension of the decomposition methods from standard linear model to additive model. We claim that our new approaches are more flexible and accurate in modeling the underlying relationship and calculating the attribution values. We use simulation examples to demonstrate the superior performance of our new approaches over traditional methods. We further illustrate the value of our proposed approaches using a real advertising campaign dataset. △ Less

Submitted 17 October, 2017; originally announced October 2017.

A computational continuum model of poroelastic beds

Authors: Uǧis Lācis, Giuseppe Antonio Zampogna, Shervin Bagheri

Abstract: Despite the ubiquity of fluid flows interacting with porous and elastic materials, we lack a validated non-empirical macroscale method for characterizing the flow over and through a poroelastic medium. We propose a computational tool to describe such configurations by deriving and validating a continuum model for the poroelastic bed and its interface with the above free fluid. We show that, using… ▽ More Despite the ubiquity of fluid flows interacting with porous and elastic materials, we lack a validated non-empirical macroscale method for characterizing the flow over and through a poroelastic medium. We propose a computational tool to describe such configurations by deriving and validating a continuum model for the poroelastic bed and its interface with the above free fluid. We show that, using stress continuity condition and slip velocity condition at the interface, the effective model captures the effects of small changes in the microstructure anisotropy correctly and predicts the overall behaviour in a physically consistent and controllable manner. Moreover, we show that the performance of the effective model is accurate by validating with fully microscopic resolved simulations. The proposed computational tool can be used in investigations in a wide range of fields, including mechanical engineering, bio-engineering and geophysics. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: Main paper: 28 pages, 14 figures; Supplementary appendicies: 7 pages; under consideration for publication in Roy. Soc. Proc. A

Passive control of a falling sphere by elliptic-shaped appendages

Authors: Uǧis Lācis, Stefano Olivieri, Andrea Mazzino, Shervin Bagheri

Abstract: The majority of investigations characterizing the motion of single or multiple particles in fluid flows consider canonical body shapes, such as spheres, cylinders, discs, etc. However, protrusions on bodies -- being either as surface imperfections or appendages that serve a function -- are ubiquitous in both nature and applications. In this work, we characterize how the dynamics of a sphere with a… ▽ More The majority of investigations characterizing the motion of single or multiple particles in fluid flows consider canonical body shapes, such as spheres, cylinders, discs, etc. However, protrusions on bodies -- being either as surface imperfections or appendages that serve a function -- are ubiquitous in both nature and applications. In this work, we characterize how the dynamics of a sphere with an axis-symmetric wake is modified in the presence of thin three-dimensional elliptic-shaped protrusions. By investigating a wide range of three-dimensional appendages with different aspect ratios and lengths, we clearly show that the sphere with an appendage may robustly undergo an inverted-pendulum-like (IPL) instability. This means that the position of the appendage placed behind the sphere and aligned with the free-stream direction is unstable, in a similar way that an inverted pendulum is unstable under gravity. Due to this instability, non-trivial forces are generated on the body, leading to turn and drift, if the body is free to fall under gravity. Moreover, we identify the aspect ratio and length of the appendage that induces the largest side force on the sphere, and therefore also the largest drift for a freely falling body. Finally, we explain the physical mechanisms behind these observations in the context of the IPL instability, i.e., the balance between surface area of the appendage exposed to reversed flow in the wake and the surface area of the appendage exposed to fast free-stream flow. △ Less

Submitted 26 January, 2017; v1 submitted 12 October, 2016; originally announced October 2016.

Comments: 16 pages, 13 figures, 2 tables, under consideration for publication in Phys. Rev. Fluids; revision

Journal ref: Phys. Rev. Fluids 2, 033901 (2017)

A framework for computing effective boundary conditions at the interface between free fluid and a porous medium

Authors: Uǧis Lācis, Shervin Bagheri

Abstract: Interfacial boundary conditions determined from empirical or ad-hoc models remain the standard approach to model fluid flows over porous media, even in situations where the topology of the porous medium is known. We propose a non-empirical and accurate method to compute the effective boundary conditions at the interface between a porous surface and an overlying flow. Using multiscale expansion (ho… ▽ More Interfacial boundary conditions determined from empirical or ad-hoc models remain the standard approach to model fluid flows over porous media, even in situations where the topology of the porous medium is known. We propose a non-empirical and accurate method to compute the effective boundary conditions at the interface between a porous surface and an overlying flow. Using multiscale expansion (homogenization) approach, we derive a tensorial generalized version of the empirical condition suggested by Beavers & Joseph (1967). The components of the tensors determining the effective slip velocity at the interface are obtained by solving a set of Stokes equations in a small computational domain near the interface containing both free flow and porous medium. Using the lid-driven cavity flow with a porous bed, we demonstrate that the derived boundary condition is accurate and robust by comparing an effective model to direct numerical simulations. Finally, we provide an open source code that solves the microscale problems and computes the velocity boundary condition without free parameters over any porous bed. △ Less

Submitted 8 November, 2016; v1 submitted 11 April, 2016; originally announced April 2016.

Comments: 23 pages, 6 figures, 2 tables, under consideration for publication in Journal of Fluid Mechanics

Solving the G-problems in less than 500 iterations: Improved efficient constrained optimization by surrogate modeling and adaptive parameter control

Authors: Samineh Bagheri, Wolfgang Konen, Michael Emmerich, Thomas Bäck

Abstract: Constrained optimization of high-dimensional numerical problems plays an important role in many scientific and industrial applications. Function evaluations in many industrial applications are severely limited and no analytical information about objective function and constraint functions is available. For such expensive black-box optimization tasks, the constraint optimization algorithm COBRA was… ▽ More Constrained optimization of high-dimensional numerical problems plays an important role in many scientific and industrial applications. Function evaluations in many industrial applications are severely limited and no analytical information about objective function and constraint functions is available. For such expensive black-box optimization tasks, the constraint optimization algorithm COBRA was proposed, making use of RBF surrogate modeling for both the objective and the constraint functions. COBRA has shown remarkable success in solving reliably complex benchmark problems in less than 500 function evaluations. Unfortunately, COBRA requires careful adjustment of parameters in order to do so. In this work we present a new self-adjusting algorithm SACOBRA, which is based on COBRA and capable to achieve high-quality results with very few function evaluations and no parameter tuning. It is shown with the help of performance profiles on a set of benchmark problems (G-problems, MOPTA08) that SACOBRA consistently outperforms any COBRA algorithm with fixed parameter setting. We analyze the importance of the several new elements in SACOBRA and find that each element of SACOBRA plays a role to boost up the overall optimization performance. We discuss the reasons behind and get in this way a better understanding of high-quality RBF surrogate modeling. △ Less

Submitted 31 December, 2015; originally announced December 2015.

A stable fluid-structure-interaction solver for low-density rigid bodies using the immersed boundary projection method

Authors: Uǧis Lācis, Kunihiko Taira, Shervin Bagheri

Abstract: Dispersion of low-density rigid particles with complex geometries is ubiquitous in both natural and industrial environments. We show that while explicit methods for coupling the incompressible Navier-Stokes equations and Newton's equations of motion are often sufficient to solve for the motion of cylindrical particles with low density ratios, for more complex particles - such as a body with a prot… ▽ More Dispersion of low-density rigid particles with complex geometries is ubiquitous in both natural and industrial environments. We show that while explicit methods for coupling the incompressible Navier-Stokes equations and Newton's equations of motion are often sufficient to solve for the motion of cylindrical particles with low density ratios, for more complex particles - such as a body with a protrusion - they become unstable. We present an implicit formulation of the coupling between rigid body dynamics and fluid dynamics within the framework of the immersed boundary projection method. Similarly to previous work on this method, the resulting matrix equation in the present approach is solved using a block-LU decomposition. Each step of the block-LU decomposition is modified to incorporate the rigid body dynamics. We show that our method achieves second-order accuracy in space and first-order in time (third-order for practical settings), only with a small additional computational cost to the original method. Our implicit coupling yields stable solution for density ratios as low as $10^{-4}$. We also consider the influence of fictitious fluid located inside the rigid bodies on the accuracy and stability of our method. △ Less

Submitted 5 November, 2015; v1 submitted 16 October, 2015; originally announced October 2015.

Comments: Article has undergone final revision stage, minor adjustments in text, fixed citations

Eberlein-Smulian compactness and Kolmogorov extension theorems; a model theoretic approach

Authors: Seyed-Mohammad Bagheri, Karim Khanaki

Abstract: This paper has two parts. First, we complete the proof of the Kolmogorov extension theorem for unbounded random variables using compactness theorem of integral logic which was proved for bounded case in [8]. Second, we give a proof of the Eberlein-Smulian compactness theorem by Ramsey's theorem and point out the correspondence between this theorem and a result in Shelah's classification theory. This paper has two parts. First, we complete the proof of the Kolmogorov extension theorem for unbounded random variables using compactness theorem of integral logic which was proved for bounded case in [8]. Second, we give a proof of the Eberlein-Smulian compactness theorem by Ramsey's theorem and point out the correspondence between this theorem and a result in Shelah's classification theory. △ Less

Submitted 20 September, 2015; v1 submitted 10 September, 2015; originally announced September 2015.

Comments: This paper has been withdrawn by the author due to a crucial sign error

MSC Class: 03B48; 03C20; 03C45; 46B50