ES2952720B2 - Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical field - Google Patents
Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical fieldInfo
- Publication number
- ES2952720B2 ES2952720B2 ES202330342A ES202330342A ES2952720B2 ES 2952720 B2 ES2952720 B2 ES 2952720B2 ES 202330342 A ES202330342 A ES 202330342A ES 202330342 A ES202330342 A ES 202330342A ES 2952720 B2 ES2952720 B2 ES 2952720B2
- Authority
- ES
- Spain
- Prior art keywords
- medium
- refractive index
- size distribution
- particle size
- carrying
- 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.)
- Active
Links
- 239000002245 particle Substances 0.000 title claims description 32
- 238000000034 method Methods 0.000 title claims description 30
- 230000003287 optical effect Effects 0.000 title claims description 13
- 230000008033 biological extinction Effects 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000005855 radiation Effects 0.000 claims description 4
- 230000001427 coherent effect Effects 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000004590 computer program Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 9
- 239000000243 solution Substances 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 8
- 230000005672 electromagnetic field Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000000149 argon plasma sintering Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000843 powder Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000000443 aerosol Substances 0.000 description 1
- 239000012615 aggregate Substances 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000000839 emulsion Substances 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000013101 initial test Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000004005 microsphere Substances 0.000 description 1
- 239000004810 polytetrafluoroethylene Substances 0.000 description 1
- 229920001343 polytetrafluoroethylene Polymers 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000012798 spherical particle Substances 0.000 description 1
- 239000000725 suspension Substances 0.000 description 1
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 description 1
- 239000012085 test solution Substances 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
- G01N15/0205—Investigating 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)
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 complejo Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical field
SECTOR DE LA TÉCNICA TECHNICAL 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 coherent detection of forward scattered light.
ESTADO DE LA TÉCNICA STATE OF THE ART
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, sample delivery system installation.
(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 light scattered forward by the sample to compare it with the intensity without sample and from this 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 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) Inverse analysis: determining 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. An inverse 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 calculation data, obtain 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: For modeling the extinction coefficient <t>(<v>) of the sample in the direct analysis step (ii) it is common to use the following expression:
DondeN0es la cantidad total de partículas por unidad de volumen,Des el diámetro de las partículas,ves la frecuencia del campo electromagnético empleado en la medida,f (D )es la DTP normalizada del medio yQext(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, Ds the diameter of the particles, v is the frequency of the electromagnetic field used in the measurement, f (D ) is the normalized PTD 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 the appropriate discretization of the integral, can be expressed in the form of a system of linear equations:
Sin embargo, la obtención de la soluciónfmediante la inversión de la matrizKno 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. Distintos métodos de análisis inverso aplicados a sistemas de partículas se muestran en 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. However, obtaining the solutionf by inverting the matrixK is not possible, since these systems usually present ill-conditioning, i.e., small changes in the values of their elements generate radically different solutions. For the proper solution of this type of problem, it is necessary to apply regularized inversion techniques. Different inverse analysis methods applied to particle systems are shown in 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.
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ámetroD¡según el método de Twomey tiene la siguiente forma: A well-known method for the inversion of particle size distribution from extinction measurements (step iv) is the Twomey method. This belongs to the category of iterative methods which obtain an estimate of the solution after repeated multiplication of the previous solution by a certain factor starting from an initial test solution. In particular, the iterative solution of the distribution for a diameter D¡ according to the Twomey method has the following form:
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(fv).Por otro lado, este cociente está pesado por una función (Wí7-) 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 quotient between the measured extinction (rmeas) and the calculated extinction (<t>^ ) with the Mie model (equation 1) using the solution of the previous iteration (fv). On the other hand, this quotient is weighted by a function (Wí7-) previously calculated from the Mie model (equation 1) and whose function is to regulate the contribution of each particle size proportionally to its contribution in 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. In addition, 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 accuracy in calculating the particle size distribution.
RESUMEN DE LA INVENCIÓN SUMMARY 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 estimating the particle size distribution of a heterogeneous medium from scattered signal intensity information gives the same results as estimating it from phase information and that aggregating the 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 and field modulus information to obtain the particle size distribution more accurately.
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 method and system for estimating particle size distribution by combining information on the modulus and phase of the scattered electromagnetic field. To do this, the detection of the forward scattered signal must be performed with means capable of obtaining information on the modulus 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 FIGURAS BRIEF 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 assist in a better understanding of the characteristics of the invention and to complement this description, the following figures, whose character is illustrative and not limiting, are attached as an integral part thereof:
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 procedure and that obtained by the procedure of the invention.
DESCRIPCIÓN DETALLADA DETAILED 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 method 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; - detecting a portion of forward scattered radiation by means of a detector capable of obtaining the complex optical field, i.e. capable of extracting information from both the modulus 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 such a medium (reference field), i.e. measure the electromagnetic field of the light scattered forward by the sample to compare it with the field without 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 sample and reference field modules. Calculate the phase difference A0 as the subtraction of the sample and reference field phases;
- derivar los parámetros ópticos del medio heterogéneo bajo prueba(nexp, índice de refracción yrexp, 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 (nexp, refractive index and rexp, extinction) from the experimental transmission information and the phase difference calculated above using the following formulas:
dondedes el espesor de la muestra,vla frecuencia ycla velocidad de la luz en el vacío; where d is the thickness of the sample, v is the frequency and c is the speed of light in 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: modelling the optical parameters of the heterogeneous sample, for which the Waterman-Truell theory can be used for media with large amounts of particles, which takes into account the multiple scattering of light, the effect of which is noticeable starting from small volumetric fractions of particles (>1%). This theory allows the calculation of the refractive index and the extinction with the following formula:
dondeRe(-) y¡m(-)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 ynpara una sola partícula.kmed = ^ n medes la constante de propagación del medio matriz (es decir, sin tener en cuenta el efecto de las partículas agregadas) ynmedel índice de refracción de dicho medio. El símbolo ( ) indica un promedio de las amplitudesSpesadas con la DTP del medio; whereRe(-) and¡m(-)refer to taking the real and imaginary parts of the complex number in the parentheses. (S(0/n)) are the Mie amplitudes for angles 0 and nfor a single particle.kmed = ^ n medthe propagation constant of the parent medium (i.e., without taking into account the effect of added particles) and nmethe refractive index of that medium. The symbol ( ) indicates an average of the S amplitudesweighed with the PDT of the medium;
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: Obtaining 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:
dondenexpes el índice de refracción medido experimentalmente yn(3a)lces el índice de refracción calculado usando la ecuación 5 con la DTP de la iteración p.Wfison los pesos calculados mediante la ecuación. 3 yW[¡es el análogo empleando el índice de refracción where exp is the experimentally measured refractive index and n(3a)l is the refractive index calculated using equation 5 with the DTP of iteration p.Wfi are the weights calculated using equation 3 and W[¡ is the analogue using the refractive index
- el parámetroarepresenta 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 eligeacomo aquel que produce un mínimo de la métrica ERM, donde ERM es - the parametera 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 are more similar to the experimental ones. For this purpose,a is chosen as the one that produces a minimum of the ERM metric, where ERM is
- 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 deaque se encuentran próximos a los valores óptimos de forma que puede considerarse un buen indicador para determinarasin 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 (MRE) 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 ofa that are close to the optimal values so that it can be considered a good indicator to determinea without further information than the experimental measurements. In case the value of MRE does not present a very pronounced minimum, it is possible to average the solutions in a range of values of a 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. Figures 2a and 2b show optical parameters of granulated samples made from PTFE powder and glass microspheres. The accumulated volume fraction of the spheres is 5%. Figure 2c shows the results of the DTP inversion 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 (rcaíc,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ámetrosay la posición del mínimo empleado en las inversiones. In Figures 2a and 2b the measured optical parameters (black) as well as those calculated from the light scattering model (rcaíc,ncalc) using the inverted DTP as input can be observed. The results of the inversion using both the traditional Twomey method that only takes into account the extinction (dotted line) and the proposed method combining extinction and refractive index information (solid line) are shown below together with the DTP given by the manufacturer. In Fig. 2d the ERM values obtained with various parameters and the position of the minimum used in the inversions are shown.
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%. The proposed method manages to reduce the deviation from 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 Terahertz Time Domain Spectroscopy instrument that allows to obtain the extinction coefficient and the 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 the information on the magnitude and phase of the scattered electromagnetic field in their corresponding spectral ranges of operation and therefore 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 claimed.
Claims (5)
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| ES202330342A ES2952720B2 (en) | 2023-04-28 | 2023-04-28 | Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical field |
| 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 ES2952720B2 (en) | 2023-04-28 | 2023-04-28 | Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical field |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| ES2952720A1 ES2952720A1 (en) | 2023-11-03 |
| ES2952720B2 true ES2952720B2 (en) | 2024-11-11 |
Family
ID=88558708
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| ES202330342A Active ES2952720B2 (en) | 2023-04-28 | 2023-04-28 | Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical field |
Country Status (2)
| Country | Link |
|---|---|
| ES (1) | ES2952720B2 (en) |
| WO (1) | WO2024223970A1 (en) |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5037202A (en) * | 1990-07-02 | 1991-08-06 | International Business Machines Corporation | Measurement of size and refractive index of particles using the complex forward-scattered electromagnetic field |
-
2023
- 2023-04-28 ES ES202330342A patent/ES2952720B2/en active Active
-
2024
- 2024-04-23 WO PCT/ES2024/070254 patent/WO2024223970A1/en not_active Ceased
Also Published As
| Publication number | Publication date |
|---|---|
| ES2952720A1 (en) | 2023-11-03 |
| WO2024223970A1 (en) | 2024-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nepomnyashchaya et al. | Inverse problem of laser correlation spectroscopy for analysis of polydisperse solutions of nanoparticles | |
| CN105823756B (en) | A kind of metal Terahertz to far infrared complex refractivity index joint inversion method | |
| Vega et al. | Latex particle size distribution by dynamic light scattering: novel data processing for multiangle measurements | |
| Riviere et al. | Multispectral polarized BRDF: design of a highly resolved reflectometer and development of a data inversion method. | |
| Ye et al. | Measurements of particle size distribution based on mie scattering theory and markov chain inversion algorithm. | |
| Miles et al. | Sources of error and uncertainty in the use of cavity ring down spectroscopy to measure aerosol optical properties | |
| KR20150146415A (en) | Dynamic light scattering measurement device and dynamic light scattering measurement method | |
| KR100947031B1 (en) | Apparatus and method for measuring refractive index and thickness of phase object using three wavelength light source | |
| Kassianov et al. | Simultaneous retrieval of effective refractive index and density from size distribution and light-scattering data: Weakly absorbing aerosol | |
| CN118363016A (en) | A method for calculating the boundary value of extinction coefficient based on micro-pulse laser radar | |
| ES2952720B2 (en) | Procedure and system for estimating particle size distribution in heterogeneous media from the complex optical field | |
| Macquart et al. | Understanding the radio variability of Sagittarius A | |
| CN105115940B (en) | Optical material refractive index curve measuring method and device | |
| Ma et al. | The measurement and modeling investigation on the spectral polarized BRDF of brass | |
| CN106290094B (en) | Mie applied to fugitive dust particulate matter on-line monitoring scatters quick calculation method | |
| Levesque et al. | Correction of the calibration measurement by taking into account the Spectralon spectro-polarimetric BRDF model | |
| Mroczka | Method of moments in light scattering data inversion in the particle size distribution function | |
| Wabnitz et al. | Characterization of homogeneous tissue phantoms for performance tests in diffuse optics | |
| Wang et al. | Modified degree of polarization function for rough metallic surface parameter estimation based on multispectral polarimetric measurement | |
| CN110501267B (en) | Method for correcting absorption coefficient of particulate matter measured by T-mode | |
| Liu et al. | Development of the ZJU polarized near-infrared high spectral resolution lidar | |
| Vande Hey | Theory of lidar | |
| US20120116739A1 (en) | Infrared Microspectroscopy for Intact Fibers | |
| Riviere et al. | Multi-band 3D LiDAR performance analysis for classification and segmentation | |
| Obukhov | Sound and light propagation in a weakly inhomogeneous atmosphere |
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 |
|
| FG2A | Definitive protection |
Ref document number: 2952720 Country of ref document: ES Kind code of ref document: B2 Effective date: 20241111 |