Progress in numerical modeling of non-premixed combustion

Progress is made on the numerical modeling of both laminar and turbulent non-premixed flames. Instead of solving the transport equations for the numerous species involved in the combustion process, the present study proposes reduced-order combustion models based on local flame structures.

For laminar non-premixed flames, curvature and multi-dimensional diffusion effects are found critical for the accurate prediction of sooting tendencies. A new numerical model based on modified flamelet equations is proposed. Sooting tendencies are calculated numerically using the proposed model for a wide range of species. These first numerically-computed sooting tendencies are in good agreement with experimental data. To further quantify curvature and multi-dimensional effects, a general flamelet formulation is derived mathematically. A budget analysis of the general flamelet equations is performed on an axisymmetric laminar diffusion flame. A new chemistry tabulation method based on the general flamelet formulation is proposed. This new tabulation method is applied to the same flame and demonstrates significant improvement compared to previous techniques.

For turbulent non-premixed flames, a new model to account for chemistry-turbulence interactions is proposed. %It is found that these interactions are not important for radicals and small species, but substantial for aromatic species. The validity of various existing flamelet-based chemistry tabulation methods is examined, and a new linear relaxation model is proposed for aromatic species. The proposed relaxation model is validated against full chemistry calculations. To further quantify the importance of aromatic chemistry-turbulence interactions, Large-Eddy Simulations (LES) have been performed on a turbulent sooting jet flame. %The aforementioned relaxation model is used to provide closure for the chemical source terms of transported aromatic species. The effects of turbulent unsteadiness on soot are highlighted by comparing the LES results with a separate LES using fully-tabulated chemistry. It is shown that turbulent unsteady effects are of critical importance for the accurate prediction of not only the inception locations, but also the magnitude and fluctuations of soot.

Avaliação da função endotelial microvascular, da rigidez arterial e da resposta inflamatória sistêmica em pacientes submetidos a cirurgia cardíaca com circulação extracorpórea

A Cirurgia de Revascularização do Miocárdio, realizada com o auxílio da Circulação Extracorpórea, está associada a alterações importantes na microcirculação e na produção e circulação de citocinas e marcadores inflamatórios. No presente estudo, foram avaliados 23 pacientes com indicação de Revascularização do Miocárdio, no dia do procedimento e 7 e 28 dias após a cirurgia. A microcirculação cutânea, enquanto reflexo da microcirculação coronariana, foi estudada através da hiperemia térmica e/ ou reativa pós oclusiva e da iontoforese de substâncias vasoativas por mecanismos dependentes e independentes do endotélio. A rigidez arterial foi aferida através da análise da onda de pulso digital. Foi avaliado ainda o impacto da doença e do procedimento cirúrgico sobre a produção e circulação sérica de citocinas e marcadores inflamatórios, tais como: PCR-HS, nitrito/ nitrato, IL-6, Il-7, IL-8, IL-10, IFN-γ, TNF-α e G-CSF. Foi observada uma tendência à redução da vasodilatação da microcirculação cutânea após a administração de doses acumulativas de acetilcolina (endotélio dependente) através da iontoforese de 7 e 28 dias após o procedimento cirúrgico. A hiperemia térmica foi mais pronunciada na avaliação basal do que aos 7 e 28 dias. A hiperemia reativa pós oclusiva não demonstrou alterações 7 dias após o procedimento. Aos 28 dias, houve um aumento da condutância microvascular cutânea. Quando avaliada a vasodilatação endotélio-independente (nitroprussiato de sódio), observamos aumento do fluxo microvascular cutâneo diretamente proporcional à carga/ dose aplicada, sem diferenças nos valores obtidos no basal e 7 e 28 dias após o procedimento. A rigidez arterial não apresentou alterações. A análise dos fatores inflamatórios e das citocinas demonstrou aumento marcante da IL-6 e da IL-8 após 7 dias do procedimento cirúrgico, com retorno parcial aos níveis basais da IL-6 e total da IL-8 após 28 dias. O IFN-γ, TNF-α e G-CSF apenas apresentaram níveis detectáveis na avaliação basal e IL-7 e IL-10 não demonstraram alterações significativas nos tempos avaliados. A PCR-HS demonstrou níveis mais elevados após 7 dias e retorno parcial aos níveis basais após 28 dias. O nitrito/ nitrato, após 7 dias, apresentou leve queda em sua concentração plasmática. Concluímos que a pequena diferença entre o valores obtidos entre o basal e após 7 dias do procedimento cirúrgico com a iontoforese de acetilcolina resulta em minimização do impacto endotelial e um valor constante deste dado após 28 dias, sugere recomposição fisiológica completa. Este resultado foi semelhante com a análise da hiperemia térmica e reativa pós oclusiva. As interleucinas IL-6 e IL-8, bem como a PCR-HS apresentaram comportamento correlacionável, refletindo a cinética inflamatória. A rigidez arterial não demonstrou alterações.

Condições mecânicas do sistema respiratório em adultos com fibrose cística: estudo pela espirometria, pletismografia e técnica de oscilações forçadas

A fibrose cística (FC) é a doença autossômica recessiva mais comum na população branca que leva à redução na expectativa de vida. A doença pulmonar é a maior causa de morbidade e mortalidade. A relevância do presente estudo se dá diante de alguns fatores: aumento drástico da sobrevida média nos últimos 60 anos na FC, a fisiopatologia pulmonar não é bem compreendida, ausência de estudos reportados na literatura, até o momento, utilizando a técnica de oscilações forçadas (TOF) exclusivamente em adultos com FC. Assim sendo os objetivos deste estudo são: analisar as alterações da mecânica respiratória em adultos com FC através da espirometria, pletismografia e TOF; correlacionar os resultados da TOF aos espirométricos e pletismográficos e avaliar a sensibilidade e especificidade da TOF nestes indivíduos. É um estudo de corte transversal descritivo, no qual foram analisados dois grupos de indivíduos: controle (n=23) e FC (n=27). Os resultados foram expressos através média desvio-padrão. As técnicas funcionais respiratórias foram realizadas na seguinte sequência: TOF, espirometria, pletismografia. Na pletismografia foram avaliados os parâmetros: CPT (capacidade pulmonar total), CRF (capacidade residual funcional) e VR (volume residual), CRF/CPT e VR/CPT, resistência (Rva) e condutância específica das vias aéreas (SGva). Na espirometria: volume expiratório forçado no primeiro segundo (VEF1), capacidade vital forçada (CVF), fluxo expiratório entre 25% e 75% (FEF25%-75%) da CVF (FEF25%-75%) e razões VEF1/CVF (%) e FEF/CVF (%). Na TOF: propriedades resistivas do sistema respiratório- R0 (resistência no intercepto), Rm (resistência média) e S (inclinação da reta de resistência) e propriedades reativas: Cdin,sr (complacência dinâmica do sistema respiratório), Xm (reatância média), frequência de ressonância (fr); e o módulo da impedância em 4 Hz (׀Zrs4Hz׀). Na espirometria o distúrbio ventilatório obstrutivo (DVO) com CVF reduzida foi predominante, com marcante redução do FEF25%-75% no grupo FC (p<0,0001) em relação ao controle. Na pletismografia: destacou-se a elevação de VR, na presença de CPT normal e elevação da Rva e redução da SGva no grupo FC. Alterações da TOF ocorridas no grupo FC em relação ao controle: aumento de R0 e Rm (p<0,0001) e fr (p<0,0002), relacionados à obstrução das vias aéreas; redução de S (p<0,0006), Xm (p<0,0001) associadas à não-homogeneidade do sistema respiratório e Cdin,sr (p<0,0001), relacionada à redução da complacência pulmonar; aumento do módulo da impedância em 4 Hz (׀Zrs4Hz׀) (representando a carga mecânica total do sistema respiratório) resultante da interação das demais alterações da TOF citadas. Os parâmetros da TOF apresentaram correlações muito boas com a espirometria e moderadas com a pletismografia. Rm foi o único parâmetro que não se relacionou com nenhuma destas técnicas. A sensibilidade e especificidade da TOF em adultos com FC apresentaram valores elevados, sobretudo nos parâmetros reativos, em especial, Xm (85,2% e 73,9% respectivamente e área sob a curva de 0,86).

Boundary-layer transition on a slender cone in hypervelocity flow with real gas effects

The laminar to turbulent transition process in boundary layer flows in thermochemical nonequilibrium at high enthalpy is measured and characterized. Experiments are performed in the T5 Hypervelocity Reflected Shock Tunnel at Caltech, using a 1 m length 5-degree half angle axisymmetric cone instrumented with 80 fast-response annular thermocouples, complemented by boundary layer stability computations using the STABL software suite. A new mixing tank is added to the shock tube fill apparatus for premixed freestream gas experiments, and a new cleaning procedure results in more consistent transition measurements. Transition location is nondimensionalized using a scaling with the boundary layer thickness, which is correlated with the acoustic properties of the boundary layer, and compared with parabolized stability equation (PSE) analysis. In these nondimensionalized terms, transition delay with increasing CO₂ concentration is observed: tests in 100% and 50% CO₂, by mass, transition up to 25% and 15% later, respectively, than air experiments. These results are consistent with previous work indicating that CO₂ molecules at elevated temperatures absorb acoustic instabilities in the MHz range, which is the expected frequency of the Mack second-mode instability at these conditions, and also consistent with predictions from PSE analysis. A strong unit Reynolds number effect is observed, which is believed to arise from tunnel noise. N_Tr for air from 5.4 to 13.2 is computed, substantially higher than previously reported for noisy facilities. Time- and spatially-resolved heat transfer traces are used to track the propagation of turbulent spots, and convection rates at 90%, 76%, and 63% of the boundary layer edge velocity, respectively, are observed for the leading edge, centroid, and trailing edge of the spots. A model constructed with these spot propagation parameters is used to infer spot generation rates from measured transition onset to completion distance. Finally, a novel method to control transition location with boundary layer gas injection is investigated. An appropriate porous-metal injector section for the cone is designed and fabricated, and the efficacy of injected CO₂ for delaying transition is gauged at various mass flow rates, and compared with both no injection and chemically inert argon injection cases. While CO₂ injection seems to delay transition, and argon injection seems to promote it, the experimental results are inconclusive and matching computations do not predict a reduction in N factor from any CO₂ injection condition computed.

The influence of freestream turbulence, freestream nuclei populations and a drag-reducing polymer on cavitation inception on two axisymmetric bodies

The influence upon the basic viscous flow about two axisymmetric bodies of (i) freestream turbulence level and (ii) the injection of small amounts of a drag-reducing polymer (Polyox WSR 301) into the test model boundary layer was investigated by the schlieren flow visualization technique. The changes in the type and occurrence of cavitation inception caused by the subsequent modifications in the viscous flow were studied. A nuclei counter using the holographic technique was built to monitor freestream nuclei populations and a few preliminary tests investigating the consequences of different populations on cavitation inception were carried out.

Both test models were observed to have a laminar separation over their respective test Reynolds number ranges. The separation on one test model was found to be insensitive to freestream turbulence levels of up to 3.75 percent. The second model was found to be very susceptible having its critical velocity reduced from 30 feet per second at a 0.04 percent turbulence level to 10 feet per second at a 3.75 percent turbulence level. Cavitation tests on both models at the lowest turbulence level showed the value of the incipient cavitation number and the type of cavitation were controlled by the presence of the laminar separation. Cavitation tests on the second model at 0.65 percent turbulence level showed no change in the inception index, but the appearance of the developed cavitation was altered.

The presence of Polyox in the boundary layer resulted in a cavitation suppression comparable to that found by other investigators. The elimination of the normally occurring laminar separation on these bodies by a polymer-induced instability in the laminar boundary layer was found to be responsible for the suppression of inception.

Freestream nuclei populations at test conditions were measured and it was found that if there were many freestream gas bubbles the normally present laminar separation was elminated and travelling bubble type cavitation occurred - the value of the inception index then depended upon the nuclei population. In cases where the laminar separation was present it was found that the value of the inception index was insensitive to the free stream nuclei populations.

Remoção de metais pesados de efluentes sintéticos aquosos por nanofiltração

O objetivo deste trabalho foi estudar a remoção de metais pesados de efluentes sintéticos aquosos através de nanofiltração, com utilização de membrana de poliamida. A membrana utilizada constituída de poliamida apresentou uma boa permeabilidade hidráulica, constatada pela linearidade de fluxo com a variação de pressão de trabalho. Os metais estudados foram cobre, cádmio, zinco, chumbo e níquel com especiação de nitrato e sulfato. O desempenho do sistema foi verificado pela variação dos parâmetros operacionais (pressão e fluxo), associados também com a variação da salinidade e da concentração do metal. Em uma etapa inicial estudou-se a remoção de metais (nitrato e sulfato) isoladamente em uma mesma concentração molar, na pressão de 10 bar, comparando-se as rejeições com e sem agente complexante (EDTA). Em uma segunda etapa, estudou-se o desempenho da membrana, com misturas dos metais (nitrato) e mistura de metais (sulfato), também com e sem agente complexante, os resultados mostraram que a presença de agente complexante melhorou a rejeição dos metais. Nestas misturas estudaram-se os efeitos das concentrações (0,001; 0,0005 e 0,0001mol/L) e das pressões (5; 7; 10 e 12,5 bar) no sistema de remoção, constatou-se que o aumento da concentração e diminuição da pressão afeta a remoção. Em uma terceira etapa, estudou-se a influência da salinidade pela mistura de dois metais de mesma especiação (sulfato), comparando a remoção na pressão de 10 bar. Os resultados atestaram uma rejeição maior que 94 % para todos metais, na pressão de 10 bar, indicando um excelente desempenho e fluxo adequado, daí viabilidade de processo para todas as concentrações testadas. A pertinência do teste se justifica pela adequação do permeado ao atendimento aos padrões ambientais de concentração de metais e confirmam a eficácia do sistema de nanofiltração na remoção de metais pesados

Characterization and modeling of premixed turbulent n-heptane flames in the thin reaction zone regime

n-heptane/air premixed turbulent flames in the high-Karlovitz portion of the thin reaction zone regime are characterized and modeled in this thesis using Direct Numerical Simulations (DNS) with detailed chemistry. In order to perform these simulations, a time-integration scheme that can efficiently handle the stiffness of the equations solved is developed first. A first simulation with unity Lewis number is considered in order to assess the effect of turbulence on the flame in the absence of differential diffusion. A second simulation with non-unity Lewis numbers is considered to study how turbulence affects differential diffusion. In the absence of differential diffusion, minimal departure from the 1D unstretched flame structure (species vs. temperature profiles) is observed. In the non-unity Lewis number case, the flame structure lies between that of 1D unstretched flames with "laminar" non-unity Lewis numbers and unity Lewis number. This is attributed to effective Lewis numbers resulting from intense turbulent mixing and a first model is proposed. The reaction zone is shown to be thin for both flames, yet large chemical source term fluctuations are observed. The fuel consumption rate is found to be only weakly correlated with stretch, although local extinctions in the non-unity Lewis number case are well correlated with high curvature. These results explain the apparent turbulent flame speeds. Other variables that better correlate with this fuel burning rate are identified through a coordinate transformation. It is shown that the unity Lewis number turbulent flames can be accurately described by a set of 1D (in progress variable space) flamelet equations parameterized by the dissipation rate of the progress variable. In the non-unity Lewis number flames, the flamelet equations suggest a dependence on a second parameter, the diffusion of the progress variable. A new tabulation approach is proposed for the simulation of such flames with these dimensionally-reduced manifolds.

MEMS for Diabetic Retinopathy

As the worldwide prevalence of diabetes mellitus continues to increase, diabetic retinopathy remains the leading cause of visual impairment and blindness in many developed countries. Between 32 to 40 percent of about 246 million people with diabetes develop diabetic retinopathy. Approximately 4.1 million American adults 40 years and older are affected by diabetic retinopathy. This glucose-induced microvascular disease progressively damages the tiny blood vessels that nourish the retina, the light-sensitive tissue at the back of the eye, leading to retinal ischemia (i.e., inadequate blood flow), retinal hypoxia (i.e., oxygen deprivation), and retinal nerve cell degeneration or death. It is a most serious sight-threatening complication of diabetes, resulting in significant irreversible vision loss, and even total blindness.

Unfortunately, although current treatments of diabetic retinopathy (i.e., laser therapy, vitrectomy surgery and anti-VEGF therapy) can reduce vision loss, they only slow down but cannot stop the degradation of the retina. Patients require repeated treatment to protect their sight. The current treatments also have significant drawbacks. Laser therapy is focused on preserving the macula, the area of the retina that is responsible for sharp, clear, central vision, by sacrificing the peripheral retina since there is only limited oxygen supply. Therefore, laser therapy results in a constricted peripheral visual field, reduced color vision, delayed dark adaptation, and weakened night vision. Vitrectomy surgery increases the risk of neovascular glaucoma, another devastating ocular disease, characterized by the proliferation of fibrovascular tissue in the anterior chamber angle. Anti-VEGF agents have potential adverse effects, and currently there is insufficient evidence to recommend their routine use.

In this work, for the first time, a paradigm shift in the treatment of diabetic retinopathy is proposed: providing localized, supplemental oxygen to the ischemic tissue via an implantable MEMS device. The retinal architecture (e.g., thickness, cell densities, layered structure, etc.) of the rabbit eye exposed to ischemic hypoxic injuries was well preserved after targeted oxygen delivery to the hypoxic tissue, showing that the use of an external source of oxygen could improve the retinal oxygenation and prevent the progression of the ischemic cascade.

The proposed MEMS device transports oxygen from an oxygen-rich space to the oxygen-deficient vitreous, the gel-like fluid that fills the inside of the eye, and then to the ischemic retina. This oxygen transport process is purely passive and completely driven by the gradient of oxygen partial pressure (pO₂). Two types of devices were designed. For the first type, the oxygen-rich space is underneath the conjunctiva, a membrane covering the sclera (white part of the eye), beneath the eyelids and highly permeable to oxygen in the atmosphere when the eye is open. Therefore, sub-conjunctival pO₂ is very high during the daytime. For the second type, the oxygen-rich space is inside the device since pure oxygen is needle-injected into the device on a regular basis.

To prevent too fast or too slow permeation of oxygen through the device that is made of parylene and silicone (two widely used biocompatible polymers in medical devices), the material properties of the hybrid parylene/silicone were investigated, including mechanical behaviors, permeation rates, and adhesive forces. Then the thicknesses of parylene and silicone became important design parameters that were fine-tuned to reach the optimal oxygen permeation rate.

The passive MEMS oxygen transporter devices were designed, built, and tested in both bench-top artificial eye models and in-vitro porcine cadaver eyes. The 3D unsteady saccade-induced laminar flow of water inside the eye model was modeled by computational fluid dynamics to study the convective transport of oxygen inside the eye induced by saccade (rapid eye movement). The saccade-enhanced transport effect was also demonstrated experimentally. Acute in-vivo animal experiments were performed in rabbits and dogs to verify the surgical procedure and the device functionality. Various hypotheses were confirmed both experimentally and computationally, suggesting that both the two types of devices are very promising to cure diabetic retinopathy. The chronic implantation of devices in ischemic dog eyes is still underway.

The proposed MEMS oxygen transporter devices can be also applied to treat other ocular and systemic diseases accompanied by retinal ischemia, such as central retinal artery occlusion, carotid artery disease, and some form of glaucoma.