Skip to main content Accessibility help

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

  • Petr V. Kungurtsev (a1) and Matthew P. Juniper (a1)


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.


Corresponding author

Email address for correspondence:


Hide All
Beltman, W. M.1998 Viscothermal wave propagation including acousto–elastic interaction. PhD thesis, University of Twente, The Netherlands.
Beltman, W. M. 1999 Viscothermal wave propagation including acousto–elastic interaction. Part I. Theory. J. Sound Vib. 227 (3), 555586.
Berggren, M., Bernland, A. & Noreland, D. 2018 Acoustic boundary layers as boundary conditions. J. Comput. Phys. 371, 633650.
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.
Carslaw, H. S. & Jaeger, J. C. 1986 Conduction of Heat in Solids. Clarendon Press.
Christensen, R. 2017 Topology optimization of thermoviscous acoustics in tubes and slits with hearing aid applications. In COMSOL Conference 2017 Rotterdam, at Rotterdam. COMSOL.
Chu, B. T. 1965 On the energy transfer to small disturbances in fluid flow (part I). Acta Mechanica 1 (3), 215234.
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.
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.
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.
Gunzburger, M. 2002 Perspectives in Flow Control and Optimization. Society for Industrial and Applied Mathematics.
Hoath, S. D.(Ed.) 2016 Fundamentals of Inkjet Printing, Wiley-VCH Verlag GmbH & Co. KGaA.
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.
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.
Kampinga, W. R., Wijnant, Y. H. & de Boer, A. 2010 Performance of several viscothermal acoustic finite elements. Acta Acust. united Acust. 96, 115124.
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.
Logg, A., Mardal, K.-A. & Wells, G. N. 2012 Automated Solution of Differential Equations by the Finite Element Method. Springer.
Luchini, P. & Bottaro, A. 2014 Adjoint equations in stability analysis. Annu. Rev. Fluid Mech. 46, 493517.
Moser, M. 1980 Damping of structure-borne sound by the viscosity of a layer between 2 plates. Acustica 46 (2), 210217.
Müller, B. 1998 Low-Mach-number asymptotics of the Navier–Stokes equations. J. Engng Maths 34 (1), 97109.
Myers, M. K. 1980 On the acoustic boundary condition in the presence of flow. J. Sound Vib. 71 (3), 429434.
Rienstra, S. W. & Hirschberg, A. 2013 An introduction to acoustics. Eindhoven University of Technology.
Rognes, M. E. & Logg, A. 2013 Automated goal-oriented error control I: stationary variational problems. SIAM J. Sci. Comput. 35 (3), C173C193.
Samareh, J. A. 2001 Survey of shape parameterization techniques for high-fidelity multidisciplinary shape optimization. AIAA J. 39 (5), 877884.
Schmidt, S. & Schulz, V. 2010 Shape derivatives for general objective functions and the incompressible Navier–Stokes equations. Control Cybern. 39 (3), 677713.
Seccombe, D. et al. 1997 Ink-cooled thermal inkjet printhead. US Patent 5657061.
Sokolowski, J. & Zolesio, J.-P. 1992 Introduction to Shape Optimization. Springer.
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/
Svanberg, K. 1987 The method of moving asymptotes – a new method for structural optimization. Intl J. Numer. Meth. Engng 24 (2), 359373.
Tijdeman, H. 1975 On the propagation of sound waves in cylindrical tubes. J. Sound Vib. 39 (1), 133.10.1016/S0022-460X(75)80206-9
MathJax is a JavaScript display engine for mathematics. For more information see

JFM classification


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed