GithubHelp home page GithubHelp logo

matteobernardini / streams Goto Github PK

View Code? Open in Web Editor NEW
79.0 12.0 37.0 355 KB

Supersonic TuRbulEnt Accelerated navier stokes Solver

License: GNU General Public License v3.0

Makefile 0.88% Fortran 99.06% C 0.07%

streams's People

Contributors

davidem88 avatar francescosalvadore avatar matteobernardini avatar

Stargazers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

Watchers

 avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar  avatar

streams's Issues

boundary condition

For swbli case, the boundary condition specified at outlet is non reflecting boundary condition but which non reflecting boundary condition is specified is not mentioned anywhere, Can you please give information on which kind of non reflecting bc is implemented at outlet?

w(1,i,j,k)

Hello,
does w(1,i,j,k) anywhere in the code represent the instantaneous density at particular i,j,k?

Thank you

dtsave variable

Does the variable dtsave = 5 in input.dat file signify that field files are written after every 5 seconds?
What information is being obtained in plotxy_000.dat file?

Segmentation fault while running plot3d files

Hello,
I tried to run supersonic channel again in CPU and compiled using intel(2020) compiler.
compile

Compilation was successfully done with no errors or warnings but still when I tried to run the code with enable_plot3d and enable_vtk as 1, it showed segmentation error. When I disabled enable_plot3d and enable_vtk and tried to run it again, it showed segmentation error while writing .bin files.
segmentation
So, I wanted to know if there are any changes that are supposed to be done in code when writing plot3d and vtk files

stat.bin file not found

Hello,
I tried to run the code, and it successfully ran for 158910 iterations, but when I again resubmitted the job it threw an error as shown in this file.
sbli_supersonic.log
But it had stat.bak file and no stat.bin file was generated in the last run, so is it possible to make the code read stat.bak instead of stat.bin and continue the iterations.

Time step

Hello,
I am trying to understand at what actual time steps are field files getting generated. The time variable mentioned in output files are at intervals of 5 as specified by us using "dtsave", but what does the digit signify. And is the time step constant or does it vary?
Thank you

cf files

The properties in the cf files are written for which loation of x,y,z?

Streams2

Hello,
I have run few cases using STREAmS 1, is it possible to continue those runs using STREAmS 2.

plotxy_nnn

Is the data obtained in plotxy_nnn.dat files the mean data or instantaneous data at the end of a particular run?

cf files

Hello,
Can you please provide insight on what each column in cfz_nnn_nnnn and wavplot_nnn_nnnn.dat files represent. These files are generated but there is not much information on what data it is trying to write.

Deflection Angle

Hello,
There are few doubts I had related to the variable deflec_shock in input file.

According to the explanation in "preparing input.dat file", deflec_shock is the shock angle (beeta) whose default value given was 8deg in the input.dat file, but according to the paper https://doi.org/10.1016/j.cpc.2021.107906 the 8deg given is mentioned as Incidence angle of shock generator (theeta). So, I wanted to know what exactly the term deflec_shock is specifying, is it the shock angle (beeta) or deflection angle(theeta).

And if its the deflection angle what is the maximum deflection angle that we can feed into the input file.

Field files and Shock Impingement location (Shock wave Boundary layer Interaction case)

Hello,

  1. I wanted to know at what intervals are field files getting generated, if I want them to be written for every 10,000 iterations is it possible to do so? or is there any other logic behind how these files are being written?
  2. If I wanted to run the SBLI case for different Mach Numbers or for different shock deflection Angles will I have to change the variable "xshock_imp"?

Thank you.

Time series data

Hello,
I am trying to do low frequency unsteadiness study for sbli case and I wanted to know how to obtain time series data.
Is it possible for us to change few parameters in input.dat and get time-series data or should we ourselves make changes in code to get that data?

Thank You

turbulence kinetic energy terms

Hello,
Are the turbulent kinetic Energy terms such as advection, Turbulent transport, production being written in any of the output files, or should we separately make the solver print these terms.

Mean Wall Pressure plot

Hello,
For the case of Shock wave - Boundary layer interaction,

  1. I wanted to know if there is any file that is generated with pressure information or should we calculate it using the convective variable vector that is obtained to obtain the mean wall pressure plot.
  2. Is P1 same as Pinf here?
  3. How to obtain the proper value of Interaction length(L)

free stream velocity

Which value of T(free stream temperature) should be considered to calculate free stream velocity.

Time series data

Is there a possibility of adding probes, to get time series data at a particular location for low frequency unsteadiness study?

Thank you

clarification on database_bl.dat for SBLI example

If I want to simulate a complete new case on Shock impingement on flat plate with input conditions which are different from example test case given in the repository, do I need to change the database_bl.dat file. If yes, then what are the changes that I am supposed to do in the file database_bl.dat.

Thank you

Appending plotxy_nnn files

I generated plotxy files and while appending it, there were zones where no data was printed.
plotxy_whitespace
plotxy_zoomed
How to get rid of these white spaces that are generated between each plotxy file?
Thank you

mean properties

Can you also let me know how can I look at mean flow properties such as mean flow pressure, density and velocity.
I tried to get the slice data and average it for mean flow properties but it has not worked. Is there any way to obtain this data.

Thank You

field files

Hello,
I am trying to understand the unsteady mechanism in the Interaction Zone, in order to do this I will need the field files to be written at a constant time interval, is it possible for me to print field files at a constant time interval?
Thank you

Reattachment shock

Hello,
The paper that is referred here talks about impinging shock and reflected shock but nothing about reattachment shock, neither did I see a reattachment shock in schlieren image that I generated after running the case of SWBLI.
So, will the reattachment shock be captured, if I do necessary changes in input file for the case of SWBLI?

Ref_temp

The reference temperature in input.dat file of Shock-Wave boundary layer interaction, is it the freestream temperature?

ghost nodes

Hello,
Can you please clarify if in the code j=1 represents wall or j = 4 (because of three ghost nodes) represents wall?

Thank you

probe data

Hello,
I wanted time series data at certain points, so I added few lines at 2 places. I have attached screenshots of the same. Can you please confirm if it is correct, because those edits which I considered to be correct are not fetching me required results.

  1. These few lines I have added in Solver.F90, I was tryin to open 7 new files at 7 different locations.
    solver
  2. I added the below lines in printres.F90 to get the data with 7 columns (iteration number, time step, rho, rho-u, rho-v, rho-w, rho-e)
    at 7 different locations among which first 4 locations are at wall
    printres

Thank You

dt_limitations

What does each column in file dt_limitations represent?

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    ๐Ÿ–– Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. ๐Ÿ“Š๐Ÿ“ˆ๐ŸŽ‰

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google โค๏ธ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.