This paper presents a numerical investigation of plaque growth in a diseased artery using the two-way fluid–structural interaction (FSI) technique. An axis-asymmetric 45% stenosis model is used as the base model to start the plaque growth approximation. The blood is modelled as a non-Newtonian fluid described by the Casson model. The artery tissue is assumed to be a nonlinear material. The two-way FSI simulation is carried out in a way that mimics the unsteady blood flow through a diseased artery by using a pulsatile flow condition. After each flow velocity cycle, the numerical results are extracted and used to modify the stenosis geometry based upon critical wall shear stress (WSS) values and an accepted relationship between the concentration of low density lipoprotein and WSS. The simulation procedure is repeated until the growth-updated stenosis morphology reaches 79% severity. The behaviour of the flow velocity is analysed at each growth stage, together with the WSS, to determine the change of plaque morphology due to growth. The effects of WSS and pressure on the artery wall at the final stage (79% severity) of the plaque growth model are also compared with results from the authors’ previous work, to demonstrate the importance of the morphology change in plaque growth modelling.

Selective withdrawal of a two-layer fluid is considered. The fluid layers are weakly compressible, miscible and viscous and therefore flow rotationally. The lower, denser fluid flows with constant velocity out through one or more drain holes in the bottom of a rectangular tank. The drain is opened impulsively and the subsequent draw-down produces waves in the interface which travel outward to the edges of the tank and are reflected back with a $18{0}^{\circ } $ change of phase. The points on the interface that have the highest absolute gradient form regions of high vorticity in the tank, enabling mixing of the fluids. An inviscid linearized interface is computed and compared to contour plots of density for the viscous solution. The two agree closely at early times in the withdrawal process, but as time increases, nonlinear and viscous effects take over. The time at which the lighter fluid starts to flow out of the tank is dependent on the number of drains, their width, and the fluid flow rate and density, and is investigated here.

We develop a computational method for solving an optimal control problem governed by a switched impulsive dynamical system with time delay. At each time instant, only one subsystem is active. We propose a computational method for solving this optimal control problem where the time spent by the state in each subsystem is treated as a new parameter. These parameters and the jump strengths of the impulses are decision parameters to be optimized. The gradient formula of the cost function is derived in terms of solving a number of delay differential equations forward in time. Based on this, the optimal control problem can be solved as an optimization problem.

Temperature distributions recorded by thermocouples in a solid body (slab) subject to surface heating are used in a mathematical model of two-dimensional heat conduction. The corresponding Dirichlet problem for a holomorphic function (complex potential), involving temperature and a heat stream function, is solved in a strip. The Zhukovskii function is reconstructed through singular integrals, involving an auxiliary complex variable. The complex potential is mapped onto an auxiliary half-plane. The flow net (orthogonal isotherms and heat lines) of heat conduction is compared with the known Carslaw–Jaeger solution and shows a puzzling topology of three regimes of energy fluxes for temperature boundary conditions common in passive thermal insulation. The simplest regime is realized if cooling of a shaded zone is mild and heat flows in a slightly distorted “resistor model” flow tube. The second regime emerges when cooling is stronger and two disconnected separatrices demarcate the back-flow of heat from a relatively hot segment of the slab surface to the atmosphere through relatively cold parts of this surface. The third topological regime is characterized by a single separatrix with a critical point inside the slab, where the thermal gradient is nil. In this regime the back-suction of heat into the atmosphere is most intensive. The closed-form solutions obtained can be used in assessment of efficiency of thermal protection of buildings.

A theoretical investigation of the unsteady flow of a Newtonian fluid through a channel is presented using an alternative boundary condition to the standard no-slip condition, namely the Navier boundary condition, independently proposed over a hundred years ago by both Navier and Maxwell. This boundary condition contains an extra parameter called the slip length, and the most general case of a constant but different slip length on each channel wall is studied. An analytical solution for the velocity distribution through the channel is obtained via a Fourier series, and is used as a benchmark for numerical simulations performed utilizing a finite element analysis modified with a penalty method to implement the slip boundary condition. Comparison between the analytical and numerical solution shows excellent agreement for all combinations of slip lengths considered.