Levantamento de plantas daninhas aquáticas no reservatório de Alagados, Ponta Grossa-PR

O monitoramento de macrófitas da represa Alagados um reservatório artificial de água que abastece parcialmente a população do município de Ponta Grossa-PR se faz necessário, uma vez que resíduos urbanos, bem como dejetos de áreas rurais, são descartados direta ou indiretamente nesse ambiente. Os objetivos deste trabalho foram realizar o inventário de plantas aquáticas coletadas entre agosto de 2007 e julho de 2008 no entorno do lago, como medida de monitoramento da área; indicar a frequência das espécies e apresentar um diagnóstico sobre o crescimento populacional destas; e identificar áreas onde a colonização por macrófitas é intensa, podendo tornar-se problema ambiental. Para realizar esse inventário, foram efetuadas 10 coletas em 48 pontos de amostragem ao longo das margens do lago da represa. Exsicatas foram depositadas no Herbário da Universidade Estadual de Ponta Grossa (HUPG). Foram amostradas 54 espécies, distribuídas em 16 famílias de angiospermas, 8 monocotiledôneas e 8 dicotiledôneas, e outras duas famílias de pteridófitas.

Monitoramento de fitoplâncton e microcistina no reservatório da UHE Americana

Este trabalho foi realizado na UHE Americana, pertencente à Companhia Paulista de Força e Luz, e faz parte de um projeto de pesquisa e desenvolvimento realizado em conjunto com a Faculdade de Ciências Agronômicas - UNESP, de Botucatu. As amostragens de água foram realizadas nos meses de fevereiro, abril, junho e outubro de 2004. As características analisadas foram: temperatura da água, pH, oxigênio dissolvido, condutividade, nitrogênio total, nitrito, nitrato, amônia, fósforo total, fosfato, fosfato inorgânico, juntamente com análise qualitativa e quantitativa da comunidade fitoplanctônica e a toxicidade. O reservatório apresentou valores elevados de fósforo total, variando de 18 a 509 µg L-1; fosfato, de 4 a 463 µg L-1; nitrogênio total, de 0,99 a 17,25 mg L-1; e nitrato, de 0,26 a 15,29 mg L-1. Para a comunidade fitoplanctônica foram encontrados 103 táxons em todo o período amostrado; a maior riqueza foi encontrada no ponto P06, e a maior pobreza de táxons, nos pontos localizadas no corpo central do reservatório (P02, P03, P04 e P05). A maior concentração de cianofícea ocorreu em abril de 2004: 5.375.175 ind. L-1. As espécies que apresentaram as maiores densidades foram Microcystis aeruginosa, Anabaena spiroides, Microcystis sp. e Pseudoanabaena mucicola; a maior densidade foi apresentada por Anabaena spiroides, com 4.178.084 ind. L-1. Nos meses de junho e outubro a classe Cryptophyceae teve uma grande contribuição para a densidade total. Apesar da grande densidade de cianobactérias, os valores de toxicidade ficaram abaixo do limite permitido pela Portaria nº 1.469.

Estudo da movimentação de plantas aquáticas imersas presentes no reservatório da UHE Eng. Souza Dias - Jupiá

O objetivo deste estudo foi desenvolver um método para determinação de pontos de reprodução e de suas respectivas importâncias na disseminação das espécies E. densa, E. najas e C. demersum no reservatório de Jupiá. Foram monitorados dez locais nos quais se encontrou registro de altas infestações dessas espécies. Dois desses locais situam-se no rio Paraná, nas lagoas denominadas de Ferradura e Pernilongo; os outros oito sítios estão localizados no rio Tietê, nas lagoas marginais Testemunha, Barrenta, Vírgula e Flórida, e no leito do rio nos pontos Acima da Ponte dos Barrageiros, Abaixo da Ponte, Baía ao Lado da Ponte e em frente da Praia de Itapura. em cada um desses sítios de reprodução foram liberados e monitorados dez blocos de plantas aquáticas imersas por mês. Os blocos possuíam volume de 0,14 cm³ e receberam uma etiqueta plástica de identificação externa e uma bóia de cor laranja, com o objetivo de serem facilmente localizadas ou visualizadas à distância. Após a soltura, cada bloco foi georreferenciado por meio de um aparelho GPS, sendo o seu deslocamento avaliado a cada 15 dias. Os sítios Lagoa Vírgula, Lagoa Testemunha e Lagoa Barrenta destacaram-se quanto ao número de blocos (30, 20 e 19 blocos, respectivamente), que percorreram distâncias superiores a 500 m e, conseqüentemente, saíram de seus respectivos locais de reprodução. A grande maioria dos blocos permaneceu dentro da Lagoa Flórida; somente 12 blocos saíram e alcançaram o leito do rio Tietê. O sítio denominado Leito Acima da Ponte foi o que mais se destacou dentre os três localizados no leito do rio Tietê, pois forneceu 18 blocos de plantas, enquanto os sítios Abaixo da Ponte e Baía ao Lado da Ponte forneceram 15 e 14 blocos, respectivamente. O sítio Praia de Itapura foi o que apresentou menor importância em relação ao fornecimento de plantas imersas no rio Tietê, sendo localizadas apenas 11 bóias com deslocamento superior a 500 m. Os sítios de reprodução localizados no rio Paraná também contribuíram para o fornecimento de plantas imersas ao reservatório de Jupiá; nove blocos saíram da Lagoa Pernilongo e apresentaram deslocamentos superiores a 500 m. Entretanto, os sítios Lagoa Barrenta e Lagoa Ferradura foram aqueles que mais contribuíram com a chegada de plantas na tomada de água da UHE Eng. Souza Dias, uma vez que foram recolhidos, respectivamente, seis e quatro blocos de plantas provenientes desses locais de dispersão. A metodologia utilizada mostrou-se eficaz quanto à avaliação da movimentação de plantas aquáticas imersas dentro do reservatório de Jupiá.

Levantamento da infestação de plantas aquáticas em Porto Primavera antes do enchimento final do reservatório

O objetivo deste trabalho foi identificar as plantas aquáticas e os níveis de infestação de cada espécie em Porto Primavera antes do enchimento final do reservatório. Foram avaliados todos os focos de vegetação aquática no reservatório (72 pontos), sendo os pontos demarcados com um aparelho de GPS. As espécies foram identificadas e estimouse visualmente (tamanho da área) a distribuição proporcional das plantas no foco de infestação. Após a identificação, foram encontradas 18 espécies de plantas aquáticas vegetando no reservatório, das quais foram determinadas a frequência de espécie de planta aquática e a distribuição dentro do sistema de geração de energia. As espécies encontradas no reservatório foram: Eichhornia crassipes, Eichhornia azurea, Pistia stratiotes, Paspalum repens, Cyperus brevifolius, Paspalum conspersum, Echinochloa polystachya, Egeria densa, Egeria najas, Polygonum hidropiperoides, Polygonum lapathifolium, Alternanthera philoxeroides, Eleocharis sellowiana, Nymphaea ampla, Pontederia cordata, Salvinia auriculata, Salvinia rotundifolia e Typha angustifolia. As maiores frequências relativas foram observadas em: E. azurea (36,11%), E. crassipes (16,67%), P. stratiotes (13,89%), S. auriculata (13,89%), C. brevifolius (11,11%) e P. lapathifolium (6,94%).

Seleção de sementes de milho durante a simulação da semeadura com disco perfurado horizontal

Teve-se o objetivo de verificar a possível seleção de sementes de milho - arredondadas e achatadas - no processo de semeadura com disco perfurado horizontal, bem como as conseqüências na distribuição longitudinal. Os registros dos dados deram-se pela avaliação visual de duas pessoas, com contador manual, no final da esteira carpetada em movimento de 1,54 m s-1, com 4,60 m de perspectiva para avaliação. Assim, um sensor registrava o número de orifícios do disco horizontal que devia passar pelo tubo condutor, e os espaçamentos falhos e múltiplos eram computados por pessoas diferentes. Os espaçamentos aceitáveis foram obtidos pela subtração dos espaçamentos falhos e múltiplos do total contabilizado. A medição do comprimento e espessura das sementes foi realizada com auxílio de paquímetro. Não houve alterações significativas nas variáveis estudadas, até a passagem de 55 mil orifícios do disco horizontal pelo tubo condutor; portanto, não houve seleção de sementes de milho durante a simulação do processo de semeadura.

Modelagem e simulação da deposição de parafinas em escoamento turbulento

The flow assurance has become one of the topics of greatest interest in the oil industry, mainly due to production and transportation of oil in regions with extreme temperature and pressure. In these operations the wax deposition is a commonly problem in flow of paraffinic oils, causing the rising costs of the process, due to increased energy cost of pumping, decreased production, increased pressure on the line and risk of blockage of the pipeline. In order to describe the behavior of the wax deposition phenomena in turbulent flow of paraffinic oils, under different operations conditions, in this work we developed a simulator with easy interface. For that we divided de work in four steps: (i) properties estimation (physical, thermals, of transport and thermodynamics) of n-alkanes and paraffinic mixtures by using correlations; (ii) obtainment of the solubility curve and determination the wax appearance temperature, by calculating the solid-liquid equilibrium of parafinnic systems; (iii) modelling wax deposition process, comprising momentum, mass and heat transfer; (iv) development of graphic interface in MATLAB® environment for to allow the understanding of simulation in different flow conditions as well as understand the matter of the variables (inlet temperature, external temperature, wax appearance temperature, oil composition, and time) on the behavior of the deposition process. The results showed that the simulator developed, called DepoSim, is able to calculate the profile of temperature, thickness of the deposit, and the amount of wax deposited in a simple and fast way, and also with consistent results and applicable to the operation

Desenvolvimento de sistemas de pastas leves para aplicação em cimentação de poços petrolíferos com baixo gradiente de fratura

All around the world, naturally occurring hydrocarbon deposits, consisting of oil and gas contained within rocks called reservoir rocks , generally sandstone or carbonate exists. These deposits are in varying conditions of pressure and depth from a few hundred to several thousand meters. In general, shallow reservoirs have greater tendency to fracture, since they have low fracture gradient, ie fractures are formed even with relatively low hydrostatic columns of fluid. These low fracture gradient areas are particularly common in onshore areas, like the Rio Grande do Norte basin. During a well drilling, one of the most favorable phases for the occurrence of fractures is during cementing, since the cement slurry used can have greater densities than the maximum allowed by the rock structure. Furthermore, in areas which are already naturally fractured, the use of regular cement slurries causes fluid loss into the formation, which may give rise to failures cementations and formation damages. Commercially, there are alternatives to the development of lightweight cement slurries, but these fail either because of their enormous cost, or because the cement properties were not good enough for most general applications, being restricted to each transaction for which the cement paste was made, or both reasons. In this work a statistical design was made to determine the influence of three variables, defined as the calcium chloride concentration, vermiculite concentration and nanosilica concentration in the various properties of the cement. The use of vermiculite, a low density ore present in large amounts in northeastern Brazil, as extensor for cementing slurries, enabled the production of stable cements, with high water/cement ratio, excellent rheological properties and low densities, which were set at 12.5 lb / gal, despite the fact that lower densities could be achieved. It is also seen that the calcium chloride is very useful as gelling and thickening agent, and their use in combination with nanosilica has a great effect on gel strength of the cement. Hydrothermal Stability studies showed that the pastes were stable in these conditions, and mechanical resistance tests showed values of the order of up to 10 MPa

Estudo do equilíbrio líquido-vapor do sistema água+etanol+líquido iônico visando a separação do álcool anidro

Anhydrous ethanol is used in chemical, pharmaceutical and fuel industries. However, current processes for obtaining it involve high cost, high energy demand and use of toxic and pollutant solvents. This problem occurs due to the formation of an azeotropic mixture of ethanol + water, which does not allow the complete separation by conventional methods such as simple distillation. As an alternative to currently used processes, this study proposes the use of ionic liquids as solvents in extractive distillation. These are organic salts which are liquids at low temperatures (under 373,15 K). They exhibit characteristics such as low volatility (almost zero/ low vapor ), thermal stability and low corrosiveness, which make them interesting for applications such as catalysts and as entrainers. In this work, experimental data for the vapor pressure of pure ethanol and water in the pressure range of 20 to 101 kPa were obtained as well as for vapor-liquid equilibrium (VLE) of the system ethanol + water at atmospheric pressure; and equilibrium data of ethanol + water + 2-HDEAA (2- hydroxydiethanolamine acetate) at strategic points in the diagram. The device used for these experiments was the Fischer ebulliometer, together with density measurements to determine phase compositions. The experimental data were consistent with literature data and presented thermodynamic consistency, thus the methodology was properly validated. The results were favorable, with the increase of ethanol concentration in the vapor phase, but the increase was not shown to be pronounced. The predictive model COSMO-SAC (COnductor-like Screening MOdels Segment Activity Coefficient) proposed by Lin & Sandler (2002) was studied for calculations to predict vapor-liquid equilibrium of systems ethanol + water + ionic liquids at atmospheric pressure. This is an alternative for predicting phase equilibrium, especially for substances of recent interest, such as ionic liquids. This is so because no experimental data nor any parameters of functional groups (as in the UNIFAC method) are needed

Modelagem e simulação de um separador bifásico com inversão de fases induzida

The present work has the main goal to study the modeling and simulation of a biphasic separator with induced phase inversion, the MDIF, with the utilization of the finite differences method for the resolution of the partial differencial equations which describe the transport of contaminant s mass fraction inside the equipment s settling chamber. With this aim, was developed the deterministic differential model AMADDA, wich was admensionalizated and then semidiscretizated with the method of lines. The integration of the resultant system of ordinary differential equations was realized by means of a modified algorithm of the Adam-Bashfort- Moulton method, and the sthocastic optimization routine of Basin-Hopping was used in the model s parameter estimation procedure . With the aim to establish a comparative referential for the results obtained with the model AMADDA, were used experimental data presented in previous works of the MDIF s research group. The experimental data and those obtained with the model was assessed regarding its normality by means of the Shapiro-Wilk s test, and validated against the experimental results with the Student s t test and the Kruskal-Wallis s test, depending on the result. The results showed satisfactory performance of the model AMADDA in the evaluation of the MDIF s separation efficiency, being possible to determinate that at 1% significance level the calculated results are equivalent to those determinated experimentally in the reference works

Modelagem e simulação da solubilização da parafina em tubos circulares e em regime laminar

O óleo produzido nos novos campos de petróleo está cada vez mais parafínico e viscoso, com isso, à medida que o óleo é escoado, parafinas são depositadas sobre as paredes internas do tubo, e ao longo do tempo, tendem a reduzir drasticamente a área transversal ao escoamento. Visando estudar o processo de solubilização da parafina em dutos, esse trabalho objetiva desenvolver modelos matemáticos que represente o processo, com base nos fenômenos envolvidos no mesmo tais como transferência de massa, transferência de energia e equilíbrio sólido-líquido, implementando-os em um ambiente de desenvolvimento VBA (Visual Basic) for Excel ®. O presente trabalho foi realizado em quatro etapas: i) modelagem dos fenômenos de transferência de calor e massa, ii) modelagem da rotina dos coeficientes de atividade através do modelo UNIFAC e modelagem do sistema de equilíbrio sólido-líquido; iii) modelagem matemática do processo de solubilização e cálculo da espessura da parafina ao longo do tempo; iv) implementação dos modelos em um ambiente de desenvolvimento VBA for Excel® e criação de um simulador com uma interface gráfica, para simular o processo de solubilização da parafina depositada em dutos e sua otimização. O simulador conseguiu produzir soluções bastante adequadas, mantendo continuidade das equações diferenciáveis do balanço de energia e de massa, com uma interpretação física viável, sem a presença de dissipação de oscilações nos perfis de temperatura e massa. Além disso, esse simulador visa permitir a simulação nas diversas condições de escoamento, bem como compreender a importância das variáveis (vazão, temperatura de entrada, temperatura externa, cadeia carbônica do solvente). Através dos resultados foram possíveis verificar os perfis de temperatura, fração molar e o de solubilização

Recuperação e purificação de quitosanases usando adsorção em leito expandido com streamline DEAE com modelagem e simulação usando redes neurais

Expanded Bed Adsorption (EBA) is an integrative process that combines concepts of chromatography and fluidization of solids. The many parameters involved and their synergistic effects complicate the optimization of the process. Fortunately, some mathematical tools have been developed in order to guide the investigation of the EBA system. In this work the application of experimental design, phenomenological modeling and artificial neural networks (ANN) in understanding chitosanases adsorption on ion exchange resin Streamline® DEAE have been investigated. The strain Paenibacillus ehimensis NRRL B-23118 was used for chitosanase production. EBA experiments were carried out using a column of 2.6 cm inner diameter with 30.0 cm in height that was coupled to a peristaltic pump. At the bottom of the column there was a distributor of glass beads having a height of 3.0 cm. Assays for residence time distribution (RTD) revelead a high degree of mixing, however, the Richardson-Zaki coefficients showed that the column was on the threshold of stability. Isotherm models fitted the adsorption equilibrium data in the presence of lyotropic salts. The results of experiment design indicated that the ionic strength and superficial velocity are important to the recovery and purity of chitosanases. The molecular mass of the two chitosanases were approximately 23 kDa and 52 kDa as estimated by SDS-PAGE. The phenomenological modeling was aimed to describe the operations in batch and column chromatography. The simulations were performed in Microsoft Visual Studio. The kinetic rate constant model set to kinetic curves efficiently under conditions of initial enzyme activity 0.232, 0.142 e 0.079 UA/mL. The simulated breakthrough curves showed some differences with experimental data, especially regarding the slope. Sensitivity tests of the model on the surface velocity, axial dispersion and initial concentration showed agreement with the literature. The neural network was constructed in MATLAB and Neural Network Toolbox. The cross-validation was used to improve the ability of generalization. The parameters of ANN were improved to obtain the settings 6-6 (enzyme activity) and 9-6 (total protein), as well as tansig transfer function and Levenberg-Marquardt training algorithm. The neural Carlos Eduardo de Araújo Padilha dezembro/2013 9 networks simulations, including all the steps of cycle, showed good agreement with experimental data, with a correlation coefficient of approximately 0.974. The effects of input variables on profiles of the stages of loading, washing and elution were consistent with the literature

Estudo de sistemas microemulsionados utilizando água do mar na recuperação avançada de petróleo

Petroleum exists in the nature in certain underground formations where it is adsorbed into the rocks pores. For the conventional recovery methods usually only 30% of the oil is extracted and this can be credited, basically, to three aspects: high viscosity of the oil, geology of the formation and high interfacial tensions between the reservoir s fluids. The enhanced recovery methods use the injection of a fluid or fluids mixture in a reservoir to act in points where the conventional process didn't reach the recovery rates. Microemulsion flooding, considered an enhanced method, has the purpose to desorb the oil from the rock formation and to attain an efficient displacement of the oil emulsion. With this in mind, this work was accomplished with two main objectives: the study of the parameters effect that influence a microemulsified system (surfactant and cosurfactant types, C/S rate and salinity) and the evaluation of displacement efficiency with the microemulsions that showed stability in the rich aqueous area. For the analyzed parameters it was chose the microemulsions composition used in the recovery stage: 25% water, 5% kerosene, 46.7% of butanol as cosurfactant and 23.3% of BC or SCO cosurfactant. The core plugs of Assu and Botucatu sandstones were appraised in porosity and permeability tests and then submitted to the steps of saturation with seawater and oil, conventional recovery with water and enhanced recovery with the selected microemulsions. The Botucatu sandstone presented better recovery parameters, and the microemulsion composed with BS surfactant had larger recovery efficiency (26.88%)

Determinação de dados de equilíbrio de fases para sistemas aquosos com eletrólitos

In the present work are established initially the fundamental relationships of thermodynamics that govern the equilibrium between phases, the models used for the description of the behavior non ideal of the liquid and vapor phases in conditions of low pressures. This work seeks the determination of vapor-liquid equilibrium (VLE) data for a series of multicomponents mixtures of saturated aliphatic hydrocarbons, prepared synthetically starting from substances with analytical degree and the development of a new dynamic cell with circulation of the vapor phase. The apparatus and experimental procedures developed are described and applied for the determination of VLE data. VLE isobarics data were obtained through a Fischer s ebulliometer of circulation of both phases, for the systems pentane + dodecane, heptane + dodecane and decane + dodecane. Using the two new dynamic cells especially projected, of easy operation and low cost, with circulation of the vapor phase, data for the systems heptane + decane + dodecane, acetone + water, tween 20 + dodecane, phenol + water and distillation curves of a gasoline without addictive were measured. Compositions of the equilibrium phases were found by densimetry, chromatography, and total organic carbon analyzer. Calibration curves of density versus composition were prepared from synthetic mixtures and the behavior excess volumes were evaluated. The VLE data obtained experimentally for the hydrocarbon and aqueous systems were submitted to the test of thermodynamic consistency, as well as the obtained from the literature data for another binary systems, mainly in the bank DDB (Dortmund Data Bank), where the Gibbs-Duhem equation is used obtaining a satisfactory data base. The results of the thermodynamic consistency tests for the binary and ternary systems were evaluated in terms of deviations for applications such as model development. Later, those groups of data (tested and approved) were used in the KijPoly program for the determination of the binary kij parameters of the cubic equations of state original Peng-Robinson and with the expanded alpha function. These obtained parameters can be applied for simulation of the reservoirs petroleum conditions and of the several distillation processes found in the petrochemistry industry, through simulators. The two designed dynamic cells used equipments of national technology for the determination of VLE data were well succeed, demonstrating efficiency and low cost. Multicomponents systems, mixtures of components of different molecular weights and also diluted solutions may be studied in these developed VLE cells

Estudo da secagem do pedúnculo do caju em sistemas convencional e solar: modelagem e simulação do processo

This work aims to study the drying of cashew-nut pulp with different lay-out of dryers using conventional and solar energy. It concerns with the use of exceeding of the regional raw material and the suitable knowledge for the applicability of the drying systems as pathway for food conservation. Besides, it used renewable sources as solar energy to dry these agroindustrial products. Runs were carried out using a conventional tray-dryer with temperature, air velocity control and cashew slice thickness of 55°C, 65°C, 75°C; 3.0; 4.5, 6.0 m s-1; 1.0; 1.5 and 2.0 cm, respectively, in order to compare the studied systems. To evaluate the conventional tray-dryer, it was used a diffusional model of 2nd Fick´s law, where the drying curves were quite well fitted to an infinite flat plate design. For the drying runs where the room temperature had no control, it was developed a phenomenological-mathematical model for the solar dryer with indirect radiation under natural and forced convection based on material and energy balances of the system. Besides, it was carried out assays in the in natura as well as dehydrated, statistic analysis of the experimental drying data, sensorial analysis of the final dry product and a simplified economical analysis of the systems studied

Viabilidade técnica de uma rota não convencional para a produção de leite de cabra em pó em cooperativas do Rio Grande do Norte

Nowadays, Brazil has both the greatest goat herd and the greatest goat milk production of South America. The state of Rio Grande do Norte, located in northeast of Brazil, has an average year production of three thousand cubic meters of goat milk in natura. Part of this milk production is homemade and it comes from small farms, which unite in rural cooperatives created to encourage the production and implementation of industrial processes for preservation and processing of milk. Results presented by literature and obtained from preliminary essays in this thesis show that non conventional dryer of spouted bed with inert particles is able to produce powder milk from in natura milk (cattle or goat), with the same quality of spray dryer, however, operating at low cost. The method of drying in spouted bed consists of injecting milk emulsion on the bed of inert particles gushed by hot air. This emulsion covers the particles with a thin film, which dries and is reduced to powder during the circulation of inerts inside the bed. The powder is dragged by exhaustion air and separated in the cyclone. The friction among particles resulted from the particles circulation, encourages high taxes of shear in the thin film of emulsion, breaking the cohesive forces and making this process possible. Studying the drying process and the powder goat milk production in one unit of spouted bed with inert particles, seeing the development of a low cost technological route for powder milk production is the aim of this thesis. The powder milk produced by this route must attend the local demand of food industries which need an intermediate product to be used as a food ingredient (ice-cream, milk candy). In order to reach this aim, this thesis approaches the aspects related to physical, thermodynamics and physic-chemicals characteristics of goat milk, whose complete data are still inexistent in the literature. The properties of materials are of great importance to the project of any process which involves the operations of transportation of movement, heat and mass quantity, such as the dryers which operate in fluid dynamically active regime, like the spouted bed. It was obtained new data related to the goat milk properties in function of concentration of solids and temperature. It is also important to mention the study developed about the kinetic of solids retention in the bed of inert particles during the drying of goat milk. It was found more adequate processes conditions to the proposed technological route to be implemented in small and micro-industries, with simplifications in the system of milk injection as well as in the form of operation of the dryer. Important data were obtained for a posterior stage of this research which involves the v modeling, simulation, control and optimization of the process. The results obtained in this thesis, in relation to process performance as well as to the quality of produced powder milk validate the proposal of using the spouted bed dryer in the production of powder goat milk