Informe de la semanas 1-19/9/2008

Norberto: estas semanas han sido de trabajo intenso con el reactor de polimerización. Probamos las siguientes configuraciones:

1. Reactor agitado con agua.
2. Reactor agitado con fluido no newtoniano (ley de potencia).
3. Reactor agitado con fluido no newtoniano con dos sistemas de referencia (uno fijo al sistema de laboratorio donde se ubican las entradas y salidas, y otro movil para las paletas).

Los dos primeros funcionaron bien, usamos algunas UDF's, el ultimo no esta andando correctamente. Con Mixing Planes no lo pudimos hacer andar y con Sliding Meshes funciona parcialmente. Los resultados de todo esto Ignacio los uso para armar su presentación para Montreal.
Parte tambien usamos en la reunión con P. Tanguy. De la reunión sacamos en claro dos cosas:

1. Vamos a trabajar en un modelo con mu_0 promedio y ley de potencia y viscosidad promedio tambien. Con este modelo mediremos tiempos de residencia de partículas, una vez obtenida una distribución a la salida aplicaremos los modelos de reactor batch de Diana & Cia.

2. El segundo modelo es igual al anterior en general pero contempla el seguimiento de reactores batch dentro del reactor, los cuales afectarán esencialmente la viscosidad, este modelo será acoplado, probablemente incluyamos el código de MATLAB de Diana como una UDF en Fluent.

Respecto de las materias todo va en orden, este jueves tenemos el primer examen de Algoritmos y Estructuras de Datos con Mario.

Saludos... Santiago.

Informe de la semana 25-29/8/2008

Norberto: respecto de la semana pasada lo más importante tuvo que ver con HIPS. Las materias avanzando, en FEM Alberto termino con la introducción y Victor Fachinotti desarrolló un programar sencillo de FEM en MATLAB, muy interesante.
Respecto de HIPS, Ignacio Gamba y Diana se reunieron con Mario el lunes por la mañana un momento, justo yo no estaba. Simplemente le mostraron una tesis de maestría de una chica que Ignacio conoció en su viaje.
Trata de un caso similar al que estamos manejando, ella modeló un reactor de HIPS en Fluent, en 3D e incluyendo la química. Imprimí y leí la tesis completa. Esta muy buena, luego la discutí con Mario ya que hace uso de marcos de referencia en movimiento.
Luego discutimos algunas cosas con Ignacio y nos armamos una geometría en 3D de un tacho con un agitador como para ir haciendo las primeras pruebas. Ahora lo estamos corriendo en la máquina de el.
Ignacio tiene un quad-core que sin demasiadas complicaciones salió laburando en paralelo con Fluent, asi que estamos haciendo uso de los cuatro cores de la maquina. En estos días voy a tomar tiempos para ver que tal escalea.
La idea ahora es reunirnos con Diana y Mario nuevamente y mostrarles lo que hemos hecho, vamos a aprovechar el laburo de esta chica para agregar un modelo sencillo de quimica y la variación de la viscosidad con la tasa de corte y el peso molecular de lo que se va produciendo (poliestireno).
Asi van las cosas por acá. Saludos...
Santiago.

Informe de la semanas 11-22/8/2008

Norberto: seguimos con las corridas de la cavidad y la esfera. Comencás las materias del segundo cuatrimestre y me inscribí para el cargo de JTP. Este es el detalle.

1) HIPS: Nacho Gamba nos mostró las ecuaciones que tienen cargadas en el modelo hecho en MATLAB, básicamente son las ecuaciones para reactor batch, incluyen la cinética pero responden a un reactor puntual. Yo esperaba ver ecuaciones de transporte de especies y de interacción quimica-fluidodinámica y viceversa pero vamos a tener que esperar. Quedamos con el en ver el modelo en MATLAB funcionando.
2) Las materias comenzaron, la semana pasada con Mario, Algoritmos y Estructuras de Datos y esta semana Elementos Finitos con Alberto Cardona. Todo en orden, voy llevando la lectura al día. La semana proxima nos darán los prácticos.
3) Respecto de la esfera: sigo corriendo con LES, me faltan cinco puntos más. Si bien el modelado de la estela es bastante más realista que con los otros métodos, hasta ahora no veo mejoras en el coeficiente de drag. Una vez que termine las corridas te mando las gráficas de resumen. La cavidad va bien, estuve comparando los graficos del las FIGS. 6-7 del paper que me pasaste y hay cierto acuerdo. Tendrías que decirme cuales son las magnitudes que queremos comparar (de las que se discuten en el paper) o cuales usualemente se comparan. Tenemos tambien los calculos que hizo Mario con PETSC-FEM.
4) JTP: ayer me anoté para el cargo. Totalmente de acuerdo con el tema de la dedicación. Igualmente con el tema de que no es una oportunidad para desaprovechar, no sabemos si se presentará otra. Especialmente por lo que me decía Mario, el hecho de que muchas de las cátedras ya están completas. Veremos la semana que viene si la gente de la FICH me comenta algo.

Asi van las cosas por acá. Muchos Saludos... Santiago.

Informe de la semana 4-8/8/2008

Norberto: la semana pasada estuve trabajando con LES, las corridas son largas, ya estoy obteniendo los primeros resultados, estoy buscando constancia en el valor del CD. Recibí tu mail respecto de las corridas que se deberán hacer a saber:

1.- malla
2.- paso de tiempo
3.- solver : simple, simplec,piso, coupled
4.- modelado de la turbulencia

Una vez que capte bien lo que esta pasando con LES ire haciendo las corridas necesarias.
Respecto del multifásico Damian me había comentado algo y habiamos quedado de hacerlo, asi que en estos nos ponemos de acuerdo y le damos, te comento luego.
Respecto de HIPS, en un momento tendremos una reunión con Diana e Ignacio al respecto, ellos nos mostrarán los modelos químicos que vamos a cargar en Fluent, te iré informando de los resultados.
Con Mario en Calculo Paralelo, me quedó 100 en el segundo parcial, en unos días nos comunican la calificación final del curso.
Saludos... Santiago.

Informe de la semana 28/7-01/8/2008

Norberto: aca estamos probando este nuevo método para informar. De esta manera quedarán todos los informes registrados y ordenados.
El miércoles de la semana pasada rendimos Fluidos, el examen fue resoluble aunque un poco corto el tiempo para los tres ejercicios, el lunes nos comunicarán las notas finales.
Mañana tenemos el repaso final con Mario por el curso de Calculo Paralelo.
Respecto de la modelación de la esfera, como te comenté refiné la malla sin embargo no obtuve resultados mucho mejores, ahora estoy preparando todo para correr con LES según lo que charlamos por chat y con Mario.
Respecto de trámites, estuve hablando con la gente de Tesoreria, resulta que pensaron que yo había OPTADO por Banco Patagonia y por lo tanto nunca me iban a abrir la caja en Francés, ya hable con ellos y con los de becas para aclarar el tema, como siempre fue una falta de comunicación entre ellos. Saludos...
Santiago.

Informe de la semana 21-25/7/08

Norberto: la semana pasada estuve trabajando en los prácticos de Fluidos y la preparación para el examen, que es mañana. Corrí tambien el caso de la esfera, tal como habíamos quedado hice más corridas para el rango de la crisis del drag, esto es en el intervalo de Re: 1e5-1e6. Las corridas con k-e standard, realizable y RNG dieron lo mismo, ni se entera de la crisis. Para correr con LES el caso tiene que ser unsteady así que no pude probar, finalmente probe con lo que según el manual es el modelo más elaborado de turbulencia en Fluent, el Reynolds Stress Method (RSM), con este se ve la crisis de drag y mediante un postprocesador se observa el cambio en la posición del punto de separación. Ahora tengo que pulir algunas cosas ya que los valores de drag no dan similares y el valor del Reynolds para el punto de crisis no es del todo el esperado.
Estuve hablando el tema con Mario y quedamos en que extendería el dominio un poco mas hacia los lados para evitar problemas de reflexión por la condición de contorno slip, tambien voy a refinar la malla un poco más en los alrededores de la esfera para modelar un poco mejor la capa limite. Ni bien rindamos, etc. me pongo con eso. Saludos...

Informe de la semana 7-11/7/2008

Norberto: tal como te comenté por chat la semana pasada estuve centrado en terminar los prácticos de Mario y prepararme para el segundo parcial. Tuve que venirme el 9 de julio pero quedo todo entregado.
Te paso un link donde te podes bajar el informe del primer tutorial de Fluent. Preferí profundizar en este tutorial antes que avanzar mucho en los demas sin concretar, además para los tutorials posteriores me esta faltando bastante más teoría. Esta hecho el tutorial completo, ampliado respecto de lo que aparece en la explicación original, le agregué además algunos tips de teoría relacionada, que tiene que ver con las decisiones que toman a la hora de elegir las opciones de los solvers, etc.
Igualmente estuve viendo otros tutorials posteriores para poder volcar lo que esta en el informe. Ahora tengo un panorama del uso del soft, los menúes, las herramientas de postproceso, los solvers, condiciones de contorno, etc. Tendría que mirar algo de Gambit como para manejar la generación de geometrias y mallas.
Respecto del reactor de butenos, debido a que el tiempo es escaso y para tratar de hacer la comparación en las condiciones más parecidas posibles decidimos con Damián importar el reactor a Fluent y correrlo para comparar, sin embargo la opción de importación de archivos CFX no anda y la gente del soporte después de pedirme que haga algunas pruebas no me ha respondido nada más, me parece que se les quemaron los papeles.
Seguiré con los prácticos que quedan de fluidos, cuando tengas los tests para correr en Fluent avisame nomás, le pondremos todo el esfuerzo. Acá va el link:

http://santiagomarquezd.googlepages.com/doctorado

El trabajo es: 1-Introducción: Flujo turbulento y ecuación del calor

te agradecería que en la medida que puedas me hagas un feedback bien a fondo así sacamos lo máximo del ejemplo.

Saludos... Santiago.