Método numérico de Matriz Resposta acoplado a um esquema de reconstrução espacial analítica para cálculos unidimensionais de transporte de nêutrons na formulação de ordenadas discretas multigrupo de energia com fonte fixa

Um método de Matriz Resposta (MR) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas multigrupo de transporte de nêutrons com fonte fixa e em geometria unidimensional na formulação de ordenadas discretas (SN). Portanto, o método multigrupo MR com esquema iterativo de inversão nodal parcial (NBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações multigrupo SN, afora os erros de arredondamento da aritmética finita computacional. É também desenvolvido um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em cada grupo de energia em um intervalo qualquer do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.

Fish assemblages associated with three types of artificial reefs: density of assemblages and possible impacts on adjacent fish abundance

We evaluated the effectiveness of wooden artificial reefs (ARs) as fish habitat. Three types of ARs, made of cedar logs, broadleaf tree logs, and PVC pipes, respectively, were deployed in triplicate at 8-m depth off Maizuru, Kyoto Prefecture, Sea of Japan, in May 2004. Fish assemblages associated with each of the nine ARs were observed by using SCUBA twice a month for four years. Fish assemblages in the adjacent habitat were also monitored for two years before and four years after reef deployment. In the surveyed areas (ca. 10 m2) associated with each of the cedar, broadleaf, and PVC ARs, the average number of fish species was 4.14, 3.49, and 3.00, and the average number of individuals was 40.7, 27.9, and 20.3, respectively. The estimated biomass was also more greater when associated with the cedar ARs than with other ARs. Visual censuses of the habitat adjacent to the ARs revealed that the number of fish species and the density of individuals were not affected by the deployment of the ARs. Our results support the superiority of cedar as an AR material and indicate that deployment of wooden ARs causes no reduction of fish abundance in adjacent natural reefs.

Evaluation of the capture efficiency and size selectivity of four pot types in the prospective fishery for North Pacific giant octopus (Enteroctopus dofleini)

Over 230 metric tons of octopus is harvested as bycatch annually in Alaskan trawl, long-line, and pot fisheries. An expanding market has fostered interest in the development of a directed fishery for North Pacific giant octopus (Enteroctopus dofleini). To investigate the potential for fishery development we examined the efficacy of four different pot types for capture of this species. During two surveys in Kachemak Bay, Alaska, strings of 16 –20 sablefish, Korean hair crab, shrimp, and Kodiak wooden lair pots were set at depths ranging between 62 and 390 meters. Catch per-unit-of-ef for t estimates were highest for sablefish and lair pots. Sablefish pots caught significantly heavier North Pacific giant octopuses but also produced the highest bycatch of commercially important species, such as halibut (Hippoglossus stenolepis), Pacific cod (Gadus macrocephalus), and Tanner crab (Chionoecetes bairdi).

Modelagem computacional da estabilização de sistemas subcríticos segundo o modelo unidimensional de difusão de nêutrons monoenergéticos

Um método numérico espectronodal (END) livre de erros de truncamento espacial é desenvolvido para problemas unidimensionais de difusão de nêutrons monoenergéticos em duas versões. Na versão de problemas de autovalor, o método gera soluções numéricas para o perfil do fluxo escalar e para o fator de multiplicação efetivo (k), que coincidem com a solução analítica dominante, afora os erros da aritmética finita computacional. Na versão de fonte fxa, o método também gera soluções numéricas analíticas para o problema de fonte fixa correspondente, onde a fonte de fissão, com dependência espacial, é obtida analiticamente, a partir da reconstrução espacial do fluxo escalar gerado pelo método END para problemas de autovalor. Alguns experimentos numéricos são apresentados para dois problemas modelos a fim de ilustrar a precisão do método.

Analysis of multivariate stochastic signals sampled by on-line particle analyzers: Application to the quantitative assessment of occupational exposure to NOAA in multisource industrial scenarios (MSIS)

In multisource industrial scenarios (MSIS) coexist NOAA generating activities with other productive sources of airborne particles, such as parallel processes of manufacturing or electrical and diesel machinery. A distinctive characteristic of MSIS is the spatially complex distribution of aerosol sources, as well as their potential differences in dynamics, due to the feasibility of multi-task configuration at a given time. Thus, the background signal is expected to challenge the aerosol analyzers at a probably wide range of concentrations and size distributions, depending of the multisource configuration at a given time. Monitoring and prediction by using statistical analysis of time series captured by on-line particle analyzers in industrial scenarios, have been proven to be feasible in predicting PNC evolution provided a given quality of net signals (difference between signal at source and background). However the analysis and modelling of non-consistent time series, influenced by low levels of SNR (Signal-Noise Ratio) could build a misleading basis for decision making. In this context, this work explores the use of stochastic models based on ARIMA methodology to monitor and predict exposure values (PNC). The study was carried out in a MSIS where an case study focused on the manufacture of perforated tablets of nano-TiO2 by cold pressing was performed

Análise espectral das equações de transporte de nêutrons na formulação de ordenadas discretas em meios multiplicativos

É presentada nesta dissertação uma análise espectral das equações de transporte de nêutrons, independente do tempo, em geometria unidimensional e bidimensional, na formulação de ordenadas discretas (SN), utilizando o modelo de uma velocidade e multigrupo, considerando meios onde ocorrem o fenômeno da fissão nuclear. Esta análise espectral constitui-se na resolução de problemas de autovalores e respectivos autovetores, e reproduz a expressão para a solução geral analítica local das equações SN (para geometria unidimensional) ou das equações nodais integradas transversalmente (geometria retangular bidimensional) dentro de cada região homogeneizada do domínio espacial. Com a solução geral local determinada, métodos numéricos, tais como os métodos de matriz de resposta SN, podem ser derivados. Os resultados numéricos são gerados por programas de computadores implementados em MatLab, versão 2012, a fim de verificar a natureza dos autovalores e autovetores correspondentes no espaço real ou complexo.

Análise dinâmica não-linear de pisos mistos (aço-concreto) submetidos a ações humanas rítmicas.

Atualmente, o crescimento dos problemas de vibrações excessivas sobre pisos mistos (aço-concreto) tem conduzido à necessidade de desenvolvimento de critérios específicos para projetos estruturais submetidos à ação de atividades humanas rítmicas. Com base no desenvolvimento desta dissertação de mestrado, objetiva-se, principalmente, verificar a influência das ligações estruturais (ligações viga-viga), sobre a resposta dinâmica não-linear de pisos mistos (aço-concreto) de edificações, quando submetidos a cargas dinâmicas humanas rítmicas. Deste modo, o carregamento dinâmico empregado para a simulação das atividades humanas sobre o modelo estrutural investigado foi obtido através de testes experimentais com indivíduos praticando atividades rítmicas e não rítmicas. O modelo analisado nesta dissertação corresponde a um piso misto (aço-concreto) com uma área total de 1600m2 e consiste de um ambiente onde serão desenvolvidas atividades de ginástica aeróbica. O sistema estrutural é constituído por lajes de concreto armado apoiadas sobre vigas de aço, simulando o comportamento de um sistema estrutural misto (aço-concreto) com interação total. A metodologia de análise desenvolvida emprega técnicas usuais de discretização presentes no método dos elementos finitos, com base no emprego do programa ANSYS. A modelagem do sistema contempla ligações estruturais do tipo rígidas, semirrígidas e flexíveis. Os valores das acelerações de pico foram comparados com os limites recomendados por normas de projeto, baseando-se em critérios de conforto humano. As conclusões alcançadas ao longo deste trabalho de pesquisa revelam que as ligações estruturais do tipo viga-viga não apresentam influência significativa, no que diz respeito a resposta dinâmica não-linear da estrutura. Por outro lado, as acelerações de pico obtidas com base na análise dinâmica não-linear apresentam valores elevados indicando que o piso misto (aço-concreto) investigado apresenta problemas de vibração excessiva inerentes ao conforto humano.

Análise dinâmica e controle de vibrações de passarelas de pedestres submetidas ao caminhar humano.

Passarelas de pedestres mistas (aço-concreto) e de aço são frequentemente submetidas a ações dinâmicas de magnitude variável, devido à travessia de pedestres sobre a laje de concreto. Estas ações dinâmicas podem produzir vibrações excessivas e dependendo de sua magnitude e intensidade, estes efeitos adversos podem comprometer a confiabilidade e a resposta do sistema estrutural e, também, podem levar a uma redução da expectativa de vida útil da passarela. Por outro lado, a experiência e o conhecimento dos engenheiros estruturais em conjunto com o uso de novos materiais e tecnologias construtivas têm produzido projetos de passarelas mistas (aço-concreto) bastante arrojados. Uma consequência direta desta tendência de projeto é um aumento considerável das vibrações estruturais. Com base neste cenário, esta dissertação visa investigar o comportamento dinâmico de três passarelas de pedestres mistas (aço-concreto) localizadas no Rio de Janeiro, submetidas ao caminhar humano. Estes sistemas estruturais são constituídos por uma estrutura principal de aço e laje em concreto e são destinados à travessia de pedestres. Deste modo, foram desenvolvidos modelos numérico-computacionais, adotando-se as técnicas tradicionais de refinamento presentes em simulações do método de elementos finitos, com base no uso do software ANSYS. Estes modelos numéricos permitiram uma completa avaliação dinâmica das passarelas investigadas, especialmente em termos de conforto humano. As respostas dinâmicas foram obtidas em termos de acelerações de pico e comparadas com valores limites propostas por diversos autores e normas de projeto. Os valores de aceleração de pico e aceleração rms encontrados na presente investigação indicaram que as passarelas analisadas apresentaram problemas relacionados com o conforto humano. Assim sendo, considerando-se que foi detectado que estas estruturas poderiam atingir níveis elevados de vibração que possam vir a comprometer o conforto dos usuários, foi verificado que uma estratégia para o controle estrutural era necessária, a fim de reduzir as vibrações excessivas nas passarelas. Finalmente, uma investigação foi realizada com base em alternativas de controle estrutural objetivando atenuar vibrações excessivas, a partir do emprego de sistemas de atenuadores dinâmicos sintonizados (ADS).

The Fisheries for Mangrove Cockles,Anadaraspp., from Mexico to Peru, With Descriptions of Their Habitats and Biology, the Fishermen’s Lives, and the Effects of Shrimp Farming

This paper provides the first description of the mangrove cockle, Anadara spp., fisheries throughout their Latin American range along the Pacific coast from Mexico to Peru. Two species, A. tuberculosa and A. grandis, are found over the entire range, while A. similis occurs from El Salvador to Peru. Anadara tuberculosa is by far the most abundant, while A. grandis has declined in abundance during recent decades. Anadara tuberculosa and A. similis occur in level mud sediments in mangrove swamps, comprised mostly of Rhizophora mangle, which line the main-lands and islands of lagoons, whereas A. grandis inhabits intertidal mud flats along the edges of the same mangrove swamps. All harvested cockles are sexually mature. Gametogenesis of the three species occurs year round, and juvenile cockles grow rap-idly. Cockle densities at sizes at least 16–42 mm long ranged from 7 to 24/m2 in Mexico. Macrofaunal associates of cockles include crustaceans, gastropods, and finfishes. The mangrove swamps are in nearly pristine condition in every country except Honduras, Ecuador, and Peru, where shrimp farms constructed in the 1980’s and 1990’s have destroyed some mangrove zones. In addition, Hurricane Mitch destroyed some Honduran mangrove swamps in 1998. About 15,000 fishermen, including men, women, and children, harvest the cockles. Ecuador has the largest tabulated number of fishermen, 5,055, while Peru has the fewest, 75. Colombia has a large number, perhaps exceeding that in Ecuador, but a detailed census of them has never been made. The fishermen are poor and live a meager existence; they do not earn sufficient money to purchase adequate food to allow their full health and growth potential. They travel almost daily from their villages to the harvesting areas in wooden canoes and fiberglass boats at low tide when they can walk into the mangrove swamps to harvest cockles for about 4 h. Harvest rates, which vary among countries owing to differences in cockle abundances, range from about 50 cockles/fisherman/day in El Salvador and Honduras to 500–1,000/ fisherman/day in Mexico. The fishermen return to their villages and sell the cockles to dealers, who sell them mainly whole to market outlets within their countries, but there is some exporting to adjacent countries. An important food in most countries, the cockles are eaten in seviche, raw on the half-shell, and cooked with rice. The cockles are under heavy harvesting pressure, except in Mexico, but stocks are not yet being depleted because they are harvested at sizes which have already spawned. Also some spawning stocks lie within dense mangrove stands which the fishermen cannot reach. Consumers fortunately desire the largest cockles, spurning the smallest. Cockles are important to the people, and efforts to reduce the harvests to prevent overfishing would lead to severe economic suffering in the fishing communities. Pro-grams to conserve and improve cockle habitats may be the most judicious actions to take. Preserving the mangrove swamps intact, increasing their sizes where possible, and controlling cockle predators would lead to an increase in cockle abundance and harvests. Fishes that prey on juvenile cockles might be seined along the edges of swamps before the tide rises and they swim into the swamps to feed. Transplanting mangrove seedlings to suitable areas might increase the size of those habitats. The numbers of fishermen may increase in the future, because most adults now have several children. If new fishermen are tempted to harvest small, immature cockles and stocks are not increased, minimum size rules for harvestable cockles could be implemented and enforced to ensure adequate spawning.

Ensaios de tratabilidade de águas residuárias da indústria de piso de madeira por processos oxidativos avançados: ozônio e UV-C.

A indústria moveleira em geral e de piso de madeira é caracterizada por processos a seco. Em decorrência da limpeza e lavagem de maquinário e superfícies gera águas residuárias em baixos volumes porém altamente contaminadas com substâncias recalcitrantes de difícil tratamento. Com o objetivo de oferecer uma alternativa tecnológica de fácil operação (com uso mínimo de reagentes químicos), capaz de tratar tais águas, reduzir toxicidade e aumentar biodegradabilidade do efluente, viabilizando desta forma uma etapa subsequente de tratamento biológico, o presente trabalho conduziu ensaios de tratabilidade baseados em processos oxidativos avançados com águas residuárias provenientes de lavagens em uma indústria de piso de madeira em Nybro, Suécia, com DQO inicial de aproximadamente 45.000 mg/L e pH de 2,3-2,6. Após tratamento primário in situ com redução de 89% de DQO, os seguintes processos geradores de radicais OH foram testados em escala de laboratório: ozonização sem ajuste de pH (Etapa I); ozonização com ajuste de pH (Etapa II); e ozonização + UV-C com ajuste de pH (Etapa III). Os experimentos foram realizados em semi-batelada, utilizando-se um aparato com dois reatores conectados (de O3 e de UV), entre os quais foi mantido fluxo de recirculação em contracorrente com o gás. Para o desenho experimental, utilizou-se a metodologia de planejamento fatorial do tipo Delineamento Composto Central Rotacional (DCCR) utilizando-se como variáveis independentes, (i) vazão de recirculação entre os reatores (VR) com níveis variando de 0,3 a 3L/min; (ii) concentração de O3 na mistura gasosa injetada ([O3]) ou modo de dosagem cujos níveis variaram de 40 a 115 g/Nm3 e; (iii) pH inicial (Etapas II e III) com níveis que variaram de 3 a 11. A quantidade total de O3 fornecida até o final de cada ensaio foi fixada em 2 gramas de O3 para cada grama de DQO inicial (2:1). A eficiência de cada tratamento foi avaliada em função da redução de DQO, da remoção de COT e do consumo de O3. Para ozonização sem ajuste de pH foram alcançadas reduções de DQO e COT de 65 e 31% respectivamente. No tratamento com ajuste de pH sem UV, foram alcançadas reduções de DQO e COT de 85% e 43% com pHinicial = 7,0 enquanto que no tratamento com UV foi alcançada redução de DQO e COT de 93% e 56% respectivamente, com pHi = 9,38. Os resultados de respirometria com os efluentes dos tratamentos indicaram que a ozonização destas águas em pH baixo pode resultar na geração de subprodutos menos biodegradáveis. Entretanto, tais produtos não são gerados quando a ozonização é realizada em pH mais elevado (7,0), mesmo sem aplicação de UV, com doses de O3:DQOinicial menores que 1:1. O tratamento por ozonização pode, portanto, ser aplicado como uma etapa intermediária de tratamento das águas residuárias em questão, inclusive precedendo o tratamento biológico.

Gross energy requirement in fishing operations

Energy is a key input into the fish harvesting process. Efficient use of energy helps in reducing operational costs and environmental impact, while increasing profits. Energy optimisation is an important aspect of responsible fishing as enunciated in the Code of Conduct for Responsible Fisheries. Gross Energy Requirement (GER) is the sum of all non-renewable energy resources consumed in making available a product or service and is expressed in energy units per physical unit of product or service delivered. GER is a measure of intensity of non-renewable resource use and it reflects the amount of depletion of earth’s inherited store of non-renewable energy in order to create and make available a product or service. In this study, GER in fish harvesting up to the point of landing is estimated in selected fish harvesting systems in the small-mechanised sectors of Indian fisheries and compared with reported results from selected non mechanised and motorised fishing systems to reflect the situation during 1997-1998. Among the fish harvesting systems studied, GER t fish-1 ranged from 5.54 and 5.91 GJ, respectively, for wooden and steel purse seiners powered by 156 hp engines; 6.40 GJ for wooden purse seiner with 235 hp engine; 25.18 GJ for mechanised gillnet/line fishing vessel with 89 hp engines; to 31.40 and 36.97 GJ, respectively, for wooden and steel trawlers powered by 99-106 hp engines.

Uso de web services e softwares livres na disseminação de informações georreferenciadas sobre produtos da floresta amazônica.

Georreferenced information has been increasingly required for the planning and decision-making in different sectors of society. New ways of dissemination of data, such as the Open Geospatial Consortium (OGC) web services, have contributed to the ease of access to this information. Even with all the technological advances in the area of data distribution, there is still low availability of georreferenced data about the Amazon. The goal of the present work is the development of a spatial data infrastructure (SDI), that is, an environment of sharing and use of georreferenced data based on the technology of web services, metadata and interfaces that allow the user easy access to these data. The present work discussess the OGC patterns, the most relevant georeferrenced data servers, the main web clients, and the revolution in the dissemination of georeferrenced data which geobrowsers and web clients offered to regular users. Data to be released for the case study come from the project Exploitation of Non-wooden Forest Products-PFNM-in progress at the National Institute of Research in the Amazon-INPA-as well as from inventories of NGOs and other government bodies. Besides contributing to the enhancement of PFNM, this project aims at encouraging the use of GIS in the state of Amazonas offering tech support for the deployment of geographic databases and sharing between agencies, optimizing the resources applied in this area through the use of free software and integration of diffuse information currently available.

Modelagem do comportamento dinâmico e análise de fadiga de pontes rodoviárias mistas (aço-concreto).

Considerando-se os diversos carregamentos que solicitam a estrutura de uma ponte rodoviária, ao longo de sua vida útil, alguns possuem um comportamento essencialmente dinâmico, ou seja, variam com o tempo, diferentemente do que é considerado na prática corrente de projeto desse tipo de obra de arte. Em geral, os projetistas desse tipo de estrutura têm tratado carregamentos como o de vento, tráfego de veículos e de pedestres como ações de natureza estática, ignorando seu perfil cíclico. Tal consideração, em diversas situações de projeto, tende a minorar os efeitos das ações dinâmicas sobre o sistema estrutural. Além disso, estruturas submetidas a solicitações cíclicas, sob uma carga inferior à caga máxima suportada pelo material, estão sujeitas ao fenômeno da fadiga. A consideração adequada de todos estes aspectos mostra-se fundamentalmente importante para correta avaliação dos níveis de esforços solicitantes do sistema estrutural e, bem como, para a identificação de fenômenos importantes como o da fadiga que pode vir a provocar, por exemplo, a ruptura de componentes estruturais sem aviso prévio ou motivo aparente. Para tal, nesta dissertação as técnicas para a contagem de ciclos de tensão e a aplicação das regras de dano acumulado foram analisadas através de curvas do tipo S-N, associadas a diversas normas de projeto. A ponte rodoviária mista (aço-concreto) investigada neste estudo é constituída por seis vigas de aço longitudinais com enrijecedores transversais, oito transversinas e por um tabuleiro de concreto armado. O modelo numérico-computacional, desenvolvido para a análise dinâmica da ponte, foi elaborado com base em técnicas usuais de discretização através do método dos elementos finitos. As mesas e almas das vigas, assim como os enrijecedores , foram modelados por elementos de casca e laje de concreto armado, por elementos sólidos. O carregamento dinâmico avaliado no presente estudo diz respeito ao tráfego de veículos, cuja representação se dá a partir de sistemas "massa-mola-amortecedor". Os comboios formados são adotados como sendo semi-infinitos, deslocando-se com velocidade constante sobre a ponte. As conclusões da presente investigação versam acerca da vida útil de serviço dos elementos estruturais de pontes mistas (aço-concreto).

Condições de contorno tipo albedo para cálculos globais de reatores nucleares com o modelo de dois grupos de energia na formulação de transporte SN em geometria bidimensional cartesiana

Os eventos de fissão nuclear, resultados da interação dos nêutrons com os núcleos dos átomos do meio hospedeiro multiplicativo, não estão presentes em algumas regiões dos reatores nucleares, e.g., moderador, refletor, e meios estruturais. Nesses domínios espaciais não há geração de potência nuclear térmica e, além disso, comprometem a eficiência computacional dos cálculos globais de reatores nucleares. Propomos nesta tese uma estratégia visando a aumentar a eficiência computacional dessas simulações eliminando os cálculos numéricos explícitos no interior das regiões não-multiplicativas (baffle e refletor) em torno do núcleo ativo. Apresentamos algumas modelagens e discutimos a eficiência da aplicação dessas condições de contorno aproximadas tipo albedo para uma e duas regiões nãomultiplicativas, na formulação de ordenadas discretas (SN) para problemas de autovalor a dois grupos de energia em geometria bidimensional cartesiana. A denominação Albedo, palavra de origem latina para alvura, foi originalmente definida como a fração da luz incidente que é refletida difusamente por uma superfície. Esta denominação latina permaneceu como o termo científico usual em astronomia e, nesta tese, este conceito é estendido para reflexão de nêutrons. Estas condições de contorno tipo albedo SN não-convencional substituem aproximadamente as regiões de baffle e refletor no em torno do núcleo ativo do reator, desprezando os termos de fuga transversal no interior dessas regiões. Se o problema, em particular, não possui termos de fuga transversal, i.e., trata-se de um problema unidimensional, então as condições de contorno albedo, como propostas nesta tese, são exatas. Por eficiência computacional entende-se a análise da precisão dos resultados numéricos em comparação com o tempo de execução computacional de cada simulação de um dado problema-modelo. Resultados numéricos considerando dois problemas-modelo com de simetria são considerados para ilustrar esta análise de eficiência.

Estimativa de fontes externas uniformes de nêutrons em regiões combustíveis para uma distribuição prescrita de potência

Projetos de reatores nucleares foram classificados em quatro gerações (Gen) pelo Departamento de Energia dos Estados Unidos da América (DOE), quando o DOE introduziu o conceito de reatores de geração IV (Gen IV). Reatores Gen IV são um conjunto de projetos de reator nuclear, em sua maioria teóricos, atualmente sendo pesquisados. Entre os projetos Gen IV, incluem-se os projetos dos ADS (Accelerator Driven Systems), que são sistemas subcríticos estabilizados por fontes externas estacionárias de nêutrons. Estas fontes externas de nêutrons são normalmente geradas a partir da colisão de prótons com alta energia contra os núcleos de metais pesados presentes no núcleo do reator, fenômeno que é conhecido na literatura como spallation, e os prótons são acelerados num acelerador de partículas que é alimentado com parte da energia gerada pelo reator. A criticalidade de um sistema mantido por reações de fissão em cadeia depende do balanço entre a produção de nêutrons por fissão e a remoção por fuga pelos contornos e absorção de nêutrons. Um sistema está subcrítico quando a remoção por fuga e absorção ultrapassa a produção por fissão e, portanto, tende ao desligamento. Entretanto, qualquer sistema subcrítico pode ser estabilizado pela inclusão de fontes estacionárias de nêutrons em seu interior. O objetivo central deste trabalho é determinar as intensidades dessas fontes uniformes e isotrópicas de nêutrons, que se deve inserir em todas as regiões combustíveis do sistema, para que o mesmo estabilize-se gerando uma distribuição prescrita de potência elétrica. Diante do exposto, foi desenvolvido neste trabalho um aplicativo computacional em linguagem Java que estima as intensidades dessas fontes estacionárias de nêutrons, que devem ser inseridas em cada região combustível para que estabilizem o sistema subcrítico com uma dada distribuição de potência definida pelo usuário. Para atingir este objetivo, o modelo matemático adotado foi a equação unidimensional de transporte de nêutrons monoenergéticos na formulação de ordenadas discretas (SN) e o convencional método de malha fina diamond difference (DD) foi utilizado para resolver numericamente os problemas SN físicos e adjuntos. Resultados numéricos para dois problemas-modelos típicos são apresentados para ilustrar a acurácia e eficiência da metodologia proposta.