Temperatura_A=TA, Temperatura_B=TB, Conductividad=k, Nodos=nx, Volúmenes=nvx, Delta=delta) print(T) # # Definimos la fuente # Su = np.zeros(ivx) Su[:] = q # # Definimos el esquema de disccretización # dif_scheme = tDiffusion1D(malla, Su, dt=dt, Gamma=k) # # Definimos la ecuación a resolver # poisson = PDE(rod, T) # # Preparamos el sistema lineal y creamos la matriz # poisson.setNumericalScheme(dif_scheme) # # Solución analítica # x, _, _ = malla.coordinatesMeshFVM() xa = np.linspace(0, longitud, 100) Ta = analyticSol(xa) #
printInfo(Longitud=longitud, Temperatura_A=TA, Temperatura_B=TB, Conductividad=k, Nodos=nx, Volúmenes=nvx, Delta=delta, dt=dt) # # Definimos la fuente # S = np.zeros(ivx) # Por ahora no hay fuente. # # Definimos el esquema de disccretización # dif_scheme = tDiffusion1D(malla, S, Gamma=k, dt=dt) # # Definimos la ecuación a resolver # laplace = PDE(barra, T) # # Preparamos el sistema lineal y creamos la matriz # laplace.setNumericalScheme(dif_scheme) # # Preparamos la visualización con VisCoFlow # axis_par = [{ 'title': 'Solución Numérica', 'xlabel': 'x [cm]', 'ylabel': 'T [$^o$C]'