Hostname: page-component-586b7cd67f-t8hqh Total loading time: 0 Render date: 2024-11-24T14:39:27.008Z Has data issue: false hasContentIssue false

Adjoint-based shape optimization of the microchannels in an inkjet printhead

Published online by Cambridge University Press:  17 May 2019

Petr V. Kungurtsev
Affiliation:
Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, UK
Matthew P. Juniper*
Affiliation:
Department of Engineering, University of Cambridge, Cambridge CB2 1PZ, UK
*
Email address for correspondence: [email protected]

Abstract

In drop-on-demand inkjet printheads, ink is pumped steadily through small channels, each of which contains an actuator and a nozzle. When an actuator pulses, a droplet is forced through the nozzle, after which acoustic oscillations reverberate within the channel. Manufacturers would like to damp the residual reverberations, without increasing the pressure drop required to drive the steady flow. In this paper we use gradient-based optimization to show that this can be achieved by constricting the channel where the acoustic velocity is largest and enlarging the channel where the acoustic velocity is smallest. This increases the viscothermal dissipation of the acoustics without changing the viscous dissipation of the steady flow. We separate the compressible Navier–Stokes equations into equations for a steady flow with no oscillations and equations for oscillations with no steady flow. We define two objective functions: the viscous dissipation of the steady flow and the dissipation of the oscillations. We then derive the adjoints for both sets of equations, and obtain expressions for the gradient of each objective function with respect to boundary deformations in Hadamard form. We combine these with a gradient-based optimization algorithm, incorporating constraints such as the shapes of the actuator and nozzle. This algorithm quickly converges to a design that has the same viscous dissipation for the steady flow but a 50 % larger decay rate for the oscillating flow. We show that this design is nearly optimal. It is a shape that inkjet manufacturers, using physical insight and trial and error, have probably not yet considered. We also show how the adjoint fields provide physical insight into the mechanisms affecting each objective function. The main requirements of this method are that the steady flow Mach number and oscillating flow Mach number are small, and that dissipation is dominated by thermoviscous mechanisms. These requirements are often satisfied in microfluidics, so the method in this paper could be applied to many other applications.

Type
JFM Papers
Copyright
© 2019 Cambridge University Press 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

Beltman, W. M.1998 Viscothermal wave propagation including acousto–elastic interaction. PhD thesis, University of Twente, The Netherlands.Google Scholar
Beltman, W. M. 1999 Viscothermal wave propagation including acousto–elastic interaction. Part I. Theory. J. Sound Vib. 227 (3), 555586.Google Scholar
Berggren, M., Bernland, A. & Noreland, D. 2018 Acoustic boundary layers as boundary conditions. J. Comput. Phys. 371, 633650.Google Scholar
Bogy, D. B. & Talke, F. E. 1984 Experimental and theoretical study of wave propagation phenomena in drop-on-demand ink jet devices. IBM J. Res. Dev. 28 (3), 314321.Google Scholar
Carslaw, H. S. & Jaeger, J. C. 1986 Conduction of Heat in Solids. Clarendon Press.Google Scholar
Christensen, R. 2017 Topology optimization of thermoviscous acoustics in tubes and slits with hearing aid applications. In COMSOL Conference 2017 Rotterdam, at Rotterdam. COMSOL.Google Scholar
Chu, B. T. 1965 On the energy transfer to small disturbances in fluid flow (part I). Acta Mechanica 1 (3), 215234.Google Scholar
Culick, F., Heitor, M. V. & Whitelaw, J. H. 2012 Equations for unsteady motions in combustion chambers. In Unsteady Combustion, vol. 306. chap. 3, Springer Science & Business Media, 3‐1–3‐14.Google Scholar
Delfour, M. & Zolésio, J. 2011 Shape and tangential differential calculus. In Shapes and Geometries, 2nd edn. chap. 9, pp. 457518. Society for Industrial and Applied Mathematics.Google Scholar
Geuzaine, C. & Remacle, J.-F. 2009 Gmsh: a 3D finite element mesh generator with built-in pre- and post-processing facilities. Intl J. Numer. Meth. Engng 79 (11), 13091331.Google Scholar
Gunzburger, M. 2002 Perspectives in Flow Control and Optimization. Society for Industrial and Applied Mathematics.Google Scholar
Hoath, S. D.(Ed.) 2016 Fundamentals of Inkjet Printing, Wiley-VCH Verlag GmbH & Co. KGaA.Google Scholar
Homentcovschi, D., Miles, R. N., Loeppert, P. V. & Zuckerwar, A. J. 2014 A microacoustic analysis including viscosity and thermal conductivity to model the effect of the protective cap on the acoustic response of a MEMS microphone. Microsystem Technol. 20 (2), 265272.Google Scholar
Homentcovschi, D., Murray, B. T. & Miles, R. N. 2010 An analytical formula and FEM simulations for the viscous damping of a periodic perforated MEMS microstructure outside the lubrication approximation. Microfluid. Nanofluid. 9 (4–5), 865879.Google Scholar
Kampinga, W. R., Wijnant, Y. H. & de Boer, A. 2010 Performance of several viscothermal acoustic finite elements. Acta Acust. united Acust. 96, 115124.Google Scholar
Kampinga, W. R., Wijnant, Y. H. & de Boer, A. 2011 An efficient finite element model for viscothermal acoustics. Acta Acust. united Acust. 97 (4), 618631.Google Scholar
Logg, A., Mardal, K.-A. & Wells, G. N. 2012 Automated Solution of Differential Equations by the Finite Element Method. Springer.Google Scholar
Luchini, P. & Bottaro, A. 2014 Adjoint equations in stability analysis. Annu. Rev. Fluid Mech. 46, 493517.Google Scholar
Moser, M. 1980 Damping of structure-borne sound by the viscosity of a layer between 2 plates. Acustica 46 (2), 210217.Google Scholar
Müller, B. 1998 Low-Mach-number asymptotics of the Navier–Stokes equations. J. Engng Maths 34 (1), 97109.Google Scholar
Myers, M. K. 1980 On the acoustic boundary condition in the presence of flow. J. Sound Vib. 71 (3), 429434.Google Scholar
Rienstra, S. W. & Hirschberg, A. 2013 An introduction to acoustics. Eindhoven University of Technology.Google Scholar
Rognes, M. E. & Logg, A. 2013 Automated goal-oriented error control I: stationary variational problems. SIAM J. Sci. Comput. 35 (3), C173C193.Google Scholar
Samareh, J. A. 2001 Survey of shape parameterization techniques for high-fidelity multidisciplinary shape optimization. AIAA J. 39 (5), 877884.Google Scholar
Schmidt, S. & Schulz, V. 2010 Shape derivatives for general objective functions and the incompressible Navier–Stokes equations. Control Cybern. 39 (3), 677713.Google Scholar
Seccombe, D. et al. 1997 Ink-cooled thermal inkjet printhead. US Patent 5657061.Google Scholar
Sokolowski, J. & Zolesio, J.-P. 1992 Introduction to Shape Optimization. Springer.Google Scholar
Sonntag, M., Schmidt, S. & Gauger, N. R. 2016 Shape derivatives for the compressible Navier–Stokes equations in variational form. J. Comput. Appl. Maths 296, 334351.10.1016/j.cam.2015.09.010Google Scholar
Svanberg, K. 1987 The method of moving asymptotes – a new method for structural optimization. Intl J. Numer. Meth. Engng 24 (2), 359373.Google Scholar
Tijdeman, H. 1975 On the propagation of sound waves in cylindrical tubes. J. Sound Vib. 39 (1), 133.10.1016/S0022-460X(75)80206-9Google Scholar