RESUMEN
We consider the variational foam model, where the goal is to minimize the total surface area of a collection of bubbles subject to the constraint that the volume of each bubble is prescribed. We apply sharp interface methods to develop an efficient computational method for this problem. In addition to simulating time dynamics, we also report on stationary states of this flow for ≤ 21 bubbles in two dimensions and ≤ 17 bubbles in three dimensions. For small numbers of bubbles, we recover known analytical results, which we briefly discuss. In two dimensions, we also recover previous numerical results, computed using other methods. Particular attention is given to locally optimal foam configurations and heterogeneous foams, where the volumes of the bubbles are not equal. Configurational transitions are reported for the quasi-stationary flow where the volume of one of the bubbles is varied and, for each volume, the stationary state is computed. The results from these numerical experiments are described and accompanied by many figures and videos.
Asunto(s)
Modelos Teóricos , Algoritmos , Análisis Numérico Asistido por Computador , TermodinámicaRESUMEN
We consider Steklov eigenvalues of reflection-symmetric, nearly circular, planar domains. Treating such domains as perturbations of the disc, we obtain a second-order formal asymptotic estimate in the domain perturbation parameter. We conclude with a discussion of implications for isoperimetric inequalities. Namely, our results corroborate the results of Weinstock and Brock that state, respectively, that the disc is the maximizer for the area and perimeter constrained problems. They also support the result of Hersch, Payne and Schiffer that the product of the first two eigenvalues is maximal among all open planar sets of equal perimeter. In addition, our results imply that the disc is not the maximizer of the area constrained problems for higher even numbered Steklov eigenvalues, as suggested by previous numerical results.
RESUMEN
OBJECTIVE: Deep brain stimulation (DBS) is a growing treatment option for movement and psychiatric disorders. As DBS technology moves toward directional leads with increased numbers of smaller electrode contacts, trial-and-error methods of manual DBS programming are becoming too time-consuming for clinical feasibility. We propose an algorithm to automate DBS programming in near real-time for a wide range of DBS lead designs. APPROACH: Magnetic resonance imaging and diffusion tensor imaging are used to build finite element models that include anisotropic conductivity. The algorithm maximizes activation of target tissue and utilizes the Hessian matrix of the electric potential to approximate activation of neurons in all directions. We demonstrate our algorithm's ability in an example programming case that targets the subthalamic nucleus (STN) for the treatment of Parkinson's disease for three lead designs: the Medtronic 3389 (four cylindrical contacts), the direct STNAcute (two cylindrical contacts, six directional contacts), and the Medtronic-Sapiens lead (40 directional contacts). MAIN RESULTS: The optimization algorithm returns patient-specific contact configurations in near real-time-less than 10 s for even the most complex leads. When the lead was placed centrally in the target STN, the directional leads were able to activate over 50% of the region, whereas the Medtronic 3389 could activate only 40%. When the lead was placed 2 mm lateral to the target, the directional leads performed as well as they did in the central position, but the Medtronic 3389 activated only 2.9% of the STN. SIGNIFICANCE: This DBS programming algorithm can be applied to cylindrical electrodes as well as novel directional leads that are too complex with modern technology to be manually programmed. This algorithm may reduce clinical programming time and encourage the use of directional leads, since they activate a larger volume of the target area than cylindrical electrodes in central and off-target lead placements.