Derivative-free optimization and filter methods to solve nonlinear constrained problems

In real optimization problems, usually the analytical expression of the objective function is not known, nor its derivatives, or they are complex. In these cases it becomes essential to use optimization methods where the calculation of the derivatives, or the verification of their existence, is not necessary: the Direct Search Methods or Derivative-free Methods are one solution. When the problem has constraints, penalty functions are often used. Unfortunately the choice of the penalty parameters is, frequently, very difficult, because most strategies for choosing it are heuristics strategies. As an alternative to penalty function appeared the filter methods. A filter algorithm introduces a function that aggregates the constrained violations and constructs a biobjective problem. In this problem the step is accepted if it either reduces the objective function or the constrained violation. This implies that the filter methods are less parameter dependent than a penalty function. In this work, we present a new direct search method, based on simplex methods, for general constrained optimization that combines the features of the simplex method and filter methods. This method does not compute or approximate any derivatives, penalty constants or Lagrange multipliers. The basic idea of simplex filter algorithm is to construct an initial simplex and use the simplex to drive the search. We illustrate the behavior of our algorithm through some examples. The proposed methods were implemented in Java.

Analysis of diffusion process in fractured reservoirs using fractional derivative approach

The fractal geometry is used to model of a naturally fractured reservoir and the concept of fractional derivative is applied to the diffusion equation to incorporate the history of fluid flow in naturally fractured reservoirs. The resulting fractally fractional diffusion (FFD) equation is solved analytically in the Laplace space for three outer boundary conditions. The analytical solutions are used to analyze the response of a naturally fractured reservoir considering the anomalous behavior of oil production. Several synthetic examples are provided to illustrate the methodology proposed in this work and to explain the diffusion process in fractally fractured systems.

Descending pain modulation in a Kaolin-induced hydrocephalus rat model

Pain transmission at the spinal cord is modulated by descending actions that arise from supraspinal areas which collectively form the endogenous pain control system. Two key areas involved of the endogenous pain control system have a circunventricular location, namely the periaqueductal grey (PAG) and the locus coeruleus (LC). The PAG plays a crucial role in descending pain modulation as it conveys the input from higher brain centers to the spinal cord. As to the LC, it is involved in descending pain inhibition by direct noradrenergic projections to the spinal cord. In the context of neurological defects, several diseases may affect the structure and function of the brain. Hydrocephalus is a congenital or acquired disease characterized by an enlargement of the ventricles which leads to a distortion of the adjacent tissues, including the PAG and LC. Usually, patients suffering from hydrocephalus present dysfunctions in learning and memory and also motor deficits. It remains to be evaluated if lesions of the periventricular brain areas involved in pain control during hydrocephalus may affect descending pain control and, herein, affect pain responses. The studies included in the present thesis used an experimental model of hydrocephalus (the rat injected in the cisterna magna with kaolin) to study descending modulation of pain, focusing on the two circumventricular regions referred above (the PAG and the LC). In order to evaluate the effects of kaolin injection into the cisterna magna, we measured the degree of ventricular dilatation in sections encompassing the PAG by standard cytoarquitectonic stanings (thionin staining). For the LC, immunodetection of the noradrenaline-synthetizing enzyme tyrosine hydroxylase (TH) was performed, due to the noradrenergic nature of the LC neurons. In general, rats with kaolin-induced hydrocephalus presented a higher dilatation of the 4th ventricle, along with a tendency to a higher area of the PAG. Due to the validated role of detection the c-fos protooncogene as a marker of neuronal activation, we also studied neuronal activation in the several subnuclei which compose the PAG, namely the dorsomedial, dorsolateral, lateral and ventrolateral (VLPAG) parts. A decrease in the numbers of neurons immunoreactive for Fos protein (the product of activation of the c-fos protooncogene) was detected in rats injected with kaolin, whereas the remaining PAG subnuclei did not present changes in Fos-immunoreactive nuclei. Increases in the levels of TH in the LC, namely at the rostral parts of the nucleus, were detected in hydrocephalic animals. The following pain-related parameters were measured, namely 1) pain behavioural responses in a validated pain inflammatory test (the formalin test) and 2) the nociceptive activation of spinal cord neurons. A decrease in behavioral responses was detected in rats with kaolin-induced hydrocephalus was detected, namely in the second phase of the test (inflammatory phase). This is the phase of the formalin test in which the motor behaviour is less important, which is important since a semi-quantitative analysis of the motor performance of rats injected with kaolin indicates that these animals may present some motor impairments. Collectively, the results of the behavioral studies indicate that rats with kaolin-induced hydrocephalus exhibit hypoalgesia. A decrease in Fos expression was detected at the superficial dorsal layers of the spinal cord in rats with kaolin-induced hydrocephalus, further indicating that hydrocephalus decreases nociceptive responses. It remains to be ascertained if this is due to alterations in the PAG and LC in the rats with kaolin-induced hydrocephalus, which may affect descending pain modulation. It remains to be evaluated what are the mechanisms underlying the increased pain inhibition at the spinal dorsal horn in the hydrocephalus rats. Regarding the VLPAG, the decrease in neuronal activity may impair descending modulation. Since the LC has higher levels of TH in rats with kaolininduced hydrocephalus, which also appears to increase the noradrenergic innervation in the spinal dorsal horn, it is possible that an increase in the release of noradrenaline at the spinal cord accounts for pain inhibition. Our studies also determine the need to study in detail patients with hydrocephalus namely in what concerns their thresholds to pain and to perform imaging studies focused on the structure and function of pain control areas in the brain.

Strange Dynamics in a Fractional Derivative of Complex-Order Network of Chaotic Oscillators

We study the peculiar dynamical features of a fractional derivative of complex-order network. The network is composed of two unidirectional rings of cells, coupled through a "buffer" cell. The network has a Z3 × Z5 cyclic symmetry group. The complex derivative Dα±jβ, with α, β ∈ R+ is a generalization of the concept of integer order derivative, where α = 1, β = 0. Each cell is modeled by the Chen oscillator. Numerical simulations of the coupled cell system associated with the network expose patterns such as equilibria, periodic orbits, relaxation oscillations, quasiperiodic motion, and chaos, in one or in two rings of cells. In addition, fixing β = 0.8, we perceive differences in the qualitative behavior of the system, as the parameter c ∈ [13, 24] of the Chen oscillator and/or the real part of the fractional derivative, α ∈ {0.5, 0.6, 0.7, 0.8, 0.9, 1.0}, are varied. Some patterns produced by the coupled system are constrained by the network architecture, but other features are only understood in the light of the internal dynamics of each cell, in this case, the Chen oscillator. What is more important, architecture and/or internal dynamics?