skip to main content

Title: Fluidic Topology Optimization with an Anisotropic Mixture Model
Fluidic devices are crucial components in many industrial applications involving fluid mechanics. Computational design of a high-performance fluidic system faces multifaceted challenges regarding its geometric representation and physical accuracy. We present a novel topology optimization method to design fluidic devices in a Stokes flow context. Our approach is featured by its capability in accommodating a broad spectrum of boundary conditions at the solid-fluid interface. Our key contribution is an anisotropic and differentiable constitutive model that unifies the representation of different phases and boundary conditions in a Stokes model, enabling a topology optimization method that can synthesize novel structures with accurate boundary conditions from a background grid discretization. We demonstrate the efficacy of our approach by conducting several fluidic system design tasks with over four million design parameters.  more » « less
Award ID(s):
2106768 2008584 1763638
Author(s) / Creator(s):
; ; ; ; ; ;
Date Published:
Journal Name:
ACM Transactions on Graphics
Page Range / eLocation ID:
1 to 14
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Simulating the time evolution of Partial Differential Equations (PDEs) of large-scale systems is crucial in many scientific and engineering domains such as fluid dynamics, weather forecasting and their inverse optimization problems. However, both classical solvers and recent deep learning-based surrogate models are typically extremely computationally intensive, because of their local evolution: they need to update the state of each discretized cell at each time step during inference. Here we develop Latent Evolution of PDEs (LE-PDE), a simple, fast and scalable method to accelerate the simulation and inverse optimization of PDEs. LE-PDE learns a compact, global representation of the system and efficiently evolves it fully in the latent space with learned latent evolution models. LE-PDE achieves speedup by having a much smaller latent dimension to update during long rollout as compared to updating in the input space. We introduce new learning objectives to effectively learn such latent dynamics to ensure long-term stability. We further introduce techniques for speeding-up inverse optimization of boundary conditions for PDEs via backpropagation through time in latent space, and an annealing technique to address the non-differentiability and sparse interaction of boundary conditions. We test our method in a 1D benchmark of nonlinear PDEs, 2D Navier-Stokes flows into turbulent phase and an inverse optimization of boundary conditions in 2D Navier-Stokes flow. Compared to state-of-the-art deep learning-based surrogate models and other strong baselines, we demonstrate up to 128x reduction in the dimensions to update, and up to 15x improvement in speed, while achieving competitive accuracy. 
    more » « less
  2. In this paper, we consider a coupled chemotaxis-fluid system that models self-organized collective behavior of oxytactic bacteria in a sessile drop. This model describes the biological chemotaxis phenomenon in the fluid environment and couples a convective chemotaxis system for the oxygen-consuming and oxytactic bacteria with the incompressible Navier–Stokes equations subject to a gravitational force, which is proportional to the relative surplus of the cell density compared to the water density. We develop a new positivity preserving and high-resolution method for the studied chemotaxis-fluid system. Our method is based on the diffuse-domain approach, which we use to derive a new chemotaxis-fluid diffuse-domain (cf-DD) model for simulating bioconvection in complex geometries. The drop domain is imbedded into a larger rectangular domain, and the original boundary is replaced by a diffuse interface with finite thickness. The original chemotaxis-fluid system is reformulated on the larger domain with additional source terms that approximate the boundary conditions on the physical interface. We show that the cf-DD model converges to the chemotaxis-fluid model asymptotically as the width of the diffuse interface shrinks to zero. We numerically solve the resulting cf-DD system by a second-order hybrid finite-volume finite-difference method and demonstrate the performance of the proposed approach on a number of numerical experiments that showcase several interesting chemotactic phenomena in sessile drops of different shapes, where the bacterial patterns depend on the droplet geometries. 
    more » « less
  3. null (Ed.)
    Abstract Aerosol jet printing (AJP) is a direct-write additive manufacturing technique, which has emerged as a high-resolution method for the fabrication of a broad spectrum of electronic devices. Despite the advantages and critical applications of AJP in the printed-electronics industry, AJP process is intrinsically unstable, complex, and prone to unexpected gradual drifts, which adversely affect the morphology and consequently the functional performance of a printed electronic device. Therefore, in situ process monitoring and control in AJP is an inevitable need. In this respect, in addition to experimental characterization of the AJP process, physical models would be required to explain the underlying aerodynamic phenomena in AJP. The goal of this research work is to establish a physics-based computational platform for prediction of aerosol flow regimes and ultimately, physics-driven control of the AJP process. In pursuit of this goal, the objective is to forward a three-dimensional (3D) compressible, turbulent, multiphase computational fluid dynamics (CFD) model to investigate the aerodynamics behind: (i) aerosol generation, (ii) aerosol transport, and (iii) aerosol deposition on a moving free surface in the AJP process. The complex geometries of the deposition head as well as the pneumatic atomizer were modeled in the ansys-fluent environment, based on patented designs in addition to accurate measurements, obtained from 3D X-ray micro-computed tomography (μ-CT) imaging. The entire volume of the constructed geometries was subsequently meshed using a mixture of smooth and soft quadrilateral elements, with consideration of layers of inflation to obtain an accurate solution near the walls. A combined approach, based on the density-based and pressure-based Navier–Stokes formation, was adopted to obtain steady-state solutions and to bring the conservation imbalances below a specified linearization tolerance (i.e., 10−6). Turbulence was modeled using the realizable k-ε viscous model with scalable wall functions. A coupled two-phase flow model was, in addition, set up to track a large number of injected particles. The boundary conditions of the CFD model were defined based on experimental sensor data, recorded from the AJP control system. The accuracy of the model was validated using a factorial experiment, composed of AJ-deposition of a silver nanoparticle ink on a polyimide substrate. The outcomes of this study pave the way for the implementation of physics-driven in situ monitoring and control of AJP. 
    more » « less
  4. Lakhtakia, Akhlesh ; Martín-Palma, Raúl J. ; Knez, Mato (Ed.)
    In this study, the implementation and performance of bipennate topology fluidic artificial muscle (FAM) bundles operating under varying boundary conditions is investigated and quantified experimentally. Soft actuators are of great interest to design engineers due to their inherent flexibility and potential to improve safety in human robot interactions. McKibben fluidic artificial muscles are soft actuators which exhibit high force to weight ratios and dynamically replicate natural muscle movement. These features, in addition to their low fabrication cost, set McKibben FAMs apart as attractive components for an actuation system. Previous studies have shown that there are significant advantages in force and contraction outputs when using bipennate topology FAM bundles as compared to the conventional parallel topology1 . In this study, we will experimentally explore the effects of two possible boundary conditions imposed on FAMs within a bipennate topology. One boundary condition is to pin the muscle fiber ends with fixed pin spacings while the other is biologically inspired and constrains the muscle fibers to remain in contact. This paper will outline design considerations for building a test platform for bipennate fluidic artificial muscle bundles with varying boundary conditions and present experimental results quantifying muscle displacement and force output. These metrics are used to analyze the tradespace between the two boundary conditions and the effect of varying pennation angles. 
    more » « less
  5. Abstract The problem of parameter identification appears in many physical applications. A parameter of particular interest in cancer treatment is permeability, which modulates the fluidic streamlines in the tumor microenvironment. Most of the existing permeability identification techniques are invasive and not feasible to identify the permeability with minimal interference with the porous structure in their working conditions. In this paper, a theoretical framework utilizing partial differential equation (PDE)-constrained optimization strategies is established to identify a spatially distributed permeability of a porous structure from its modulated external velocity field measured around the structure. In particular, the flow around and through the porous media are governed by the steady-state Navier–Stokes–Darcy model. The performance of our approach is validated via numerical and experimental tests for the permeability of a 3D printed porous surrogate in a micro-fluidic chip based on the sampled optical velocity measurement. Both numerical and experimental results show a high precision of the permeability estimation. 
    more » « less