The numerical method for the optimal supporting position and related optimal control for the catalytic reaction system

This paper considers the numerical approximation for the optimal supporting position and related optimal control of a catalytic reaction system with some control and state constraints, which is governed by a nonlinear partial differential equations with given initial and boundary conditions. By the Galerkin finite element method, the original problem is projected into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then the control parameterization method is applied to approximate the control and reduce the original system to an optimal parameter selection problem, in which both the position and related control are taken as decision variables to be optimized. This problem can be solved as a nonlinear optimization problem by a particle swarm optimization algorithm. The numerical simulations are given to illustrate the effectiveness of the proposed numerical approximation method.     

Optimal production control of a service-oriented manufacturing system with customer balking behavior

The great challenge for service-oriented manufacturing (SOM) is how to cope with customer behavior while making decision on production planning and scheduling. In this paper, we consider a single-stage manufacturing system for SOM with impatient customers. In order to represent customer balking behavior caused by backlog, we employ a balking function, which is an arbitrary non-decreasing function of the backlog for characterizing the customer’s response to the backlog. The objective is to find the optimal production policy that minimizes the system cost. The problem is formulated as a Markov decision process. The optimal production policy is proved to be a base-stock policy. The effects of system parameters on the optimal base-stock level are analytically investigated, and the impact of customer balking behavior on the system is illustrated by numerical example in which linear balking function is employed. Numerical example shows that customer balking has a significant impact on the optimal control and the performance measures of the system under the optimal policy.     

An optimal control model of diabetes mellitus

Engineering optimal control theory is applied to equations describing insulin and glucose interactions. The nature of the optimal controller is established. It is shown how the results can be utilized in a closed loop feedback control system.     

A finite volume method solution for the bidomain equations and their application to modelling cardiac ischaemia

This paper presents an implementation of the finite volume method with the aim of studying subendocardial ischaemia during the ST segment. In this implementation, based on hexahedral finite volumes, each quadrilateral sub-face is split into two triangles to improve the accuracy of the numerical integration in complex geometries and when fibre rotation is included. The numerical method is validated against previously published solutions obtained from slab and cylindrical models of the left ventricle with subendocardial ischaemia and no fibre rotation. Epicardial potential distributions are then obtained for a half-ellipsoid model of the left ventricle. In this case it is shown that for isotropic cardiac tissue the degree of subendocardial ischaemia does not affect the epicardial potential distribution, which is consistent with previous findings from analytical studies in simpler geometries. The paper also considers the behaviour of various preconditioners for solving numerically the resulting system of algebraic equations resulting from the implementation of the finite volume method. It is observed that each geometry considered has its own optimal preconditioner.     

A Theory of Fluid Flow in Compliant Tubes

Starting with the Navier-Stokes equations, a system of equations is obtained to describe quasi-one-dimensional behavior of fluid in a compliant tube. The nonlinear terms which cannot be shown to be small in the original equations are retained, and the resulting equations are nonlinear. A functional pressure-area relationship is postulated and the final set of equations are quasi-linear and hyperbolic, with two independent and two dependent variables. A method of numerical solution of the set of equations is indicated, and the application to cases of interest is discussed.     

Optimal HIV treatment by maximising immune response

We present an optimal control model of drug treatment of the human immunodeficiency virus (HIV). Our model is based upon ordinary differential equations that describe the interaction between HIV and the specific immune response as measured by levels of natural killer cells. We establish stability results for the model. We approach the treatment problem by posing it as an optimal control problem in which we maximise the benefit based on levels of healthy CD4+ T cells and immune response cells, less the systemic cost of chemotherapy. We completely characterise the optimal control and compute a numerical solution of the optimality system via analytic continuation.Research supported by the Natural Science and Engineering Research Council (NSERC) and the Mathematics of Information Technology and Complex Systems (MITACS) of Canada     

Wavefront propagation in an activation model of the anisotropic cardiac tissue: asymptotic analysis and numerical simulations

In this paper we present a macroscopic model of the excitation process in the myocardium. The composite and anisotropic structure of the cardiac tissue is represented by a bidomain, i.e. a set of two coupled anisotropic media. The model is characterized by a non linear system of two partial differential equations of parabolic and elliptic type. A singular perturbation analysis is carried out to investigate the cardiac potential field and the structure of the moving excitation wavefront. As a consequence the cardiac current sources are approximated by an oblique dipole layer structure and the motion of the wavefront is described by eikonal equations. Finally numerical simulations are carried out in order to analyze some complex phenomena related to the spreading of the wavefront, like the front-front or front-wall collision. The results yielded by the excitation model and the eikonal equations are compared.     

Theoretical approach to blood ejection from the human left ventricle

An ejection dynamics mathematical model of human left ventricle (LV) based on physiological data of human heart is proposed in this study. The mathematical equations were expressed in terms of vorticity-stream function equations in a prolate spheroidal coordinate system. These equations combined with specified boundary conditions were numerically solved by using an alternating-direction-implicit (ADI) algorithm with second order accuracy. The unsteady aspects of the ejection process were subsequently introduced into the numerical simulation. The numerical results have shown that the present ellipsoidal model could be available to simulate the ejection process of the human LV. Such a model combined with cardiac muscle mechanics could be studied further to determine altered left ventricular function in cardiac diseases.     

Optimal control of a rabies epidemic model with a birth pulse

A system of ordinary differential equations describes the population dynamics of a rabies epidemic in raccoons. The model accounts for the dynamics of a vaccine, including loss of vaccine due to animal consumption and loss from factors other than racoon uptake. A control method to reduce the spread of disease is introduced through temporal distribution of vaccine packets. This work incorporates the effect of the seasonal birth pulse in the racoon population and the attendant increase in new-borns which are susceptible to the diseases, analysing the impact of the timing and length of this pulse on the optimal distribution of vaccine packets. The optimization criterion is to minimize the number of infected raccoons while minimizing the cost of distributing the vaccine. Using an optimal control setting, numerical results illustrate strategies for distributing the vaccine depending on the timing of the infection outbreak with respect to the birth pulse.     

Computational optimal control for the time fractional convection-diffusion-reaction system

This paper proposes a numerical approximation method for computational optimal control of a time fractional convection-diffusion-reaction system. The proposed method involves discretizing the spatial domain by finite element method, approximating the admissible controls by control parameterization, and then obtaining an optimal parameter selection problem which can be solved by numerical optimization algorithms such as sequential quadratic programming. Specifically, an implicit finite difference method is employed to solve the time fractional system, and the sensitivity method for gradient computation in integer order optimal control problems is adjusted to the fractional order case. Simulation results demonstrate the validity and accuracy of the proposed numerical approximation method.     

Optimal control for micro‐algae on a raceway model

We apply numerical optimal control methods to an existing algae growth model with the aim to determine the best performance of the model under known conditions using a variety of decision variables. We transform the system of differential algebraic equations in the existing model to a system of ordinary differential equations which introduces dynamics for average light intensity and chlorophyll. In addition, we allow for variable nitrogen concentration of the inflow as well as variable initial nitrogen concentration of the raceway. Our main focus is on optimizing of the production of lipids. We calculate both open and closed loop optimal controllers and test their robustness. Finally, we also consider raceway depth as a decision variable. © 2017 American Institute of Chemical Engineers Biotechnol. Prog., 34:107–119, 2018     

Predictive algorithms for neuromuscular control of human locomotion.

The problem of quantifying muscular activity of the human body can be formulated as an optimal control problem. The current methods used with large-scale biomechanical systems are non-derivative techniques. These methods are costly, as they require numerous integrations of the equations of motion. Additionally, the convergence is slow, making them impractical for use with large systems. We apply an efficient numerical algorithm to the biomechanical optimal control problem. Using direct collocation with a trapezoidal discretization, the equations of motion are converted into a set of algebraic constraint equations. An augmented Lagrangian formulation is used for the optimization problem to handle both equality and inequality constraints. The resulting min-max problem is solved with a generalized Newton method. In contrast to the prevalent optimal control implementations, we calculate analytical first- and second-derivative information and obtain local quadratic convergence. To demonstrate the efficacy of the method, we solve a steady-state pedaling problem with 7 segments and 18 independent muscle groups. The computed muscle activations compare well with experimental EMG data. The computational effort is significantly reduced and solution times are a fraction of those of the non-derivative techniques.     

Optimal Foraging and Predator–Prey Dynamics

A system consisting of a population of predators and two types of prey is considered. The dynamics of the system is described by differential equations with controls. The controls model how predators forage on each of the two types of prey. The choice of these controls is based on the standard assumption in the theory of optimal foraging which requires that each predator maximizes the net rate of energy intake during foraging. Since this choice depends on the densities of populations involved, this allows us to link the optimal behavior of an individual with the dynamics of the whole system. Simple qualitative analysis and some simulations show the qualitative behavior of such a system. The effect of the optimal diet choice on the stability of the system is discussed.     

The central core model of the renal medulla: A particular solution

A simple particular solution of a set of differential equations that model the behavior of the human kidney is obtained. This solution is extended by means of an expansion procedure, and numerical results are presented in graphical form. These results may have application to kidney malfunction and also provide values with which to compare numerical results of the basic equations.     

Shape optimization in unsteady blood flow: a numerical study of non-Newtonian effects

This paper presents a numerical study of non-Newtonian effects on the solution of shape optimization problems involving unsteady pulsatile blood flow. We consider an idealized two dimensional arterial graft geometry. Our computations are based on the Navier-Stokes equations generalized to non-Newtonian fluid, with the modified Cross model employed to account for the shear-thinning behavior of blood. Using a gradient-based optimization algorithm, we compare the optimal shapes obtained using both the Newtonian and generalized Newtonian constitutive equations. Depending on the shear rate prevalent in the domain, substantial differences in the flow as well as in the computed optimal shape are observed when the Newtonian constitutive equation is replaced by the modified Cross model. By varying a geometric parameter in our test case, we investigate the influence of the shear rate on the solution.     

Numerical simulation for the propagation of nonlinear pulsatile waves in arteries.

The behavior of nonlinear pulsatile flow of incompressible blood contained in an elastic tube is examined. The theory takes into account the nonlinear convective terms of the Navier-Stokes equations. The motion of the arterial wall is characterized by a set of linearized differential equations. The region bounded by the flexible arterial wall is mapped into a fixed area in which numerical discretization takes place. The finite element method (Galerkin weighted residual approach) is used for the solution of this nonlinear system. The results obtained are pressure distribution, velocity profile, flow rate and wall displacements along the elastic tube (20 cm long).     

Coupled Attitude-Orbit Dynamics and Control for an Electric Sail in a Heliocentric Transfer Mission

The paper discusses the coupled attitude-orbit dynamics and control of an electric-sail-based spacecraft in a heliocentric transfer mission. The mathematical model characterizing the propulsive thrust is first described as a function of the orbital radius and the sail angle. Since the solar wind dynamic pressure acceleration is induced by the sail attitude, the orbital and attitude dynamics of electric sails are coupled, and are discussed together. Based on the coupled equations, the flight control is investigated, wherein the orbital control is studied in an optimal framework via a hybrid optimization method and the attitude controller is designed based on feedback linearization control. To verify the effectiveness of the proposed control strategy, a transfer problem from Earth to Mars is considered. The numerical results show that the proposed strategy can control the coupled system very well, and a small control torque can control both the attitude and orbit. The study in this paper will contribute to the theory study and application of electric sail.     

An optimal shape problem related to the realistic design of river fishways

A river fishway is a hydraulic structure that facilitates fish in overcoming obstacles (dams, waterfalls, etc.) to their spawning and other migrations in rivers. In this work we present a mathematical formulation of an optimal design problem for a vertical slot fishway, where the state system is given by the 2D shallow water equations fixing the height and velocity of water, the design variables are the geometry of the slots, and the objective function is determined by the existence of rest areas for fish and of a water velocity suitable for fish swimming capability. We also derive an expression for the gradient of the objective function via the adjoint system. From the numerical point of view, we present a characteristic-Galerkin method for solving the shallow water equations, and a direct search algorithm for the computation of the optimal design variables. Finally, we give numerical results obtained for a standard ten pools channel.     

Shape optimization in unsteady blood flow: A numerical study of non-Newtonian effects

This paper presents a numerical study of non-Newtonian effects on the solution of shape optimization problems involving unsteady pulsatile blood flow. We consider an idealized two dimensional arterial graft geometry. Our computations are based on the Navier–Stokes equations generalized to non-Newtonian fluid, with the modified Cross model employed to account for the shear-thinning behavior of blood. Using a gradient-based optimization algorithm, we compare the optimal shapes obtained using both the Newtonian and generalized Newtonian constitutive equations. Depending on the shear rate prevalent in the domain, substantial differences in the flow as well as in the computed optimal shape are observed when the Newtonian constitutive equation is replaced by the modified Cross model. By varying a geometric parameter in our test case, we investigate the influence of the shear rate on the solution.