Ciencia computacional

Soy relativamente nuevo en el modelado de CFD y estoy haciendo un modelo VOF de una caja rectangular ( m) con una entrada ( ) en el lado largo.l×w×h=120×80×8l×w×h=120×80×8l \times w \times h = 120\times 80 \times8w×h=10×8w×h=10×8w\times h = 10\times 8 La salida de presión es la cara superior de …

8 fluid-dynamics  mesh-generation 

La ecuación de advección necesita ser discretizada para ser utilizada para el método de Crank-Nicolson. ¿Alguien puede mostrarme cómo hacer eso?

8 pde  advection  crank-nicolson 

En mis cursos sobre análisis numérico, me enseñaron que la principal y principal motivación para preacondicionar sistemas lineales de ecuaciones es aumentar la tasa de convergencia de los solucionadores iterativos para ese LSE. Pero, ¿hay algún efecto en la precisión de la solución calculada? Puedo recordar un resultado sobre la …

8 linear-algebra  condition-number 

Según Wolfram Alpha y el sistema de álgebra computacional Sage , la siguiente identidad se mantiene: cos(arctan(l1−l2d))=11+(l1−l2)2d2−−−−−−−−√cos⁡(arctan⁡(l1−l2d))=11+(l1−l2)2d2 \cos\left(\arctan\left(\frac{l_1-l_2}{d}\right)\right) = \frac{1}{\sqrt{1 + \frac{(l_1-l_2)^2}{d^2}}} Sin embargo, cuando traté de verificarlo con un ejemplo arbitrario en NumPy, noté una diferencia bastante grande en los valores reales calculados por ambos lados de la identidad. …

8 python  numerics 

Supongamos que queremos resolver una ley de conservación hiperbólica . Realmente me gusta usar Lax-Wendroff, que dicetut+ f( u )X= 0ut+f(u)x=0u_t+f(u)_x=0 tun + 1j= unortej- Δ tΔ x( g( unortej + 1, unortej) - g( unortej, unortej - 1) )ujn+1=ujn−ΔtΔx(g(uj+1n,ujn)−g(ujn,uj−1n))u_j^{n+1} = u_j^n -\frac{\Delta t}{\Delta x}(g(u_{j+1}^n,u_j^n)-g(u_j^n,u_{j-1}^n)) dónde sol( v , w …

8 fluid-dynamics  hyperbolic-pde 

En "Factorizaciones LU, QR y Cholesky usando capacidades vectoriales de GPU", de Vasily Volkov y James Demmel , hay una forma interesante de interpretar las latencias, tamaños de línea y tamaños de página de un caché desde un gráfico como este: Aparentemente, el tamaño de la línea es donde el …

8 performance  gpu  benchmarking 

Tengo un código de modelado basado en el filtro Kalman que he desarrollado para una aplicación de mapeo ionosférico regional casi en tiempo real. El código asimila datos de diferentes sensores en un mapa (descrito por un conjunto de funciones básicas) utilizando un filtro Kalman. Estoy tratando de escalar esto …

8 linear-algebra  sparse-matrix 

Dada la ecuación de Poisson 3D y el lado derecho y el dominio, soy libre de imponer cualquier condición de límite (BC) en la función , ¿o tienen que ser de alguna manera consistentes con el lado derecho? En particular, si impongo BC periódica, ¿habrá exactamente una solución para cualquier …

8 pde  boundary-conditions  elliptic-pde 

No estoy realmente seguro de cómo explicar este problema claramente, así que tengan paciencia conmigo. Tengo una base de 3 vectores de unidades ortonormales y una posición, una matriz de transformación 4x4 estándar en gráficos de computadora. También tengo varios puntos (compensaciones) en ese espacio que me transformo en espacio …

8 optimization  regression  basis-set  least-squares 

¿Existe un algoritmo para reorganizar una matriz en forma de diagonal de bloque, dado que la matriz es de naturaleza diagonal de bloque pero aleatoria con una elección de base imprudente? En particular, ¿hay algún módulo de Python para esto?

8 matrix  python 

¿Cuál es una buena manera de extraer la diagonal de una matriz simétrica que ya es casi diagonal, pero donde no tiene los elementos de la matriz (solo la capacidad de aplicarla a los vectores)? Otras restricciones son, (1) aplicar la matriz n-por-n n veces para construir explícitamente la diagonal …

8 linear-algebra  preconditioning 

Estoy buscando encontrar una solución al siguiente problema, pero tengo problemas para formularlo con sensatez y luego encontrar un algoritmo adecuado para resolverlo. Considere una lista de artículos colocados en una bolsa de compras: 1, 2, 3, 4 ... Cada artículo puede ser parte de una o más promociones: A, …

8 optimization  algorithms 

Hola, si alguien puede ofrecer información sobre cómo resolver mi problema, ¡sería genial! Estoy buscando calcular el área de intersección entre un volumen 3D y un plano 2D. Volumen 3D: definido por 6 puntos (siempre será una cuña 3D que se asemeja a un triángulo grueso) Plano 2D: definido por …

8 computational-geometry 

Estoy leyendo un libro sobre métodos numéricos y el cuadrado de la norma discreta L2L2L^2 se define como ||x||22=h∑1Nx2i||x||22=h∑1Nxi2||x||^2_2=h\sum_1^Nx^2_i Cada punto obtiene un "peso", que es hhh , por lo tanto, es como un promedio sobre los cuadrados de los valores en todos los puntos. De hecho, esto proviene de …

8 pde 

Necesito aprender y utilizar un paquete de solución automatizada de volumen finito para un proyecto en el que estoy trabajando y lo he reducido a estos dos paquetes. Me preguntaba si alguien tiene experiencia con ambos paquetes y quizás podría comentar una comparación. Considerando el objetivo similar de ambos paquetes, …

8 openfoam  finite-volume