ES2952720A1 - Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding) - Google Patents

Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding) Download PDF

Info

Publication number
ES2952720A1
ES2952720A1 ES202330342A ES202330342A ES2952720A1 ES 2952720 A1 ES2952720 A1 ES 2952720A1 ES 202330342 A ES202330342 A ES 202330342A ES 202330342 A ES202330342 A ES 202330342A ES 2952720 A1 ES2952720 A1 ES 2952720A1
Authority
ES
Spain
Prior art keywords
medium
size distribution
refractive index
particle size
procedure
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
ES202330342A
Other languages
Spanish (es)
Inventor
Chorro Miguel Angel Báez
Rodriguez Borja Vidal
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Universidad Politecnica de Valencia
Original Assignee
Universidad Politecnica de Valencia
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Universidad Politecnica de Valencia filed Critical Universidad Politecnica de Valencia
Priority to ES202330342A priority Critical patent/ES2952720A1/en
Publication of ES2952720A1 publication Critical patent/ES2952720A1/en
Priority to PCT/ES2024/070254 priority patent/WO2024223970A1/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/02Investigating particle size or size distribution
    • G01N15/0205Investigating particle size or size distribution by optical means

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field. Procedure and system to estimate the particle size distribution by combining information on the modulus and phase of the scattered electromagnetic field. The detection of the forward spread signal is carried out with means capable of obtaining the module information and the phase of the electric field. From these measurements it is possible to directly derive the extinction coefficient and the refractive index. The size distribution is then calculated using an iterative Twomey method, where the factor for each iteration is modified so that it includes information about the extinction coefficient and refractive index. (Machine-translation by Google Translate, not legally binding)

Description

DESCRIPCIÓNDESCRIPTION

Procedimiento y sistema para la estimación de la distribución del tamaño de partículas en medios heterogéneos a partir del campo óptico complejoProcedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field

SECTOR DE LA TÉCNICATECHNIQUE SECTOR

La invención pertenece al sector del estudio de sistemas de partículas y en especial en agregados, aerosoles, emulsiones, polvos y suspensiones.The invention belongs to the sector of the study of particle systems and especially aggregates, aerosols, emulsions, powders and suspensions.

Más en particular, la presente invención se refiere a un procedimiento para la estimación de la distribución del tamaño de múltiples partículas a partir de la detección coherente de la luz esparcida hacia adelante.More particularly, the present invention relates to a method for estimating the size distribution of multiple particles from the coherent detection of forward scattered light.

ESTADO DE LA TÉCNICASTATE OF THE TECHNIQUE

Existe la necesidad de comprender la naturaleza de los sistemas de partículas en numerosos campos, desde la farmacología a la ciencia de materiales. Con diferencia, la característica más importante de las muestras de partículas es la distribución del tamaño. Esta característica determina la calidad, la seguridad y el rendimiento del producto final.There is a need to understand the nature of particle systems in numerous fields, from pharmacology to materials science. By far the most important characteristic of particle samples is the size distribution. This characteristic determines the quality, safety and performance of the final product.

El proceso de estimación de la distribución del tamaño de partículas (DTP) generalmente se organiza en los siguientes pasos:The particle size distribution (PSD) estimation process is generally organized into the following steps:

(i) Preparación: procesamiento de muestras, configuración de instrumentos, calibración y alineamiento, instalación del sistema de entrega de muestras.(i) Preparation: sample processing, instrument setup, calibration and alignment, installation of sample delivery system.

(ii) Medición: iluminar la muestra con una señal y medir la intensidad de la luz esparcida hacia adelante por la muestra para compararla con la intensidad sin muestra y de ahí derivar el coeficiente de extinción de la muestra.(ii) Measurement: Illuminate the sample with a signal and measure the intensity of the forward scattered light from the sample to compare it with the intensity without sample and from there derive the extinction coefficient of the sample.

(iii) Análisis directo: determinar el efecto del esparcimiento de la luz provocado por las partículas de la muestra bajo prueba sobre el coeficiente de extinción, asumiendo parámetros específicos para las partículas (forma, índice de refracción, concentración, etc.). (iii) Direct analysis: determine the effect of light scattering caused by the particles of the sample under test on the extinction coefficient, assuming specific parameters for the particles (shape, refractive index, concentration, etc.).

(iv) Análisis inverso: determinar las características desconocidas de la muestra utilizando los resultados del análisis directo y los datos de medición correspondientes. Un paso esencial de esta etapa es el análisis del contenido de información de las ecuaciones a invertir. Un análisis inverso generalmente se complementa con información a priori sobre la solución de interés.(iv) Reverse analysis: determine the unknown characteristics of the sample using the results of the forward analysis and the corresponding measurement data. An essential step of this stage is the analysis of the information content of the equations to be inverted. A reverse analysis is usually supplemented with a priori information about the solution of interest.

(v) Estadísticas y análisis de datos: comprender los datos de cálculo, obtener la distribución del tamaño de las partículas utilizando una variedad de modelos (número/volumen/masa).(v) Statistics and data analysis: understand the calculation data, obtain the particle size distribution using a variety of models (number/volume/mass).

Para la modelización del coeficiente de extinción t ( v ) de la muestra en el paso de análisis directo (ii) es común utilizar la siguiente expresión:To model the extinction coefficient t (v) of the sample in the direct analysis step (ii), it is common to use the following expression:

Figure imgf000003_0001
Figure imgf000003_0001

00

Donde N0 es la cantidad total de partículas por unidad de volumen, D es el diámetro de las partículas, v es la frecuencia del campo electromagnético empleado en la medida, f (D) es la DTP normalizada del medio y Qext(v,D) es el factor de eficiencia de extinción para una sola partícula esférica, dado por la teoría de Mie.Where N0 is the total number of particles per unit volume, D is the diameter of the particles, v is the frequency of the electromagnetic field used in the measurement, f ( D) is the normalized DTP of the medium and Qext ( v,D) is the extinction efficiency factor for a single spherical particle, given by Mie theory.

Para obtener la DTP es necesario invertir la ecuación anterior que, tras la adecuada discretización de la integral, puede expresarse en forma de sistema de ecuaciones lineales:To obtain the DTP it is necessary to invert the previous equation which, after adequate discretization of the integral, can be expressed in the form of a system of linear equations:

Figure imgf000003_0002
Figure imgf000003_0002

Sin embargo, la obtención de la solución f mediante la inversión de la matriz K no es posible, ya que estos sistemas suelen presentar un mal condicionamiento, es decir, pequeños cambios en los valores de sus elementos generan soluciones radicalmente diferentes. Para la adecuada solución de este tipo de problema es necesario aplicar técnicas de inversión regularizada.However, obtaining the solution f by inverting the matrix K is not possible, since these systems usually present poor conditioning, that is, small changes in the values of their elements generate radically different solutions. For the adequate solution of this type of problem, it is necessary to apply regularized investment techniques.

Un método conocido para la inversión de la distribución de tamaño de partículas a partir de medidas de extinción, (paso iv) es el método de Twomey. Éste corresponde a la categoría de métodos iterativos que obtienen una estimación de la solución tras la repetida multiplicación de la solución previa por un cierto factor partiendo desde una solución de prueba inicial. En concreto, la solución iterativa de la distribución para un diámetro Dj según el método de Twomey tiene la siguiente forma:A known method for inversion of particle size distribution from extinction measurements (step iv) is the Twomey method. This corresponds to the category of iterative methods that obtain an estimate of the solution after repeatedly multiplying the previous solution by a certain factor starting from an initial test solution. In Specifically, the iterative solution of the distribution for a diameter Dj according to the Twomey method has the following form:

Figure imgf000004_0001
Figure imgf000004_0001

El factor multiplicativo en cada iteración depende, por un lado, del cociente entre la extinción medida (rmeas) y la calculada (t^ ) con el modelo de Mie (ecuación 1) empleando la solución de la iteración previa (fp). Por otro lado, este cociente está pesado por una función (Wtj ) calculada previamente a partir del modelo de Mie (ecuación 1) y cuya función es regular la contribución de cada tamaño de partícula proporcionalmente a su contribución en el valor de los parámetros ópticos. El factor iterativo finalmente se obtiene como el producto de las cantidades anteriores para cada frecuencia, v¿.The multiplicative factor in each iteration depends, on the one hand, on the ratio between the measured extinction (rmeas) and that calculated ( t ^ ) with the Mie model (equation 1) using the solution from the previous iteration ( fp). On the other hand, this quotient is weighted by a function ( Wtj ) previously calculated from the Mie model (equation 1) and whose function is to regulate the contribution of each particle size proportionally to its contribution to the value of the optical parameters. The iterative factor is finally obtained as the product of the previous quantities for each frequency, v¿.

Una ventaja de este tipo de métodos iterativos es que aseguran que la solución tenga un valor positivo. Además, permiten trabajar con cantidades experimentales cuya descripción matemática no esté relacionada linealmente con la distribución de tamaños de partícula.An advantage of this type of iterative methods is that they ensure that the solution has a positive value. Furthermore, they allow working with experimental quantities whose mathematical description is not linearly related to the particle size distribution.

Puesto que conocer la distribución del tamaño de partículas es crítico en muchas aplicaciones, es deseable disponer de un método que consiga una mayor precisión en el cálculo de dicha distribución.Since knowing the particle size distribution is critical in many applications, it is desirable to have a method that achieves greater precision in the calculation of said distribution.

RESUMEN DE LA INVENCIÓNSUMMARY OF THE INVENTION

La presente invención parte de la superación del prejuicio técnico de que la estimación de la distribución de tamaños de partícula de un medio heterogéneo a partir de la información sobre intensidad de la señal dispersada ofrece los mismos resultados que la estimación a partir de la información de fase y que la agregación de estas dos es redundante. Por el contrario, se ha encontrado que una modificación del método de Twomey para la resolución del problema inverso a partir de la luz esparcida hacia delante permite incorporar tanto la información de fase como la del módulo del campo para obtener la distribución de tamaños de partícula de forma más precisa.The present invention is based on overcoming the technical prejudice that the estimation of the particle size distribution of a heterogeneous medium from the information on the intensity of the dispersed signal offers the same results as the estimation from the phase information. and that the aggregation of these two is redundant. On the contrary, it has been found that a modification of the Twomey method for solving the inverse problem from forward scattered light allows incorporating both phase information and field modulus information to obtain the particle size distribution of more precise way.

La invención proporciona un procedimiento y sistema para estimar la distribución del tamaño de las partículas combinando pues la información del módulo y fase del campo electromagnético esparcido. Para ello, la detección de la señal esparcida hacia delante se debe realizar con medios capaces de obtener la información del módulo y fase del campo eléctrico. A partir de estas medidas es posible derivar directamente el coeficiente de extinción y el índice de refracción.The invention provides a procedure and system to estimate the size distribution of particles, thus combining the information of the module and phase of the field. scattered electromagnetic. To do this, the detection of the signal spread forward must be carried out with means capable of obtaining information on the module and phase of the electric field. From these measurements it is possible to directly derive the extinction coefficient and the refractive index.

BREVE DESCRIPCIÓN DE LAS FIGURASBRIEF DESCRIPTION OF THE FIGURES

Con objeto de ayudar a una mejor comprensión de las características de la invención y para complementar esta descripción, se acompañan como parte integrante de la misma las siguientes figuras, cuyo carácter es ilustrativo y no limitativo:In order to help a better understanding of the characteristics of the invention and to complement this description, the following figures are attached as an integral part of it, the nature of which is illustrative and not limiting:

La Figura 1 muestra los elementos esenciales de la invención.Figure 1 shows the essential elements of the invention.

Las Figuras 2a-2d muestran resultados experimentales comparando la distribución obtenida mediante el procedimiento del estado de la técnica y aquella obtenida mediante el procedimiento de la invención.Figures 2a-2d show experimental results comparing the distribution obtained by the state of the art method and that obtained by the method of the invention.

DESCRIPCIÓN DETALLADADETAILED DESCRIPTION

El procedimiento para la estimación de la distribución del tamaño de múltiples partículas de la invención comprende los siguientes pasos:The procedure for estimating the size distribution of multiple particles of the invention comprises the following steps:

- Iluminar la muestra con una fuente de radiación que incluya distintas frecuencias (Fig. 1) ;- Illuminate the sample with a radiation source that includes different frequencies (Fig. 1);

- detectar una porción de radiación esparcida hacia adelante mediante un detector capaz de obtener el campo óptico complejo, es decir, capaz de extraer la información tanto del módulo como de la fase del campo electromagnético;- detect a portion of radiation spread forward by means of a detector capable of obtaining the complex optical field, that is, capable of extracting information from both the module and the phase of the electromagnetic field;

- medir tanto el efecto del medio heterogéneo en el campo complejo (campo muestra) como el campo en ausencia de dicho medio (campo referencia), es decir, medir el campo electromagnético de la luz esparcida hacia adelante por la muestra para compararla con el campo sin muestra- measure both the effect of the heterogeneous medium on the complex field (sample field) and the field in the absence of said medium (reference field), that is, measure the electromagnetic field of the light scattered forward by the sample to compare it with the field no sample

- calcular la transmisión, T, a través de la muestra como el valor absoluto del cociente de los módulos de los campos muestra y referencia. Calcular la diferencia de fase A0 como la resta de las fases de los campos muestra y referencia; - calculate the transmission, T, through the sample as the absolute value of the quotient of the modules of the sample and reference fields. Calculate the phase difference A0 as the subtraction of the phases of the sample and reference fields;

- derivar los parámetros ópticos del medio heterogéneo bajo prueba ( n exp, índice de refracción y rexp, extinción) a partir de la información experimental de transmisión y de la diferencia de fase calculadas anteriormente mediante las siguientes fórmulas:- derive the optical parameters of the heterogeneous medium under test ( n exp, refractive index and rexp, extinction) from the experimental transmission information and the phase difference calculated previously using the following formulas:

Figure imgf000006_0003
Figure imgf000006_0003

donde d es el espesor de la muestra, v la frecuencia y c la velocidad de la luz en el vacío;where d is the thickness of the sample, v the frequency and c the speed of light in a vacuum;

modelizar los parámetros ópticos de la muestra heterogénea, para lo que puede emplearse la teoría de Waterman-Truell para medios con grandes cantidades de partículas, que tiene en cuenta el esparcimiento múltiple de la luz, cuyo efecto es notable a partir de pequeñas fracciones volumétrica de partículas (>1%). Esta teoría permite calcular el índice de refracción y la extinción con la siguiente fórmula:model the optical parameters of the heterogeneous sample, for which the Waterman-Truell theory can be used for media with large quantities of particles, which takes into account the multiple scattering of light, the effect of which is notable from small volumetric fractions of particles (>1%). This theory allows the refractive index and extinction to be calculated with the following formula:

Figure imgf000006_0004
Figure imgf000006_0004

donde Re(-) y Im(-) hacen referencia a tomar la parte real e imaginaria del número complejo en el paréntesis. (S(0/n)) son las amplitudes de Mie para ángulos 0 y n para una sola partícula. kmed — - n med es la constante de propagación del medio matriz (es decir, sin tener en cuenta el efecto de las partículas agregadas) y nmed el índice de refracción de dicho medio. El símbolo {) indica un promedio de las amplitudes S pesadas con la DTP del medio;where Re ( -) and Im ( -) refer to taking the real and imaginary part of the complex number in the parentheses. ( S ( 0/n)) are the Mie amplitudes for angles 0 and n for a single particle. kmed — - n med is the propagation constant of the matrix medium (i.e., without taking into account the effect of added particles) and nmed is the refractive index of said medium. The symbol {) indicates an average of the heavy S amplitudes with the DTP of the middle;

obtener la distribución del tamaño de las partículas a partir de los parámetros ópticos mediante un nuevo algoritmo basado en la modificación del método de Twomey para incorporar al factor la contribución de las medidas de índice de refracción además de las del coeficiente de extinción:obtain the particle size distribution from the optical parameters using a new algorithm based on the modification of the Twomey method to incorporate into the factor the contribution of the refractive index measurements in addition to those of the extinction coefficient:

Figure imgf000006_0001
Figure imgf000006_0001

donde nexp es el índice de refracción medido experimentalmente

Figure imgf000006_0002
es el índice de refracción calculado usando la ecuación 5 con la DTP de la iteración p. Wij son los pesos calculados mediante la ecuación. 3 y Wi] es el análogo empleando el índice de refracciónwhere nexp is the experimentally measured refractive index
Figure imgf000006_0002
is the refractive index calculated using equation 5 with the DTP of iteration p. Wij are the weights calculated using the equation. 3 and Wi] is the analogue using the refractive index

W n _ ncalcjv i>Dj ) ™ ij l t n calc(vi,DJ) ( ) - el parámetro a representa un equilibrio entre la importancia que tiene la extinción o el índice de refracción en la inversión. Este parámetro se elegirá como aquel que produzca una distribución de tamaños cuyos parámetros ópticos calculados se asemejan más a los experimentales. Para ello se elige a como aquel que produce un mínimo de la métrica ERM, donde ERM es W n _ ncalcjv i>Dj ) ™ ij ltn calc ( vi,DJ) ( ) - the parameter a represents a balance between the importance of extinction or refractive index in the inversion. This parameter will be chosen as the one that produces a size distribution whose calculated optical parameters most closely resemble the experimental ones. To do this, a is chosen as the one that produces a minimum of the ERM metric, where ERM is

Figure imgf000007_0001
Figure imgf000007_0001

- esta métrica se basa en determinar el error relativo medio (ERM) que es la suma de la media de los residuos relativos de la extinción y el índice de refracción. Se ha determinado que los valores mínimos de esta métrica corresponden con valores de a que se encuentran próximos a los valores óptimos de forma que puede considerarse un buen indicador para determinar a sin más información que las medidas experimentales. En caso de que el valor de ERM no presente un mínimo muy pronunciado, es posible promediar las soluciones en un rango de valores de a alrededor del mínimo. Se ha observado que un rango 0.25 proporciona una estimación robusta de la verdadera DTP en la mayoría de los casos.- This metric is based on determining the mean relative error (ERM), which is the sum of the mean of the relative residuals of the extinction and the refractive index. It has been determined that the minimum values of this metric correspond to values of a that are close to the optimal values so that it can be considered a good indicator to determine a without more information than experimental measurements. In case the ERM value does not present a very pronounced minimum, it is possible to average the solutions in a range of values from to around the minimum. It has been observed that a range of 0.25 provides a robust estimate of the true DTP in most cases.

Resultados experimentales:Experimental results:

La técnica propuesta para la obtención del tamaño de partículas ha sido validada experimentalmente. En las figuras 2a y 2b se muestran parámetros ópticos de muestras granuladas hechas a partir de polvo de PTFE y microesferas de vidrio. La fracción de volumen acumulada de las esferas es de 5%. La figura 2c muestra los resultados de inversión de la DTP utilizando el método de Twomey original (línea punteada) y utilizando el método propuesto (línea continua) junto a la DTP dada por el fabricante. En la figura 2d se muestra el valor de la métrica utilizada para determinar el valor del parámetro a.The proposed technique for obtaining particle size has been experimentally validated. Optical parameters of granulated samples made from PTFE powder and glass microspheres are shown in Figures 2a and 2b. The cumulative volume fraction of the spheres is 5%. Figure 2c shows the DTP inversion results using the original Twomey method (dotted line) and using the proposed method (solid line) together with the DTP given by the manufacturer. Figure 2d shows the value of the metric used to determine the value of parameter a.

En las Figuras 2a y 2b se pueden observar los parámetros ópticos medidos (negro) así como los calculados a partir del modelo de esparcimiento de la luz (Tcatc, ncalc) usando como entrada la DTP invertida. Los resultados de la inversión empleando tanto el método de Twomey tradicional que sólo tiene en cuenta la extinción (línea punteada) como el método propuesto aunando información de extinción e índice de refracción (línea continua) se muestran en la inferior junto a la DTP dada por el fabricante. En la Fig. 2d se muestran los valores de ERM obtenidos con varios parámetros a y la posición del mínimo empleado en las inversiones. El método propuesto consigue disminuir la desviación respecto a la DTP tabulada consiguiendo un error menor que el método de Twomey. Concretamente, la mejora resulta de un 32.6%.In Figures 2a and 2b you can see the measured optical parameters (black) as well as those calculated from the light scattering model ( Tcatc, ncalc) using the inverted DTP as input. The investment results using both the Twomey method traditional that only takes into account the extinction (dotted line) as the proposed method combining extinction information and refractive index (solid line) are shown in the bottom next to the DTP given by the manufacturer. Fig. 2d shows the ERM values obtained with various parameters a and the position of the minimum used in the inversions. The proposed method manages to reduce the deviation with respect to the tabulated DTP, achieving a smaller error than the Twomey method. Specifically, the improvement is 32.6%.

Esta validación experimental se ha efectuado en el infrarrojo lejano empleando un instrumento de Espectroscopía de Terahercios en Dominio Temporal que permite obtener el coeficiente de extinción y el índice de refracción de muestras en la banda de 0.2 a 1.6 THz. Un experto entenderá que sistemas que empleen antenas y detectores de RF o microondas, interferómetros espectrales de baja coherencia e interferómetros con fuentes coherentes sintonizables también son capaces de proporcionar la información sobre módulo y fase del campo electromagnético esparcido en sus correspondientes rangos espectrales de operación y por lo tanto son también aptos para su uso en el presente procedimiento.This experimental validation has been carried out in the far infrared using a Time Domain Terahertz Spectroscopy instrument that allows obtaining the extinction coefficient and refractive index of samples in the band from 0.2 to 1.6 THz. An expert will understand that systems that employ RF or microwave antennas and detectors, low-coherence spectral interferometers and interferometers with tunable coherent sources are also capable of providing information on the modulus and phase of the electromagnetic field spread in their corresponding spectral ranges of operation and therefore Therefore they are also suitable for use in the present procedure.

A la vista de esta descripción y figuras, el experto en la materia podrá entender que la invención ha sido descrita según algunas realizaciones preferentes de la misma, pero que múltiples variaciones pueden ser introducidas en dichas realizaciones preferentes, sin exceder el objeto de la invención tal y como ha sido reivindicada. In view of this description and figures, the person skilled in the art will be able to understand that the invention has been described according to some preferred embodiments thereof, but that multiple variations can be introduced in said preferred embodiments, without exceeding the object of the invention as such. and how it has been claimed.

Claims (5)

REIVINDICACIONES 1. Procedimiento para el cálculo de la estimación de la distribución del tamaño de partículas en medios heterogéneos, donde el procedimiento comprende los pasos de: - a) iluminar el medio con una fuente de radiación a distintas frecuencias;1. Procedure for calculating the estimation of the particle size distribution in heterogeneous media, where the procedure includes the steps of: - a) illuminating the medium with a radiation source at different frequencies; - b) medir la luz esparcida hacia adelante por el medio mediante un detector capaz de obtener el campo óptico complejo para compararla con la señal sin medio;- b) measure the light scattered forward through the medium using a detector capable of obtaining the complex optical field to compare it with the signal without medium; - c) calcular la transmisión, T, a través del medio como el valor absoluto del cociente de los módulos de los campos con medio y sin medio y calcular la diferencia de fase A<p como la resta de las fases de dichos campos muestra y referencia- c) calculate the transmission, T, through the medium as the absolute value of the quotient of the modules of the fields with medium and without medium and calculate the phase difference A<p as the subtraction of the phases of said sample fields and reference - d) derivar el índice de refracción nexp y el coeficiente de extinción rexp, del medio heterogéneo a partir de la información experimental de transmisión y de la diferencia de fase calculadas anteriormente mediante las siguientes fórmulas:- d) derive the refractive index nexp and the extinction coefficient rexp, of the heterogeneous medium from the experimental transmission information and the phase difference calculated previously using the following formulas:
Figure imgf000009_0001
Figure imgf000009_0001
donde d es el espesor de la muestra, v la frecuencia y c la velocidad de la luz en el vacío;where d is the thickness of the sample, v the frequency and c the speed of light in a vacuum; - e) modelizar los parámetros ópticos del medio heterogéneo empleando la teoría de Waterman-Truell para medios con grandes cantidades de partículas para obtener los valores del índice de refracción ncaic y el coeficiente de extinción r caíc;- e) model the optical parameters of the heterogeneous medium using the Waterman-Truell theory for media with large amounts of particles to obtain the values of the refractive index ncaic and the extinction coefficient r caíc; - f) obtener la distribución del tamaño de las partículas mediante el método iterativo de Twomey, donde el factor para cada iteración está modificado según la siguiente expresión:- f) obtain the particle size distribution using the Twomey iterative method, where the factor for each iteration is modified according to the following expression:
Figure imgf000009_0002
Figure imgf000009_0002
donde nexp es el índice de refracción medido experimentalmente y n f y es el índice de refracción calculado con la DTP de la iteración p.W-j son los pesos calculados mediante la ecuación. 3 y W-j es el análogo empleando el índice de refracciónwhere nexp is the experimentally measured refractive index and nfy is the refractive index calculated with the DTP of the pW-j iteration are the weights calculated using the equation. 3 and Wj is the analogue using the refractive index
Figure imgf000009_0003
Figure imgf000009_0003
y donde el parámetro a se elige como aquel que produce un mínimo de la métrica ERM, donde ERM esand where the parameter a is chosen as the one that produces a minimum of the ERM metric, where ERM is
Figure imgf000010_0001
Figure imgf000010_0001
2. Sistema para llevar a cabo el procedimiento de la reivindicación 1 que comprende: una fuente de radiación con distintas frecuencias, un detector capaz de obtener el campo óptico complejo y un procesador provisto de medios de programa para realizar los pasos c-f.2. System for carrying out the procedure of claim 1 comprising: a radiation source with different frequencies, a detector capable of obtaining the complex optical field and a processor provided with program means to carry out steps c-f. 3. Sistema según la reivindicación 2 caracterizado por que el detector es un espectroscopio de terahercios en dominio temporal, un detector de RF o microondas, un interferómetro espectral de baja coherencia o un interferómetro con fuentes coherentes sintonizables.3. System according to claim 2 characterized in that the detector is a time domain terahertz spectroscope, an RF or microwave detector, a low coherence spectral interferometer or an interferometer with tunable coherent sources. 4. Producto de programa de ordenador que comprende instrucciones para llevar a cabo los pasos c-f de la reivindicación 1.4. Computer program product comprising instructions for carrying out steps c-f of claim 1. 5. Medio legible por ordenador que comprende almacenado el programa capaz de llevar a cabo los pasos c-f de la reivindicación 1. 5. Computer-readable medium comprising stored the program capable of carrying out steps c-f of claim 1.
ES202330342A 2023-04-28 2023-04-28 Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding) Pending ES2952720A1 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
ES202330342A ES2952720A1 (en) 2023-04-28 2023-04-28 Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding)
PCT/ES2024/070254 WO2024223970A1 (en) 2023-04-28 2024-04-23 Method and system for estimating particle size distribution in heterogeneous media based on the complex optical field

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
ES202330342A ES2952720A1 (en) 2023-04-28 2023-04-28 Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding)

Publications (1)

Publication Number Publication Date
ES2952720A1 true ES2952720A1 (en) 2023-11-03

Family

ID=88558708

Family Applications (1)

Application Number Title Priority Date Filing Date
ES202330342A Pending ES2952720A1 (en) 2023-04-28 2023-04-28 Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding)

Country Status (2)

Country Link
ES (1) ES2952720A1 (en)
WO (1) WO2024223970A1 (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0464337A2 (en) * 1990-07-02 1992-01-08 International Business Machines Corporation Measurement of size and refractive index of particles using the complex forward-scattered electromagnetic field

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0464337A2 (en) * 1990-07-02 1992-01-08 International Business Machines Corporation Measurement of size and refractive index of particles using the complex forward-scattered electromagnetic field

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SAWYER N B E ET AL. Amplitude and phase microscopy for sizing of spherical particles. Applied Optics 1 Aug. 2003 Opt. Soc. America USA. , 01/08/2003, Vol. 42, Páginas 4488 - 4498 ISSN 0003-6935 (print), *
SWIRNIAK GRZEGORZ ET AL. Forward and inverse analysis for particle size distribution measurements of disperse samples: A review. MEASUREMENT, 20211019 INSTITUTE OF MEASUREMENT AND CONTROL. LONDON, GB. , 19/10/2021, Vol. 187, ISSN 0263-2241, (DOI: doi:10.1016/j.measurement.2021.110256) Apartados 1, 2.3, 4.3.2 *
TAUBENBLATT M A ET AL. Measurement of the size and refractive index of a small particle using the complex forward-scattered electromagnetic field. Applied Optics 20 Nov. 1991 USA. , 20/11/1991, Vol. 30, Páginas 4972 - 4979 ISSN 0003-6935 (print), *
ZHOU J ET AL. Particle size and refractive index retrieval from the backscattering spectrum of white light using the Twomey iterative method: Simulation and experiment. Applied Optics 20060910 OSA - The Optical Society usa. , 10/09/2006, Vol. 45, Páginas 6876 - 6885 ISSN 1559-128X (print) ISSN 1539-4522 (electronic), (DOI: doi:10.1364/AO.45.006876) Apartado 2 *

Also Published As

Publication number Publication date
WO2024223970A1 (en) 2024-10-31

Similar Documents

Publication Publication Date Title
Onasch et al. Single scattering albedo monitor for airborne particulates
Koda et al. Are peculiar velocity surveys competitive as a cosmological probe?
Hill Review of optical scintillation methods of measuring the refractive-index spectrum, inner scale and surface fluxes
Bellagamba et al. Optimal filtering of optical and weak lensing data to search for galaxy clusters: application to the COSMOS field
Zhang et al. Nondestructive measurement of hemoglobin in blood bags based on multi-pathlength VIS-NIR spectroscopy
Ye et al. Measurements of Particle Size Distribution Based on Mie Scattering Theory and Markov Chain Inversion Algorithm.
Magazù et al. FTIR spectroscopy studies on the bioprotective effectiveness of trehalose on human hemoglobin aqueous solutions under 50 Hz electromagnetic field exposure
Tang et al. Retrieval of spheroid particle size distribution from spectral extinction data in the independent mode using PCA approach
CN107478555A (en) Gas particles thing measuring method and device
KR20090105545A (en) Apparatus for measuring the Refractive Index and Thickness of Phase Object using Three-different-wavelength Light Sources and the Method therewith
ES2952720A1 (en) Procedure and system for estimating the particle size distribution in heterogeneous media from the complex optical field (Machine-translation by Google Translate, not legally binding)
Brenner et al. Optical coherence tomography images simulated with an analytical solution of Maxwell’s equations for cylinder scattering
Yang et al. Blood cell counting and classification by nonflowing laser light scattering method
Ma et al. The Measurement and modeling investigation on the spectral polarized BRDF of brass
RU2610942C1 (en) Method for optical measurement of calculating concentration of dispersed particles in liquid environments and device for its implementation
Sviridenkov et al. Parameterization of asymmetry factor of atmospheric aerosol scattering phase function
CN108872152B (en) Particle refractive index measuring method, computer device and computer readable storage medium
Buccini et al. Weakly constrained Lucy–Richardson with applications to inversion of light scattering data
Lastovka Light mixing spectroscopy and the spectrum of light scattered by thermal fluctuations in liquids.
Ustinov On inverse reconstruction problems of the erythrocyte size distribution in laser diffractometry
Chen et al. Optical aerosol sizing method without prior refractive index
Yurkin Symmetry relations for the Mueller scattering matrix integrated over the azimuthal angle
Huang et al. Dual-prism based terahertz time-domain attenuated total reflection spectroscopy and its application to characterise the hydration state of L-threonine in solution
Giacomelli et al. Review of the application of T-matrix calculations for determining the structure of cell nuclei with angle-resolved light scattering measurements
Gu et al. Mass concentration calculation with the pulse height distribution of aerosols and system calibration

Legal Events

Date Code Title Description
BA2A Patent application published

Ref document number: 2952720

Country of ref document: ES

Kind code of ref document: A1

Effective date: 20231103