FACULTAD DE INGENIERÍA QUÍMICA
Temas Selectos
de Ingeniería
Química:
Agitación y CFD
Manual de Prácticas
Flujo Turbulento en una Tubería
Agitación y CFD
Csrlos Rubio Atoche
Flujo Turbulento en una Tubería
Especificación del Problema
Figura 1 Diagrama del problema que se va a resolver
Considere un fluido que circula a través de una tubería de sección circular, de radio constante,
como se muestra en la figura 1. El diámetro (D) de la tubería es de 0.2 metros y la longitud
(L) de 8 metros. Considere que la velocidad de entrada (Uz) es constante e igual a 1 m/s en
toda la sección transversal de la tubería. La presión a la salida de la tubería es de 1 atm.
Considerando que se desea simular un flujo laminar, se ha fijado un valor del número de
Reynolds igual a 100. Para lograr esto se considerará la densidad del fluido ρ = 1 kg / m3 y
el coeficiente de viscosidad µ = 2 x 10 -5 kg / (m·s). Evidentemente esto no corresponde a
ningún fluido real, pero para fines didácticos es muy adecuado, ya que nos permite establecer
la dependencia entre el régimen de flujo y el perfil de velocidades establecido. De igual modo
vamos a observar la dependencia de los resultados obtenidos con la calidad de la malla
empleada para la simulación.
Para lograr esto, se resolverá el problema planteado utilizando FLUENT en ANSYS.
Como evidencia deben presentarse los siguientes resultados.
Vectores de velocidad
Contorno de la magnitud de velocidad
Contornos de presión
Perfil de velocidad en diferentes puntos, incluyendo a la salida
Variación de la velocidad a lo largo del eje del tubo
El coeficiente de fricción de pared (skin friction) a lo largo de la pared del tubo.
Una vez obtenidos estos resultados es fundamental compararlos con la solución
analítica de flujo completamente desarrollado, para validar los resultados obtenidos.
Flujo Turbulento
El flujo turbulento exhibe fluctuaciones a pequeña escala en el tiempo. Por lo general,
no es posible resolver estas fluctuaciones en un cálculo de CFD. Por lo tanto, las variables
de flujo como la velocidad, la presión, etc. tienen un promedio de tiempo.
Desafortunadamente, las ecuaciones de gobierno promediadas en el tiempo no están cerradas.
(es decir, contienen cantidades fluctuantes que deben modelarse utilizando un modelo de
Agitación y CFD
Csrlos Rubio Atoche
turbulencia). Actualmente no hay un modelo de turbulencia disponible que sea válido para
todos los tipos de flujos, por lo que es necesario elegir y ajustar un modelo en función de las
clases particulares de flujos.
Este ejercicio, se realizará con las variantes del modelo k-ε. Pero en el mundo real,
ande con mucha precaución: debe evaluar la validez de sus cálculos utilizando un modelo de
turbulencia con mucho cuidado (lo cual, ejem, significa que no hay forma de evitar estudiar
los conceptos de dinámica de fluidos y métodos numéricos con mucho cuidado). FLUENT
no debe usarse como una caja negra. Los modelos k-ε consisten en dos ecuaciones
diferenciales: una para la energía cinética turbulenta k y la disipación turbulenta ε. Estas dos
ecuaciones deben resolverse junto con las ecuaciones de continuidad, cantidad de momento
y energía promediadas en el tiempo. Por lo tanto, los cálculos de flujo turbulento son mucho
más difíciles y lentos que los cálculos de flujo laminar. Este es un ejercicio para despertar el
interés por los cálculos de flujo turbulento..
𝝏(𝒗𝒛 )
=𝟎
𝝏𝒓
Figura 2. Diagrama del problema que se va a resolver
Geometría
Para diseñar la geometría es preciso seleccionar el software sobre el que se llevará a
cabo la simulación.
Selección de proyecto de flujo de fluido (FLUENT)
En el lado izquierdo de la ventana del banco de trabajo, verá una caja de herramientas
llena de varios sistemas de análisis. A la derecha, verá un espacio de trabajo vacío. Este es el
lugar donde organizará su proyecto. En la parte inferior de la ventana, verá mensajes de
ANSYS.
Agitación y CFD
Csrlos Rubio Atoche
Tipo de análisis
En la pestaña Project Schematic de la ventana del Workbench, haga clic derecho
en Geometry y seleccione Properties.
El menú de propiedades aparecerá a la derecha de la ventana de Workbench. En
Advance Geometry Options, cambie en Analysis Type a un análisis 2D.
Iniciar Design Modeler
En la opción Project Schematic, haga doble clic en Geometry para comenzar
a construir la geometría.
En este punto, se abrirá una nueva ventana, ANSYS Design Modeler. Se le pedirá que
seleccione la unidad de longitud de trabajo. Use la unidad de medida predeterminada (metros)
y haga clic en Aceptar OK.
Definir el boceto de la geometría.
Comience creando un boceto en el plano XY. En Tree Outline, seleccione
XYPlane, luego haga clic en Sketching justo antes de Details View. Esto abrirá las
cajas de herramientas de dibujo Sketching Toolboxes.
Haga clic sobre el eje +Z en la esquina inferior derecha de la ventana Graphics para
presentar el plano XY.
En las cajas de herramientas Dibujo (Sketching), seleccione Rectangle. En la
ventana Gráficos, cree un rectángulo aproximado haciendo clic una vez en el origen y luego
haciendo clic una vez en algún lugar del plano XY positivo. (Asegúrese de ver una letra P
en el origen antes de hacer clic. La P implica que el cursor está directamente sobre el punto
de intersección de los ejes coordenados).
Figura 3. Rectángulo diseñado para generar la geometría de trabajo.
Agitación y CFD
Csrlos Rubio Atoche
Dimensiones
En este punto, el rectángulo se podrá dimensionar correctamente. En Sketching de
las cajas de herramientas, seleccione la pestaña Dimensions, use las herramientas de
dimensionamiento predeterminadas. Acote la geometría como se muestra en la siguiente
imagen.
Figura 4. Rectángulo con sus dimensiones
Debajo de la tabla Details View (ubicada en la esquina inferior izquierda),
establezca la altura del rectángulo V1 = 0.1 m y H2 = 8 m.
Creación de la superficie del cuerpo.
Para crear la superficie del cuerpo, de primero (Haga clic) Concepto> Superficie a
partir de bocetos
Esto creará una nueva superficie SurfaceSK1. En la pestaña vista de detalles
(Details View), seleccione el boceto1 (Sketch1) como objeto base (Base Objects)
seleccionando una de las líneas del boceto y haciendo clic en Aplicar (Apply). Luego
seleccione el grosor para que sea 0.1m y haga clic en el botón Generate para generar la
superficie.
Guardar (Saving)
En este punto, puede cerrar el Design Modeler y volver a la página del proyecto
en Workbench.
Guarde el proyecto haciendo clic en el botón "Guardar como ..." , que se
encuentra en la parte superior de la página del proyecto Workbench. Guarde el proyecto como
"Tubo Laminar" en su directorio de trabajo. Cuando guarde en ANSYS, se creará un archivo
y una carpeta. Por ejemplo, si guarda como "Tubo Laminar", aparecerá un archivo "Tubo
Laminar" y una carpeta llamada "Tubo Laminar_files". Para volver a abrir los archivos
ANSYS en el futuro, necesitará tanto el archivo ".wbpj" como la carpeta. Si no tiene
AMBOS, no podrá acceder a su proyecto.
Agitación y CFD
Csrlos Rubio Atoche
Mallado
En esta sección se generará la malla necesaria para realizar la simulación que nos
ocupa. La geometría estará mallada con 500 elementos. Es decir, la tubería se dividirá en 100
elementos en la dirección axial y 5 elementos en la dirección radial.
Iniciar el mallado (Mesher)
Para comenzar el proceso de mallado, vaya a la página del proyecto de Workbench,
luego haga un doble clic en el apartado Malla (Mesh).
Malla predeterminada
En esta sección se generará la malla predeterminada. Esto se puede llevar a cabo de
dos maneras. La primera es (clic con el botón derecho) Malla> Generar malla ((Right
Click) Mesh > Generate Mesh).
La segunda manera en la que se puede generar la malla de forma predeterminada es
haciendo clic en Malla> Generar malla (Mesh > Generate Mesh), en los botones
debajo del menú de comandos.
Cualquiera de los métodos deberá dar los mismos resultados. La malla predeterminada
que se genere debe ser semejante a la siguiente imagen.
Figura 5. Rectángulo con la malla predeterminada.
Tenga en cuenta que en Workbench generalmente hay al menos dos formas de
implementar acciones como se ha mostrado anteriormente. Por simplicidad, el método de
"menú" para implementar acciones se utilizará como referencia en el resto de las indicaciones
para la práctica.
Malla de cara mapeada
Como se puede ver en la figura 5, la malla predeterminada tiene elementos irregulares.
Existe un particular interés en crear un estilo de malla de cuadrícula que se pueda asignar a
un dominio rectangular. Este estilo de mallado se denomina Mallado de caras mapeadas
(Mapped Face Meshing). Para incorporar este estilo de malla (Hacer clic) Control de
malla> Malla de cara mapeada ((Click) Mesh Control > Mapped Face Meshing )
en el botón Mesh control.
Ahora, la malla de cara mapeada (Mapped Face Meshing) debe aplicarse a la
geometría de la tubería. Para hacerlo, primero haga clic en el cuerpo de la tubería, esta
Agitación y CFD
Csrlos Rubio Atoche
operación debe resaltarlo en verde. A continuación, (Haga clic) Aplicar (Apply) en la tabla
Detalles de Malla de cara mapeada (Details of Mapped Face Meshing ) a la
izquierda normalmente.
Ahora, genere la malla utilizando cualquiera de los métodos de la sección "Malla
predeterminada" anterior. Debe obtener por resultado una malla comparable a la siguiente
imagen.
Figura 6. Rectángulo con la malla de cara mapeada predeterminada.
Tamaño de los bordes (número de divisiones)
La malla deseada tiene un número específico de divisiones (Edge Sizing) a lo largo
de la dirección radial y axial (horizontal y vertical respectivamente). Para obtener la cantidad
requerida de divisiones, se debe usar Edge Sizing. Las divisiones a lo largo de la dirección
axial se especificarán primero. Ahora, se debe insertar un tamaño de borde. Primero, haga
clic en Control de malla>Tamaño (Mesh Control > Sizing).
Ahora, se debe especificar la geometría, sobre la cual se aplicará la instrucción, y el
número de divisiones requeridas. Para esto, primero haga clic en Filtro de selección de
bordes, Edge Selection Filter, . Luego mantenga presionado el botón "Control" y
haga clic en el borde inferior y superior del rectángulo para seleccionarlos. Ambos lados
deben resaltar en verde. A continuación, presione Aplicar (Apply) debajo de la tabla
Detalles de tamaño (Details of Sizing).
Ahora, cambie la forma de establecer el tamaño (Type) a Número de divisiones, en el
menú Detalles de tamaño (Details of Sizing). Establezca Número de divisiones
(Number of Divisions) en 100 como se muestra a continuación.
Siga el mismo procedimiento que para el tamaño del borde en la dirección radial,
excepto que seleccione el lado izquierdo y derecho en lugar de la parte superior e inferior y
establezca el Número de divisiones (Number of Divisions) en 30. Luego, genere la
malla utilizando cualquiera de los métodos de la "Malla predeterminada "sección anterior.
Debe obtener la siguiente malla.
Figura 7. Rectángulo con la malla resultado de especificar bordes.
También vamos a necesitar sesgarlo. Esto se debe a que queremos divisiones más
pequeñas cuanto más nos acerquemos a la pared..
Agitación y CFD
Csrlos Rubio Atoche
A continuación, aplicaremos un tamaño de borde con sesgo a la entrada, el extremo
izquierdo de la tubería. Haga clic en Mesh Control > Sizing. Con la herramienta de
selección de bordes, resalte la entrada (extremo izquierdo) de la tubería y haga clic en Apply
junto a Geometría. Como en el tutorial de Flujo de tubería laminar, cambie el Tipo a Número
de divisiones e ingrese 30. Cambie el comportamiento Behavior a Hard. Ahora,
apliquemos un sesgo al tamaño del borde. En Bias Type, seleccione la segunda opción, - -
- ----. Ingrese un Bias Factor de 10.
Ahora nos gustaría aplicar un tamaño de borde a la salida, el extremo derecho de la
tubería. Una vez más, utilizaremos 30 divisiones, con un factor de sesgo de 10 y con las
divisiones más pequeñas en la parte superior, cerca de la pared. Esta vez, al seleccionar Tipo
de sesgo, elija la primera opción, ---- - - -. Esto pondrá las divisiones más pequeñas en la
parte superior. Aparte de esto, el procedimiento es el mismo que para la entrada.
Haga clic derecho en Malla y seleccione "Generar malla". El factor de polarización
genera una malla más fina cerca de la pared de la tubería. Esto se hace para calcular la
pequeña fluctuación en la propiedad del fluido cerca de la pared.
Como resultado, en la malla de arriba hay 540 elementos, cuando debería haber solo
500. Las estadísticas de malla se pueden encontrar haciendo clic en Malla Mesh en el menú
Árbol y luego expandiendo Estadísticas Statistics en la tabla Detalles de Malla Details
of Mesh. Para obtener la malla deseada de 500 elementos, el comportamiento de la malla
debe cambiarse de suave a duro (Soft Hard) para ambos tamaños de borde (Edge
Sizing's). Para llevar a cabo esta operación, primero, expanda Malla en el menú de árbol,
luego haga clic en Tamaño del borde Edge Sizing y luego cambie el “comportamiento”
Behavior a Duro Hard en la tabla de la izquierda en Detalles del tamaño del borde
Details of Edge Sizing.
Posteriormente, establezca el “comportamiento” Behavior a Duro Hard para el
Tamaño de borde 2 Edge Sizing 2. A continuación, genere la malla utilizando cualquiera
de los métodos de la sección "Malla predeterminada" anterior. Luego debe obtener la
siguiente malla de 3000 elementos.
Agitación y CFD
Csrlos Rubio Atoche
Figura 8. Rectángulo con la malla de 500 elementos, resultado de especificar bordes y
el comportamiento de estos.
Nombrado de los bordes
Aquí, los bordes de la geometría recibirán nombres para que se puedan asignar
posteriormente las condiciones frontera en Fluent. El lado izquierdo de la tubería se llamará
"Entrada o Inlet" y el lado derecho se llamará "Salida o Outlet". El lado superior del
rectángulo se llamará "Pared o Pipe Wall" y el lado inferior del rectángulo se llamará "Eje o
Ceterline" como se muestra en la imagen a continuación. Los nombres en inglés facilitan la
asignación de las condiciones frontera.
Figura 9. Rectángulo con los nombres en inglés de cada frontera del sistema
Para crear una selección con nombre primero (Haga clic) en el filtro de selección de bordes,
Edge Selection Filter, . Luego haga clic en el lado izquierdo del rectángulo y
debería resaltarse el borde de verde. Luego, haga clic derecho en el lado izquierdo del
rectángulo y elija Crear selección con nombre Create Named Selection.
Seleccione el borde izquierdo y haga clic derecho y seleccione Crear selección con nombre
Create Named Selection. Teclee entrada o Inlet y haga clic en OK.
Ahora, cree selecciones con nombre para los tres lados restantes y asígneles un nombre de
acuerdo al diagrama.
Guardar, salir y actualizar.
Primero guarde el proyecto. A continuación, cierre la ventana Mesher. Luego, vaya a
la página del proyecto Workbench Workbench Project Page y haga clic en el botón
Actualizar proyecto Update Project .
Configuración del modelado (Physics Setup)
En este punto deben de tener marcas de verificación a la derecha de las celdas
Geometría y Malla. Como siguiente paso, los datos de la malla y la geometría deben leerse
en FLUENT. Para leer los datos (clic con el botón derecho) Configuración> Actualizar
(Right Click) Setup > Refresh en la página del proyecto de Workbench
Workbench Project Page. Si la opción de actualización no está disponible,
simplemente omita este paso.
Agitación y CFD
Csrlos Rubio Atoche
Después de hacer clic en Actualizar, debe aparecer un signo de interrogación a la
derecha de la celda Configuración. Esto indica que el proceso de configuración aún no se ha
completado.
Inicio de Fluent
Haga doble clic en setup (double Click) Setup en la página del proyecto
Workbench Workbench Project Page, que abrirá FLUENT. Cuando aparezca
FLUENT Launcher FLUENT Launcher, cambie las opciones a "Doble precisión" y
luego haga clic en Aceptar. La opción Doble precisión se utiliza para seleccionar el
solucionador de doble precisión. En el solucionador de doble precisión, cada número de
punto flotante se representa usando 64 bits en contraste con el solucionador de precisión
simple que usa 32 bits. Los bits adicionales aumentan no solo la precisión, sino también el
rango de magnitudes que se pueden representar. La desventaja de usar doble precisión es que
requiere más memoria.
Espere un poco mientras se inicia la interfaz de FLUENT. Es en este punto donde
especificaremos las ecuaciones de gobierno y las condiciones frontera para nuestro problema.
En el lado izquierdo de la interfaz FLUENTE, vemos varios elementos enumerados en
Configuración del problema. Trabajaremos de arriba a abajo de los elementos de
Configuración del problema para configurar la física de nuestro problema. En el lado
derecho, tenemos el panel de Gráficos y, debajo de eso, el panel de Comando.
Comprobar y mostrar la malla
Primero, se comprobará la malla para verificar que se haya importado correctamente
desde Workbench. Para obtener las estadísticas sobre la malla (clic) Malla> Información>
Tamaño (Click) Mesh > Info > Size.
La malla que se creó anteriormente tiene 3000 elementos (30 radiales x 100 axiales).
Tenga en cuenta que los elementos en FLUENT se llaman celdas. El resultado indica que
hay 3000 celdas, lo cual es una buena señal. A continuación, se le pedirá a FLUENT que
verifique la malla en busca de errores. Para llevar a cabo el procedimiento de verificación de
malla (Hacer clic) Malla> Verificar.
No debería ver errores en el Panel de comandos Command Pane. Ahora que se
ha verificado la malla, se discutirán las opciones de visualización de la malla. Para que
aparezcan las opciones de visualización (Haga clic) General> Malla> Mostrar (Click)
General > Mesh > Display.
El paso anterior debería hacer que se abra la ventana para poder observar la malla.
Tenga en cuenta que ahora aparecen los lados con los nombres definidos en los pasos de
mallado.
Agitación y CFD
Csrlos Rubio Atoche
Debe tener todas las superficies que se muestran en la instantánea anterior. Al hacer
clic en el nombre de una superficie en el menú mostrar malla Mesh Display, se alternará
entre seleccionar y deseleccionar. Al hacer clic en Mostrar, se mostrarán todas las entidades
de superficie actualmente seleccionadas en el panel de gráficos. Anule la selección de todas
las superficies y luego seleccione cada una a la vez para ver a qué parte del dominio o límite
corresponde la entidad de superficie particular (necesitará acercar / alejar y trasladar el
modelo a medida que lo hace).
Asegúrese de que estén seleccionados los 5 elementos en Superficies Surfaces. El
botón al lado de Superficies selecciona todos los límites, mientras que el botón
deselecciona todos los límites a la vez. Una vez que se han seleccionado los 5 límites, haga
clic en mostrar Display, luego cierre la ventana mostrar malla Mesh Display. El
rectángulo largo y delgado que se muestra en la ventana gráfica corresponde a nuestro
dominio de solución. Algunas de las operaciones disponibles en la ventana de gráficos para
interrogar la geometría y la malla son:
Traslación: El modelo se puede trasladar en cualquier dirección manteniendo
presionado el botón izquierdo del ratón Left Mouse Button y luego moviendo el ratón
en la dirección deseada.
Acercar (Zoom in): Mantenga presionado el botón central Middle Mouse
Button del ratón y arrastre un cuadro desde la esquina superior izquierda Upper Left
Hand Corner a la esquina inferior derecha Lower Right Hand Corner sobre el área
que desea ampliar.
Alejar (Zoom out): Mantenga presionado el botón central del ratón Middle
Mouse Button y arrastre un cuadro desde la esquina inferior derecha Lower Right
Hand Corner hasta la esquina superior izquierda Upper Left Hand Corner.
Use estas operaciones para acercar e interrogar la malla.
Definir propiedades del solucionador (solver)
En esta sección, se especificarán las diversas propiedades del solucionador para
obtener la solución adecuada para el flujo laminar en la tubería. Primero, se debe especificar
la naturaleza aximétrica del sistema, en la geometría. En General > Solver > 2D
Space, seleccione Axisymmetric.
A continuación, se especificarán los parámetros del modelo viscoso Viscous
Model que deberá considerarse para simular el comportamiento del fluido. Para abrir las
Opciones de modelo viscoso Models > Viscous - Laminar > Edit... De forma
predeterminada, las opciones del Modelo viscoso están configuradas en laminar, por lo que
deberá cambiarse el modelo viscoso para resolver el problema.
Agitación y CFD
Csrlos Rubio Atoche
Haga clic en Editar ... y elija el modelo k-epsilon (2eqn). Observe que la ventana
se expande y se muestran opciones adicionales al elegir el modelo de turbulencia k-epsilon.
En el apartado Near-Wall Treatment, elija Enhanced Wall Treatment. Esta
opción utiliza una función combinada para ir entre un modelo de dos capas y las funciones
de muro estándar. Si la malla cerca de la pared es lo suficientemente fina, se usa el modelo
de dos capas. De lo contrario, se utilizan funciones de pared estándar. Alternativamente,
puede utilizar las funciones de pared estándar Standard Wall Functions; esto
funcionará bien cuando 30 < y + < 100.
Acto seguido, se especificarán los parámetros del modelo de energía. Para abrir las
Opciones de modelo de energía Models > Energy-Off > Edit... Para un flujo
incompresible, la ecuación de energía se desacopla de las ecuaciones de continuidad y
momento. Necesitamos resolver la ecuación de energía solo si estamos interesados en
determinar la distribución de temperatura. No trataremos con la temperatura en este ejemplo.
Deje la ecuación de energía Energy Equation configurada en apagado y haga clic en
Cancelar para salir del menú.
Definir las propiedades de los materiales
Ahora, se especificarán las propiedades del fluido que se está modelando. Las
propiedades del fluido se definieron en la sección Especificación del Problema. Para crear
un nuevo fluido (Click) Materials > Fluid > Create/Edit ...
En el menú Create/Edit Materials, establezca la densidad en 1 kg / m3
(constante) y la viscosidad en 2e-5 kg / (ms) (constante).
Haga clic en Change/Create y cierre la ventana.
Definir las condiciones frontera
En este punto, se especificarán las condiciones frontera para las cuatro bordes de la
geometría.
Presión de Operación
Configuración del problema> Condiciones frontera> Condiciones de funcionamiento
... Problem Setup > Boundary conditions > Operating Conditions... Recuerde que
para todos los flujos, FLUENT utiliza internamente la presión manométrica. Cada vez que
se necesita una presión absoluta, se genera agregando la presión de funcionamiento a la
presión manométrica. Usaremos el valor predeterminado de 1 atm (101,325 Pa) como la
presión de funcionamiento Operating Pressure.
Haga clic en Cancel para dejar el valor predeterminado en su lugar.
La condición límite de la frontera para la entrada se especificarán a continuación.
Condición límite de entrada (INLET)
Agitación y CFD
Csrlos Rubio Atoche
Para comenzar el proceso (haga clic) Boundary Conditions > inlet > Edit....
Tenga en cuenta que el tipo de condición límite Boundary Condition Type
debería haberse configurado automáticamente en velocity-inlet si el borde fue nombrado
como INLET, en caso contrario estará configurado como pared WALL. Ahora, se
especificará la velocidad en la entrada. En el menú de velocidad de entrada, configure el
Método de especificación de velocidad Velocity Specification Method en
Componentes Components, y establezca la velocidad axial Axial-Velocity (m/s)
(m/s) en 1 m/s. Luego, haga clic en Aceptar para cerrar el menú de especificación de la
velocidad de entrada. Seleccione Intensity and Hydraulic Diameter junto a
Turbulence Specification Method posteriormente indique 1% para Turbulence
Intensity y 0.2 m para Hydraulic Diameter, Finalmente presione sobre OK para fijar
las condiciones frontera de la entrada.
Condición límite de salida (OUTLET)
Primero, seleccione la salida en el menú condiciones frontera Boundary
Conditions. Como se puede apreciar, la condición Type debe haberse configurado
automáticamente en presión de salida pressure-outlet. Si la condición no está
configurada en presión de salida, configúrelo bajo esta condición. Presione sobre Edit. El
valor por defecto de la presión manométrica Gauge Pressure es 0. La presión (absoluta)
en la salida es de 1 atm. Dado que la presión de funcionamiento se establece en 1 atm, la
presión manométrica de salida = presión absoluta de salida - presión de funcionamiento = 0.
Elija la salida en Zona. El tipo de este límite es la salida de presión. Haz clic en Editar. El
valor predeterminado de la presión manométrica es 0. Haga clic en Cancelar para dejar los
valores predeterminados en su lugar. Bajo las condiciones de la salida,
Condición límite de la línea central (AXIS)
Seleccione la línea central en el menú condiciones frontera. Como se puede apreciar,
la condición Type se ha configurado automáticamente a pared Wall, lo cual no representa
nuestro caso, por lo que no es correcto. Cambie la condición a Axis.
Cuando aparezcan los cuadros de diálogo, haga clic en Sí Yes para cambiar el tipo
de límite. Luego haga clic en OK para aceptar eje o "centerline" como el nombre de la zona.
Condición límite de la pared de la tubería
Primero, seleccione pared o Pipe_wall en el menú condiciones frontera Boundary
Conditions.
Como podrá apreciar, la condición Type debería haberse configurado
automáticamente en pared wall. Si la condición no está establecida en pared, configúrelo en
pared. Ahora, no se necesitan más cambios para la condición límite de la pared del tubo
(Pared o pipe_wall).
Guardar
Agitación y CFD
Csrlos Rubio Atoche
Para guardar su trabajo haga clic en (Click)File > Save.
Solución numérica
Esquema de segundo orden
Se utilizará un esquema de discretización de segundo orden para aproximar la
solución. Para implementar el esquema de segundo orden, haga clic en Métodos de solución
Solution Methods, luego haga clic en las ecuaciones de Momentum Turbulence
Kinetic Energy and Turbulence Dissipation Rate y seleccione Second
Order Upwind.
Establecer el valor inicial.
Aquí, la zona de flujo se inicializará a los valores en la entrada. Para llevar a cabo la
inicialización, haga clic en Solution Initialization, luego haga clic en Standard
Initialization, Compute from, y seleccione la entrada inlet.
Luego, haga clic en el botón Inicializar, . Esto completa el proceso de
inicialización.
Establecer los criterios de convergencia.
FLUENT reporta un residuo para cada ecuación de gobierno que se resuelve. El
residual es una medida de qué tan bien la solución actual satisface la forma discreta de cada
ecuación de gobierno. Repetiremos la solución hasta que el residual para cada ecuación caiga
por debajo de 1e-6. Para especificar los criterios residuales (Click) Monitors >
Residuals > Edit....
A continuación, cambie el residual bajo el criterio de convergencia Convergence
Criterion para continuity, x-velocity, y-velocity, así como k-e, todo a 1e-6.
Por último, haga clic en OK para cerrar el menú Residual Monitors.
Cálculo
Antes de ejecutar el cálculo, se debe establecer el número máximo de iteraciones.
Para especificar el número máximo de iteraciones, haga clic en Run Calculation y luego
establezca el Número de iteraciones en 500.
Para evitar perder datos, guarde el proyecto. Ahora, haga clic en Calculate dos
veces para ejecutar el cálculo. Los residuos para cada iteración se imprimen y se trazan en la
ventana gráfica a medida que se calculan.
Agitación y CFD
Csrlos Rubio Atoche
Los residuos caen por debajo del criterio de convergencia especificado de 1e-6 en
aproximadamente 150 iteraciones, como se muestra a continuación. El número real de pasos
de convergencia puede variar ligeramente.
En este punto, guarde el proyecto una vez más.
Resultados
Una vez que nuestro sistema a convergido debemos graficar y tabular (en su caso) los
siguientes valores.
Vectores de Velocidad
Los vectores de velocidad nos permiten visualizar el comportamiento del fluido al
interior de la geometría que estamos estudiando. De igual modo, nos dan una idea de la
calidad de los resultados obtenidos. A través de un balance de materia es posible establecer
la certeza de los mismos.
A partir de ellos podemos establecer las líneas de flujo.
Línea de Velocidad Central
El graficar los vectores de velocidad en el eje del tubo nos permite determinar y evaluar
el desarrollo de la corriente del flujo dentro de la geometría de estudio.
Coeficiente de rozamiento en superficie
El coeficiente de rozamiento es un parámetro adimensional definido como la relación
entre el esfuerzo cortante en la pared y la presión dinámica de referencia:
𝜏𝜔
𝐶𝑓 ≡
1⁄ 𝜌 𝑣 2
2 𝑟𝑒𝑓 𝑟𝑒𝑓
Donde 𝜏𝜔 es el esfuerzo cortante en la pared, y 𝜌𝑟𝑒𝑓 así como 𝑣𝑟𝑒𝑓 son los valores de
referencia definidos para la densidad y la velocidad.
Perfiles de Velocidad
Los perfiles de velocidad nos ayudan a evaluar la evolución y el desarrollo del flujo.
En el caso del flujo laminar el perfil de velocidades corresponde a un desarrollo parabólico.
y+
El flujo turbulento se ve significativamente afectado por la presencia de paredes. El
modelo de turbulencia k-épsilon es principalmente válido lejos de las paredes y se requiere
un tratamiento especial para hacerlo válido cerca de las paredes. El modelo cercano a la pared
es sensible a la resolución de la cuadrícula que se evalúa en la unidad de pared y +.
Agitación y CFD
Csrlos Rubio Atoche
Revisaremos los detalles por ahora y utilizaremos la siguiente regla general: seleccione la
resolución cercana a la pared de manera que y +> 30 o <5 para la celda adyacente a la pared
cuando use la opción Enhanced Wall Treatment.
Verificación validación
Es muy importante que se tome el tiempo de verificar su solución. Esta sección lo guía
a través de algunos controles sobre la solución.
Verificar las condiciones frontera
Se deben verificar las condiciones frontera como son velocidad en la interfase,
velocidad en ejes, los coeficientes de rozamiento.
Verificar el balance de materia
En la barra de menú, haga clic en Informe> Informes de resultados Report >
Result Reports. Asegúrese de que Fluxes esté resaltado y haga clic en Configurar ....
Verifique que se haya seleccionado Caudal másico Mass Flow Rate antes de seleccionar
todos los límites, excepto interior-surface_body.
Verificar los errores numéricos
Se deben verificar los errores numéricos, resultado de los cálculos realizados dentro de
la malla.
Agitación y CFD
Csrlos Rubio Atoche
