Orientadores: Carlos Alberto Mariotoni, Caio Glauco Sanchez / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Civil / Made available in DSpace on 2018-07-25T01:21:05Z (GMT). No. of bitstreams: 1 Goulart_EduardoAntonio_M.pdf: 4150112 bytes, checksum: e2e80be604c26e414ff20623c3219692 (MD5) Previous issue date: 1999 / Resumo: o crescimento das grandes cidades, da industrialização e do consumo, tem forçado a civilização a conviver com um processo mais intenso de degradação do meio ambiente, ameaçando os recursos naturais e energéticos, além da maior produção de resíduos sólidos, líquidos e gasosos. Isto tem trazido preocupação com relação à sustentabilidade do processo de crescimento das economias mundiais e ao futuro dos recursos esgotáveis. Dessa forma, torna-se essencial a utilização de combustíveis obtidos a partir de fontes renováveis, como por exemplo os resíduos sólidos, líquidos e gasosos, gerados pela atividade industrial e pelas cidades. Uma das metas do planejamento energético de um país é minimizar a produção de descartes, e incentivar a reinserção dos resíduos e sub-produtos nos diversos processos produtivos. Além disso deve-se minimizar os impactos desses resíduos no meio ambiente. Como exemplo pode-se citar a produção de polímeros, que nas últimas décadas, com a evolução da pesquisa, tem provocado um crescimento na geração de resíduos, que devido à sua composição, permanecem intactos nos aterros por longos períodos de tempo sem que haja sensível degradação. Por outro lado, grande parte destes materiais possuem características bastante adequadas ao reaproveitamento ou reciclagem e, que podem apresentar razoável potencial energético. Entre os polímeros mais comuns estão as borrachas, principalmente as empregadas nos pneus automotivos, que tem gerado, nos países com maior frota, grandes problemas ambientais, ainda não completamente solucionados. Procurando contribuir para a solução destas questões, apresenta-se a proposta de um reator que utiliza pneus usados para a obtenção de óleo combustível através de sua gaseificação. A eficiência e a viabilidade da utilização deste processo foi avaliada através do seu balanço energético e da caracterização do óleo obtido como energético / Abstract: The growth of the big cities, industrialization and goods consumptions, have forced the civilization to live with more intensive environmental degradation process, threatening both the nature and energy resources, besides the bigger production of gas, liquid and solid wastes. That brings a greater worry related to the sustentability of the growing economy world process as well as to the future of the renewable resources. Therefore, it is important to utilize the fuels obtained from renewable sources, like the gas, liquid and solid wastes generated by both the industrial activities and the cities. Any energy planning aims to minimize the wastes production as well as to stimulate both the waste and the sub-products re-utilization in the several industrial processes. Besides that it is important to minimize their environmental impacts. As example it is mentioned the production of polymers, with have contributed to increase the wastes generation what stay untouched for many years at the landfills due to their compositions. Therefore most of those wastes can be re-utilized or recycled, and can have a reasonable energy potential. Among the more usual polymers are the tires rubber which generate, in the developed countries, great non-solved environmental problems. In order to contribute to the solution of those problems, it is presented in this work, a proposal of a reactor which uses scrap tires to get fuel oil using the gasification process. The efficiency and the viability of using this process have been analyzed by means of the energy balance and the characterization of the combustion oil as energy vector / Mestrado / Saneamento / Mestre em Engenharia Civil

Orientador: Milton Mori / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-27T17:58:36Z (GMT). No. of bitstreams: 1 Pelissari_DanielCicero_M.pdf: 3433960 bytes, checksum: 82461e75c8ef3d10fb44cd8d35de3eb9 (MD5) Previous issue date: 2015 / Resumo: A aplicação de fluidodinâmica computacional (CFD) em estudos de otimização e projeto de novos equipamentos de processos industriais vem aumentando significativamente, uma vez que apresenta custo reduzido e possibilidade de avaliar equipamentos complexos e de extremas condições de operação. Dentre os processos mais estudados via CFD está o processo de craqueamento catalítico fluidizado (FCC), onde as frações pesadas do petróleo de baixo valor são convertidas em produtos de maior valor agregado, sendo uma das aplicações de fluidização gás-sólido mais importante na indústria de petróleo. O presente trabalho avaliou especificamente a zona de injeção do FCC, na qual a matéria-prima, alimentada por bicos injetores, se mistura a sólidos quentes (catalisador) e a vapor de fluidização. A performance desses dispersores de carga para garantir uma boa distribuição das gotículas de gasóleo com o catalisador é a chave para melhorar a eficiência do riser de FCC. Desta forma, o principal objetivo deste trabalho foi avaliar o efeito de diferentes designs, ângulos (30°, 45° e 60°) e configurações dos injetores sobre o escoamento gás-sólido e o desempenho do riser. Para tal, simulou-se um escoamento gás-sólido reativo tridimensional baseado em uma abordagem Eulerian-Eulerian. Nas simulações foram utilizados o modelo cinético de 12-lumps de Wu et al. (2009), modelo de turbulência k-? e modelo de arraste de Gidaspow. Foi observado que o design, o ângulo e a configuração dos bicos injetores exercem uma forte influência sobre a fluidodinâmica e a performance do riser, sendo o ângulo a variável que apresentou maior influência. Pôde-se observar que o design de bico tipo multi-orifícios (Caso 3) foi o que apresentou os melhores resultados, sendo a partir deste avaliados os ângulos, onde notou-se que o aumento do ângulo de 30° para 60° melhorou a mistura entre as fases e o rendimento. A análise dos arranjos foi realizada utilizando o design de bico do Caso 3 e o ângulo de 45°, e observou-se que o arranjo com bicos intercalados (Arranjo 2) apresentou uma mistura mais homogênea entre as fases e, consequentemente, uma melhor conversão e rendimento de produtos desejados. Em geral, os resultados obtidos no presente trabalho salientam a importância da utilização de geometrias mais detalhadas para os bicos, uma vez que influenciam a mistura entre as fases, a qual afeta o desempenho do riser / Abstract: The Computational Fluid Dynamic (CFD) application in industrial process optimization and new equipments design studies has increased significantly, once it presents low cost and the possibility of evaluating complex and extreme operating conditions equipments . Among the most widely studied processes via CFD is the fluidized catalytic cracking process (FCC), where the oil heavy fractions of low-value are converted into higher value-added products and which is one of the most important gas-solid fluidization applications in the oil industry . The present study specifically evaluated the FCC injection zone, in which the feedstock fed by nozzles, is mixed with hot solids (catalyst) and fluidization steam. The nozzles performance to guarantee a good gas oil droplets distribution with the catalyst is the key to improve the efficiency of FCC riser. Thus, the study main objective was to evaluate the different nozzles designs, angles (30 °, 45 ° and 60 °) and arrangements effect on the gas-solid flow and the riser performance. For this purpose, it was simulated a three-dimensional reactive gas-solid flow based on an Eulerian-Eulerian approach. In the simulations it was used the 12 lumps kinetic model by Wu et al. (2009), turbulence model k- ? and drag model Gidaspow. It was observed that the nozzle design, angle and configuration have a strong influence on fluid dynamics and on the riser performance, and the angle was the variable with the greatest influence. It can be observed that the nozzle design of multi-orifice type (Case 3) showed the best results, and that¿s why it was used to evaluate the angle, in which was noted that the angle increase of 30 ° to 60 ° improved phases mixing and the yield. The arrangement analysis was performed using Case 3 nozzle design and the design angle of 45 °, and it was observed that the arrangement with intercalated nozzles (Arrangement 2) showed a more homogeneous phases mixture and therefore a better conversion and desired product yield. In general, the results obtained in this work highlighted the importance of using more detailed geometries for the nozzles, since they influence the mixing of the phases, which affects the riser performance / Mestrado / Engenharia Química / Mestre em Engenharia Química

Orientador: Milton Dall'Aglio Sobrinho / Banca: Tsunao Matsumoto / Banca: Gilson Arimura Arima / Resumo: O reator aeróbio de leito fluidizado com circulação apresenta características que permitem efetuar o tratamento em apenas uma etapa. O uso deste reator torna-se atraente devido à possibilidade de mantê-lo com uma alta concentração de biomassa ativa, propiciando seu porte compacto, conseqüência do reduzido período de detenção hidráulica. A eficiência do biofilme no tratamento depende das condições hidrodinâmicas do reator, uma vez que sua espessura pode ser alterada de acordo com a velocidade de circulação, e a transferência de oxigênio através de suas camadas depende da vazão de ar injetada no sistema. O trabalho propôs-se a investigar a influência da relação entre diâmetros dos tubos concêntricos em quatro configurações, mantendo o diâmetro externo constante e igual a 250 mm e alterando o diâmetro interno entre 100, 125, 150 e 200 mm, proporcionando relações Ae/Ai iguais a 5,25; 3,00; 1,78 e; 0,56, respectivamente. O tubo interno possui 2,22 m e é posicionado 0,05m acima do fundo do reator. A altura estática da coluna de fluido foi 2,35 m. Os ensaios hidrodinâmicos - velocidade do líquido, fração volumétrica de ar, tempo de mistura - e os de transferência de oxigênio foram realizados para vazões de ar até 2500 l/h em meio bifásico e em meio trifásico com concentrações de areia iguais a 30, 70, 100 e 150 g/l. O diâmetro interno de 125 mm foi o mais eficiente na indução de velocidade de circulação do líquido. A fração volumétrica de ar aumentou com o aumento do diâmetro interno para mesmas velocidades superficiais de ar, com exceção do menor diâmetro, o qual apresentou comportamento diferenciado dos demais. Os reator com diâmetro interno de 150 mm apresentou maiores valores de KLa para as mesmas vazões de ar. Observou-se que a adição de sólidos causou um decréscimo da velocidade do líquido e da fração... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Aerobic internal-loop airlift reactor has some characteristics that allow wastewater treatment in a single stage. The attractive to use this kind of reactor is mainly because of its capacity of maintaining a high active biomass concentration into it, which induces its small structure, in a consequence of its small hydraulic detention time. The biofilm efficiency in the treatment depends on reactor's hydrodynamic conditions. The structure of the biofilm and the oxygen transfer into its layers is affected by the liquid velocity induced by the air flow rate. This work proposed an investigation about the influence of the diameter of the concentric tubes, using four geometric configurations. The model had a 250 mm external tube and the inner draft tube with diameters of 100, 125, 150 and 200 mm, which means AD/AR equals to 5.25, 3.00, 1.78 and 0.56, respectively. The draft tube's height was 2,22 m and was positioned 0,05 m above the bottom of the reactor. The static fluid height was 2,35 m. The hydrodynamic tests - liquid velocities, gas hold up, mixture time - and the oxygen transfer tests were conducted for two and three-phase systems, with variable air flow rate up to 2500 l/h. The solid phase concentrations used were 30, 70, 100 and 150 g/l. The 125 mm reactor was the most efficient in inducing the liquid circulation velocity. The gas hold up increased with the increasing of the draft tube diameter, for same superficial gas velocities, excepting the 100 mm, which showed a different behavior. The 150 mm reactor presented the highest values for KLa for the same air rate. It was noticed that the solids addiction decreased the liquid velocity, gas hold up and oxygen transfer. Hydrodynamics equations found in literature showed a limited applicability, evidencing the necessity of relations capable to reach the most different configuration for this kind of reactor. A mass transfer mathematic...(Complete abstract click electronic acces below) / Mestre

Orientador: Marco Aurelio Cremasco / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-09T06:13:25Z (GMT). No. of bitstreams: 1 Castilho_GuilhermeJosede_M.pdf: 55720465 bytes, checksum: ff5e31cd4371ec0283beb9eec7df3f35 (MD5) Previous issue date: 2007 / Resumo: Regimes gás-sólido têm inúmeras aplicações industriais, sendo empregados, por exemplo, em sistemas de transporte pneumático das partículas e em leito fluidizado. Nas últimas duas décadas notou-se a aplicação crescente do leito fluidizado circulante, principalmente para facilidades de geração de energia, como caldeira e reações catalíticas. No processo de craqueamento catalítico fluidizado (FCC), dentre os vários equipamentos que compõem uma unidade conversora, é no riser onde praticamente ocorrem todas as reações de craqueamento. Medidas de grandezas como concentração de sólidos são fundamentais para a compreensão da fluidodinâmica nestes reatores e é pré-requisito para o projeto técnico bem sucedido de sistemas de transporte pneumático. Dentre as várias técnicas utilizadas para medidas de grandezas básicas em sistemas particulados, a sonda de fibra óptica é relativamente simples, de alta precisão e utilizada em várias situações que envolvem o sistema fluido-partícula. Nesta Dissertação, a oncentração volumétrica dos sólidos é medida por meio de uma sonda de fibra óptica. A concentração volumétrica dos sólidos é tomada radialmente e em diversas posições axiais ao longo do riser sob diferentes condições de operação, na Unidade Multipropósito do Laboratório de Processos em Meios Porosos da Faculdade de Engenharia Química da Universidade Estadual de Campinas. Resultados mostram que a concentração volumétrica de sólidos é maior na parede em comparação ao centro do riser, para todas posições axiais, exceto para o ponto mais alto do reator, onde há a influência da região de saída do escoamento. Sobre a distribuição axial de de concentração, verificou-se que o escoamento é mais concentrado na parte inferior do equipamento e tende a ficar menos concentrado ao longo do reator. Através dos espectros de sinais obtidos pela sonda, pôde-se analisar as flutuações através da distribuição de probabilidade dos sinais e da distribuição do desvio padrão do sinal médio. Distribuições de probabilidade permitem identificar a presença de clusters em determinados pontos radiais, revelando que para as condições de mais concentradas, as flutuações são as maiores. As distribuições de desvio padrão apresentam picos, onde sugere maior interação entre as partículas / Abstract: Gas-solid systems have innumerable industrial applications, being used, for example, in the pneumatic conveying of particles and in fluidized bed. In the last two decades it was noticed increasing application of the circulating fluidized bed, mainly for easinesses of energy generation, as boilers and catalytic reactions. In the process of fluidized catalytic cracking (FCC), amongst some equipment that composes a converting unit, it is in riser where the cracking reactions practically occur at all. Measures as solid concentration are basic for the understanding of the fluidodynamic in these reactors and are prerequisite for the successful design of pneumatic transport systems. Amongst the several techniques used for measures in particulates systems, the optical-fiber probe is relatively simple, of high precision and used in some situations that involve fluid-particle system. In this Dissertation, the volumetric concentration of solids is measured by means of an optical-fiber probe. The volumetric concentration of solids is taken radially and in diverse axial positions throughout riser under different operation conditions, in the Multipurpose Unit of the Laboratory of Processes in Porous Means of the College of Chemical Engineering of the State University of Campinas. Results show that the volumetric solid concentration is bigger in the wall in comparison to the center of the riser, for all axial positions, except forthe highest point of the reactor, where it has the influence of the exit region of the flow. On the axial distribution of concentration, it was verified that the flow is more concentrated in the inferior part of the equipment. Through the electrical signs specters gotten by the probe, it could be analyzed the fluctuations through the distribution of probability of the signals and the distribution of the standard deviation of the average signal. The distribution of probability allows identifying the presence of clusters in determined radial points, disclosing that for the most concentrated conditions, the fluctuations are the greater. The distributions of the standard deviation present peaks, where it suggests greater interaction between particles / Mestrado / Engenharia de Processos / Mestre em Engenharia Química

Orientador: Milton Mori / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-09-11T21:07:59Z (GMT). No. of bitstreams: 1 Bastos_JaciCarloSchrammCamara_M.pdf: 4360168 bytes, checksum: 100c57a7c93213cf37fcdab4ee05e24f (MD5) Previous issue date: 2005 / Resumo: A presente pesquisa apresenta a modelagem matemática e a simulação de uma operação de fluidização rápida em um longo duto cilíndrico vertical, operação importante em vários processos industriais, sendo que sua principal aplicação está vinculada ao craqueamento catalítico do petróleo para a conversão em gasolina. Tem como objetivos a obtenção de uma compreensão contínua e o conhecimento do desenvolvimento do escoamento gás-sólido, bem como a soma de experiências às pesquisas em escala industrial motivados pela possível predição do desempenho deste tipo de escoamento. O modelo tridimensional, turbulento e bifásico usado para a predição do escoamento gás-sólido, consiste num conjunto de equações de conservação da massa e momento para cada uma das fases, formuladas seguindo a aproximação Euleriana-Euleriana. As variáveis fluidodinâmicas foram estimadas pela solução do modelo, com o emprego de correlações empíricas da literatura e disponibilizadas pelo código computacional de CFD, para garantir o fechamento do modelo e sua solução numérica. Desde de que a predição da dinâmica do escoamento complexo, em dutos com alto fluxo ascendente de sólidos não é possível por meio somente de equações fundamentais, a maioria dos modelos requerem entradas empíricas, as quais somente são adquiridas com a experimentação. Estes dados foram obtidos dos estudos de PÄRSSINEN e ZHU (2001). A geometria e a malha numérica estrutural do duto vertical foram geradas pelo software (CAD) ICEM, subdividido em DDN (geometria) e Hexa (malha). A adaptação do modelo matemático para a geração do modelo numérico foi alcançada com o uso do simulador comercial CFX 5.7. Os resultados obtidos foram avaliados com respeito à teoria apresentada ao longo da dissertação, sendo finalmente feita uma comparação entre as predições numéricas com o modelo e dados experimentais da literatura / Abstract: The present research presents the mathematical modeling and the simulation of an operation of fast fluidization in a long vertical cylindrical duct line, which is an important operation in many industrial processes, where its main application is tied with the catalytic cracking of oil for gasoline synthesis. It has as objective the attainment of a continuous understanding and knowledge of the development of the gas-solid flow, as well as adding of experiences to research of industrial scale equipments motivated by the possibility of prediction of the performance of this type of flow. The three-dimensional, turbulent and two-phase model used for the prediction of the gas-solid flow, consists of a set of conservation equations of mass and momentum for each phase, which was formulated following the Eulerian-Eulerian approach. The fluid dynamics variables was estimated by the solution of the model with the use of correlations found in the literature and available in the computational CFD code, in order to guarantee the closure of the model and its numerical solution. Since the prediction of the dynamics of the complex flow in ducts with high ascending solid flow is not possible by solely using the basic equations, the majority of the models require the setting empirical, parameters which are acquired only with experimentation. These data were obtained from the studies of PÄRSSINEN and ZHU (2001). The geometry and the structural numerical mesh of the vertical duct was generated by the software (CAD) ICEM, subdivided in DDN (geometry) and Hexa (meshing). The adaptation of the mathematical model for the numerical model generation was reached with the use of the commercial simulator CFX 5.7. The results were evaluated with respect to the theory presented along this dissertation, where comparisons between the numerical predictions with the model and experimental data from the literature were performed. / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Milton Mori / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-10T23:15:40Z (GMT). No. of bitstreams: 1 Rosa_LeonardoMachadoda_D.pdf: 1367270 bytes, checksum: 1819e68c66a949c1c8fae024725df552 (MD5) Previous issue date: 2008 / Resumo: A aplicação de Fluidodinâmica Computacional (CFD) em operações industriais vem aumentando sensivelmente, em razão desta auxiliar no projeto de novos equipamentos, bem como a otimização dos já existentes. O presente trabalho apresenta uma série de simulações tridimensionais, bifásicas e adiabáticas, utilizando-se o modelo de dois-fluidos para o cálculo do escoamento gás-sólido em risers de diferentes tamanhos ¿ um com as dimensões de 10m de altura e 76mm de diâmetro interno, e outro com 18m de altura e 50mm de diâmetro interno. Utilizou-se o modelo k-epsilon para tratar a fase gasosa como turbulenta. Foram considerados os modelos cinéticos de 4 lumps e 10 lumps para o cálculo das taxas de craqueamento catalítico presentes no reator. Dados experimentais para comparação do escoamento foram obtidos em literatura, e os resultados obtidos mostram que, para a correta predição dos fenômenos presentes no interior de um riser, deve-se levar em consideração também as reações de craqueamento, visto que a alteração da composição química da fase gasosa interfere diretamente no comportamento dinâmico do escoamento / Abstract: The application of Computational Fluid Dynamics (CFD) in industrial operations has been increased not only in the project of new equipments but also in the optimization of the existing ones. The present work shows a set of tridimensional two-phase and adiabatic conditions, using the two-fluids model for the calculation of gas-solids flow in a high-flux risers of two different geometries, one with 10m of height and 76mm of internal diameter, and the second with 18m of height and 50mm of internal diameter. The k-epsilon model was applied to consider the influence of turbulence over the gas phase. The four-lump and ten-lump schemes kinetic models were used to predict the catalytic cracking reactions in the two different reactors. The riser with 10m of height and 76 mm of internal diameter had two entrances, one for particles inlet and one for the organic materials. The size of mesh for this case was 584.000 control volumes. The second riser studied with geometry of 18m of height and 50mm of internal iameter had one entrance for vapor, one for organic materials, and a third for particles inlet. The mesh generated over this geometry was considered 650.000 control volumes. The numerical results of gas-solids dispersion flow in a high-flux were compared with the xperimental data from the work of P¨arssinen and Zhu (2001), the four lump scheme for he mass concentration of chemical species was compared with the work of (ALI et al., 997), and the ten lump scheme was compares with the work of NAYAK et al. (2005). The catalytic cracking reactions models show an important role in the prediction of all henomena present in the reactor. The simulated results show that the kinetics models are important not only because of the variation of the gas phase and the temperature profiles inside the reactor but also because these parameters interfere directly in the dynamic behavior of the gas-solid flows / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientador: Milton Mori / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Química / Made available in DSpace on 2018-08-26T06:44:29Z (GMT). No. of bitstreams: 1 AlvarezCastro_HelverCrispiniano_D.pdf: 5280150 bytes, checksum: 96fb29e2849d11c06890186cb8af2f51 (MD5) Previous issue date: 2014 / Resumo: Conversão de matéria-prima e rendimento são estudados através de um modelo 3D que simula o principal reator do processo de craqueamento catalítico fluidizado (FCC). Fluidodinâmica computacional (CFD) é usada como técnica com abordagem Euleriana-Euleriana para prever o comportamento do riser da FCC. O modelo considera 12 grupos com desativação do catalisador dada por coque e envenenamento por nitritos alcalinas e adsorção aromático policíclico para estimar o comportamento cinético, o qual partindo de uma determinada matéria prima, se craquea a diversos produtos. Composições diferentes para matérias-primas foram consideradas. Para validação o modelo foi comparado com dados experimentais em condições de operação industrial. O modelo é capaz de representar com precisão a geração de produtos para as diferentes condições de operação consideradas. Todas as condições consideradas foram resolvidos usando o solver Ansys CFX 14.0. As diferentes variáveis do processo de operação e efeitos hidrodinâmicos do riser industrial de craqueamento catalítico fluidizado (FCC) foram avaliadas. As previsões do modelo são apresentados e comparadas; conversão e produtos são apresentados. Recomendações para criar condições para obter maior rendimento do produto no processo industrial foram dadas / Abstract: Feedstock conversion and yield products are studied through a 3D model simulating the main reactor of the fluid catalytic cracking (FCC) process. Computational Fluid Dynamic (CFD) is used with Eulerian-Eulerian approach to predict the fluid catalytic cracking behavior. The model considers 12 lumps with catalyst deactivation by coke and poisoning by alkaline nitrides and polycyclic aromatic adsorption to estimate the kinetic behavior which, starting from a given feedstock, produces several cracking products. Different feedstock compositions are considered. The model is compared with sampling data at industrial operation conditions. The simulation model is able to represent accurately the products behavior for the different operating conditions considered. All the conditions considered were solved using a solver Ansys CFX 14.0. The different operation process variables and hydrodynamic effects of the industrial riser of a fluid catalytic cracking (FCC) are evaluated. Predictions from the model are shown and comparison with experimental conversion and yields products are presented; recommendations are drawn to establish the conditions to obtain higher product yields in the industrial process / Doutorado / Engenharia de Processos / Doutor em Engenharia Química