Formação de biofilmes por Enterococcus faecium sensível e multirresistente aos antimicrobianos: características à proteômica

Enterococcus faecium tem se destacado no cenário das infecções hospitalares, particularmente, amostras portadoras de características de multirresistência. Apesar de serem responsabilizados como agentes etiológicos em diferentes quadros clínicos, fatores associados a patogênese das infecções ainda não estão esclarecidos. Entretanto, sabe-se que a capacidade de formação de biofilmes pode ser responsabilizada como um dos atributos capazes de promover o papel patogênico desses microrganismos. Assim sendo, este estudo se propôs a investigar a capacidade de formação de biofilmes por duas amostras de E. faecium: SS-1274 (derivada da amostra tipo da espécie) e CL-6729 (multirresistente e pertencente a um complexo clonal globalmente disperso). As amostras foram caracterizadas fenotipica e genotipicamente para confirmação da identificação, determinação da susceptibilidade a 17 antimicrobianos, caracterização da concentração mínima inibitória para ampicilina, gentamicina e vancomicina, determinação do genótipo vanA e detecção de genes associados a expressão dos genes asa1, cylA, esp, gelE e hyl. A análise quantitativa da formação por 24h e 72h dos biofilmes foi realizada por metodologia do cristal violeta, em placas de microtitulação de poliestireno. Foram construídas curvas de formação pela avaliação da DO570nm das preparações coradas por cristal violeta em períodos de tempo de 2h a 74h. A influência de subCMIs de ampicilina, gentamicina e vancomicina na formação de biofilmes de E. faecium foi também caracterizada em ensaios de quantificação da biomassa. A presença de DNA e proteínas foi avaliada em ensaios de destacamento e por fluorimetria. A arquitetura, distribuição espacial e reação aos fluorocromos Syto9 e SYPRO Ruby Protein foram evidenciadas por microscopia confocal de varredura a laser. Análises proteômicas através da avaliação por SDS-PAGE e por ESI-Q-TOF também foram empregadas para avaliação de biofilmes versus crescimento planctônico. Nossos resultados demonstraram que a amostra CL-6729 apresentou uma maior biomassa nos tempos analisados. Apesar da menor quantidade de biomassa da amostra SS-1274, o perfil da curva de formação foi semelhante a CL-6729. As análises da matriz por ensaios quantitativos e microscopia confocal revelaram que proteína parece ser um importante constituinte para ambas as amostras, entretanto biofilmes formados na presença de da CMI e durante 24h, parecem sofrer alterações na constituição. Os resultados relativos às espectrometrias de massa sugerem que células de biofilmes de E. faecium podem também estar metabolicamente ativas, devido a identificação de um considerável número de proteínas relacionadas ao metabolismo e divisão celular, similarmente às células planctônicas. No entanto, investigações complementares são necessárias para quantificar as diferenças relacionadas a sua expressão. Foi observado que ambas as amostras independente das variações fenotípicas e genotípicas são capazes de formar biofilmes maduros exibindo constituição e arquitetura similares. Entretanto, nossos resultados sugeriram que cada amostra responde as diferentes situações de acordo com seus determinantes de virulência e resistência.

Produção de plantas e metabólitos de Cleome dendroides Schult. & Schult. f. (Cleomaceae) utilizando diferentes sistemas de cultivo in vitro

Cleome dendroides é uma espécie endêmica da Mata Atlântica dos estados do Rio de Janeiro e Espírito Santo, bioma alterado por intensa atividade antrópica, o que constitui uma ameaça à preservação de suas populações. Não existem estudos dos pontos de vista fisiológico, biotecnológico, fitoquímico ou farmacológico sobre a espécie. Considerando-se o perfil fitoquímico e o potencial medicinal do gênero, torna-se relevante definirem-se protocolos para a produção de plantas e metabólitos de C. dendroides utilizando diferentes sistemas de cultivo in vitro. No presente trabalho, foram realizados estudos sobre a germinação in vivo da espécie, avaliando-se a influência do substrato, temperatura e luz. Não se observou qualquer tipo de dormência, sendo as temperaturas de 20, 25 e 20-30C, em areia ou vermiculita, apropriadas para a germinação in vivo. Definiu-se também uma metodologia eficiente de germinação sob condições in vitro, e as plântulas obtidas foram utilizadas como fonte de explantes para os estudos de propagação in vitro. A resposta morfogênica foi avaliada considerando-se a origem e tipo do explante, tipos e concentrações de reguladores de crescimento e número de subculturas. A metodologia empregada mostrou-se eficiente para a produção de brotos e manutenção de estoques de plantas in vitro que serviram como fonte de explantes. A melhor condição para a propagação in vitro foi definida em meio solidificado contendo BA, independentemente do tipo de explante e da origem. Os brotos obtidos foram alongados, enraizados e aclimatizados. Também foi estabelecida a cultura de raízes e a regeneração de brotos a partir destas culturas. Avaliou-se o efeito da origem do explante, assim como dos tipos e concentrações de fitorreguladores sobre a proliferação de raízes e a regeneração de brotos. O fitorregulador AIB propiciou maior multiplicação das raízes, enquanto BA mostrou-se eficiente na regeneração de brotos a partir das raízes recém-formadas. Foi estabelecido ainda um protocolo de cultura de calos e de suspensões celulares. Avaliou-se o efeito da origem e do tipo de explante, dos tipos e das concentrações de fitorreguladores sobre a calogênese. A combinação de PIC com KIN foi a mais eficiente para a indução de calos em explantes de plântulas obtidas a partir de germinação in vitro, produzindo calos que foram mantidos por pelo menos dois anos. As suspensões celulares também foram estabelecidas em meio contendo PIC + KIN, mantendo uma produção de biomassa de cerca de cinco vezes o peso fresco inicial por três sucessivas subculturas. Análises histoquímicas e fitoquímicas revelaram a presença de alcaloides nos calos e nas suspensões celulares. Foram realizadas análises fitoquímicas de plantas de campo, plantas aclimatizadas, plantas mantidas em estoque in vitro e culturas de raízes, as quais indicaram a presença de derivados de glicosinolatos. Os resultados demonstraram a viabilidade de produção de material vegetal de C. dendroides por meio de métodos biotecnológicos e a produção in vitro de metabólitos de importância medicinal

Uma Simulação computacional do passeio aleatório simples

Em 1828 foi observado um fenômeno no microscópio em que se visualizava minúsculos grãos de pólen mergulhados em um líquido em repouso que mexiam-se de forma aleatória, desenhando um movimento desordenado. A questão era compreender este movimento. Após cerca de 80 anos, Einstein (1905) desenvolveu uma formulação matemática para explicar este fenômeno, tratado por movimento Browniano, teoria cada vez mais desenvolvida em muitas das áreas do conhecimento, inclusive recentemente em modelagem computacional. Objetiva-se pontuar os pressupostos básicos inerentes ao passeio aleatório simples considerando experimentos com e sem problema de valor de contorno para melhor compreensão ao no uso de algoritmos aplicados a problemas computacionais. Foram explicitadas as ferramentas necessárias para aplicação de modelos de simulação do passeio aleatório simples nas três primeiras dimensões do espaço. O interesse foi direcionado tanto para o passeio aleatório simples como para possíveis aplicações para o problema da ruína do jogador e a disseminação de vírus em rede de computadores. Foram desenvolvidos algoritmos do passeio aleatório simples unidimensional sem e com o problema do valor de contorno na plataforma R. Similarmente, implementados para os espaços bidimensionais e tridimensionais,possibilitando futuras aplicações para o problema da disseminação de vírus em rede de computadores e como motivação ao estudo da Equação do Calor, embora necessita um maior embasamento em conceitos da Física e Probabilidade para dar continuidade a tal aplicação.

Avaliação do reúso agrícola de águas oleosas da indústria de petróleo tratadas por eletrofloculação e osmose inversa: efeitos na germinação, no desenvolvimento de plântulas de girassol e em atributos do solo

A presente tese de doutorado teve como objetivo estabelecer parâmetros para avaliar a viabilidade do reúso agrícola de águas oleosas da indústria de petróleo, como as águas de produção (AP) de campos de exploração onshore do semiárido nordestino, na irrigação de culturas de girassol destinadas à produção de biodiesel. A AP foi produzida sinteticamente e tratada utilizando-se as técnicas de eletrofloculação (EF) e de osmose inversa (OI). Foram analisados os efeitos da AP não tratada, tratada por EF e por EF combinada com osmose inversa (EF+OI) na germinação, desenvolvimento e produção de biomassa de plântulas de girassol e também nos atributos de um solo característico do semiárido nordestino. Na melhor condição operacional do tratamento da AP por EF (28,6 A m-2 durante 4 min.) foram obtidas eficiências de remoção de óleos e graxas (O&G), demanda química de oxigênio (DQO), cor e turbidez superiores a 95%. O tratamento EF+OI promoveu a remoção do excesso de salinidade e de ferro oriundo da etapa de EF, enquadrando-se esses parâmetros dentro de níveis de referência recomendados para água de irrigação. Níveis de O&G e DQO superiores, respectivamente, a 337 mg L-1 e 1.321 mg O2 L-1 na AP bruta produziram efeitos tóxicos, reduzindo-se o índice de velocidade de germinação (IVG) e o percentual de plântulas normais do girassol. Por outro lado, os efluentes tratados por EF e EF+OI produziram efeitos similares no percentual de sementes germinadas, no IVG, no percentual de plântulas normais e na produção de biomassa do girassol. O uso da AP tratada por EF, com ou sem diluição, contribuiu significativamente para o aumento da salinidade e dos teores de sódio do solo, diferentemente da tratada por EF+OI, que produziu efeitos similares ao do controle (água destilada)

Avaliação de modelos de permeabilidade em meios porosos não consolidados

As simulações computacionais tem sido amplamente empregadas no estudo do escoamento darciano e não-darciano em meios porosos consolidados e não-consolidados. Neste trabalho, através de uma nova formulação para a equação de Forchheimer, foram identificadas duas novas propriedades denominados fator de comportamento do fluido, que atua paralelamente a permeabilidade, e permeabilidade equivalente global, resultado da relação anterior. Este comportamento foi estudado e validado através da implementação de um aparato experimental e um código computacional baseado no modelo de regressão-linear que, além disso, demonstrou que o escoamento, ainda que em regime não darciano, comporta-se linearmente como a equação de Darcy, ainda que o coeficiente angular desta diminuiu de acordo com a faixa do número de Reynolds atingida, sendo esta dependente do tipo de leito empregado. Ainda neste trabalho, foi implementado o método de otimização R2W para estimar os parâmetros da equação de Kozeny-Carman a partir de dados experimentais obtidos por Dias et al, a fim de simular o escoamento darciano em meios porosos. Por fim, foi alcançada excelente concordância entre os dados simulados pelo método R2W / equação de Kozeny-Carman e os dados reais.

Experimentos numéricos para deposição de parafinas com modelo Multisólido

Este trabalho objetiva a construção de estruturas robustas e computacionalmente eficientes para a solução do problema de deposição de parafinas do ponto de vista do equilíbrio sólido-líquido. São avaliados diversos modelos termodinâmicos para a fase líquida: equação de estado de Peng-Robinson e os modelos de coeficiente de atividade de Solução Ideal, Wilson, UNIQUAC e UNIFAC. A fase sólida é caracterizada pelo modelo Multisólido. A previsão de formação de fase sólida é inicialmente prevista por um teste de estabilidade termodinâmica. Posteriormente, o sistema de equações não lineares que caracteriza o equilíbrio termodinâmico e as equações de balanço material é resolvido por três abordagens numéricas: método de Newton multivariável, método de Broyden e método Newton-Armijo. Diversos experimentos numéricos foram conduzidos de modo a avaliar os tempos de computação e a robustez frente a diversos cenários de estimativas iniciais dos métodos numéricos para os diferentes modelos e diferentes misturas. Os resultados indicam para a possibilidade de construção de arcabouços computacionais eficientes e robustos, que podem ser empregados acoplados a simuladores de escoamento em dutos, por exemplo.

Comunidades macrofaunais bênticas associadas às algas de deriva e à grama marinha Halodule wrightii Ascherson na Ilha do Japonês, Cabo Frio, RJ

As extensas pradarias submersas formadas pelas gramas marinhas são importantes habitats da costa, onde ocorrem interações ecológicas entre diversas espécies da vegetação subaquática, invertebrados bentônicos e peixes. As gramas marinhas e algas de deriva são conhecidas como macrófitas marinhas e, por ocuparem o mesmo tipo de substrato, são normalmente encontradas juntas, proporcionando oxigênio, alimento, proteção, abrigo além de sítios de reprodução e pastagem para os animais associados a essas pradarias. Amostras de algas de deriva e de H. wrightii foram coletadas, ao longo de transectos fixos de 50 m paralelos à Ilha do Japonês, a fim de analisar a existência de relações positivas entre as espécies de macrófitas marinhas e sua macrofauna associada, comparar as duas comunidades e avaliar a estruturação da comunidade macrofaunal bêntica do local. Os transectos foram alocados de acordo com a posição do banco de grama marinha. Observou-se que a densidade de eixos e a biomassa de H. wrightii não explicam a variação da biomassa, riqueza de espécies e diversidade (Índice de Simpson) das algas de deriva. A grande movimentação das algas de deriva ao longo do banco de grama marinha faz com que elas se homogeneízem e ocupem diferentes lugares ao acaso na pradaria, muitos desses locais com baixa biomassa de H. wrightii devido à grande variabilidade na distribuição dessa espécie no local de estudo. Os descritores ecológicos da grama marinha também não tiveram relações positivas com sua macrofauna bêntica associada. A comunidade macrofaunal associada às gramas marinhas foi mais densa, rica e diversa do que a comunidade macrofaunal associada às algas de deriva. Os moluscos Anomalocardia flexuosa, Cerithium atratum, Ostrea sp, Tellina lineata e Divalinga quadrissulcata dominaram o ambiente de gramas marinhas. A maior complexidade estrutural das algas de deriva forneceu um habitat protegido mais atrativo para os crustáceos como, Pagurus criniticornis, Cymadusa filosa e Batea catharinensis. A malacofauna associada às algas não foi abundante, mas um novo registro foi a ocorrência do bivalve invasor Lithopaga aristatus, perfurando uma concha de Ostrea sp. As relações entre os descritores da biomassa algal foram comprovadas para a maioria dos descritores de sua fauna associada. As relações das macrófitas marinhas com a macrofauna total associada seguiram o mesmo padrão das relações das algas de deriva. As análises de agrupamento e ordenação mostraram que as comunidades macrofaunais bênticas do local são estruturadas de acordo com os táxons dos organismos associados mais dominantes influenciados pelo tipo de vegetação basibionte (algas de deriva ou grama marinha). Destaca-se com o presente estudo a importância de medidas de maior proteção no local para a preservação e manutenção do ecossistema da Ilha do Japonês, RJ, Brasil

Produção in vitro, análise fitoquímica e avaliação da atividade antineoplásica de metabólitos de Hovenia dulcis Thunb. (Rhamnaceae) obtidos por diferentes estratégias biotecnológicas

Hovenia dulcis Thunberg, natural da Ásia Oriental, é cultivada no Brasil onde é conhecida como uva-do-japão. A espécie possui várias indicações na medicina popular e alguns estudos apontam o seu potencial antineoplásico, tripanocida e hepatoprotetor. Metabólitos secundários são substâncias não essenciais para a sobrevivência celular, mas que fornecem vantagens adaptativas aos vegetais, sendo atribuído, para algumas delas, atividades biológicas importantes. Substâncias de interesse medicinal têm sido obtidas por técnicas da cultura de tecidos vegetais, como a calogênese e a cultura de células em suspensão, que permitem a síntese de matéria-prima de forma contínua e homogênea, independentemente de fatores ambientais e sazonais. O presente estudo objetivou o estabelecimento de culturas in vitro de H. dulcis, visando à produção de metabólitos de interesse, com vistas à avaliação do seu potencial antineoplásico sobre células K562. Foram testados protocolos para o estabelecimento de diferentes sistemas, como culturas de calos, de células em suspensão (CCS) e compact callus clusters (CCC) e ainda a avaliação do uso de elicitores na otimização de metabólitos produzidos in vitro. Foi verificado que a adição dos fitorreguladores KIN e TDZ, substituindo o BAP, não foi capaz de induzir a formação de calos friáveis, bem como a manutenção das culturas em ausência de luz. O uso do nitrato de prata promoveu a friabilidade de calos em todas as concentrações testadas, considerando-se 2,0 mg.L-1 a melhor concentração. Foram alcançadas taxas de 100% de formação de CCS tanto na presença, quanto em ausência de AgNO3. O maior acúmulo de biomassa foi verificado na concentração mais baixa de PIC (0,625 mg.L-1). A análise dos espectros de RMN indicou a presença de (+)-dihidromiricetina, (+)-galocatequina, hovenitina II, hovenosideo G, hodulosideo III, hodulosideo IV, hodulosideo I e hovenidulciosideo B1 nas culturas de calos friáveis. No estabelecimento de culturas CCC, observou-se a formação de calos compactos verdes em todas as concentrações de ANA testadas. O aumento da velocidade de rotação para 135 rpm aumentou a dispersão das células com consequente formação dos agregados celulares desejados. A seleção de linhagens celulares demonstrou ser um método eficiente na uniformização do tamanho desses agregados e tal uniformidade se manteve estável por mais de cinco subcultivos em 100% das culturas. Uma fração rica em saponinas foi obtida a partir dos agregados celulares, correspondendo a 1,46% da massa seca. A análise por RMN sugeriu a presença das saponinas Hovenosideo G e dos hovenidulciosideos A2 e B2. O uso de elicitores em cultura de calos mostrou-se adequado à produção de metabólitos secundários, sem alterações morfológicas nos mesmos. A elicitação alterou o perfil cromatográfico analisado por HPLC. Na elicitação com 5,0 mg.L-1 de extrato de levedura foi verificado um aumento de quase três vezes (12,280 3,396 equivalentes de quercetina/mg de extrato) na síntese de flavonoides. Finalmente, os estudos de ação antitumoral in vitro demonstraram citotoxicidade dos extratos de calos não elicitados de H. dulcis sobre linhagem de leucemia mieloide crônica (IC50 de 74,05 μg.mL-1.) e inibição do crescimento de tais células (K562), sugerindo o potencial antineoplásico para um produto biotecnológio (calo) desta espécie.

Um estudo de caso envolvendo a retroanálise de tensões residuais de cravação.

As tensões residuais de cravação ocorrem quando a estaca e o solo não são totalmente descomprimidos após dissipação completa da energia transferida à estaca durante sua instalação. Este assunto tem merecido a atenção de vários pesquisadores, uma vez que a presença das tensões residuais na estaca causa uma alteração no seu comportamento quando carregada. A presente dissertação reuniu estudos anteriores que estabeleceram os principais fatores que influenciam as tensões residuais de cravação. Procurou-se executar uma série de simulações objetivando melhor entender o desenvolvimento das tensões residuais com a porcentagem de resistência de ponta, em continuidade aos estudos de Costa (1994). Costa (1994) observou que a porcentagem de carga na ponta é um fator relevante na avaliação das tensões residuais de cravação. Uma análise paramétrica efetuada confirmou estudos anteriores, verificando que a razão entre a carga residual na ponta em relação à capacidade de carga global cresce, a medida que a porcentagem de carga na ponta aumenta, chegando a um valor máximo para, em seguida, diminuir. Este comportamento é similar ao da curva de compactação do solo, quando, ao aumentar a umidade, o peso específico seco aumenta, até um valor máximo, correspondente à umidade ótima, para em seguida reduzir. Verificouse, ainda, semelhantemente ao aumento da energia de compactação, que vai transladando a curva para cima e para à esquerda do gráfico umidade x peso específico seco, que o aumento do comprimento da estaca apresenta um comportamento similar. O aumento do comprimento leva a curva para cima e para a esquerda. Finalmente, selecionou-se um caso de obra com condições propícias para o desenvolvimento de altas tensões residuais de cravação. Trata-se de um caso de estacas metálicas longas, embutidas em solo residual jovem, de elevada resistência, bem documentado com extensa instrumentação, por ocasião da instalação. A retroanálise de cinco estacas do banco de dados mostrou que as cargas residuais previstas, num programa de simulação de cravação, se aproximaram muito dos valores experimentais. Foi observado também que a profundidade do ponto neutro previsto e medido apresentou uma excelente concordância. O resultado mais relevante desta pesquisa foi quando os valores de porcentagem de ponta foram introduzidos no eixo das abscissas e os valores da razão entre a carga residual na ponta e a capacidade de carga global, no eixo das ordenadas, e se observou o aspecto da curva semelhante à de compactação. A dissertação ilustra ser possível e simples a previsão das tensões residuais de cravação através de uma análise pela equação da onda.

Avaliação de um biorreator piloto de lodo ativado no controle de odores.

Apesar da importância das estações de tratamento de efluentes industriais (ETEIs) na conservação dos ecossistemas, estas podem ser geradoras de gases com maus odores contendo compostos orgânicos voláteis - COVs. Os maus odores têm sido motivos de protestos e reclamações por parte da população circunvizinha às fontes emissoras. Em virtude da conscientização ambiental, e dos impactos sobre a saúde do homem, o objetivo geral deste trabalho é avaliar a eficiência de um biorreator aeróbio piloto no controle de gases odoríferos emitidos em estação de tratamento de efluentes de indústrias de alimentos. Foi desenvolvido um sistema de difusão de ar odorífero em um reator aeróbio de lodo ativado, hermeticamente fechado, operado no regime de batelada sequencial, durante os dias 14, 21, 23, e 30 do mês de julho.Foram realizadas análises dos parâmetros físico-químicos do lodo ativado utilizado no reator aeróbio piloto, como determinação dos sólidos, DBO5, DQO, OD, pH, temperatura e IVL. A atividade da biomassa do lodo ativado foi avaliada por meio do teste de Respirometria. A eficiência do reator quanto à redução da DQO dos gases foi analisada por meio da absorção dos gases em solução de dicromato de potássio. Para avaliação da mensuração dos odores utilizou-se o método de cromatografia gasosa e espectrometria de massa, quantificando amostras de gases odorantes adsorvidas em tubos de carvão ativado, na entrada e na saída dobiorreator. Os resultados obtidos confirmaram o potencial do sistema de difusão em lodos ativados para o tratamento de gases odoríferos em ETEI, com eficiência de remoção dos COV`s variando de 97,3% a 98,9%.

Implementação paralelizada de métodos de resolução de sistemas algébricos na simulação de reservatórios de gás

Desde a década de 1960, devido à pertinência para a indústria petrolífera, a simulação numérica de reservatórios de petróleo tornou-se uma ferramenta usual e uma intensa área de pesquisa. O principal objetivo da modelagem computacional e do uso de métodos numéricos, para a simulação de reservatórios de petróleo, é o de possibilitar um melhor gerenciamento do campo produtor, de maneira que haja uma maximização na recuperação de hidrocarbonetos. Este trabalho tem como objetivo principal paralelizar, empregando a interface de programação de aplicativo OpenMP (Open Multi-Processing), o método numérico utilizado na resolução do sistema algébrico resultante da discretização da equação que descreve o escoamento monofásico em um reservatório de gás, em termos da variável pressão. O conjunto de equações governantes é formado pela equação da continuidade, por uma expressão para o balanço da quantidade de movimento e por uma equação de estado. A Equação da Difusividade Hidráulica (EDH), para a variável pressão, é obtida a partir deste conjunto de equações fundamentais, sendo então discretizada pela utilização do Método de Diferenças Finitas, com a escolha por uma formulação implícita. Diferentes testes numéricos são realizados a fim de estudar a eficiência computacional das versões paralelizadas dos métodos iterativos de Jacobi, Gauss-Seidel, Sobre-relaxação Sucessiva, Gradientes Conjugados (CG), Gradiente Biconjugado (BiCG) e Gradiente Biconjugado Estabilizado (BiCGStab), visando a uma futura aplicação dos mesmos na simulação de reservatórios de gás. Ressalta-se que a presença de heterogeneidades na rocha reservatório e/ou às não-linearidades presentes na EDH para o escoamento de gás aumentam a necessidade de métodos eficientes do ponto de vista de custo computacional, como é o caso de estratégias usando OpenMP.

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.