Jato transversal de gotas: simulações por ALE/FEM e efeitos interfaciais.

Um código computacional para escoamentos bifásicos incorporando metodologia híbrida entre oMétodo dos Elementos Finitos e a descrição Lagrangeana-Euleriana Arbitrária do movimento é usado para simular a dinâmica de um jato transversal de gotas na zona primária de quebra. Os corpos dispersos são descritos por meio de um método do tipo front-tracking que produz interfaces de espessura zero através de malhas formadas pela união de elementos adjacentes em ambas as fases e de técnicas de refinamento adaptativo. Condições de contorno periódicas são implementadas de modo variacionalmente consistente para todos os campos envolvidos nas simulações apresentadas e uma versão modificada do campo de pressão é adicionada à formulação do tipo um-fluido usada na equação da quantidade de movimento linear. Simulações numéricas diretas em três dimensões são executadas para diferentes configurações de líquidos imiscí veis compatíveis com resultados experimentais encontrados na literatura. Análises da hidrodinâmica do jato transversal de gotas nessas configurações considerando trajetórias, variação de formato de gota, espectro de pequenas perturbações, além de aspectos complementares relativos à qualidade de malha são apresentados e discutidos.

Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Structural and thermodynamic properties of different phases of supercooled liquid water.

Computer simulation results are reported for a realistic polarizable potential model of water in the supercooled region. Three states, corresponding to the low density amorphous ice, high density amorphous ice, and very high density amorphous ice phases are chosen for the analyses. These states are located close to the liquid-liquid coexistence lines already shown to exist for the considered model. Thermodynamic and structural quantities are calculated, in order to characterize the properties of the three phases. The results point out the increasing relevance of the interstitial neighbors, which clearly appear in going from the low to the very high density amorphous phases. The interstitial neighbors are found to be, at the same time, also distant neighbors along the hydrogen bonded network of the molecules. The role of these interstitial neighbors has been discussed in connection with the interpretation of recent neutron scattering measurements. The structural properties of the systems are characterized by looking at the angular distribution of neighboring molecules, volume and face area distribution of the Voronoi polyhedra, and order parameters. The cumulative analysis of all the corresponding results confirms the assumption that a close similarity between the structural arrangement of molecules in the three explored amorphous phases and that of the ice polymorphs I(h), III, and VI exists.

发布/订阅系统可伸缩性支持技术研究

在发布/订阅系统中，订阅者向系统注册其对某一个或某一类事件的兴趣，发布者将生成的事件发布到系统中，发布/订阅系统负责将事件与订阅进行匹配，并向对事件感兴趣的订阅者发出事件通知。发布者和订阅者在时间、空间和控制流上完全解耦合，能很好地满足大型分布式应用对松耦合交互的需求。发布/订阅系统在金融服务、新闻服务、物流等领域有着很大的实际应用价值。随着传感器技术的发展和普适应用的增长， 发布/订阅系统应用于传感器网络和 RFID应用等普适环境下的研究方兴未艾。在普适环境下，分布式大规模多种类事件给发布/订阅系统带来了更多的系统效率和可伸缩性问题。 来自传感器或RFID设备的原子事件数据被采集获取后即被发送到功能更为强大的代理服务器上处理。当传感器节点很多、传感数据规模很大时，大量收集到的数据容易造成网络拥塞以及某些代理服务器负载过重，从而成为性能瓶颈。因此，需要设计有效的代理部署算法，选取代理服务器的最佳位置，构建数据收集网络的拓扑结构，以高效地利用代理服务器的处理能力和网络带宽，提高数据的采集效率，从而支持发布/订阅系统更高效地处理大量事件和订阅。 在基于内容的发布/订阅系统中，基于内容的订阅是高度多样化的，不同的事件可能满足大量的感兴趣的订阅者。 系统中， 大规模的订阅信息和事件的到达，可能造成某些代理服务器负载过重，网络消息流量也过大，使得事件不能及时送达订阅者。因此，需要设计一个有效的消息路由策略，将发布者发布的事件消息高效地传递给订阅者，并且将事件匹配的责任分派到网络中不同的代理服务器上，使得网络中各代理服务器的负载分布相对均衡。 本文围绕发布/订阅系统的可伸缩性，对原子事件数据收集中的代理服务器部署问题，以及原子事件/订阅的路由问题展开研究。本文针对一般部署问题，即存在禁止部署区域的情形，提出了基于栅格化和爬山法的代理部署算法；对于无禁止区域的部署问题， 提出了更优的基于 Voronoi 图和Steiner 树的代理部署算法。本文实现了上述两个算法，实验表明，算法能获得在部署代价、通信开销和负载均衡度方面均较优的部署方案。在原子事件/订阅的路由上，本文提出了基于事件空间划分的渠道路由算法。算法将事件空间划分为多个事件渠道，每个事件渠道分配给单个代理服务器负责。代理服务器对到达的事件或订阅，根据其所属的渠道，决定转发或在本地处理。算法还通过实时监控代理服务器的负载，在负载较重时，分割渠道，在负载较轻时，合并渠道，从而达到动态负载均衡的效果。本文在OncePubSub系统上实现了上述算法，实验表明，算法能够有效地降低事件和订阅的传递开销，减小事件处理延迟，并能均衡代理服务器的负载。

煤层气数值试井非结构网格生成方法研究

分析了数值试井中的非结构网格的生成方法,利用Delaunay三角剖分的优良特性,设计了"翼边"数据结构,实现了满足数值计算前处理要求的Delaunay网格生成算法,并对网格质量进行了控制。实例表明,生成的网格具有边界约束一致性和较好的质量,可以满足煤层气数值试井计算的要求。

乡镇级人口统计数据空间化研究

Population data which collected and saved according to administrative region is a kind of statistical data. As a traditional method of spatial data expression, average distribution in every administrative region brings population data on a low spatial and temporal precision. Now, an accurate population data with high spatial resolution is becoming more and more important in regional planning, environment protection, policy making and rural-urban development. Spatial distribution of population data is becoming more important in GIS study area. In this article, the author reviewed the progress of research on spatial distribution of population. Under the support of GIS, correlative geographical theories and Grid data model, Remote Sensing data, terrain data, traffic data, river data, resident data, and social economic statistic were applied to calculate the spatial distribution of population in Fujian province, which includes following parts: (1) Simulating of boundary at township level. Based on access cost index, land use data, traffic data, river data, DEM, and correlative social economic statistic data, the access cost surface in study area was constructed. Supported by the lowest cost path query and weighted Voronoi diagram, DVT model (Demarcation of Villages and Towns) was established to simulate the boundary at township level in Fujian province. (2) Modeling of population spatial distribution. Based on the knowledge in geography, seven impact factors, such as land use, altitude, slope, residential area, railway, road, and river were chosen as the parameters in this study. Under the support of GIS, the relations of population distribution to these impact factors were analyzed quantificationally, and the coefficients of population density on pixel scale were calculated. Last, the model of population spatial distribution at township level was established through multiplicative fusion of population density coefficients and simulated boundary of towns. (3) Error test and analysis of population spatial distribution base on modeling. The author not only analyzed the numerical character of modeling error, but also its spatial distribution. The reasons of error were discussed.

Motion Planning with Six Degrees of Freedom

The motion planning problem is of central importance to the fields of robotics, spatial planning, and automated design. In robotics we are interested in the automatic synthesis of robot motions, given high-level specifications of tasks and geometric models of the robot and obstacles. The Mover's problem is to find a continuous, collision-free path for a moving object through an environment containing obstacles. We present an implemented algorithm for the classical formulation of the three-dimensional Mover's problem: given an arbitrary rigid polyhedral moving object P with three translational and three rotational degrees of freedom, find a continuous, collision-free path taking P from some initial configuration to a desired goal configuration. This thesis describes the first known implementation of a complete algorithm (at a given resolution) for the full six degree of freedom Movers' problem. The algorithm transforms the six degree of freedom planning problem into a point navigation problem in a six-dimensional configuration space (called C-Space). The C-Space obstacles, which characterize the physically unachievable configurations, are directly represented by six-dimensional manifolds whose boundaries are five dimensional C-surfaces. By characterizing these surfaces and their intersections, collision-free paths may be found by the closure of three operators which (i) slide along 5-dimensional intersections of level C-Space obstacles; (ii) slide along 1- to 4-dimensional intersections of level C-surfaces; and (iii) jump between 6 dimensional obstacles. Implementing the point navigation operators requires solving fundamental representational and algorithmic questions: we will derive new structural properties of the C-Space constraints and shoe how to construct and represent C-Surfaces and their intersection manifolds. A definition and new theoretical results are presented for a six-dimensional C-Space extension of the generalized Voronoi diagram, called the C-Voronoi diagram, whose structure we relate to the C-surface intersection manifolds. The representations and algorithms we develop impact many geometric planning problems, and extend to Cartesian manipulators with six degrees of freedom.

Calculations of the minimal perimeter for N deformable cells of equal area confined in a circle

Cox, S.J. (2006) Calculations of the minimal perimeter for N deformable cells of equal area confined in a circle. Philosophical Magazine Letters. 86:569-578.

Approximately Uniform Random Sampling in Sensor Networks

Recent work in sensor databases has focused extensively on distributed query problems, notably distributed computation of aggregates. Existing methods for computing aggregates broadcast queries to all sensors and use in-network aggregation of responses to minimize messaging costs. In this work, we focus on uniform random sampling across nodes, which can serve both as an alternative building block for aggregation and as an integral component of many other useful randomized algorithms. Prior to our work, the best existing proposals for uniform random sampling of sensors involve contacting all nodes in the network. We propose a practical method which is only approximately uniform, but contacts a number of sensors proportional to the diameter of the network instead of its size. The approximation achieved is tunably close to exact uniform sampling, and only relies on well-known existing primitives, namely geographic routing, distributed computation of Voronoi regions and von Neumann's rejection method. Ultimately, our sampling algorithm has the same worst-case asymptotic cost as routing a point-to-point message, and thus it is asymptotically optimal among request/reply-based sampling methods. We provide experimental results demonstrating the effectiveness of our algorithm on both synthetic and real sensor topologies.

Subset Warping: Rubber Sheeting with Cuts

Image warping, often referred to as "rubber sheeting" represents the deformation of a domain image space into a range image space. In this paper, a technique is described which extends the definition of a rubber-sheet transformation to allow a polygonal region to be warped into one or more subsets of itself, where the subsets may be multiply connected. To do this, it constructs a set of "slits" in the domain image, which correspond to discontinuities in the range image, using a technique based on generalized Voronoi diagrams. The concept of medial axis is extended to describe inner and outer medial contours of a polygon. Polygonal regions are decomposed into annular subregions, and path homotopies are introduced to describe the annular subregions. These constructions motivate the definition of a ladder, which guides the construction of grid point pairs necessary to effect the warp itself.

Limited range coverage problems

Tal como o título indica, esta tese estuda problemas de cobertura com alcance limitado. Dado um conjunto de antenas (ou qualquer outro dispositivo sem fios capaz de receber ou transmitir sinais), o objectivo deste trabalho é calcular o alcance mínimo das antenas de modo a que estas cubram completamente um caminho entre dois pontos numa região. Um caminho que apresente estas características é um itinerário seguro. A definição de cobertura é variável e depende da aplicação a que se destina. No caso de situações críticas como o controlo de fogos ou cenários militares, a definição de cobertura recorre à utilização de mais do que uma antena para aumentar a eficácia deste tipo de vigilância. No entanto, o alcance das antenas deverá ser minimizado de modo a manter a vigilância activa o maior tempo possível. Consequentemente, esta tese está centrada na resolução deste problema de optimização e na obtenção de uma solução particular para cada caso. Embora este problema de optimização tenha sido investigado como um problema de cobertura, é possível estabelecer um paralelismo entre problemas de cobertura e problemas de iluminação e vigilância, que são habitualmente designados como problemas da Galeria de Arte. Para converter um problema de cobertura num de iluminação basta considerar um conjunto de luzes em vez de um conjunto de antenas e submetê-lo a restrições idênticas. O principal tema do conjunto de problemas da Galeria de Arte abordado nesta tese é a 1-boa iluminação. Diz-se que um objecto está 1-bem iluminado por um conjunto de luzes se o invólucro convexo destas contém o objecto, tornando assim este conceito num tipo de iluminação de qualidade. O objectivo desta parte do trabalho é então minimizar o alcance das luzes de modo a manter uma iluminação de qualidade. São também apresentadas duas variantes da 1-boa iluminação: a iluminação ortogonal e a boa !-iluminação. Esta última tem aplicações em problemas de profundidade e visualização de dados, temas que são frequentemente abordados em estatística. A resolução destes problemas usando o diagrama de Voronoi Envolvente (uma variante do diagrama de Voronoi adaptada a problemas de boa iluminação) é também proposta nesta tese.

Área metropolitana de Lisboa: descontinuidades, desenho e planeamento

O conceito de descontinuidade é abordado em função de valores associáveis à cidade num contexto de mudanças de paradigmas relacionados com o planeamento do território. O impacto que as descontinuidades têm na forma, estrutura e fronteiras urbanas é analisado recorrendo à Área Metropolitana de Lisboa como estudo de caso. Para abordar as descontinuidades, recorre-se à análise de tendências territoriais recentes, malhas e tecidos construídos e planos mais directamente relacionados com o desenho urbano. O trabalho empírico culmina com o aprofundamento do estudo de um território do interior da Península de Setúbal. Com os resultados desta investigação pretende-se articular tendências territoriais, territórios construídos e instrumentos de planeamento na definição do conceito de descontinuidade. Apresenta-se ainda como conclusão um modelo de intervenção com vista à resolução de descontinuidades. Esta conclusão pretende ser também um contributo, sob a forma de resposta, para a questão levantada por Federico Oliva, quando se interroga, a propósito da cidade de Milão, sobre o que resta dos planos na Cidade.

Development of numerical methodologies for parameter identification and shape optimization in metal forming simulations

Por parte da indústria de estampagem tem-se verificado um interesse crescente em simulações numéricas de processos de conformação de chapa, incluindo também métodos de engenharia inversa. Este facto ocorre principalmente porque as técnicas de tentativa-erro, muito usadas no passado, não são mais competitivas a nível económico. O uso de códigos de simulação é, atualmente, uma prática corrente em ambiente industrial, pois os resultados tipicamente obtidos através de códigos com base no Método dos Elementos Finitos (MEF) são bem aceites pelas comunidades industriais e científicas Na tentativa de obter campos de tensão e de deformação precisos, uma análise eficiente com o MEF necessita de dados de entrada corretos, como geometrias, malhas, leis de comportamento não-lineares, carregamentos, leis de atrito, etc.. Com o objetivo de ultrapassar estas dificuldades podem ser considerados os problemas inversos. No trabalho apresentado, os seguintes problemas inversos, em Mecânica computacional, são apresentados e analisados: (i) problemas de identificação de parâmetros, que se referem à determinação de parâmetros de entrada que serão posteriormente usados em modelos constitutivos nas simulações numéricas e (ii) problemas de definição geométrica inicial de chapas e ferramentas, nos quais o objetivo é determinar a forma inicial de uma chapa ou de uma ferramenta tendo em vista a obtenção de uma determinada geometria após um processo de conformação. São introduzidas e implementadas novas estratégias de otimização, as quais conduzem a parâmetros de modelos constitutivos mais precisos. O objetivo destas estratégias é tirar vantagem das potencialidades de cada algoritmo e melhorar a eficiência geral dos métodos clássicos de otimização, os quais são baseados em processos de apenas um estágio. Algoritmos determinísticos, algoritmos inspirados em processos evolucionários ou mesmo a combinação destes dois são usados nas estratégias propostas. Estratégias de cascata, paralelas e híbridas são apresentadas em detalhe, sendo que as estratégias híbridas consistem na combinação de estratégias em cascata e paralelas. São apresentados e analisados dois métodos distintos para a avaliação da função objetivo em processos de identificação de parâmetros. Os métodos considerados são uma análise com um ponto único ou uma análise com elementos finitos. A avaliação com base num único ponto caracteriza uma quantidade infinitesimal de material sujeito a uma determinada história de deformação. Por outro lado, na análise através de elementos finitos, o modelo constitutivo é implementado e considerado para cada ponto de integração. Problemas inversos são apresentados e descritos, como por exemplo, a definição geométrica de chapas e ferramentas. Considerando o caso da otimização da forma inicial de uma chapa metálica a definição da forma inicial de uma chapa para a conformação de um elemento de cárter é considerado como problema em estudo. Ainda neste âmbito, um estudo sobre a influência da definição geométrica inicial da chapa no processo de otimização é efetuado. Este estudo é realizado considerando a formulação de NURBS na definição da face superior da chapa metálica, face cuja geometria será alterada durante o processo de conformação plástica. No caso dos processos de otimização de ferramentas, um processo de forjamento a dois estágios é apresentado. Com o objetivo de obter um cilindro perfeito após o forjamento, dois métodos distintos são considerados. No primeiro, a forma inicial do cilindro é otimizada e no outro a forma da ferramenta do primeiro estágio de conformação é otimizada. Para parametrizar a superfície livre do cilindro são utilizados diferentes métodos. Para a definição da ferramenta são também utilizados diferentes parametrizações. As estratégias de otimização propostas neste trabalho resolvem eficientemente problemas de otimização para a indústria de conformação metálica.

Potential of fisheries restocking off the Algarve coast using aquaculture produced marine fish

Tese de dout., Ciências do Mar, da Terra e do Ambiente (Ecologia Marinha), Faculdade de Ciências e Tecnologia, Univ. do Algarve, 2012

Biologia pesqueira e dinâmica populacional de Diplodus vulgaris (geoffr.) e Spondyliosoma cantharus (L.) (pisces, sparidae) na costa sudoeste de Portugal

As espécies da família Sparidae ocupam um lugar de destaque na actividade pesqueira da Costa Sudoeste de Portugal, sendo das principais subsidiárias da pesca artesanal nesta região. Com o presente trabalho pretendeu-se clarificar aspectos básicos da biologia de dois esparídeos, Diplodus vulgaris (safia) e Spondyliosoma cantharus (choupa), nas áreas da ecologia alimentar, reprodução, idade, crescimento e mortalidade; estimar os parâmetros de selectividade das artes de pesca mais importantes sobre estas espécies e avaliar preliminarmente o estado da pescaria destes dois esparídeos. Para atingir estes objectivos cumpriu-se um plano de amostragem de frequências de comprimento nas lotas de Sagres e Sines, de Julho de 1992 a Março de 1994; estudos laboratoriais com amostras recolhidas em Sagres de Dezembro de 1992 a Março de 1994 e um estudo de selectividade de aparelho de anzol e redes de emalhar realizado em 1997/98. Os hábitos alimentares foram estudados a partir da análise de conteúdos estomacais, utilizando-se 3 métodos de avaliação da composição da dieta, 5 índices combinados da bibliografia e 3 propostos. Da análise comparativa entre os vários índices concluiu-se que o índice de importância relativa (IRI) e o índice de alimentação ponderado (IPO2) serão os mais adequados para descrever as dietas destas espécies. Registaram-se elevados coeficientes de vacuidade, motivados em grande parte pela natureza da arte de pesca. A dieta de D. vulgaris é composta essencialmente por ofiurídeos, poliquetas, anfípodes e equinóides, enquanto que S. cantharus consome principalmente poliquetas, anfípodes e hidrozoários. As duas dietas são significativamente diferentes, principalmente pelo maior consumo de equinodermes no caso de D. vulgaris. A estrutura e dimensões da boca nestes esparídeos são distintas, podendo influenciar o seu comportamento trófico. À semelhança de grande parte dos esparídeos, D. vulgaris e S. cantharus são espécies bentívoras, que adoptam uma estratégia alimentar generalista, com fortes afinidades com o alimento disponível no meio. Da análise do ciclo reprodutivo verificou-se que as épocas de postura são extensas: D. vulgaris: Dezembro a Março com maior intensidade em Janeiro/Fevereiro; S. cantharus: Fevereiro a Abril, com máxima intensidade em Março. Não existiram diferenças significativas na proporção entre os sexos (M:F = 1.01) ao longo do ano e por tamanho em D. vulgaris, enquanto que em S. cantharus as fêmeas foram mais abundantes ao longo do ano (M:F=0.57) e nas classes de comprimento inferiores. Os tamanhos de 1ª maturação (L50%) para o conjunto dos sexos e indivíduos indeterminados foram de: D. vulgaris, 18.33 cm, não existindo diferenças significativas entre machos e fêmeas; S. cantharus, 20.10 cm, existindo diferenças significativas entre machos (22.41 cm, TL) e fêmeas (19.98 cm, TL). Os tamanhos de 1ª maturação estimados foram consideravelmente superiores aos tamanhos mínimo legais de desembarque (TML) em Portugal, para D. vulgaris (TML=15.0 cm) e ligeiramente inferiores no caso de S. cantharus (TML=23 cm). A fecundidade média absoluta (Fa) e a fecundidade relativa (Fr) foram de: D. vulgaris - Fa = 131127 ovos; Fr = 526 ovos/g; S. cantharus - Fa = 61396 ovos; Fr = 346 ovos/g. As relações entre a fecundidade absoluta e o comprimento total (TL) e o peso somático (SW) foram as seguintes: D. vulgaris: Fa = 0.1853TL4.1903; Fa = 57.715SW1.4067; S. cantharus: Fa = 436.27TL1.5747; Fa = 2979.7SW0.585. A estratégia reprodutiva destas espécies é caracterizada por hermafroditismo, rudimentar com eventual protândria parcial em D. vulgaris e protogínia em S. cantharus. A determinação da idade foi efectuada pela análise de estruturas duras, otólitos para D. vulgaris e otólitos e escamas para S. cantharus. As estimações do comprimento à idade ajustaram-se bem ao modelo de Von Bertalanffy. De entre os métodos de estimação dos parâmetros de crescimento, escolheu-se o ajuste não linear aplicado a todos os pares de comprimentos à idade, sendo esta corrigida pela data de captura. Para além do maior rigor desta análise, os resultados produzidos foram considerados os mais aproximados da realidade: D. vulgaris (otólitos): L=28.1 cm K=0.30 ano-1 e t0= -1.618 anos (validade: 12.5-30.5 cm; 1-10 anos); S. cantharus (escamas): L=35 cm, K=0.32 ano-1 e t0=-0.481 anos (validade: 14.3-33.5 cm; 1-9 anos). A análise de distribuição de frequências de comprimento, apresentou um L mais próximo da realidade: L=39.6 cm, K=0.32 ano-1 e t0=-0.481 anos (D. vulgaris); L=40.0 cm, K=0.24 ano-1 e t0=-0.646 anos (S. cantharus). As relações peso-comprimento de D. vulgaris e S. cantharus para a Costa Sudoeste foram do tipo potencial W=a Lb, sendo definidas pelos seguintes parâmetros: a=0.0223 e b= 2.895, r2=0.89 (D. vulgaris); a=0.0106 e b= 3.085, r2= 0.89 (S. cantharus). Os valores de mortalidade natural (M), estimados por métodos indirectos, para D. vulgaris e S. cantharus foram de 0.39 e 0.30 ano-1, respectivamente. O valor de M de esparídeos pode ser estimado preliminarmente mediante a utilização de uma regressão multilinear que integra as três variáveis (M, K e L): M = -0.162 + 1.714K + 0.00273L (r2=0.77). A análise das curvas de captura permitiu obter os seguintes valores de mortalidade total (Z): D. vulgaris: Z=0.642 (otólitos) e 0.727 ano-1 (frequências de comprimento); S. cantharus: Z=0.676 (escamas) e 0.576 ano-1 (frequências de comprimento). As artes de pesca mais importantes para D vulgaris e S. cantharus e para a região considerada são, por ordem de grandeza, o aparelho de anzol e redes de emalhar, e em menor escala a arte do cerco. Os esparídeos diversos e dentro destes a choupa, que deveria ser discriminada nas estatísticas oficiais, constituíram os únicos esparídeos que apresentaram uma evolução com tendência negativa de 1987 a 1998. A estrutura demográfica de desembarques de safia é caracterizada por cerca de 86.9 % indivíduos entre 20.5 a 27.5 cm (2 a 6 anos), enquanto que para a choupa, 85.0% dos indivíduos desembarcados estão compreendidos entre 21.5 e 27.5 cm (2 a 6 anos). No estudo de selectividade foram realizadas 40 pescas experimentais em que se testaram 4 tamanhos de anzol (nº11, 12, 13 e 15) e de malha (80, 70, 60 e 50mm). D. vulgaris e S. cantharus são das espécies mais abundantes na pescaria com aparelho de anzol, tendo uma representação inferior quer em termos relativos, quer absolutos nas capturas das redes de emalhar. Existiu um decréscimo geral das taxas de captura com o tamanho do anzol e uma maior eficiência para os tamanhos intermédios da malha nas redes de emalhar. Paralelamente, existiu para ambas as espécies uma sobreposição de comprimentos médios com o tamanho dos anzóis e uma separação clara de comprimentos médios de ambas as espécies com a malhagem. Adoptaram-se, com bom ajuste, uma curva de selectividade do tipo logística para o aparelho de anzol e normal para as redes de emalhar, para ambas as espécies. Atendendo a uma eventual revisão dos tamanhos mínimos legais de desembarque, em função dos tamanhos de primeira maturação e à estrutura das capturas, as malhas de 70 e 80mm são as mais adequadas na pescaria de D. vulgaris e S. cantharus, respectivamente. O anzol número 13 será o mais apropriado a D. vulgaris, dada a pouca praticabilidade do anzol 15. Em relação a S. cantharus o anzol nº 11 apresenta a menor proporção de indivíduos com tamanho ilegal, sendo então o mais indicado à pescaria desta espécie. Para o aparelho de anzol e de acordo com os modelos de rendimento por recruta a pescaria de safia processa-se de forma moderada, a níveis inferiores a uma taxa de exploração máxima (Emax) suportável pelo "stock". A pescaria de choupa apresenta indícios de uma exploração intensa, próxima dos valores máximos comportados pelo "stock". Uma alteração nos valores actuais da taxa de exploração e/ou de comprimento de 1ª captura terá de ter em consideração o carácter multiespecifico desta pescaria. Visando a gestão destes recursos pesqueiros, dever-se-á controlar principalmente o nível do esforço de pesca praticado e rever os tamanhos mínimos de desembarque, de modo a evitar que esta pescaria entre numa fase de sobreexploração. O estabelecimento de reservas marinhas e de recifes artificiais poderá ser um contributo para uma exploração pesqueira sustentada e para a conservação da biodiversidade.