COMSOL Uncertainty Quantification Module 6.3

SECTORES

Uncertainty Quantification Module se puede utilizar con productos de toda la suite de productos COMSOL para analizar incertidumbres en simulaciones de ingeniería electromagnética, estructural, acústica, de flujo de fluidos, térmica y química. Puede combinarlo con el módulo de importación de CAD, el módulo de diseño o cualquiera de los productos LiveLink™ para CAD.

CARACTERÍSTICAS

Parámetros de entrada y cantidades de interés

Cuando se arranca un estudio de cuantificación de incertidumbre se define un conjunto de parámetros de interés en términos de una solución del modelo COMSOL®. De esta forma, las cantidades de interés son funciones de los parámetros de entrada.

En el caso de un análisis estructural, las cantidades de interés pueden ser el desplazamiento máximo, la tensión o el ángulo de deflexión. Para un análisis de transferencia de calor o CFD, las cantidades de interés pueden ser la temperatura máxima, la pérdida total de calor o la tasa de flujo total del fluido. Para una simulación electromagnética, pueden ser resistencia, capacitancia o inductancia. Dado que Uncertainty Quantification Module es aplicable a cualquier modelo físico calculado con el software COMSOL Multiphysics®, así como a cualquier expresión matemática de varias magnitudes de campo resueltas, las opciones que pueden ser valor de interés son infinitas.

Los ajustes para los parámetros de entrada y los valores de interés.

Cualquier entrada de modelo incierta, ya sea una configuración física, una dimensión geométrica, una propiedad de material o una configuración de discretización, puede tratarse como un parámetro de entrada, y cualquier salida del modelo puede utilizarse para definir las cantidades de interés.

Screening

El tipo de estudio Screening, MOAT implementa un método de cribado global ligero que proporciona una medida cualitativa de la importancia de cada parámetro de entrada. El método se basa puramente en las muestras, utiliza el método de Morris one-at-a-time (MOAT), y requiere un número relativamente pequeño de evaluaciones del modelo COMSOL. Esto lo convierte en un método ideal cuando el número de parámetros de entrada es demasiado grande para permitir estudios de cuantificación de incertidumbre más costosos desde el punto de vista computacional.

Se realiza un análisis Screening de la media y desviación estándar MOAT para un valor de interés, con siete parámetros de entrada. Un diagrama de dispersión MOAT se muestra como resultado.

Para cada valor de interés, este método MOAT calcula la media MOAT y la desviación estándar MOAT para cada parámetro de entrada. Estos valores se presentan en un diagrama de dispersión MOAT. La clasificación de la media MOAT y las desviaciones estándar MOAT dan la importancia relativa de los parámetros de entrada. Un valor alto de la media MOAT implica que el parámetro influye significativamente en el valor de interés. Un valor alto de la desviación estándar MOAT implica que el parámetro es influyente y que está interactuando fuertemente con otros parámetros o que tiene una influencia no lineal, o ambas cosas.

Análisis de sensibilidad

El tipo de estudio Análisis de sensibilidad/Sesitivity Analysis se utiliza para calcular la sensibilidad de los valores de interés con respecto a los parámetros de entrada. Este tipo de estudio incluye dos métodos: el de Sobol y el de correlación.

El método de Sobol analiza toda la distribución de parámetros de entrada y descompone la varianza de cada cantidad de interés en una suma de contribuciones de los parámetros de entrada y sus interacciones.

Para cada parámetro de entrada, el método de Sobol calcula los índices de Sobol. El índice de Sobol de primer orden muestra la varianza de una cantidad de interés atribuida a la varianza de cada parámetro de entrada individualmente. El índice de Sobol total muestra la varianza de una cantidad de interés atribuida a la varianza de cada parámetro de entrada y su interacción con los otros parámetros de entrada. Los índices de Sobol para cada cantidad de interés y todos los parámetros se presentan en un diagrama de Sobol dedicado donde los histogramas se ordenan por índice de Sobol total. La cantidad de interés es más sensible al parámetro de entrada con el índice de Sobol total más alto. La diferencia entre el índice de Sobol total y el índice de Sobol de primer orden para un parámetro de entrada mide el efecto de la interacción entre esta entrada y otras.

Un gráfico de índice de Sobol para un valor de interés 'q' utilizando 6 parámetros de entrada donde cada par de barras del mismo color representa un parámetro de entrada. La bara izquierda de cada par representa el índice Sobol de primer orden y la de la derecha el índice Sobol total.

En comparación con el método de selección, el análisis de sensibilidad se utiliza para analizar cuantitativamente cómo las incertidumbres en los valores de interés se distribuyen entre los diferentes parámetros de entrada. Este método requiere más recursos computacionales ya que el cálculo de índices de Sobol precisos se basa en un modelo sustituto de alta calidad.

El método de correlación calcula la relación lineal y monótona entre cada parámetro de entrada y los valores de interés. Para el análisis de sensibilidad basado en el método de correlación se calculan cuatro tipos de correlaciones: correlación bivariada, bivariada clasificada, parcial o correlación parcial clasificada.

Propagación de la incertidumbre

El tipo de estudio Propagación de la incertidumbre se utiliza para analizar cómo las incertidumbres de los parámetros de entrada se propagan a cada valor de interés mediante la estimación de su función de densidad de probabilidad (PDF). La física subyacente que asigna los parámetros de entrada a los valores de interés a través de las evaluaciones del modelo de COMSOL Multiphysics® es imposible de calcular analíticamente para la mayoría de las aplicaciones.

Por esta razón, es necesario un análisis de Monte Carlo para aproximar las PDF. Similar al método de Sobol, se utiliza un modelo sustituto para reducir drásticamente el coste computacional del análisis de Monte Carlo. Para cada valor de interés, se realiza una estimación de densidad kernel (KDE) y se visualiza como un gráfico, como una aproximación de la PDF. Además, en base a este análisis, una tabla de intervalos de confianza brinda, para cada valor de interés, la media; Desviación Estándar; mínimo; máximo; y los valores límite inferior y superior correspondientes a niveles de confianza del 90%, 95% y 99%.

Un gráfico de estimación de la densidad kernel para el valor de interés q.

Análisis de fiabilidad

En comparación con otros tipos de estudios de cuantificación de la incertidumbre que investigan la incertidumbre general de los valores de interés, el método Reliability analysis, EGRA aborda una pregunta más directa. Dado un diseño nominal y algunas entradas inciertas específicas, ¿cuál es la probabilidad de que el diseño falle? La falla puede ser una falla completa del diseño, pero también puede expresarse en términos de un criterio de calidad.

Para garantizar la fiabilidad, el enfoque tradicional de modelado y simulación es utilizar márgenes de seguridad y escenarios del peor de los casos. Con un análisis de fiabilidad adecuado es posible evitar la sobreestimación y la subestimación, ya que se pueden realizar estimaciones de la probabilidad real. Se puede obtener una estimación aproximada de la tabla de intervalos de confianza de la propagación de la incertidumbre para cada valor de interés. Pero con el análisis de fiabilidad, puede definirse un criterio de fiabilidad más sofisticado basado en combinaciones de los valores de interés y los umbrales correspondientes. El método de análisis de fiabilidad global eficiente (EGRA) utilizado para el estudio de análisis de fiabilidad dirige eficientemente los recursos computacionales al estado límite que separa la falla y el éxito del diseño.

Líneas de contorno para los dos valores de interés q1 (azul) y q2 (morado) donde los niveles representan los umbrales q1 = 21 y q2 = 9.95 utilizados en el análisis de fiabilidad. Aquí se utilizan dos parámetros de entrada como ilustración (datos de eje x- e y-). La región sombreada muestra los parámetros de entrada que satisfacen el criterio de fiabilidad.

Modelos sustitutos y superficies de respuesta

Los análisis de sensibilidad calculados con el método Sobol, la propagación de la incertidumbre y el análisis de fiabilidad se basan en un análisis preciso tipo Monte Carlo. Esto a menudo requiere una gran cantidad de evaluaciones para lograr una buena precisión. Para problemas realistas en los que una evaluación del modelo de COMSOL Multiphysics® podría requerir recursos significativos y donde el análisis de cuantificación de la incertidumbre implica varios parámetros, un análisis de Monte Carlo realizado solo con evaluaciones del modelo COMSOL Multiphysics® no es factible computacionalmente. Una característica clave del Uncertainty Quantification Module es su capacidad para entrenar y utilizar un modelo sustituto, también conocido como metamodelo, para un análisis UQ particular para ahorrar recursos computacionales.

Un modelo sustituto es un modelo matemático compacto construido para representar y evaluar los valores de interés en el dominio de interés definido por los parámetros de entrada. Este modelo es completamente independiente del modelo COMSOL Multiphysics® subyacente y puede, cuando se entrena adecuadamente, usarse en lugar del modelo COMSOL Multiphysics® para predecir valores para las cantidades de interés para otros valores de los parámetros de entrada distintos a los resueltos. El proceso de construcción de un modelo sustituto suele ser adaptativo y el modelo sustituto puede aproximarse al modelo original con un alto grado de precisión. Las tolerancias definidas por el usuario le permiten aumentar la precisión de los modelos sustitutos. Un nivel más alto de precisión requiere evaluaciones adicionales de modelos de COMSOL Multiphysics®.

Una vez que se ha creado un modelo sustituto puede realizarse una verificación independiente para probar aún más su validez, y puede calcularse rápidamente los datos de la superficie de respuesta para todo el espacio de parámetros de entrada. A continuación, se puede visualizar una superficie de respuesta, donde se traza una cantidad de interés frente a dos parámetros de entrada a la vez.

Un gráfico de superficie de respuesta de un valor de interés q como función de dos parámetros de entrada. Nótese que para problemas con más de dos parámetros de entrada, la superficie de respuesta puede actualizarse fácilmente por otro par de parámetros de entrada en los ajustes del gráfico.

VERSIONES