FLUBIO—An unstructured, parallel, finite-volume based Navier–Stokes and convection-diffusion like equations solver for teaching and research purposes sciencedirect.com/science…
👍︎ 24
💬︎
📅︎ Jul 08 2021
🚨︎ report
Why is there a density term in the reaction-convection-diffusion equation?

I currently am working with the reaction-convection-diffusion equation where I am modeling the concentration of a generic system that I will begin to make more and more complex. I am able to make quite a few simplifying assumptions to begin that reduce: RCD equation: https://imgur.com/a/DUey7

to

RC equation (no axial diffusion): https://imgur.com/a/5yBSk

My first problem in understanding this is why the source term is multiplied by a density (if that's what it is, the document doesn't specify what exactly rho (p) is, but I assume it is the density of products. If I'm not wrong, every group of terms in either equation must result in units of mol/(L*s), or more generally, concentration per time. If this is true, then rho must be dimensionless since R' is probably already in these units.

This problem carries on when I nondimensionalize the equation to: https://imgur.com/a/Uf4Iv

In this way, the first term on the right side cancels out (no axial diffusion still), but the source term is still multiplied by rho(p) * length/velocity, where length/velocity is equal to the residence time (Tau).

If it's not clear, my question is why is rho(p) in this equation, I've never seen it written like this. If it stands for density, then what are its units, as it seems from the above pictures that it doesn't have any units whatsoever.

👍︎ 2
💬︎
👤︎ u/Apache215
📅︎ Feb 21 2018
🚨︎ report
Numerical modelling of the convection-diffusion equation. FV or FD?

Suppose you have a transient convection-diffusion equation with a source term. Do you model this system with the finite difference method or the finite volume method? What do you take into account when choosing between these two methods?

Edit: the problem is one-dimensional.

👍︎ 5
💬︎
📅︎ Dec 31 2015
🚨︎ report
Using OpenFOAM to run a 2D convection-diffusion simulation of a compressible steady-state cfd-online.com/Forums/ope…
👍︎ 4
💬︎
👤︎ u/foadsf
📅︎ Oct 26 2021
🚨︎ report
heat equation pde for convection?

Using energy analysis on a control volume i have followed the development of the heat equation for convection in a metal rod. DT/dt =ad^2T/dx^2

I want to apply a similar analysis on a forced convection system where water of a known flow and temperature field flow over a heated plate or flow around a heated resistor

Can anyone suggest examples or resources that demonstrate this anysis

Thanks

👍︎ 10
💬︎
📅︎ Jan 08 2022
🚨︎ report
2D linear convection is solved in MATLAB. Since it is solved numerically we can see "artificial diffusion". That is originally the wave is meant to travel as a square wave but due to approximations of the differential eqs the wave diffuses. youtube.com/watch?v=40vlT…
👍︎ 21
💬︎
👤︎ u/no1sunk
📅︎ Mar 23 2021
🚨︎ report
Why does the simple Galerken approach for a convection diffusion problem fail? I know about the reasoning which relates to finite difference and central differencing but that doesn't feel satisfactory.
👍︎ 2
💬︎
📅︎ May 11 2021
🚨︎ report
What is the second term here in the heat equation? The first seems to be the accumulation term, and the last seems to be the conduction in minus the conduction out (net conduction). But what is the middle term? Convection term maybe? Thank you imgur.com/gallery/7U7JzFY
👍︎ 12
💬︎
👤︎ u/Rilandi
📅︎ Aug 18 2021
🚨︎ report
What is the significance of the grid peclet number for diffusion and convection

Hi guys,

My question is simply what I have as the subject of this post. I'm in a CFD course, and my professor asked us to discuss the relevance of a set of simulations and how their grid peclet numbers relate to the simulation. I went on Wikipedia and did some research but I'm finding it hard to understand what this number truly means.

👍︎ 2
💬︎
📅︎ Feb 20 2021
🚨︎ report
Simple question but, why is unwinding only applied to the convection term and not on the diffusion term too?

edit: I meant to say "upwinding"

👍︎ 8
💬︎
📅︎ Dec 20 2019
🚨︎ report
Help: spurious oscillations in 1-D heat diffusion equation

Hi! i'm working on solving some equaion which looks like 1-D heat diffusion equation. I tryed to solve this equation in order to check if my solver works well, but I got spurious oscillations in my solution. My code uses a Crank-Nicholson scheme for time advancing and five point scheme for espatial derivation.

The equation i'm trying to solve

I set Neumann boundary condition dY/dt=0 and solved it on a 20 points grid deffined in, with 100 time step of 1/20 length. This was my result:

Solution time evolution in a central point.

Full solution mesh.

Actualization:

I modified the code to use a ' θ scheme' ( θ =0 explicit, θ =1 implicit, θ =.5 C-N) and continued proving. I noticed pair and odd time steps give me two very different solution sets (initial values at time step 1).

Pair time steps solution set is the weirdest one.

Pair time step solution set.

Second actualization:

I solved it, it appears to me related to operations like 'A\M' (where M is a matrix). I stoped using that operations, I instead calculated inv(M) (luckly it's needed to be calculated only once) and those oscilaions disappeared.

Thanks you all guys! ^^

👍︎ 7
💬︎
👤︎ u/neutrilo
📅︎ Apr 10 2021
🚨︎ report
Can someone explain the physical interpretation of both transport by convection and transport by diffusion? Or maybe explaining the difference between the two. Google is no help
👍︎ 13
💬︎
📅︎ Jul 12 2018
🚨︎ report
solving spherical diffusion equation (FDM)

Hi everyone! i'm asking for suggestions to do it. i'm trying to solve spherical diffusion equation with FDM. Lucky i'm only interested about radial direction.

https://preview.redd.it/5jh14lgd13o61.png?width=282&format=png&auto=webp&s=0b1b0a64813dc382d0bb68ff12952b1fddf74687

At first I tried by using Runge-Kutta IV, and a 5 points derivative scheme (I already wrote a code to get 1st and 2d derivates based on 5 points). It appears to work, but it runs too slow, so that I'm thinking about try to solve it using implicit methods, like Cranck-Nicholson, but I just realized I've neverd used implicit methods in spherical coordinates.

So, do you have any suggestion to me?

👍︎ 8
💬︎
👤︎ u/neutrilo
📅︎ Mar 20 2021
🚨︎ report
If you're wondering how disease models work, it's all to do with differential equations. The basic SIR model coupled with diffusion leads to 'Travelling Wave' solutions that describe the speed at which an outbreak spreads through a population. Includes some neat PDE analysis useful for physicists. youtube.com/watch?v=uSLFu…
👍︎ 749
💬︎
📅︎ Mar 27 2020
🚨︎ report
[University Maths] Writing a code in MATLAB for solving a convection-diffusion-reaction problem using the finite difference method

I've been on this problem for ages and it's due tomorrow but I can't see to get anywhere with it. This is the imgur link because I can't format very well

https://i.imgur.com/G5vJKzF.png?1

I'm really stressing over this so please don't be too harsh on me

👍︎ 3
💬︎
👤︎ u/Lonely7717
📅︎ Dec 13 2017
🚨︎ report
Convection problem? Heat diffusion with moving plate

I'm stuck on a homework problem. It says to find the temperature profile (as a function of y) for an infinitely large plate that at t=0 starts moving and changes its temperature from T_inf to T_s. I can do the temperature profile without the plate moving, but I don't know how to solve the temperature profile when the plate is moving. Can anyone point me in the right direction?

👍︎ 2
💬︎
📅︎ Jan 31 2017
🚨︎ report
CFD Runge Kutta 4th order Advection Convection Equation PDE

Hi everyone! I'm stuck with the math of this problem, can you please identify my error? Whenever I program this, the solution explodes in RK4. This is a problem of CFD or numerical methods.

Here are the equations I worked with. I have a vector φ with an initial solution (pointy gaussian). The dimension is 1x213.

And I want to compute the solution of the PDE of advection / convection equation 1st order using upwind with Forward Euler and RK4, however, RK4 seems to be bad programmed.
In the pictures you have my equations. Can you find the error?

Initial Steps

Intermediate Steps

Final steps

Comparison between Forward Euler Upwind with RK4 Upwind ; CFL=1

Note: Forward Euler upwind with CFL=1 gives the exact solution.
The code that I put in matlab is:

n= time step; j= space step; dt=delta time; CFL= adt/dx;
K_1(n,j) = -(CFL/dt)*(phi(n,j)-phi(n,j-1));

K_2(n,j) = phi(n,j) - 0.5*CFL* (phi(n,j)-phi(n,j-1));

K_3(n,j) = phi(n,j) - 0.5*CFL* (K_2(n,j)-K_2(n,j-1));

K_4(n,j) = phi(n,j) - CFL* ( K_3(n,j) - K_3(n,j-1) );

phi(n+1,j) = phi(n,j)+ (dt/6)*(K_1(n,j) +2*K_2(n,j) +2*K_3(n,j) + K_4(n,j));

👍︎ 3
💬︎
📅︎ Sep 25 2020
🚨︎ report
[Graduate Boundary Value Problem] Please help me solve the unsteady convection diffusion problem on the upper half plane with a non-homogenous boundary and time dependent coefficient.

A graphic is available here. Many thanks

👍︎ 4
💬︎
👤︎ u/znffal
📅︎ Aug 21 2015
🚨︎ report
Tensor-driven anisotropic diffusion solves the heat equation on a manifold defined by a tensor field. At the heart of many methods for image restoration, texture enhancement, and visualization (by @gabrielpeyre) v.redd.it/javd4a8dcn651
👍︎ 91
💬︎
📅︎ Jun 23 2020
🚨︎ report
[University Partial Differential Equations] How to solve wave and heat/diffusion equation with Robin Boundary Conditions?

I've been watching the lectures over and over and I asked my professor for help, but I'm still so lost. I can solve the IVPs with the Fourier Series, but I can't understand how to derive the formula for the eigenvalues and eigenfunctions with Dirichlet and Neumann, let alone for Robin since you have to do it graphically. Can anybody help?

👍︎ 2
💬︎
👤︎ u/Quarentus
📅︎ Dec 17 2020
🚨︎ report
Diffusion-Reaction Equations with custom initial conditions
👍︎ 113
💬︎
📅︎ Jan 02 2020
🚨︎ report
Diffusion equation help needed

I need help with the following problem:

Determine the temperature distribution u(x,t) for an infinite 1D rod using the diffusion equation if it's initial temperature is given by f(x);

where

f(x)={T_0 if x<0 and 0 if x⪖0

I know I start by saying that u(x,t)=X(x)T(t) and obtaining the expressions for X(x) and T(t), or at least that's what I think, but don't know where to continue from there.

I'd really appreciate any amount of help I could get. Thanks! ​ ​

👍︎ 5
💬︎
👤︎ u/marigo3467
📅︎ Jul 12 2020
🚨︎ report
Trouble understanding the diffusion equation and context, can someone explain each part further than what the text did on a basic level?
👍︎ 3
💬︎
👤︎ u/Zoobleroni
📅︎ Sep 21 2020
🚨︎ report
Newton's solver not converging for 1D nonlinear diffusion equation.

Hello all, I am having trouble with a newton solver that I have been working on over the last week. It is in the context of modelling 2D random with proliferation walks via column averaging but that's beside the point. The PDE that describes this interaction is where D is the diffusion (migration) terma and lambda is the nonlinear (proliferation).

N_t = D * N_xx + lambda * N * (1 - N)

I have checked the equations used for the JAcobian and the f vector a dozen times to the notes in class so I'm 99% sure that's not the issue.

An implicit Euler method is used for those interested.

Newton's method fails to converge when proliferation is 'turned on' (pp > 0). The code works fine for pp = 0 (compared to the random walk data). Whenever I try pp = 0.01 i get the error('Newton failed to converge\n') message.

Any help would be greatly appreciated :D, code below

%% 1D Non-Linear Diffusion Equaltion Solver using Newton's Method
clear all
clc

% Set up variables
deltax = 0.1;
tau = 0.1; % Time step
xlims = 20; % X boundaries
xIC = 10; % Lattice Sites with the IC

pp = 0.01; % Probablity to proliferate
pm = 0.5;  % Probability for a migration


Newtol = 0.01; % newton tolerance
xarr = linspace(-xlims,xlims,2*xlims/deltax+1) + xlims; % x array

I = length(xarr); % number of lattice point
f = zeros(I,1); % the f vector used in newton iterations
J = zeros(I); % Jacobian for newton iterations
N = zeros(I,1); % Vector with averaged particle position
deltaN = N; % deltax vector for newton iterations
N(round(xlims/deltax+1)-xIC:round(xlims/deltax+1)+xIC) = 40; % Set up casona IC [-inf,0] = 1; (0,Inf] = 0


Nold = N; % used for newton iterations

curtol = Newtol + 1; % current tol, inital value must be larger than newtol
numsteps = 0; % newton steps
maxnewsteps = 10; % max newton steps
cond = 1; % convergence condition

D = pm * deltax^2/(4*tau); % diffusion constant
lam = pp / tau; % Nonlinear constant


% Set Tmax and the timestep array
Tmax = 100;

% for all the timesteps
for i = 1:(Tmax / tau)
    cond = 1; % reset convergence condition
    numsteps = 0; % reset newton step count

    % while newton has not converged
    while cond

        % Calculate the Jacobian and the f vector

        % Left Boundary
        J(1,1) = -1;
        J(1,2) = 1;
        f(1) = N(2) - N(1);

        % Right Boundary
... keep reading on reddit ➡

👍︎ 9
💬︎
📅︎ Apr 24 2020
🚨︎ report
Granular Convection and Diffusion

I just watched a video on granular convection (the "hazelnut effect"), and it brought up the question of what the implications of it there are with respect to diffusion.

Why would granular material being agitated tend to clump objects based on size (I'm looking at a case where you have a tray of objects, as opposed to a column, so gravity isn't the cause), VS diffusion spreading everything evenly?

👍︎ 2
💬︎
👤︎ u/stravant
📅︎ Apr 09 2011
🚨︎ report
Random Walk/Diffusion Equation Calculation

First time posting here but I’m completely stuck on figuring out how to either input this equation into a scientific calculator, or to solve it accurately.

I’m doing a chemistry course involving food and I’ve reached the topic of random walks and the diffusion equation.

The example in the course is measuring how long a mixture of a droplet of juice and alginate will take a develop a 1mm thick shell.

The formula given is:

D’s value is - 8*10^-10 m^2/sec

L^2 = 4Dt t = L^2 / 4D t = (10^-3 m)^2 / 4(D)

With the answer averaged at; t = 300 sec = 5 minutes

I’ve tried putting it into a scientific calculator and solving it step by step but I’m not getting the correct answer which should be averaged at 300 seconds.

👍︎ 2
💬︎
👤︎ u/j8a6m8es
📅︎ Jun 15 2020
🚨︎ report
Solutions to the 1D diffusion equation
👍︎ 5
💬︎
👤︎ u/wwwiley
📅︎ Sep 02 2020
🚨︎ report
Solutions to the diffusion equation
👍︎ 7
💬︎
👤︎ u/wwwiley
📅︎ Jul 26 2020
🚨︎ report
Where does Diffusion fall in The Radiative Transfer Equation?

Is it part of the source function? I've had a hard time breaking down the components in the RTE to specifically look at the direct and diffuse components. If anyone knows a good source (that I can access) that would be extremely helpful.

👍︎ 3
💬︎
📅︎ Aug 26 2020
🚨︎ report
The phenomenological equation for diffusion is the most fundamental one. How can one corelate it with Fick’s first law of diffusion on corelating ‘chemical potential gradient’ with ‘concentration gradient’

One can corelate it with Fick’s first law of diffusion on corelating ‘chemical potential gradient’ with ‘concentration gradient’. Can anyone please explain with/without proper equations?

👍︎ 2
💬︎
👤︎ u/adakankit
📅︎ May 05 2020
🚨︎ report
Numerical modelling of the convection-diffusion equation. FD or FV?

Suppose you have a transient convection-diffusion equation with a source term. Do you model this system with the finite difference method or the finite volume method? What do you take into account when choosing between these two methods?

Edit: the problem is one-dimensional.

👍︎ 5
💬︎
📅︎ Dec 31 2015
🚨︎ report

Please note that this site uses cookies to personalise content and adverts, to provide social media features, and to analyse web traffic. Click here for more information.