1D advection-diffusion-need to check my code; any suggestions?

In summary, the speaker is teaching themselves transport and programming and has created a program for 1-D advection-diffusion. They are looking for suggestions on how to check if their code is working properly, and are interested in finding an analytical solution or examples in Excel, Matlab, or C++. They mention that they are not familiar with Fortran or C.
  • #1
bzz77
34
0
1D advection-diffusion--need to check my code; any suggestions?

Hi everyone:

I am teaching myself transport/programming, and have created a simple program for 1-D advection-diffusion (constant velocity, diffusion coefficient, central differencing, Crank-Nicolson).

I would like to check that my code is working properly and would appreciate any advice about how to go about this. I have not been successful on Google at finding code that others have written... If anyone can direct me or provide advice about getting an analytical solution, I would be very grateful. An EXCEL-based spreadsheet or a Matlab or a c++ implementation would be useful; I am not a Fortran or C user.
 
Physics news on Phys.org
  • #2


probably the best way is to find a textbook problem and see if your solution agrees with the book's.
 

1. What is 1D advection-diffusion and why is it important in scientific research?

1D advection-diffusion is a mathematical model used to describe the transport and dispersion of particles or substances in a one-dimensional system. It is important in scientific research because it can help us understand and predict the behavior of various substances, such as pollutants in the environment or chemical reactions in biological systems.

2. Can you explain the difference between advection and diffusion in 1D advection-diffusion?

Advection refers to the movement of particles or substances due to bulk flow, such as wind or water currents. Diffusion, on the other hand, is the random movement of particles from a region of higher concentration to a region of lower concentration. In 1D advection-diffusion, both advection and diffusion play a role in the overall transport and dispersion of substances.

3. How do you check your code for accuracy in 1D advection-diffusion simulations?

To check the accuracy of your code, you can compare your simulation results to analytical solutions or experimental data. You can also perform sensitivity analyses by changing input parameters and observing the effects on the output. Additionally, debugging techniques such as step-by-step execution and printing out intermediate values can help identify any errors in the code.

4. Are there any common challenges when coding for 1D advection-diffusion simulations?

One common challenge is ensuring stability in the numerical solution. This can be achieved by using appropriate numerical schemes and ensuring that the time and space steps are small enough. Another challenge is accurately representing boundary conditions, as they can greatly affect the behavior of the system.

5. What suggestions do you have for someone new to coding 1D advection-diffusion simulations?

First, familiarize yourself with the mathematical equations and concepts involved in 1D advection-diffusion. Then, start with simple cases and gradually increase the complexity. It is also helpful to consult with experienced researchers or refer to well-established codes for guidance. And finally, always test and validate your code before using it for any important simulations.

Similar threads

I Non-Homogeneous Robin Boundary conditions and Interpretations of Signs
    • Differential Equations
Replies
2
Views
1K
A 1D Convection Diffusion Equation - Inlet Mixing Effect
    • Differential Equations
Replies
7
Views
3K
(CFD) Problem with C code for 1D linear convection equation
    • Programming and Computer Science
Replies
1
Views
943
MATLAB My Crank-Nicolson code for my diffusion equation isn't working
    • MATLAB, Maple, Mathematica, LaTeX
2
Replies
41
Views
8K
Getting Started with Application Development
    • Programming and Computer Science
Replies
6
Views
1K
I want to write my own Monte Carlo code for Neutron transport
    • Nuclear Engineering
Replies
11
Views
2K
A Chat with ChatGPT: Peerless Intelligence
    • Computing and Technology
Replies
3
Views
2K
B Understanding Physics for Coding Collision of 2 Balls
    • Mechanics
Replies
14
Views
1K
MATLAB Crank-Nicholson solution of 1D heat equation
    • MATLAB, Maple, Mathematica, LaTeX
Replies
3
Views
3K
Coding Bootcamps - Can Anyone be a Six-Figure Coder?
    • STEM Academic Advising
Replies
11
Views
1K

Hot Threads

Recent Insights

Back
Top