Showing 1–38 of 38 results for author: Bhalla, A

A consistent, volume preserving, and adaptive mesh refinement-based framework for modeling non-isothermal gas-liquid-solid flows with phase change

Authors: Ramakrishnan Thirumalaisamy, Amneet Pal Singh Bhalla

Abstract: This work expands on our recently introduced low Mach enthalpy method [1] for simulating the melting and solidification of a phase change material (PCM) alongside (or without) an ambient gas phase. The method captures PCM's volume change (shrinkage or expansion) by accounting for density change-induced flows. We present several improvements to the original work. First, we introduce consistent time… ▽ More This work expands on our recently introduced low Mach enthalpy method [1] for simulating the melting and solidification of a phase change material (PCM) alongside (or without) an ambient gas phase. The method captures PCM's volume change (shrinkage or expansion) by accounting for density change-induced flows. We present several improvements to the original work. First, we introduce consistent time integration schemes for the mass, momentum, and enthalpy equations, which enhance the method stability. Demonstrating the effectiveness of this scheme, we show that a system free of external forces and heat sources can conserve its initial mass, momentum, enthalpy, and phase composition. This allows the system to transition from a non-isothermal, non-equilibrium, phase-changing state to an isothermal, equilibrium state without exhibiting unrealistic behavior. Furthermore, we show that the low Mach enthalpy method accurately simulates thermocapillary flows without introducing spurious phase changes. We propose an analytical model to validate advanced CFD codes for simulating metal manufacturing processes like welding and 3D printing. These processes involve a heat source melting metal or alloy in an inert gas environment. Traditionally, validation relied on manipulating material properties to match complex experiments. Our model uses the Stefan problem with a density jump to provide a straightforward method for validating multiphysics simulations involving heat sources and phase changes in three-phase flows.Lastly, we demonstrate the practical utility of the method in modeling porosity defects (gas bubble trapping) during metal solidification. A field extension technique is used to accurately apply surface tension forces in a three-phase flow situation. This is where part of the bubble surface is trapped within the (moving) solidification front. △ Less

Submitted 7 July, 2024; originally announced July 2024.

Polaron-assisted dielectric relaxation processes in donor-doped BaTiO3-based ceramics

Authors: T. H. T. Rosa, M. A. Oliveira, Y. Mendez-Gonzalez, F. Guerrero, R. Guo, A. S. Bhalla, J. D. S. Guerra

Abstract: Perovskite structure materials based on the Ba1-xGdxTiO3 system, where x = 0.001, 0.002, 0.003, 0.004 and 0.005, were prepared via the Pechini chemical synthesis route. The dielectric properties have been analyzed over a wide temperature and frequency range, revealing a significant contribution of the conduction mechanisms in the dielectric response of the studied ceramics. In fact, by using the D… ▽ More Perovskite structure materials based on the Ba1-xGdxTiO3 system, where x = 0.001, 0.002, 0.003, 0.004 and 0.005, were prepared via the Pechini chemical synthesis route. The dielectric properties have been analyzed over a wide temperature and frequency range, revealing a significant contribution of the conduction mechanisms in the dielectric response of the studied ceramics. In fact, by using the Davidson-Cole formalism, the observed electrical behavior was found to be associated with relaxation processes related to intrinsic defects mobility promoted by a thermally-activated polaronic mechanism. The obtained values of the activation energy for the relaxation processes, estimated from the Arrhenius law for the mean relaxation time, revealed a decrease from 0.29 up to 0.21 eV as the Gd-doping concentration increases, which suggests the conduction process to be associated with the polaronic effects due to the coexistence of Ti4+ and Ti3+ ions in the structure. Analysis from the conductivity formalism, by using the Jonscher universal power-law, confirmed the polaron-type conduction mechanism for the dielectric dispersion, as suggested by the dielectric analysis, being the nature of the hopping mechanism governed by small polaron hopping (SPH) charge transport in the studied Ba1-xGdxTiO3 ceramics. △ Less

Submitted 13 June, 2024; v1 submitted 8 June, 2024; originally announced June 2024.

Comments: 26 pages, 9 figures, 1 table

Preventing mass loss in the standard level set method: New insights from variational analyses

Authors: Kaustubh Khedkar, Amirreza Charchi Mamaghani, Pieter Ghysels, Neelesh A. Patankar, Amneet Pal Singh Bhalla

Abstract: For decades, the computational multiphase flow community has grappled with mass loss in the level set method. Numerous solutions have been proposed, from fixing the reinitialization step to combining the level set method with other conservative schemes. However, our work reveals a more fundamental culprit: the smooth Heaviside and delta functions inherent to the standard formulation. Even if reini… ▽ More For decades, the computational multiphase flow community has grappled with mass loss in the level set method. Numerous solutions have been proposed, from fixing the reinitialization step to combining the level set method with other conservative schemes. However, our work reveals a more fundamental culprit: the smooth Heaviside and delta functions inherent to the standard formulation. Even if reinitialization is done exactly, i.e., the zero contour interface remains stationary, the use of smooth functions lead to violation of mass conservation. We propose a novel approach using variational analysis to incorporate a mass conservation constraint. This introduces a Lagrange multiplier that enforces overall mass balance. Notably, as the delta function sharpens, i.e., approaches the Dirac delta limit, the Lagrange multiplier approaches zero. However, the exact Lagrange multiplier method disrupts the signed distance property of the level set function. This motivates us to develop an approximate version of the Lagrange multiplier that preserves both overall mass and signed distance property of the level set function. Our framework even recovers existing mass-conserving level set methods, revealing some inconsistencies in prior analyses. We extend this approach to three-phase flows for fluid-structure interaction (FSI) simulations. We present variational equations in both immersed and non-immersed forms, demonstrating the convergence of the former formulation to the latter when the body delta function sharpens. Rigorous test problems confirm that the FSI dynamics produced by our simple, easy-to-implement immersed formulation with the approximate Lagrange multiplier method are accurate and match state-of-the-art solvers. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Inception Attacks: Immersive Hijacking in Virtual Reality Systems

Authors: Zhuolin Yang, Cathy Yuanchen Li, Arman Bhalla, Ben Y. Zhao, Haitao Zheng

Abstract: Today's virtual reality (VR) systems provide immersive interactions that seamlessly connect users with online services and one another. However, these immersive interfaces also introduce new vulnerabilities, making it easier for users to fall prey to new attacks. In this work, we introduce the immersive hijacking attack, where a remote attacker takes control of a user's interaction with their VR s… ▽ More Today's virtual reality (VR) systems provide immersive interactions that seamlessly connect users with online services and one another. However, these immersive interfaces also introduce new vulnerabilities, making it easier for users to fall prey to new attacks. In this work, we introduce the immersive hijacking attack, where a remote attacker takes control of a user's interaction with their VR system, by trapping them inside a malicious app that masquerades as the full VR interface. Once trapped, all of the user's interactions with apps, services and other users can be recorded and modified without their knowledge. This not only allows traditional privacy attacks but also introduces new interaction attacks, where two VR users encounter vastly different immersive experiences during their interaction. We present our implementation of the immersive hijacking attack on Meta Quest headsets and conduct IRB-approved user studies that validate its efficacy and stealthiness. Finally, we examine effectiveness and tradeoffs of various potential defenses, and propose a multifaceted defense pipeline. △ Less

Submitted 9 September, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 13 pages

A unified constraint formulation of immersed body techniques for coupled fluid-solid motion: continuous equations and numerical algorithms

Authors: Amneet Pal Singh Bhalla, Neelesh A. Patankar

Abstract: Numerical simulation of moving immersed solid bodies in fluids is now practiced routinely following pioneering work of Peskin and co-workers on immersed boundary method (IBM), Glowinski and co-workers on fictitious domain method (FDM), and others on related methods. A variety of variants of IBM and FDM approaches have been published, most of which rely on using a background mesh for the fluid equa… ▽ More Numerical simulation of moving immersed solid bodies in fluids is now practiced routinely following pioneering work of Peskin and co-workers on immersed boundary method (IBM), Glowinski and co-workers on fictitious domain method (FDM), and others on related methods. A variety of variants of IBM and FDM approaches have been published, most of which rely on using a background mesh for the fluid equations and tracking the solid body using Lagrangian points. The key idea that is common to these methods is to assume that the entire fluid-solid domain is a fluid and then to constrain the fluid within the solid domain to move in accordance with the solid governing equations. The immersed solid body can be rigid or deforming. Thus, in all these methods the fluid domain is extended into the solid domain. In this review, we provide a mathemarical perspective of various immersed methods by recasting the governing equations in an extended domain form for the fluid. The solid equations are used to impose appropriate constraints on the fluid that is extended into the solid domain. This leads to extended domain constrained fluid-solid governing equations that provide a unified framework for various immersed body techniques. The unified constrained governing equations in the strong form are independent of the temporal or spatial discretization schemes. We show that particular choices of time stepping and spatial discretization lead to different techniques reported in literature ranging from freely moving rigid to elastic self-propelling bodies. These techniques have wide ranging applications including aquatic locomotion, underwater vehicles, car aerodynamics, and organ physiology (e.g. cardiac flow, esophageal transport, respiratory flows), wave energy convertors, among others. We conclude with comments on outstanding challenges and future directions. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Exploiting Data Hierarchy as a New Modality for Contrastive Learning

Authors: Arjun Bhalla, Daniel Levenson, Jan Bernhard, Anton Abilov

Abstract: This work investigates how hierarchically structured data can help neural networks learn conceptual representations of cathedrals. The underlying WikiScenes dataset provides a spatially organized hierarchical structure of cathedral components. We propose a novel hierarchical contrastive training approach that leverages a triplet margin loss to represent the data's spatial hierarchy in the encoder'… ▽ More This work investigates how hierarchically structured data can help neural networks learn conceptual representations of cathedrals. The underlying WikiScenes dataset provides a spatially organized hierarchical structure of cathedral components. We propose a novel hierarchical contrastive training approach that leverages a triplet margin loss to represent the data's spatial hierarchy in the encoder's latent space. As such, the proposed approach investigates if the dataset structure provides valuable information for self-supervised learning. We apply t-SNE to visualize the resultant latent space and evaluate the proposed approach by comparing it with other dataset-specific contrastive learning methods using a common downstream classification task. The proposed method outperforms the comparable weakly-supervised and baseline methods. Our findings suggest that dataset structure is a valuable modality for weakly-supervised learning. △ Less

Submitted 6 January, 2024; originally announced January 2024.

Applications of machine Learning to improve the efficiency and range of microbial biosynthesis: a review of state-of-art techniques

Authors: Akshay Bhalla, Suraj Rajendran

Abstract: In the modern world, technology is at its peak. Different avenues in programming and technology have been explored for data analysis, automation, and robotics. Machine learning is key to optimize data analysis, make accurate predictions, and hasten/improve existing functions. Thus, presently, the field of machine learning in artificial intelligence is being developed and its uses in varying fields… ▽ More In the modern world, technology is at its peak. Different avenues in programming and technology have been explored for data analysis, automation, and robotics. Machine learning is key to optimize data analysis, make accurate predictions, and hasten/improve existing functions. Thus, presently, the field of machine learning in artificial intelligence is being developed and its uses in varying fields are being explored. One field in which its uses stand out is that of microbial biosynthesis. In this paper, a comprehensive overview of the differing machine learning programs used in biosynthesis is provided, alongside brief descriptions of the fields of machine learning and microbial biosynthesis separately. This information includes past trends, modern developments, future improvements, explanations of processes, and current problems they face. Thus, this paper's main contribution is to distill developments in, and provide a holistic explanation of, 2 key fields and their applicability to improve industry/research. It also highlights challenges and research directions, acting to instigate more research and development in the growing fields. Finally, the paper aims to act as a reference for academics performing research, industry professionals improving their processes, and students looking to understand the concept of machine learning in biosynthesis. △ Less

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

An effective preconditioning strategy for volume penalized incompressible/low Mach multiphase flow solvers

Authors: Ramakrishnan Thirumalaisamy, Kaustubh Khedkar, Pieter Ghysels, Amneet Pal Singh Bhalla

Abstract: The volume penalization (VP) or the Brinkman penalization (BP) method is a diffuse interface method for simulating multiphase fluid-structure interaction (FSI) problems in ocean engineering and/or phase change problems in thermal sciences. The method relies on a penalty factor (which is inversely related to body's permeability $κ$) that must be large to enforce rigid body velocity in the solid dom… ▽ More The volume penalization (VP) or the Brinkman penalization (BP) method is a diffuse interface method for simulating multiphase fluid-structure interaction (FSI) problems in ocean engineering and/or phase change problems in thermal sciences. The method relies on a penalty factor (which is inversely related to body's permeability $κ$) that must be large to enforce rigid body velocity in the solid domain. When the penalty factor is large, the discrete system of equations becomes stiff and difficult to solve numerically. In this paper, we propose a projection method-based preconditioning strategy for solving volume penalized (VP) incompressible and low-Mach Navier-Stokes equations. The projection preconditioner enables the monolithic solution of the coupled velocity-pressure system in both single phase and multiphase flow settings. In this approach, the penalty force is treated implicitly, which is allowed to take arbitrary large values without affecting the solver's convergence rate or causing numerical stiffness/instability. It is made possible by including the penalty term in the pressure Poisson equation. Solver scalability under grid refinement is demonstrated. A manufactured solution in a single phase setting is used to determine the spatial accuracy of the penalized solution. Second-order pointwise accuracy is achieved for both velocity and pressure solutions. Two multiphase fluid-structure interaction (FSI) problems from the ocean engineering literature are also simulated to evaluate the solver's robustness and performance. The proposed solver allows us to investigate the effect of $κ$ on the motion of the contact line over the surface of the immersed body. It also allows us to investigate the dynamics of the free surface of a solidifying metal △ Less

Submitted 23 June, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Memory-Efficient Solutions to Large-Graph MST Problems

Authors: Arjun Bhalla

Abstract: Minimum Spanning Trees are a well-studied subset of graph problems. While classical algorithms have existed to solve these problems for decades, new variations and application areas are constantly being discovered. When dealing with large graph problems, however, memory constraints can often be limiting, especially when using these classical methods in memory restricted environments. In this work,… ▽ More Minimum Spanning Trees are a well-studied subset of graph problems. While classical algorithms have existed to solve these problems for decades, new variations and application areas are constantly being discovered. When dealing with large graph problems, however, memory constraints can often be limiting, especially when using these classical methods in memory restricted environments. In this work, we propose an augmentation of Prim's algorithm that can be empirically shown to solve MST problems with a reduction in auxiliary memory usage of over 90%, and a margin of error of less than 0.3%. △ Less

Submitted 26 December, 2023; v1 submitted 8 May, 2023; originally announced May 2023.

A low Mach enthalpy method to model non-isothermal gas-liquid-solid flows with melting and solidification

Authors: Ramakrishnan Thirumalaisamy, Amneet Pal Singh Bhalla

Abstract: Modeling phase change problems numerically is vital for understanding many natural (e.g., ice formation, steam generation) and engineering processes (e.g., casting, welding, additive manufacturing). Almost all phase change materials (PCMs) exhibit density/volume changes during melting, solidification, boiling, or condensation, causing additional fluid flow during this transition. Most numerical wo… ▽ More Modeling phase change problems numerically is vital for understanding many natural (e.g., ice formation, steam generation) and engineering processes (e.g., casting, welding, additive manufacturing). Almost all phase change materials (PCMs) exhibit density/volume changes during melting, solidification, boiling, or condensation, causing additional fluid flow during this transition. Most numerical works consider only two phase flows (either solid-liquid or liquid-gas) for modeling phase change phenomena and some also neglect volume/density change of PCMs in the models. This paper presents a novel low Mach enthalpy method for simulating solidification and melting problems with variable thermophysical properties, including density. Additionally, this formulation allows coupling a solid-liquid PCM with a gas phase in order to simulate the free surface dynamics of PCMs undergoing melting and solidification. We revisit the two-phase Stefan problem involving a density jump between two material phases. We propose a possible means to include the kinetic energy jump in the Stefan condition while still allowing for an analytical solution. The new low Mach enthalpy method is validated against analytical solutions for a PCM undergoing a large density change during its phase transition. Additionally, a few simple sanity checks are proposed to benchmark computational fluid dynamics (CFD) algorithms that aim to capture the volume change effects of PCMs. △ Less

Submitted 27 July, 2023; v1 submitted 15 January, 2023; originally announced January 2023.

How wavelength affects the hydrodynamic performance of two accelerating mirror-symmetric slender swimmers

Authors: Zhonglu Lin, Dongfang Liang, Amneet Pal Singh Bhalla, Ahmed A. Sheikh Al-Shabab, Martin Skote, Wei Zheng, Yu Zhang

Abstract: Fish schools are capable of simultaneous linear acceleration. To reveal the underlying hydrodynamic mechanism, we numerically investigate how Reynolds number $ Re = 1000 - 2000 $, Strouhal number $ St = 0.2 - 0.7 $ and wavelength $ λ= 0.5 - 2 $ affect the mean net thrust and net propulsive efficiency of two side-by-side hydrofoils undulating in anti-phase. In total, $ 550 $ cases are simulated usi… ▽ More Fish schools are capable of simultaneous linear acceleration. To reveal the underlying hydrodynamic mechanism, we numerically investigate how Reynolds number $ Re = 1000 - 2000 $, Strouhal number $ St = 0.2 - 0.7 $ and wavelength $ λ= 0.5 - 2 $ affect the mean net thrust and net propulsive efficiency of two side-by-side hydrofoils undulating in anti-phase. In total, $ 550 $ cases are simulated using immersed boundary method. The thrust increases significantly with wavelength and Strouhal number, yet only slightly with the Reynolds number. We apply a symbolic regression algorithm to formulate this relationship. Furthermore, we find that mirror-symmetric schooling can achieve a \textit{net} thrust more than ten times that of a single swimmer, especially at low Reynolds numbers. The highest efficiency is obtained at $ St = 0.5 $ and $ λ= 1.2 $, where $ St $ is consistent with that observed in the linear-accelerating natural swimmers, \eg Crevalle jack. Six distinct flow structures are identified. The highest thrust corresponds to an asymmetric flow pattern, whereas the highest efficiency occurs when the flow is symmetric with converging vortex streets. △ Less

Submitted 22 July, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: This paper has been accepted by Physics of Fluids. This is the accepted version

How swimming style affects schooling of two fish-like wavy hydrofoils

Authors: Zhonglu Lin, Amneet Pal Singh Bhalla, Boyce E. Griffith, Zi Sheng, Hongquan Li, Dongfang Liang, Yu Zhang

Abstract: Fish swimming style and schooling behaviour are two important aspects of underwater locomotion. How swimming style affects fish schooling is investigated by a representative problem setup consisting of two tethered NACA0012 hydrofoils undulating in a free stream flow at various wavelengths $ λ^* = 0.8 - 8 $, front-back distance $ D = 0,\ 0.25,\ 0.5,\ 0.75 $, phase difference… ▽ More Fish swimming style and schooling behaviour are two important aspects of underwater locomotion. How swimming style affects fish schooling is investigated by a representative problem setup consisting of two tethered NACA0012 hydrofoils undulating in a free stream flow at various wavelengths $ λ^* = 0.8 - 8 $, front-back distance $ D = 0,\ 0.25,\ 0.5,\ 0.75 $, phase difference $ φ/π= 0,\ 0.5,\ 1,\ 1.5 $, and lateral gap distance $ G = 0.25,\ 0.3,\ 0.35 $ with fixed Reynolds number $ Re = 5000 $, Strouhal number $ St = 0.4 $, and maximum amplitude $ A_{max} = 0.1 $. In total, 336 cases were simulated by open-source software IBAMR based on immersed boundary method. The increase in wavelength generally leads to higher thrust and lateral force, consistent with a previous study regarding a single swimmer. The highest Froude efficiency is obtained at intermediate wavelength $ λ^* = 2 $ for either individual or two foils as a group. At side-by-side arrangement $ D = 0 $, the thrust force upon the two foils can be equivalent, indicating a stable formation. At staggered arrangement, $ D > 0 $, the follower can take significant advantage of the leader in locomotion performance by tuning phase difference. The follower's benefit decreases with relative distance. Various combinations of wavelength and relative distance can lead to distinct flow structures, indicating a tunable stealth capacity of the schooling swimmers. △ Less

Submitted 21 December, 2022; v1 submitted 4 September, 2022; originally announced September 2022.

Comments: 53 pages, 21 figures

An immersed peridynamics model of fluid-structure interaction accounting for material damage and failure

Authors: Keon Ho Kim, Amneet P. S. Bhalla, Boyce E. Griffith

Abstract: This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it de… ▽ More This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it describes the structural motion and resultant forces in Lagrangian form. Coupling between Eulerian and Lagrangian variables is achieved by integral transforms with Dirac delta function kernels, as in standard immersed boundary methods. The major difference between our approach and conventional immersed boundary methods is that we use peridynamics, instead of classical continuum mechanics, to determine the structural forces. We focus on non-ordinary state-based peridynamic material descriptions that allow us to use a constitutive correspondence framework that can leverage well characterized nonlinear constitutive models of soft materials. The convergence and accuracy of our approach are compared to both conventional and immersed finite element methods using widely used benchmark problems of nonlinear incompressible elasticity. We demonstrate that the immersed peridynamics method yields comparable accuracy with similar numbers of structural degrees of freedom for several choices of the size of the peridynamic horizon. We also demonstrate that the method can generate grid-converged simulations of fluid-driven material damage growth, crack formation and propagation, and rupture under large deformations. △ Less

Submitted 1 September, 2023; v1 submitted 28 July, 2022; originally announced July 2022.

A model predictive control (MPC)-integrated multiphase immersed boundary (IB) framework for simulating wave energy converters (WECs)

Authors: Kaustubh Khedkar, Amneet Pal Singh Bhalla

Abstract: In this work, we present a novel MPC-integrated multiphase IB framework that can compute the optimal energy-maximizing control force on-the-fly by dynamically interacting with a high-fidelity numerical wave tank (NWT). Due to the requirement of solving a constrained optimization problem at each time step of the IB simulation, the MPC algorithm utilizes a low-dimensional dynamical model of the devi… ▽ More In this work, we present a novel MPC-integrated multiphase IB framework that can compute the optimal energy-maximizing control force on-the-fly by dynamically interacting with a high-fidelity numerical wave tank (NWT). Due to the requirement of solving a constrained optimization problem at each time step of the IB simulation, the MPC algorithm utilizes a low-dimensional dynamical model of the device that is based on the linear potential theory (LPT). The multiphase IB solver, on the other hand, is based on the high-dimensional fictitious domain Brinkman penalization (FD/BP) method, which fully-resolves the hydrodynamic nonlinearities associated with the wave-structure interaction (WSI). A time-series forecasting auto-regressive model is implemented that predicts wave heights to estimate the future wave excitation/Froude- Krylov forces for the MPC algorithm. Moreover, we also experiment with non-linear Froude-Krylov (NLFK) forces for the first time in an MPC formulation. Under varying sea conditions, the predictions of the MPC-integrated multiphase IB solver are compared to the widely popular LPT-based solvers. Overall, six WSI/MPC solver combinations are compared for a heaving vertical cylinder. We also determine the pathway of energy transfer from the waves to the power take-off (PTO) system and verify the relationships using IB simulations. Additionally, three different sea states are simulated within the IB simulation to test the adaptive capability of MPC for WECs. MPC is demonstrated to adapt to changing sea conditions and find the optimal solution for each sea state. △ Less

Submitted 15 June, 2022; originally announced June 2022.

On immersed boundary kernel functions: a constrained quadratic minimization perspective

Authors: Amneet Pal Singh Bhalla

Abstract: In the immersed boundary (IB) approach to fluid-structure interaction modeling, the coupling between the fluid and structure variables is mediated using a regularized version of Dirac delta function. In the IB literature, the regularized delta functions, also referred to IB kernel functions, are either derived analytically from a set of postulates or computed numerically using the moving least squ… ▽ More In the immersed boundary (IB) approach to fluid-structure interaction modeling, the coupling between the fluid and structure variables is mediated using a regularized version of Dirac delta function. In the IB literature, the regularized delta functions, also referred to IB kernel functions, are either derived analytically from a set of postulates or computed numerically using the moving least squares (MLS) approach. Whereas the analytical derivations typically assume a regular Cartesian grid, the MLS method is a meshless technique that can be used to generate kernel functions on complex domains and unstructured meshes. In this note we take a viewpoint that IB kernel generation, either analytically or via MLS, is a constrained quadratic minimization problem. The extremization of a constrained quadratic function is a broader concept than kernel generation, and there are well-established numerical optimization techniques to solve this problem. For example, we show that the constrained quadratic minimization technique can be used to generate one-sided (anisotropic) IB kernels and/or to bound their values. △ Less

Submitted 23 November, 2021; v1 submitted 22 November, 2021; originally announced November 2021.

A one-sided direct forcing immersed boundary method using moving least squares

Authors: Rahul Bale, Amneet Pal Singh Bhalla, Boyce E. Griffith, Makoto Tsubokura

Abstract: This paper presents a one-sided immersed boundary (IB) method using kernel functions constructed via a moving least squares (MLS) method. The resulting kernels effectively couple structural degrees of freedom to fluid variables on only one side of the fluid-structure interface. This reduces spurious feedback forcing and internal flows that are typically observed in IB models that use isotropic ker… ▽ More This paper presents a one-sided immersed boundary (IB) method using kernel functions constructed via a moving least squares (MLS) method. The resulting kernels effectively couple structural degrees of freedom to fluid variables on only one side of the fluid-structure interface. This reduces spurious feedback forcing and internal flows that are typically observed in IB models that use isotropic kernel functions to couple the structure to fluid degrees of freedom on both sides of the interface. The method developed here extends the original MLS methodology introduced by Vanella and Balaras (J Comput Phys, 2009). Prior IB/MLS methods have used isotropic kernel functions that coupled fluid variables on both sides of the boundary to the interfacial degrees of freedom. The original IB/MLS approach converts the cubic spline weights typically employed in MLS reconstruction into an IB kernel function that satisfies particular discrete moment conditions. This paper shows that the same approach can be used to construct one-sided kernel functions (kernel functions are referred to as generating functions in the MLS literature). We also examine the performance of the new approach for a family of kernel functions introduced by Peskin. It is demonstrated that the one-sided MLS construction tends to generate non-monotone interpolation kernels with large over- and undershoots. We present two simple weight shifting strategies to construct generating functions that are positive and monotone, which enhances the stability of the resulting IB methodology. Benchmark cases are used to test the order of accuracy and verify the one-sided IB/MLS simulations in both two and three spatial dimensions. This new IB/MLS method is also used to simulate flow over the Ahmed car model, which highlights the applicability of this methodology for modeling complex engineering flows. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: In Press, Journal of Computational Physics (2021)

Critique on "Volume penalization for inhomogeneous Neumann boundary conditions modeling scalar flux in complicated geometry"

Authors: Ramakrishnan Thirumalaisamy, Nishant Nangia, Amneet Pal Singh Bhalla

Abstract: In this letter, we provide counter-examples to demonstrate that it is possible to retain second-order accuracy using Sakurai et al.'s method, even when different flux boundary conditions are imposed on multiple interfaces that do not conform to the Cartesian grid. We consider both continuous and discontinuous indicator functions in our test problems. Both indicator functions yield a similar conver… ▽ More In this letter, we provide counter-examples to demonstrate that it is possible to retain second-order accuracy using Sakurai et al.'s method, even when different flux boundary conditions are imposed on multiple interfaces that do not conform to the Cartesian grid. We consider both continuous and discontinuous indicator functions in our test problems. Both indicator functions yield a similar convergence rate for the problems considered here. We also find that the order of accuracy results for some of the cases presented in Sakurai et al. are not reproducible. This is demonstrated by re-considering the same one- and two-dimensional Poisson problems solved in Sakurai et al. in this letter. The results shown in this letter demonstrate that the spatial order of accuracy of the flux-based VP approach of Sakurai et al. is between $\mathcal{O}$(1) and $\mathcal{O}$(2), and it depends on the underlying problem/model. The spatial order of accuracy cannot simply be deduced a priori based on the imposed flux values, shapes, or grid-conformity of the interfaces, as concluded in Sakurai et al. Further analysis is required to understand the spatial convergence rate of the flux-based VP method. △ Less

Submitted 29 January, 2021; originally announced January 2021.

Handling Neumann and Robin boundary conditions in a fictitious domain volume penalization framework

Authors: Ramakrishnan Thirumalaisamy, Neelesh A. Patankar, Amneet Pal Singh Bhalla

Abstract: Sakurai et al. (J Comput Phys, 2019) presented a flux-based volume penalization (VP) approach for imposing inhomogeneous Neumann boundary conditions on embedded interfaces. The flux-based VP method modifies the diffusion coefficient of the original elliptic (Poisson) equation and uses a flux-forcing function as a source term in the equation to impose the Neumann boundary conditions. As such, the f… ▽ More Sakurai et al. (J Comput Phys, 2019) presented a flux-based volume penalization (VP) approach for imposing inhomogeneous Neumann boundary conditions on embedded interfaces. The flux-based VP method modifies the diffusion coefficient of the original elliptic (Poisson) equation and uses a flux-forcing function as a source term in the equation to impose the Neumann boundary conditions. As such, the flux-based VP method can be easily incorporated into existing fictitious domain codes. Sakurai et al. relied on an analytical construction of flux-forcing functions, which limits the practicality of the approach. Because of the analytical approach taken in the prior work, only (spatially) constant flux values on simple interfaces were considered. In this paper, we present a numerical technique for constructing flux-forcing functions for arbitrarily complex boundaries. The imposed flux values are also allowed to vary spatially in our approach. Furthermore, the flux-based VP method is extended to include (spatially varying) Robin boundary conditions, which makes the flux-based VP method even more general. We consider several two- and three-dimensional test examples to access the spatial accuracy of the numerical solutions. The method is also used to simulate flux-driven thermal convection in a concentric annular domain. We formally derive the flux-based volume penalized Poisson equation satisfying Neumann/Robin boundary condition in strong form; such a derivation was not presented in Sakurai et al., where the equation first appeared for the Neumann problem. The derivation reveals that the flux-based VP approach relies on a surface delta function to impose inhomogeneous Neumann/Robin boundary conditions. However, explicit construction of the delta function is not necessary for the flux-based VP method, which makes it different from other diffuse domain equations presented in the literature. △ Less

Submitted 9 October, 2021; v1 submitted 7 January, 2021; originally announced January 2021.

An adaptive and energy-maximizing control of wave energy converters using extremum-seeking approach

Authors: Luca Parrinello, Panagiotis Dafnakis, Giovanni Bracco, Peiman Naseradinmousavi, Giuliana Mattiazzo, Amneet Pal Singh Bhalla

Abstract: In this paper, we systematically investigate the feasibility of different extremum-seeking (ES) control schemes to improve the conversion efficiency of wave energy converters (WECs). Continuous-time and model-free ES schemes based on the sliding mode, relay, least-squares gradient, self-driving, and perturbation-based methods are used to improve the mean extracted power of a heaving point absorber… ▽ More In this paper, we systematically investigate the feasibility of different extremum-seeking (ES) control schemes to improve the conversion efficiency of wave energy converters (WECs). Continuous-time and model-free ES schemes based on the sliding mode, relay, least-squares gradient, self-driving, and perturbation-based methods are used to improve the mean extracted power of a heaving point absorber subject to regular and irregular waves. This objective is achieved by optimizing the resistive and reactive coefficients of the power take-off (PTO) mechanism using the ES approach. The optimization results are verified against analytical solutions and the extremum of reference-to-output maps. The numerical results demonstrate that except for the self-driving ES algorithm, the other four ES schemes reliably converge for the two-parameter optimization problem, whereas the former is more suitable for optimizing a single-parameter. The results also show that for an irregular sea state, the sliding mode and perturbation-based ES schemes have better convergence to the optimum, in comparison to other ES schemes considered here. The convergence of PTO coefficients towards the performance-optimal values are tested for widely different initial values, in order to avoid bias towards the extremum. We also demonstrate the adaptive capability of ES control by considering a case in which the ES controller adapts to the new extremum automatically amidst changes in the simulated wave conditions. △ Less

Submitted 7 July, 2020; originally announced July 2020.

The inertial sea wave energy converter (ISWEC) technology: device-physics, multiphase modeling and simulations

Authors: Kaustubh Khedkar, Nishant Nangia, Ramakrishnan Thirumalaisamy, Amneet Pal Singh Bhalla

Abstract: In this paper we investigate the dynamics of the inertial wave energy converter (ISWEC) device using fully-resolved computational fluid dynamics (CFD) simulations. Originally prototyped by Polytechnic University of Turin, the device consists of a floating, boat-shaped hull that is slack-moored to the sea bed. Internally, a gyroscopic power take off (PTO) unit converts the wave-induced pitch motion… ▽ More In this paper we investigate the dynamics of the inertial wave energy converter (ISWEC) device using fully-resolved computational fluid dynamics (CFD) simulations. Originally prototyped by Polytechnic University of Turin, the device consists of a floating, boat-shaped hull that is slack-moored to the sea bed. Internally, a gyroscopic power take off (PTO) unit converts the wave-induced pitch motion of the hull into electrical energy. The CFD model is based on the incompressible Navier-Stokes equations and utilizes the fictitious domain Brinkman penalization technique to couple the device physics and water wave dynamics. A numerical wave tank is used to emulate realistic sea operating conditions. A Froude scaling analysis is performed to enable two- and three-dimensional simulations for a scaled-down (1:20) ISWEC model. It is demonstrated that the scaled-down 2D model is sufficient to accurately simulate the hull's pitching motion and to predict the power generation capability of the converter. A systematic parameter study of the ISWEC is conducted, and its optimal performance in terms of power generation is determined based on the hull and gyroscope control parameters. It is demonstrated that the device achieves peak performance when the gyroscope specifications are chosen based on reactive control theory. It is shown that a proportional control of the PTO control torque is required to generate continuous gyroscope precession effects, without which the device generates no power. In an inertial reference frame, it is demonstrated that the yaw and pitch torques acting on the hull are of the same order of magnitude, informing future design investigations of the ISWEC technology. Further, an energy transfer pathway from the water waves to the hull, the hull to the gyroscope, and the gyroscope to the PTO unit is analytically described and numerically verified. △ Less

Submitted 12 May, 2020; originally announced May 2020.

Comments: Figures are compressed to comply with arXiv size requirements

A sharp interface Lagrangian-Eulerian method for rigid-body fluid-structure interaction

Authors: Ebrahim M. Kolahdouz, Amneet P. S. Bhalla, Lawrence N. Scotten, Brent A. Craven, Boyce E. Griffith

Abstract: This paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are demonstrated for a range of benchmark test cases along with large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian method,… ▽ More This paper introduces a sharp interface method to simulate fluid-structure interaction (FSI) involving rigid bodies immersed in viscous incompressible fluids. The capabilities of this methodology are demonstrated for a range of benchmark test cases along with large-scale models of biomedical FSI. The numerical approach developed herein, which we refer to as an immersed Lagrangian-Eulerian method, integrates aspects of partitioned and immersed FSI formulations by solving separate momentum equations for the fluid and solid subdomains, as in a partitioned formulation, while also using non-conforming discretizations of the dynamic fluid and structure regions, as in an immersed formulation. A Dirichlet-Neumann coupling scheme is used, in which the motion of the immersed solid is driven by fluid traction forces evaluated along the fluid-structure interface, and the motion of the fluid along that interface is constrained to match the solid velocity and thereby satisfy the no-slip condition. To develop a practical numerical method, we adopt a penalty approach that approximately imposes the no-slip condition along the fluid-structure interface. Our fluid-structure interaction scheme relies on an immersed interface method for discrete geometries, which enables the accurate determination of both velocities and stresses along complex fluid-structure interfaces. Unlike commonly used partitioned FSI methods, which can suffer from so-called added mass effect instabilities, our methodology retains stability for test cases involving extremely small, nearly equal, equal, and large solid-fluid density ratios without requiring subiterations or complex handling of the pressure. Biomedical FSI demonstration cases are also presented including the dynamics of a bileaflet mechanical heart valve in a pulse duplicator, and transport of blood clots in a patient-averaged anatomical model of the inferior vena cava. △ Less

Submitted 13 April, 2021; v1 submitted 26 March, 2020; originally announced March 2020.

Comparison of wave-structure interaction dynamics of a submerged cylindrical point absorber with three degrees of freedom using potential flow and computational fluid dynamics models

Authors: Panagiotis Dafnakis, Amneet Pal Singh Bhalla, Sergej Antonello Sirigu, Mauro Bonfanti, Giovanni Bracco, Giuliana Mattiazzo

Abstract: In this paper we compare the heave, surge, and pitch dynamics of a submerged cylindrical point absorber, simulated using potential flow and fully-resolved computational fluid dynamics (CFD) models. The potential flow model is based on the time-domain Cummins equation, whereas the CFD model uses the fictitious domain Brinkman penalization (FD/BP) technique. The submerged cylinder is tethered to the… ▽ More In this paper we compare the heave, surge, and pitch dynamics of a submerged cylindrical point absorber, simulated using potential flow and fully-resolved computational fluid dynamics (CFD) models. The potential flow model is based on the time-domain Cummins equation, whereas the CFD model uses the fictitious domain Brinkman penalization (FD/BP) technique. The submerged cylinder is tethered to the seabed using a power take-off (PTO) unit which restrains the heave, surge, and pitch motions of the converter, and absorbs energy from all three modes. It is demonstrated that the potential theory over-predicts the amplitudes of heave and surge motions, whereas it results in an insignificant pitch for a fully-submerged axisymmetric converter. It also under-estimates the slow drift of the buoy, which the CFD model is able to capture reliably. Further, we use fully-resolved CFD simulations to study the performance of a three degrees of freedom (DOF) cylindrical buoy under varying PTO coefficients, mass density of the buoy, and incoming wave heights. It is demonstrated that the PTO coefficients predicted by the linear potential theory are sub-optimal for waves of moderate and high steepness. The wave absorption efficiency improves significantly when higher than the predicted value of the PTO damping is selected. Simulations with different mass densities of the buoy show that converters with low mass densities have an increased tension in their PTO and mooring lines. Moreover, the mass density also influences the range of resonance periods of the device. Finally, simulations with different wave heights show that at higher heights, the wave absorption efficiency of the converter decreases and a large portion of available wave power remains unabsorbed. △ Less

Submitted 26 October, 2020; v1 submitted 14 August, 2019; originally announced August 2019.

Comments: Final published version

Journal ref: Physics of Fluids 32, 093307 (2020)

Simulating water-entry/exit problems using Eulerian-Lagrangian and fully-Eulerian fictitious domain methods within the open-source IBAMR library

Authors: Amneet Pal Singh Bhalla, Nishant Nangia, Panagiotis Dafnakis, Giovanni Bracco, Giuliana Mattiazzo

Abstract: In this paper we employ two implementations of the fictitious domain (FD) method to simulate water-entry and water-exit problems and demonstrate their ability to simulate practical marine engineering problems. In FD methods, the fluid momentum equation is extended within the solid domain using an additional body force that constrains the structure velocity to be that of a rigid body. Using this fo… ▽ More In this paper we employ two implementations of the fictitious domain (FD) method to simulate water-entry and water-exit problems and demonstrate their ability to simulate practical marine engineering problems. In FD methods, the fluid momentum equation is extended within the solid domain using an additional body force that constrains the structure velocity to be that of a rigid body. Using this formulation, a single set of equations is solved over the entire computational domain. The constraint force is calculated in two distinct ways: one using an Eulerian-Lagrangian framework of the immersed boundary (IB) method and another using a fully-Eulerian approach of the Brinkman penalization (BP) method. Both FSI strategies use the same multiphase flow algorithm that solves the discrete incompressible Navier-Stokes system in conservative form. A consistent transport scheme is employed to advect mass and momentum in the domain, which ensures numerical stability of high density ratio multiphase flows involved in practical marine engineering applications. Example cases of a free falling wedge (straight and inclined) and cylinder are simulated, and the numerical results are compared against benchmark cases in literature. △ Less

Submitted 3 July, 2019; v1 submitted 4 April, 2019; originally announced April 2019.

Comments: The current paper builds on arXiv:1901.07892 and re-explains some parts of it for the reader's convenience

A DLM immersed boundary method based wave-structure interaction solver for high density ratio multiphase flows

Authors: Nishant Nangia, Neelesh A. Patankar, Amneet Pal Singh Bhalla

Abstract: We present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI methodology is enabled by combining the distributed Lagrange multiplier (DLM) method of Sharma and Patankar (J Comp Phys, 2005) with a robust level set method based… ▽ More We present a robust immersed boundary (IB) method for high density ratio multiphase flows that is capable of modeling complex wave-structure interaction (WSI) problems arising in marine and coastal engineering applications. The IB/WSI methodology is enabled by combining the distributed Lagrange multiplier (DLM) method of Sharma and Patankar (J Comp Phys, 2005) with a robust level set method based multiphase flow solver. The fluid solver integrates the conservative form of the variable-coefficient incompressible Navier-Stokes equations using a hybrid preconditioner and ensures consistent transport of mass and momentum at a discrete level. The consistent transport scheme preserves the numerical stability of the method in the presence of large density ratios found in problems involving air, water, and an immersed structure. The air-water interface is captured by the level set method on an Eulerian grid, whereas the free-surface piercing immersed structure is represented on a Lagrangian mesh. The fluid-structure interaction (FSI) coupling is mediated via Peskin's regularized delta functions in an implicit manner, which obviates the need to integrate the hydrodynamic stress tensor on the complex surface of the immersed structure. The IB/WSI numerical scheme is implemented within an adaptive mesh refinement (AMR) framework, in which the Lagrangian structure and the air-water interface are embedded on the finest mesh level to capture the thin boundary layers and the vortical structures arising from WSI. We use a well-balanced force discretization for gravity force that eliminates spurious velocity currents in the hydrostatic limit due to density variation in the three phases (air, water and solid). An effective wave generation and absorption technique for a numerical wave tank is presented and used to simulate a benchmark case of water wave distortion due to a submerged structure. △ Less

Submitted 1 September, 2019; v1 submitted 21 January, 2019; originally announced January 2019.

Comments: Figures are compressed to comply with arXiv size requirements

Brownian Dynamics of Fully Confined Suspensions of Rigid Particles Without Green's Functions

Authors: Brennan Sprinkle, Aleksandar Donev, Amneet Pal Singh Bhalla, Neelesh Patankar

Abstract: We introduce a Rigid-Body Fluctuating Immersed Boundary (RB-FIB) method to perform large-scale Brownian dynamics simulations of suspensions of rigid particles in fully confined domains, without any need to explicitly construct Green's functions or mobility operators. In the RB-FIB approach, discretized fluctuating Stokes equations are solved with prescribed boundary conditions in conjunction with… ▽ More We introduce a Rigid-Body Fluctuating Immersed Boundary (RB-FIB) method to perform large-scale Brownian dynamics simulations of suspensions of rigid particles in fully confined domains, without any need to explicitly construct Green's functions or mobility operators. In the RB-FIB approach, discretized fluctuating Stokes equations are solved with prescribed boundary conditions in conjunction with a rigid-body immersed boundary method to discretize arbitrarily-shaped colloidal particles with no-slip or active-slip prescribed on their surface. We design a specialized Split--Euler--Maruyama temporal integrator that uses a combination of random finite differences to capture the stochastic drift appearing in the overdamped Langevin equation. The RB-FIB method presented in this work only solves mobility problems in each time step using a preconditioned iterative solver, and has a computational complexity that scales linearly in the number of particles and fluid grid cells. We demonstrate that the RB-FIB method correctly reproduces the Gibbs-Boltzmann equilibrium distribution, and use the method to examine the time correlation functions for two spheres tightly confined in a cuboid. We model a quasi--two-dimensional colloidal crystal confined in a narrow microchannel and hydrodynamically driven across a commensurate periodic substrate potential mimicking the effect of a corrugated wall. We observe partial and full depinning of the colloidal monolayer from the substrate potential above a certain wall speed, consistent with a transition from static to kinetic friction through propagating kink solitons. Unexpectedly, we find that particles nearest the boundaries of the domain are the first to be displaced, followed by particles in the middle of the domain. △ Less

Submitted 18 January, 2019; originally announced January 2019.

An Immersed Interface Method for Discrete Surfaces

Authors: Ebrahim M. Kolahdouz, Amneet Pal Singh Bhalla, Brent A. Craven, Boyce E. Griffith

Abstract: Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary(IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction(FSI) in such systems, but a difficulty of the IB formulation is that the pressure and viscous stress are generally discontinuous at the interface. The conventional IB method regularize… ▽ More Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary(IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction(FSI) in such systems, but a difficulty of the IB formulation is that the pressure and viscous stress are generally discontinuous at the interface. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces. The immersed interface method(IIM) is an IB-like approach to FSI that sharply imposes stress jump conditions, enabling higher-order accuracy, but prior applications of the IIM have been largely restricted to methods that rely on smooth representations of the interface geometry. This paper introduces an IIM that uses only a C0 representation of the interface,such as those provided by standard nodal Lagrangian FE methods. Verification examples for models with prescribed motion demonstrate that the method sharply resolves stress discontinuities along the IB while avoiding the need for analytic information of the interface geometry. We demonstrate that only the lowest-order jump conditions for the pressure and velocity gradient are required to realize global 2nd-order accuracy. Specifically,we show 2nd-order global convergence rate along with nearly 2nd-order local convergence in the Eulerian velocity, and between 1st-and 2nd-order global convergence rates along with 1st-order local convergence for the Eulerian pressure. We also show 2nd-order local convergence in the interfacial displacement and velocity along with 1st-order local convergence in the fluid traction. As a demonstration of the method's ability to tackle complex geometries,this approach is also used to simulate flow in an anatomical model of the inferior vena cava. △ Less

Submitted 9 June, 2022; v1 submitted 17 December, 2018; originally announced December 2018.

Comments: - Added a non-axisymmetric example (flow within eccentric rotating cylinder in Sec. 4.3) - Added a more in-depth analysis and comparison with a body-fitted approach for the application in Sec. 4.7

A robust incompressible Navier-Stokes solver for high density ratio multiphase flows

Authors: Nishant Nangia, Boyce E. Griffith, Neelesh A. Patankar, Amneet Pal Singh Bhalla

Abstract: This paper presents a robust, adaptive numerical scheme for simulating high density ratio and high shear multiphase flows on locally refined Cartesian grids that adapt to the evolving interfaces and track regions of high vorticity. The algorithm combines the interface capturing level set method with a variable-coefficient incompressible Navier-Stokes solver that is demonstrated to stably resolve m… ▽ More This paper presents a robust, adaptive numerical scheme for simulating high density ratio and high shear multiphase flows on locally refined Cartesian grids that adapt to the evolving interfaces and track regions of high vorticity. The algorithm combines the interface capturing level set method with a variable-coefficient incompressible Navier-Stokes solver that is demonstrated to stably resolve material contrast ratios of up to six orders of magnitude. The discretization approach ensures second-order pointwise accuracy for both velocity and pressure with several physical boundary treatments, including velocity and traction boundary conditions. The paper includes several test cases that demonstrate the order of accuracy and algorithmic scalability of the flow solver. To ensure the stability of the numerical scheme in the presence of high density and viscosity ratios, we employ a consistent treatment of mass and momentum transport in the conservative form of discrete equations. This consistency is achieved by solving an additional mass balance equation, which we approximate via a strong stability preserving Runga-Kutta time integrator and by employing the same mass flux (obtained from the mass equation) in the discrete momentum equation. The scheme uses higher-order total variation diminishing (TVD) and convection-boundedness criterion (CBC) satisfying limiter to avoid numerical fluctuations in the transported density field. The high-order bounded convective transport is done on a dimension-by-dimension basis, which makes the scheme simple to implement. We also demonstrate through several test cases that the lack of consistent mass and momentum transport in non-conservative formulations, which are commonly used in practice, or the use of non-CBC satisfying limiters can yield very large numerical error and very poor accuracy for convection-dominant high density ratio flows. △ Less

Submitted 1 September, 2019; v1 submitted 31 August, 2018; originally announced September 2018.

Comments: Figures are compressed to comply with arXiv size requirements

Transition in swimming direction in a model self-propelled inertial swimmer

Authors: Thomas Dombrowski, Shannon K. Jones, Georgios Katsikis, Amneet Pal Singh Bhalla, Boyce E. Griffith, Daphne Klotsa

Abstract: We propose a reciprocal, self-propelled model swimmer at intermediate Reynolds numbers ($Re$). Our swimmer consists of two unequal spheres that oscillate in antiphase generating nonlinear steady streaming (SS) flows. We show computationally that the SS flows enable the swimmer to propel itself, and also switch direction as $Re$ increases. We quantify the transition in the swimming direction by col… ▽ More We propose a reciprocal, self-propelled model swimmer at intermediate Reynolds numbers ($Re$). Our swimmer consists of two unequal spheres that oscillate in antiphase generating nonlinear steady streaming (SS) flows. We show computationally that the SS flows enable the swimmer to propel itself, and also switch direction as $Re$ increases. We quantify the transition in the swimming direction by collapsing our data on a critical $Re$ and show that the transition in swimming directions corresponds to the reversal of the SS flows. Based on our findings, we propose that SS can be an important physical mechanism for motility at intermediate $Re$. △ Less

Submitted 7 November, 2018; v1 submitted 11 January, 2018; originally announced January 2018.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. Fluids 4, 021101 (2019)

A moving control volume approach to computing hydrodynamic forces and torques on immersed bodies

Authors: Nishant Nangia, Hans Johansen, Neelesh A. Patankar, Amneet Pal Singh Bhalla

Abstract: We present a moving control volume (CV) approach to computing hydrodynamic forces and torques on complex geometries. The method requires surface and volumetric integrals over a simple and regular Cartesian box that moves with an arbitrary velocity to enclose the body at all times. The moving box is aligned with Cartesian grid faces, which makes the integral evaluation straightforward in an immerse… ▽ More We present a moving control volume (CV) approach to computing hydrodynamic forces and torques on complex geometries. The method requires surface and volumetric integrals over a simple and regular Cartesian box that moves with an arbitrary velocity to enclose the body at all times. The moving box is aligned with Cartesian grid faces, which makes the integral evaluation straightforward in an immersed boundary (IB) framework. Discontinuous and noisy derivatives of velocity and pressure at the fluid-structure interface are avoided and far-field (smooth) velocity and pressure information is used. We re-visit the approach to compute hydrodynamic forces and torques through force/torque balance equation in a Lagrangian frame that some of us took in a prior work (Bhalla et al., J Comp Phys, 2013). We prove the equivalence of the two approaches for IB methods, thanks to the use of Peskin's delta functions. Both approaches are able to suppress spurious force oscillations and are in excellent agreement, as expected theoretically. Test cases ranging from Stokes to high Reynolds number regimes are considered. We discuss regridding issues for the moving CV method in an adaptive mesh refinement (AMR) context. The proposed moving CV method is not limited to a specific IB method and can also be used, for example, with embedded boundary methods. △ Less

Submitted 15 June, 2017; v1 submitted 1 April, 2017; originally announced April 2017.

Scalable smoothing strategies for a geometric multigrid method for the immersed boundary equations

Authors: Amneet Pal Singh Bhalla, Matthew G. Knepley, Mark F. Adams, Robert D. Guy, Boyce E. Griffith

Abstract: The immersed boundary (IB) method is a widely used approach to simulating fluid-structure interaction (FSI). Although explicit versions of the IB method can suffer from severe time step size restrictions, these methods remain popular because of their simplicity and generality. In prior work (Guy et al., Adv Comput Math, 2015), some of us developed a geometric multigrid preconditioner for a stable… ▽ More The immersed boundary (IB) method is a widely used approach to simulating fluid-structure interaction (FSI). Although explicit versions of the IB method can suffer from severe time step size restrictions, these methods remain popular because of their simplicity and generality. In prior work (Guy et al., Adv Comput Math, 2015), some of us developed a geometric multigrid preconditioner for a stable semi-implicit IB method under Stokes flow conditions; however, this solver methodology used a Vanka-type smoother that presented limited opportunities for parallelization. This work extends this Stokes-IB solver methodology by developing smoothing techniques that are suitable for parallel implementation. Specifically, we demonstrate that an additive version of the Vanka smoother can yield an effective multigrid preconditioner for the Stokes-IB equations, and we introduce an efficient Schur complement-based smoother that is also shown to be effective for the Stokes-IB equations. We investigate the performance of these solvers for a broad range of material stiffnesses, both for Stokes flows and flows at nonzero Reynolds numbers, and for thick and thin structural models. We show here that linear solver performance degrades with increasing Reynolds number and material stiffness, especially for thin interface cases. Nonetheless, the proposed approaches promise to yield effective solution algorithms, especially at lower Reynolds numbers and at modest-to-high elastic stiffnesses. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Hydrodynamics of Suspensions of Passive and Active Rigid Particles: A Rigid Multiblob Approach

Authors: F. Balboa Usabiaga, B. Kallemov, B. Delmotte, A. Pal Singh Bhalla, B. E. Griffith, A. Donev

Abstract: We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally-resolved spherical blobs constrained to move as a rigid body, to arrive… ▽ More We develop a rigid multiblob method for numerically solving the mobility problem for suspensions of passive and active rigid particles of complex shape in Stokes flow in unconfined, partially confined, and fully confined geometries. As in a number of existing methods, we discretize rigid bodies using a collection of minimally-resolved spherical blobs constrained to move as a rigid body, to arrive at a potentially large linear system of equations for the unknown Lagrange multipliers and rigid-body motions. Here we develop a block-diagonal preconditioner for this linear system and show that a standard Krylov solver converges in a modest number of iterations that is essentially independent of the number of particles. For unbounded suspensions and suspensions sedimented against a single no-slip boundary, we rely on existing analytical expressions for the Rotne-Prager tensor combined with a fast multipole method or a direct summation on a Graphical Processing Unit to obtain an simple yet efficient and scalable implementation. For fully confined domains, such as periodic suspensions or suspensions confined in slit and square channels, we extend a recently-developed rigid-body immersed boundary method to suspensions of freely-moving passive or active rigid particles at zero Reynolds number. We demonstrate that the iterative solver for the coupled fluid and rigid body equations converges in a bounded number of iterations regardless of the system size. We optimize a number of parameters in the iterative solvers and apply our method to a variety of benchmark problems to carefully assess the accuracy of the rigid multiblob approach as a function of the resolution. We also model the dynamics of colloidal particles studied in recent experiments, such as passive boomerangs in a slit channel, as well as a pair of non-Brownian active nanorods sedimented against a wall. △ Less

Submitted 20 November, 2016; v1 submitted 5 February, 2016; originally announced February 2016.

Comments: Under revision in CAMCOS, Nov 2016

Journal ref: Commun. Appl. Math. Comput. Sci. 11 (2016) 217-296

Three possible mechanisms of capacitance enhancement under magnetic field: charge density gradient modulation, electron gas excitation and oscillatory magnetization- polarization coupling

Authors: Ye Wu, Ruyan Guo, Amar Bhalla

Abstract: Three mechanisms of capacitance enhance-ment by a magnetic field have been analyz-ed. Through semiclassical description of charge movement under a magnetic field, it can be shown that the charge density gradient, a term coupled with magnetic field strength, could be used to augment the magnetocapacitance. Also a magnetic field could enhance capacitance on the polarized metallic material surface du… ▽ More Three mechanisms of capacitance enhance-ment by a magnetic field have been analyz-ed. Through semiclassical description of charge movement under a magnetic field, it can be shown that the charge density gradient, a term coupled with magnetic field strength, could be used to augment the magnetocapacitance. Also a magnetic field could enhance capacitance on the polarized metallic material surface due to the electron gas excitation. Finally, a magnetic field could produce oscillation in the capacitance when relating the polarization with the magnetization through the modification of standard free energy model. By deriving these three mechanisms, it can be seen that three approaches are of potential for exploring tunable dielectric materials. △ Less

Submitted 24 November, 2015; originally announced November 2015.

Journal ref: Integrated Ferroelectrics, vol. 166, pp. 232-244, 2015

Giant magnetic pumping of photovoltage and photocurrent using dielectric lossy material

Authors: Ye Wu, Amar Bhalla, Ruyan Guo

Abstract: The enhancement of photovoltage and photocurrent of material is fundamentally significant owing to the many interesting phenomena found and the potential applications. However, vast altering of magnitude-orders of photoelectricity has been technologically challenging. Here we report two dielectric materials Li2ZnSiO4 and Li2SiO3 showing high photovoltage and photocurrent tunability. When magnetic… ▽ More The enhancement of photovoltage and photocurrent of material is fundamentally significant owing to the many interesting phenomena found and the potential applications. However, vast altering of magnitude-orders of photoelectricity has been technologically challenging. Here we report two dielectric materials Li2ZnSiO4 and Li2SiO3 showing high photovoltage and photocurrent tunability. When magnetic field increasing from 0.00015T up to 0.44T,it is found that 3850% of photovoltage tunability and 3841 % of photocurrent tunability in Li2ZnSiO4 , and 132.8% of photovoltage tunability and 132.5% of photocurrent tunability in Li2SiO3. A simple model that considers the effects of spin mixing/interaction and magnetically-tunable charge gradient is used to explain this interaction between magnetic field and the photoelectricity. This result indicates a magnetic approach can be potentially used for energy efficiency. △ Less

Submitted 24 November, 2015; originally announced November 2015.

An Immersed Boundary Method for Rigid Bodies

Authors: B. Kallemov, A. Pal Singh Bhalla, B. E. Griffith, A. Donev

Abstract: We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such pr… ▽ More We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractional-step approaches, by contrast, can impose such constraints only approximately. Imposing these constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. To solve this system efficiently, we develop a preconditioner for the constrained IB formulation that is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskin's IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We apply the method to a number of test problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations. △ Less

Submitted 23 November, 2015; v1 submitted 28 May, 2015; originally announced May 2015.

Comments: Submitted to CAMCOS

Journal ref: Commun. Appl. Math. Comput. Sci. 11 (2016) 79-141

A fully resolved active musculo-mechanical model for esophageal transport

Authors: Wenjun Kou, Amneet Pal Singh Bhalla, Boyce E. Griffith, John E. Pandolfino, Peter J. Kahrilas, Neelesh A. Patankar

Abstract: Esophageal transport is a physiological process that mechanically transports an ingested food bolus from the pharynx to the stomach via the esophagus, a multi-layered muscular tube. This process involves interactions between the bolus, the esophagus, and the neurally coordinated activation of the esophageal muscles. In this work, we use an immersed boundary (IB) approach to simulate peristaltic tr… ▽ More Esophageal transport is a physiological process that mechanically transports an ingested food bolus from the pharynx to the stomach via the esophagus, a multi-layered muscular tube. This process involves interactions between the bolus, the esophagus, and the neurally coordinated activation of the esophageal muscles. In this work, we use an immersed boundary (IB) approach to simulate peristaltic transport in the esophagus. The bolus is treated as a viscous fluid that is actively transported by the muscular esophagus, which is modeled as an actively contracting, fiber-reinforced tube. A simplified version of our model is verified by comparison to an analytic solution to the tube dilation problem. Three different complex models of the multi-layered esophagus, which differ in their activation patterns and the layouts of the mucosal layers, are then extensively tested. To our knowledge, these simulations are the first of their kind to incorporate the bolus, the multi-layered esophagus tube, and muscle activation into an integrated model. Consistent with experimental observations, our simulations capture the pressure peak generated by the muscle activation pulse that travels along the bolus tail. These fully resolved simulations provide new insights into roles of the mucosal layers during bolus transport. In addition, the information on pressure and the kinematics of the esophageal wall due to the coordination of muscle activation is provided, which may help relate clinical data from manometry and ultrasound images to the underlying esophageal motor function. △ Less

Submitted 8 January, 2015; originally announced January 2015.

Journal ref: Journal of Computational Physics, Volume 298, 1 October 2015, Pages 446-465

Separability of drag and thrust in undulatory animals and machines

Authors: Rahul Bale, Anup A. Shirgaonkar, Izaak D. Neveln, Amneet Pal Singh Bhalla, Malcolm A. MacIver, Neelesh A. Patankar

Abstract: For nearly a century, researchers have tried to understand the swimming of aquatic animals in terms of a balance between the forward thrust from swimming movements and drag on the body. Prior approaches have failed to provide a separation of these two forces for undulatory swimmers such as lamprey and eels, where most parts of the body are simultaneously generating drag and thrust. We nonetheless… ▽ More For nearly a century, researchers have tried to understand the swimming of aquatic animals in terms of a balance between the forward thrust from swimming movements and drag on the body. Prior approaches have failed to provide a separation of these two forces for undulatory swimmers such as lamprey and eels, where most parts of the body are simultaneously generating drag and thrust. We nonetheless show that this separation is possible, and delineate its fundamental basis in undulatory swimmers. Our approach unifies a vast diversity of undulatory aquatic animals (anguilliform, sub-carangiform, gymnotiform, bal- istiform, rajiform) and provides design principles for highly agile bioinspired underwater vehicles. This approach has practical utility within biology as well as engineering. It is a predictive tool for use in understanding the role of the mechanics of movement in the evolutionary emergence of morphological features relating to locomotion. For example, we demonstrate that the drag-thrust separation framework helps to predict the observed height of the ribbon fin of electric knifefish, a diverse group of neotropical fishes which are an important model system in sensory neurobiology. We also show how drag-thrust separation leads to models that can predict the swimming velocity of an organism or a robotic vehicle. △ Less

Submitted 6 February, 2014; originally announced February 2014.

Comments: 41 pages, 13 figures, 4 tables

Journal ref: Sci. Rep. 4 , 7329, 2014

A Minimally-Resolved Immersed Boundary Model for Reaction-Diffusion Problems

Authors: A. Pal Singh Bhalla, B. E. Griffith, N. A. Patankar, A. Donev

Abstract: We develop an immersed-boundary approach to modeling reaction-diffusion processes in dispersions of reactive spherical particles, from the diffusion-limited to the reaction-limited setting. We represent each reactive particle with a minimally-resolved 'blob' using many fewer degrees of freedom per particle than standard discretization approaches. More complicated or more highly resolved particle s… ▽ More We develop an immersed-boundary approach to modeling reaction-diffusion processes in dispersions of reactive spherical particles, from the diffusion-limited to the reaction-limited setting. We represent each reactive particle with a minimally-resolved 'blob' using many fewer degrees of freedom per particle than standard discretization approaches. More complicated or more highly resolved particle shapes can be built out of a collection of reactive blobs. We demonstrate numerically that the blob model can provide an accurate representation at low to moderate packing densities of the reactive particles, at a cost not much larger than solving a Poisson equation in the same domain. Unlike multipole expansion methods, our method does not require analytically-computed Green's functions, but rather, computes regularized discrete Green's functions on the fly by using a standard grid-based discretization of the Poisson equation. This allows for great flexibility in implementing different boundary conditions, coupling to fluid flow or thermal transport, and the inclusion of other effects such as temporal evolution and even nonlinearities. We develop multigrid-based preconditioners for solving the linear systems that arise when using implicit temporal discretizations or studying steady states. In the diffusion-limited case the resulting linear system is a saddle-point problem, the efficient solution of which remains a challenge for suspensions of many particles. We validate our method by comparing to published results on reaction-diffusion in ordered and disordered suspensions of reactive spheres. △ Less

Submitted 30 September, 2013; v1 submitted 13 June, 2013; originally announced June 2013.

Comments: Submitted to J.Chem.Phys

Three-dimensional polarization imaging of (Ba,Sr)TiO3:MgO composites

Authors: Patrick Irvin, Ruyan Guo, Amar Bhalla, Jeremy Levy

Abstract: The dielectric tuning and loss of (Ba,Sr)TiO3:MgO bulk composites depend strongly on the connectivity and interaction among the two phases. To investigate this relationship, the polar structure and dynamics of these composites is mapped as a function of space and time using a pair of three-dimensional probes: second-harmonic confocal scanning optical microscopy (SH-CSOM), which maps ferroelectric… ▽ More The dielectric tuning and loss of (Ba,Sr)TiO3:MgO bulk composites depend strongly on the connectivity and interaction among the two phases. To investigate this relationship, the polar structure and dynamics of these composites is mapped as a function of space and time using a pair of three-dimensional probes: second-harmonic confocal scanning optical microscopy (SH-CSOM), which maps ferroelectric polarization in three dimensions, and time-resolved scanning optical microscopy (TR-CSOM), which maps polarization dynamics along two spatial dimensions and one time dimension. SH-CSOM measurements reveal a high degree of homogeneity within the (Ba,Sr)TiO3 regions, while TR-CSOM measurements indicate that topologically connected regions respond with a spatially uniform phase. △ Less

Submitted 13 August, 2012; v1 submitted 28 September, 2004; originally announced September 2004.

Comments: 11 pages, 3 figures

Journal ref: Appl. Phys. Lett. 86, 042903 (2005)