Ciencia computacional

Estoy interesado en AMR local en mallas no estructuradas. Actualmente, estoy trabajando con la biblioteca OpenFOAM: admite AMR local completamente no estructurado: los criterios de refinamiento de celda determinan una lista de celdas que se cortan las celdas seleccionadas se refinan: se reconstruye toda la malla se crea un mapa …

8 libraries  reference-request  adaptive-mesh-refinement  openfoam 

Más o menos la pregunta. Dada una matriz general dispersa, no simétrica (tanto numérica como estructuralmente), ¿qué importancia tiene el patrón de dispersión (es decir, permutación de fila / columna de matriz / vector) para los solucionadores iterativos? Puedo ver que se vuelve importante para los solucionadores directos (LU) o …

8 linear-algebra  petsc  iterative-method 

Recientemente, vi esta pregunta: cómo medir el error de un método de diferencia finita Soy estudiante de ciencias de la simulación y desafortunadamente, para mí, no está del todo claro qué norma usar en qué contexto. Muy a menudo, usamos la norma euclidiana o la norma L2, pero ¿por qué …

8 error-estimation 

Soy nuevo en programación lineal y he formulado un programa lineal (LP) con un orden de variables y restricciones, aunque la matriz de restricciones es extremadamente escasa. 10 131013101310^{13}1013101310^{13} ¿Quería saber si un LP de esta escala es manejable o no?

8 linear-programming 

Me enfrento con el problema de encontrar el elipsoide ( es una matriz definida positiva simétrica) de volumen máximo dentro de un conjunto convexo dado como un conjunto de desigualdades lineales . Comprendí cómo se formaliza como un problema de optimización convexa como se da en "Convex Optimization, Stephen Boyd …

8 optimization  convex-optimization 

Soy nuevo en el campo de CFD. ¿Cuándo se debe optar por una cuadrícula estructurada y cuándo se debe optar por una no estructurada? (Sí, depende mucho de la geometría del problema) Más específicamente, quiero saber la diferencia en la potencia de cálculo requerida, la precisión lograda y los esfuerzos …

8 fluid-dynamics  discretization 

Quiero resolver la ecuación isotérmica de Lane-Emden [PDF, eq. 15.2.9] d2ψdξ2+2ξdψdξ=e−ψd2ψdξ2+2ξdψdξ=e−ψ\frac{d^2 \!\psi}{d \xi^2} + \frac{2}{\xi} \frac{d \psi}{d \xi} = e^{-\psi} con las condiciones iniciales ψ(ξ=0)=0dψdξ∣∣∣ξ=0=0ψ(ξ=0)=0dψdξ|ξ=0=0\psi(\xi = 0) = 0 \quad \left.\frac{d\psi}{d \xi}\right|_{\xi = 0} = 0 usando SciPyodeint() pero, como se puede ver, la ecuación es singular en el origen. …

8 python  ode 

Los diferentes métodos de álgebra lineal numérica y de optimización numérica tienen diferentes regímenes de tamaño donde son una 'buena idea', además de sus propias propiedades. Por ejemplo, para problemas de optimización muy grandes, se utilizan los métodos de gradiente, gradiente estocástico y descenso de coordenadas en lugar de los …

8 linear-algebra  optimization  nonlinear-programming 

Cualquiera conoce referencias que resumen la historia y las ideas detrás de la supercomputación, incluida la mención de desarrollos en lenguajes de programación paralelos, aplicaciones, empresas de inicio (algunas se establecieron pero no sobrevivieron), necesidades de la industria, etc.

8 parallel-computing  hpc 

Cuando pienso en resolver ecuaciones no lineales, generalmente pienso en linealizar primero, luego aplicar un preacondicionador a la matriz lineal. Se me ocurrió pensar que podría ser posible condicionar las ecuaciones no lineales ANTES de que se linealicen. ¿Existen estrategias generales para precondicionar ecuaciones no lineales que sean independientes de …

8 preconditioning 

cuáles serían buenos métodos (y / o paquetes de software) para intentar resolver un problema minimizando una función cuadrática , st 0 ≤ x i ≤ 1 , y hay más restricciones, algunas de las cuales son no lineales (y no diferenciables), por ejemplo, ∑ i x i 1 x …

8 optimization  nonlinear-programming 

Aprendí que para el solucionador BiCGStab para sistemas lineales dispersos es casi siempre necesario usar un preacondicionador. Ahora me di cuenta de que elegir una buena depende del problema. Navegando por la web descubrí que hay muchos preacondicionadores basados ​​en ILU (ILUT, MILU, etc.) Aquí es donde me confundí. Me …

8 linear-solver  fluid-dynamics  sparse-matrix  preconditioning 

Mi desafío es resolver el siguiente sistema de ecuaciones, que describe la combustión de gas en medios porosos: 1) continuidad ε∂ρg∂t+∂∂x(ρgux)=0ε∂ρg∂t+∂∂x(ρgux)=0\varepsilon \frac{\partial \rho_g}{\partial t} +\frac{\partial}{\partial x} \left(\rho_g u_x\right)=0 2) Ley de Darcy (impulso) ux=−kμ∂p∂xux=−kμ∂p∂xu_x=-\frac{k}{\mu} \frac{\partial p}{\partial x} 3) Ecuación de estado, tenga en cuenta la temperatura variable ρg=MRpRTg(x)ρg=MRpRTg(x)\rho_g=\frac{M_Rp}{RT_g(x)} 4) Ecuación …

8 finite-difference  fluid-dynamics 

Tengo una matriz hermitiana que depende de dos parámetros, digamos x e y . Cuando lo diagonalizo en dos puntos cercanos ( x 1 , y 1 ) y ( x 2 , y 2 ) obtengo dos valores propios cercanos ( ε 1 y ε 2 ) y dos …

8 eigensystem  discretization 

Recibo la advertencia en el asunto de la publicación cuando intento optimizar una función en Python con la función scipy.optimize.fmin_bfgs . La salida completa: Advertencia: el error deseado no se logra necesariamente debido a la pérdida de precisión Current function value: nan Iterations: 1 Function evaluations: 18 Gradient evaluations: 3 …

8 optimization  python