Go to the source code of this file.
Functions | |
int | callocate_arrays_global (int x, int y, int z, int max_array_length) |
int | callocate_arrays_sweeps (int max_array_length) |
int | callocate_arrays_ppm (int max_array_length) |
int | free_sweep (void) |
int | free_ppm (void) |
int | initiate_grid (int x, double xmin, double xmax, int y, double ymin, double ymax, int z, double zmin, double zmax) |
int | initiate_domain (int x, int y, int z) |
int | initiate_domain_dom (int x, int y, int z) |
int | initiate_domain_marker (int x, int y, int z) |
int | hydro_sweeps (int x, int y, int z, int direction) |
int | init_ppm (double **rho_1D, double **pre_1D, double **eng_1D, double **vx_1D, double **vy_1D, double **vz_1D, double **marker_1D, double **dx0, double **xa0, double **xa, double **dx, double **a_coef, double **ai_coef, double **b_coef, double **bi_coef, double **c_coef, double **ci_coef, double **d_x, double **da, double **ar, double **dp, double **dr, double **du, double **pl, double **p6, double **rl, double **r6, double **ul, double **u6, double **vl, double **v6, double **wl, double **w6, double **el, double **e6, double **ql, double **q6, double **deltaa, double **dv, double **dw, double **dq, double **de, double **scratch1, double **scratch2, double **scratch3, double **diffa, double **plft, double **prgh, double **ulft, double **urgh, double **rlft, double **rrgh, double **Cdtdx, double **fCdtdx, double **clft, double **crgh, double **plfti, double **prghi, double **pmid, double **pmold, double **wlft, double **wrgh, double **zlft, double **zrgh, double **umidl, double **umidr, double **umid, double **dm, double **dtbdm, double **upmid, double **xa1, double **xa2, double **xa3, double **vx_1D_old, double **e_int_1D, double **dvol, double **dvol0, double **dvol1, double **delta, double **fluxr, double **fluxu, double **fluxv, double **fluxw, double **fluxe, double **fluxq, double **dm0, double **steep, double **flat, double ***para, double **pressure_solid_1D, double **rhodown, double **rhoup, double **rhofront, double **rhoback, double **vxdown, double **vxup, double **vxfront, double **vxback, double **vydown, double **vyup, double **vyfront, double **vyback, double **vzdown, double **vzup, double **vzfront, double **vzback, int dimension) |
int | ppm_free (double **rho_1D, double **pre_1D, double **eng_1D, double **vx_1D, double **vy_1D, double **vz_1D, double **marker_1D, double **pressure_solid_1D, double **dx0, double **xa0, double **xa, double **dx, double **a_coef, double **ai_coef, double **b_coef, double **bi_coef, double **c_coef, double **ci_coef, double **d_x, double **da, double **ar, double **dp, double **dr, double **du, double **pl, double **p6, double **rl, double **r6, double **ul, double **u6, double **vl, double **v6, double **wl, double **w6, double **el, double **e6, double **ql, double **q6, double **deltaa, double **dv, double **dw, double **dq, double **de, double **scratch1, double **scratch2, double **scratch3, double **diffa, double **plft, double **prgh, double **ulft, double **urgh, double **rlft, double **rrgh, double **Cdtdx, double **fCdtdx, double **clft, double **crgh, double **plfti, double **prghi, double **pmid, double **pmold, double **wlft, double **wrgh, double **zlft, double **zrgh, double **umidl, double **umidr, double **umid, double **dm, double **dtbdm, double **upmid, double **xa1, double **xa2, double **xa3, double **vx_1D_old, double **e_int_1D, double **dvol, double **dvol0, double **dvol1, double **delta, double **fluxr, double **fluxu, double **fluxv, double **fluxw, double **fluxe, double **fluxq, double **dm0, double **steep, double **flat, double ***para, double **rhodown, double **rhoup, double **rhofront, double **rhoback, double **vxdown, double **vxup, double **vxfront, double **vxback, double **vydown, double **vyup, double **vyfront, double **vyback, double **vzdown, double **vzup, double **vzfront, double **vzback, int dimension) |
int | init_diffusion (double **temperature_1D, double **temperature_1D_future, double **mass, int **dom_1D) |
int | make_ic (int x, int y, int z) |
int | make_sod_ic (int x, int y, int z) |
int | make_kh_instabilities (int x, int y, int z) |
int | read_ic (int x, int y, int z) |
int | read_restart (int x, int y, int z) |
int | read_dom (int x, int y, int z) |
int | read_marker_file (int x, int y, int z) |
int | start_file_reader (char filename[]) |
int | calculate_pressure (int x, int y, int z, double temperature) |
long long int | ntohll (const long long int data) |
int | write_ic_tyc (int x, int y, int z) |
int | write_ic_vtk (int x, int y, int z) |
int | write_ic_amira (int x, int y, int z) |
int | write_ic_ifrit (int x, int y, int z) |
int | write_tyc (int x, int y, int z, int counter) |
int | write_vtk (int x, int y, int z, int counter) |
int | write_amira (int x, int y, int z, int counter) |
int | write_ifrit (int x, int y, int z, int counter) |
int | write_restart_tyc (int x, int y, int z) |
int | dt_calc (int x, int y, int z) |
int | dt_calc_first (int x, int y, int z) |
int | sweep_x (int x, int y, int z, int flag) |
int | sweep_y (int x, int y, int z, int flag) |
int | sweep_z (int x, int y, int z, int flag) |
int | ppm_step (int i, int j, int k, int direction, int flag, int nmin, int nmax, double *a_coef, double *ai_coef, double *b_coef, double *bi_coef, double *c_coef, double *ci_coef, double *d_x, double *diffa, double *da, double *ar, double *pl, double *p6, double *rl, double *r6, double *u6, double *ul, double *vl, double *v6, double *wl, double *w6, double *el, double *e6, double *ql, double *q6, double *dp, double *du, double *dr, double *dv, double *dw, double *dq, double *de, double *scratch1, double *scratch2, double *scratch3, double *plft, double *prgh, double *ulft, double *urgh, double *rlft, double *rrgh, double *Cdtdx, double *fCdtdx, double *steep, double *flat, double **para, double *clft, double *crgh, double *plfti, double *prghi, double *pmid, double *pmold, double *wlft, double *wrgh, double *zlft, double *zrgh, double *umidl, double *umidr, double *umid, double *dm, double *dtbdm, double *upmid, double *xa1, double *xa2, double *xa3, double *vx_1D_old, double *dvol, double *dvol0, double *dvol1, double *delta, double *fluxr, double *fluxu, double *fluxv, double *fluxw, double *fluxe, double *fluxq, double *dm0, double *e_int_1D, double *rho_1D, double *pre_1D, double *eng_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *marker_1D, double *pressure_solid_1D, double *dx0, double *xa0, double *xa, double *dx, int bound_checker, int lefter, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | set_boundary (int nmin, int nmax, int flag, int bound_checker, int lefter, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *pressure_solid_1D, double *xa0, double *dx0, double *xa, double *dx, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_zero_gradient (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_zero_gradient (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_reflecting (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_reflecting (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_small_padding_on_obstacle (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *pressure_solid_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_reflecting_on_obstacle (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *pressure_solid_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_reflecting_on_obstacle (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *pressure_solid_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_small_padding (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_small_padding (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_outflow (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_outflow (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_inflow (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_inflow (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | left_boundary_periodic (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | right_boundary_periodic (int nmin, int nmax, double *rho_1D, double *eng_1D, double *pre_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *xa0, double *dx0, double *xa, double *dx, int flag, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int viscosity_on_off, int dimension) |
int | para_coef (int nmin, int nmax, int flag1, double *a_coef, double *dx, double *ai_coef, double *b_coef, double *bi_coef, double *c_coef, double *d_x, double **para, double *ci_coef) |
int | volume (int nmin, int nmax, double *vol, double *dx, double *vol0, double *dx0) |
int | flatten (int nmin, int nmax, double *pre_1D, double *vx_1D, double *steep, double *flat) |
int | parabola (int nmin, int nmax, double *a, double *deltaa, double *a6, double *al, int flag1, double *diffa, double *da, double **para, double *ar, double *flat, double *scratch1, double *scratch2, double *scratch3) |
int | states (int nmin, int nmax, int flag, double *pre_1D, double *rho_1D, double *dx, double *Cdtdx, double *fCdtdx, double *plft, double *pl, double *dp, double *p6, double *ulft, double *ul, double *du, double *u6, double *rlft, double *rl, double *dr, double *r6, double *prgh, double *urgh, double *rrgh) |
int | riemann (int nmin, int nmax, double *clft, double *crgh, double *rlft, double *rrgh, double *plfti, double *prghi, double *pmid, double *pmold, double *plft, double *wrgh, double *prgh, double *wlft, double *zlft, double *zrgh, double *umidl, double *umidr, double *umid, double *urgh, double *ulft) |
int | evolve (int nmin, int nmax, int flag, double *rho_1D, double *dvol, double *dm, double *dtbdm, double *xa1, double *xa, double *dvol1, double *umid, double *upmid, double *pmid, double *xa2, double *dx, double *xa3, double *vx_1D_old, double *vx_1D, double *vy_1D, double *vz_1D, double *eng_1D, double *e_int_1D, double *pre_1D) |
int | remap (int nmin, int nmax, int flag, double *a_coef, double *dx, double *ai_coef, double *b_coef, double *bi_coef, double *c_coef, double *d_x, double **para, double *ci_coef, double *dr, double *r6, double *rl, double *diffa, double *da, double *ar, double *flat, double *scratch1, double *scratch2, double *scratch3, double *du, double *u6, double *ul, double *dv, double *v6, double *vl, double *w6, double *wl, double *dq, double *q6, double *ql, double *de, double *e6, double *el, double *xa, double *xa0, double *delta, double *fluxr, double *fluxu, double *fluxv, double *fluxw, double *dw, double *fluxe, double *fluxq, double *dm, double *rho_1D, double *dvol, double *dvol0, double *dm0, double *vx_1D, double *vy_1D, double *vz_1D, double *eng_1D, double *e_int_1D, double *pre_1D) |
int | advect (int nmin, int nmax, int flag, double *dx, double *vx_1D, double *marker_1D) |
int | diffusion (int x, int y, int z) |
double | alpha_heat_transfer (int x, int y, int z, int direction) |
int | copy_arrays_for_viscosity (int x, int y, int z) |
int | viscosity (int i, int j, int k, int flag, int nmin, int nmax, double *rho_1D, double *vx_1D, double *vy_1D, double *vz_1D, double *pre_1D, double *e_int_1D, double *eng_1D, int lefter, double *rhodown, double *rhoup, double *rhofront, double *rhoback, double *vxdown, double *vxup, double *vxfront, double *vxback, double *vydown, double *vyup, double *vyfront, double *vyback, double *vzdown, double *vzup, double *vzfront, double *vzback, int dimension) |
double | kinematic_viscosity (double temperature, double density) |
double | dynamic_viscosity (double temperature) |
int | wind (int x, int y, int z) |
int | upper_atmos (int x, int y, int z) |
int | velocity_field_analyser (int x, int y, int z) |
int | pressure_on_solid_calc (int x, int y, int z) |
int | boundary_velo_corrector (int x, int y, int z) |
int | insert_pressure (int x, int y, int z) |
int | reset_pressure_integrated (int x, int y, int z) |
int | prepare_the_dB_map (int x, int y, int z) |
int | pre_old_copy (int x, int y, int z) |
int advect | ( | int | nmin, |
int | nmax, | ||
int | flag, | ||
double * | dx, | ||
double * | vx_1D, | ||
double * | marker_1D | ||
) |
Advect the marker field with the velocities of the hydro-solver with a Second-order upwind scheme
Definition at line 18 of file advection.c.
References dt, spacing, x, xmax, y, ymax, z, and zmax.
Referenced by ppm_step().
double alpha_heat_transfer | ( | int | x, |
int | y, | ||
int | z, | ||
int | direction | ||
) |
Calculates the heat-transfer coefficient for the heat transfer between obstacles and gas. Here you can put your special physics for the heat transfer coefficient in. Example for flat wall.
Definition at line 25 of file alpha_heat_transfer.c.
Referenced by diffusion().
int boundary_velo_corrector | ( | int | x, |
int | y, | ||
int | z | ||
) |
example function to suppress boundary streams; !!highly dangerous, only for testing!!
Definition at line 18 of file boundary_velo_corrector.c.
References small, vx, wind_direction, x, and y.
int calculate_pressure | ( | int | x, |
int | y, | ||
int | z, | ||
double | temperature | ||
) |
int callocate_arrays_global | ( | int | x, |
int | y, | ||
int | z, | ||
int | max_array_length | ||
) |
int callocate_arrays_ppm | ( | int | max_array_length | ) |
int callocate_arrays_sweeps | ( | int | max_array_length | ) |
int copy_arrays_for_viscosity | ( | int | x, |
int | y, | ||
int | z | ||
) |
Copies the density and velocity arrays before each hydro sweep for the viscosity calculation.
Definition at line 24 of file copy_arrays_for_viscosity.c.
References rho, rho_visc, vx, vx_visc, vy, vy_visc, vz, vz_visc, x, y, and z.
Referenced by hydro_sweeps().
int diffusion | ( | int | x, |
int | y, | ||
int | z | ||
) |
In order to enable diffusion in an explicit way, this routine solves the diffusion equation in a very basic way using operator splitting. Use it with caution! This is only done within the computational domain of the obstacles.
Definition at line 25 of file diffusion.c.
References alpha_heat_transfer(), dimension, dom, dt, free_diffusion(), gasconstant, init_diffusion(), obstacle_heat_conductivity, pre, rho, specific_heat_capacity_gas, specific_heat_capacity_obstacle, volume(), x, xmax, xmin, y, ymax, ymin, z, zmax, and zmin.
Referenced by main().
int dt_calc | ( | int | x, |
int | y, | ||
int | z | ||
) |
This function calculates the hydrodynamical timestep for the simulation. It is based on the signal velocity in the medium and the courant factor, specified in the parameterfile.
Definition at line 50 of file dt_calc.c.
References courant, dimension, dom, dt, Gamma1, gasconstant, kinematic_viscosity(), max_dt(), min_dt(), olddt, pre, rho, viscosity_on_off, vx, vy, vz, x, y, z, zdx, zdy, and zdz.
Referenced by main().
int dt_calc_first | ( | int | x, |
int | y, | ||
int | z | ||
) |
This function calculates the hydrodynamical time step for the first time in the simulation It is based on the signal velocity in the medium and the courant factor, specified in the parameterfile.
Definition at line 50 of file dt_calc_first.c.
References courant, dimension, dom, dt, Gamma1, gasconstant, intial_soundspeed, kinematic_viscosity(), max_dt_first(), min_dt_first(), pre, rho, spacing, viscosity_on_off, vx, vy, vz, x, y, z, zdx, zdy, and zdz.
Referenced by main().
double dynamic_viscosity | ( | double | temperature | ) |
In order to enable viscosity in an explicit way, this routine calculates the dynamic_viscosity for a given temperature according to Sutherland's law.
Definition at line 42 of file kinematic_dyanmic_viscosity.c.
References C1_visc, dynamic_viscosity(), and S_visc.
Referenced by dynamic_viscosity(), and viscosity().
int evolve | ( | int | nmin, |
int | nmax, | ||
int | flag, | ||
double * | rho_1D, | ||
double * | dvol, | ||
double * | dm, | ||
double * | dtbdm, | ||
double * | xa1, | ||
double * | xa, | ||
double * | dvol1, | ||
double * | umid, | ||
double * | upmid, | ||
double * | pmid, | ||
double * | xa2, | ||
double * | dx, | ||
double * | xa3, | ||
double * | vx_1D_old, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | eng_1D, | ||
double * | e_int_1D, | ||
double * | pre_1D | ||
) |
Evolution of velocities and energies due to pressure acceleration and external forces.
Definition at line 32 of file evolve.c.
References dt, Gamma, grav_acc, gravity_on_off, max_evolve(), smallp, and smallr.
Referenced by ppm_step().
int flatten | ( | int | nmin, |
int | nmax, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | steep, | ||
double * | flat | ||
) |
Look for presence of a shock using pressure gradient and sign of velocity jump: shock = 1 if there is a shock in the zone, else shock = 0 Compute steepness parameter based on steepness of pressure jump IF there is a shock. The flatting is needed to eliminated post-shock oscillations.
Definition at line 47 of file flatten.c.
References max_flatten(), min_flatten(), shock_or_not, and small.
Referenced by ppm_step().
int free_ppm | ( | void | ) |
int free_sweep | ( | void | ) |
int hydro_sweeps | ( | int | x, |
int | y, | ||
int | z, | ||
int | direction | ||
) |
The forward and backward hydro sweeps the direction gives the forward [0] or backward [1] sweeps
Definition at line 18 of file hydro_sweeps.c.
References copy_arrays_for_viscosity(), dimension, sweep_x(), sweep_y(), sweep_z(), viscosity_on_off, wind(), and wind_on_off.
Referenced by main().
int init_diffusion | ( | double ** | temperature_1D, |
double ** | temperature_1D_future, | ||
double ** | mass, | ||
int ** | dom_1D | ||
) |
Definition at line 326 of file init.c.
References max_array_length.
Referenced by diffusion().
int init_ppm | ( | double ** | rho_1D, |
double ** | pre_1D, | ||
double ** | eng_1D, | ||
double ** | vx_1D, | ||
double ** | vy_1D, | ||
double ** | vz_1D, | ||
double ** | marker_1D, | ||
double ** | dx0, | ||
double ** | xa0, | ||
double ** | xa, | ||
double ** | dx, | ||
double ** | a_coef, | ||
double ** | ai_coef, | ||
double ** | b_coef, | ||
double ** | bi_coef, | ||
double ** | c_coef, | ||
double ** | ci_coef, | ||
double ** | d_x, | ||
double ** | da, | ||
double ** | ar, | ||
double ** | dp, | ||
double ** | dr, | ||
double ** | du, | ||
double ** | pl, | ||
double ** | p6, | ||
double ** | rl, | ||
double ** | r6, | ||
double ** | ul, | ||
double ** | u6, | ||
double ** | vl, | ||
double ** | v6, | ||
double ** | wl, | ||
double ** | w6, | ||
double ** | el, | ||
double ** | e6, | ||
double ** | ql, | ||
double ** | q6, | ||
double ** | deltaa, | ||
double ** | dv, | ||
double ** | dw, | ||
double ** | dq, | ||
double ** | de, | ||
double ** | scratch1, | ||
double ** | scratch2, | ||
double ** | scratch3, | ||
double ** | diffa, | ||
double ** | plft, | ||
double ** | prgh, | ||
double ** | ulft, | ||
double ** | urgh, | ||
double ** | rlft, | ||
double ** | rrgh, | ||
double ** | Cdtdx, | ||
double ** | fCdtdx, | ||
double ** | clft, | ||
double ** | crgh, | ||
double ** | plfti, | ||
double ** | prghi, | ||
double ** | pmid, | ||
double ** | pmold, | ||
double ** | wlft, | ||
double ** | wrgh, | ||
double ** | zlft, | ||
double ** | zrgh, | ||
double ** | umidl, | ||
double ** | umidr, | ||
double ** | umid, | ||
double ** | dm, | ||
double ** | dtbdm, | ||
double ** | upmid, | ||
double ** | xa1, | ||
double ** | xa2, | ||
double ** | xa3, | ||
double ** | vx_1D_old, | ||
double ** | e_int_1D, | ||
double ** | dvol, | ||
double ** | dvol0, | ||
double ** | dvol1, | ||
double ** | delta, | ||
double ** | fluxr, | ||
double ** | fluxu, | ||
double ** | fluxv, | ||
double ** | fluxw, | ||
double ** | fluxe, | ||
double ** | fluxq, | ||
double ** | dm0, | ||
double ** | steep, | ||
double ** | flat, | ||
double *** | para, | ||
double ** | pressure_solid_1D, | ||
double ** | rhodown, | ||
double ** | rhoup, | ||
double ** | rhofront, | ||
double ** | rhoback, | ||
double ** | vxdown, | ||
double ** | vxup, | ||
double ** | vxfront, | ||
double ** | vxback, | ||
double ** | vydown, | ||
double ** | vyup, | ||
double ** | vyfront, | ||
double ** | vyback, | ||
double ** | vzdown, | ||
double ** | vzup, | ||
double ** | vzfront, | ||
double ** | vzback, | ||
int | dimension | ||
) |
all quantities used in the hydro computation are initiated with zeros
Definition at line 18 of file init.c.
References max_array_length.
Referenced by sweep_x(), sweep_y(), and sweep_z().
int initiate_domain | ( | int | x, |
int | y, | ||
int | z | ||
) |
The form of the obstacle is defined in this function example of just a round thing.
Definition at line 112 of file make_ic.c.
References dom, obstacle_density, obstacle_temperature, pre, rho, vx, vy, vz, x, y, and z.
Referenced by main().
int initiate_domain_dom | ( | int | x, |
int | y, | ||
int | z | ||
) |
int initiate_domain_marker | ( | int | x, |
int | y, | ||
int | z | ||
) |
int initiate_grid | ( | int | x, |
double | xmin, | ||
double | xmax, | ||
int | y, | ||
double | ymin, | ||
double | ymax, | ||
int | z, | ||
double | zmin, | ||
double | zmax | ||
) |
In this function the grid-parameters, such as spacing and physical coordinates for each cell are calculated
Definition at line 16 of file initiate_grid.c.
References spacing, x, xmin, y, ymin, z, zdx, zdy, zdz, zmin, zxa, zxc, zya, zyc, zza, and zzc.
Referenced by main().
int insert_pressure | ( | int | x, |
int | y, | ||
int | z | ||
) |
double kinematic_viscosity | ( | double | temperature, |
double | density | ||
) |
In order to enable viscosity in an explicit way, this routine calculates the kinematic_viscosity for a given temperature according to Sutherland's law.
Definition at line 23 of file kinematic_dyanmic_viscosity.c.
References C1_visc, kinematic_viscosity(), and S_visc.
Referenced by dt_calc(), dt_calc_first(), and kinematic_viscosity().
int left_boundary_inflow | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
outflow boundary condition
Definition at line 1092 of file set_boundary.c.
References Gamma, gasconstant, inflow_density, inflow_temperature, and inflow_velocity.
Referenced by set_boundary().
int left_boundary_outflow | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
outflow boundary condition
Definition at line 973 of file set_boundary.c.
Referenced by set_boundary().
int left_boundary_periodic | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
periodic boundary conditions
Definition at line 1203 of file set_boundary.c.
Referenced by set_boundary().
int left_boundary_reflecting | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Reflecting boundary condition
Definition at line 640 of file set_boundary.c.
Referenced by set_boundary().
int left_boundary_reflecting_on_obstacle | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | pressure_solid_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Reflecting boundary condition on obstacle
Definition at line 694 of file set_boundary.c.
Referenced by set_boundary().
int left_boundary_small_padding | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Small number (e.g. 1E-50) ghost cell filling
Definition at line 865 of file set_boundary.c.
References small.
Referenced by set_boundary().
int left_boundary_small_padding_on_obstacle | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | pressure_solid_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
int left_boundary_zero_gradient | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Zero gradient boundary condition
Definition at line 475 of file set_boundary.c.
Referenced by set_boundary().
int make_ic | ( | int | x, |
int | y, | ||
int | z | ||
) |
An initial model is build in this function. If the parameter file defines strat_const_atmos a stratified atmosphere is build. If the parameter is set to zero a constant pressure/density distribution is realized.
Definition at line 19 of file make_ic.c.
References boundary::back, bound, boundary::down, boundary::front, gasconstant, grav_acc, inflow_velocity, boundary::left, pre, rho, boundary::right, small, starting_flow, strat_const_atmos, boundary::up, vx, vy, vz, x, y, ymax, and z.
Referenced by main().
int make_kh_instabilities | ( | int | x, |
int | y, | ||
int | z | ||
) |
A simple Kelvin Helmholtz Instability Setup
Definition at line 271 of file make_ic.c.
References dimension, pre, rho, vx, vy, vz, x, xmax, y, z, and zmax.
Referenced by main().
int make_sod_ic | ( | int | x, |
int | y, | ||
int | z | ||
) |
A Sod Shock Tube is generated here
Definition at line 196 of file make_ic.c.
References dimension, pre, rho, vx, vy, vz, x, y, and z.
Referenced by main().
long long int ntohll | ( | const long long int | data | ) |
converts double values to BigEndian VTK binaries are in BigEndian format as well as AMIRA MESH Files
Definition at line 21 of file big_endian_converter.c.
References x.
Referenced by write_amira(), write_ic_amira(), write_ic_vtk(), and write_vtk().
int para_coef | ( | int | nmin, |
int | nmax, | ||
int | flag1, | ||
double * | a_coef, | ||
double * | dx, | ||
double * | ai_coef, | ||
double * | b_coef, | ||
double * | bi_coef, | ||
double * | c_coef, | ||
double * | d_x, | ||
double ** | para, | ||
double * | ci_coef | ||
) |
Here the coefficients for the parabel functions are determined.
Definition at line 17 of file para_coef.c.
Referenced by ppm_step(), and remap().
int parabola | ( | int | nmin, |
int | nmax, | ||
double * | a, | ||
double * | deltaa, | ||
double * | a6, | ||
double * | al, | ||
int | flag1, | ||
double * | diffa, | ||
double * | da, | ||
double ** | para, | ||
double * | ar, | ||
double * | flat, | ||
double * | scratch1, | ||
double * | scratch2, | ||
double * | scratch3 | ||
) |
Here the picewise parabolic interpolated quantities for the PPM scheme are calculated.
Definition at line 45 of file parabola.c.
References min_par(), and signum_parabola().
Referenced by ppm_step(), and remap().
int ppm_free | ( | double ** | rho_1D, |
double ** | pre_1D, | ||
double ** | eng_1D, | ||
double ** | vx_1D, | ||
double ** | vy_1D, | ||
double ** | vz_1D, | ||
double ** | marker_1D, | ||
double ** | pressure_solid_1D, | ||
double ** | dx0, | ||
double ** | xa0, | ||
double ** | xa, | ||
double ** | dx, | ||
double ** | a_coef, | ||
double ** | ai_coef, | ||
double ** | b_coef, | ||
double ** | bi_coef, | ||
double ** | c_coef, | ||
double ** | ci_coef, | ||
double ** | d_x, | ||
double ** | da, | ||
double ** | ar, | ||
double ** | dp, | ||
double ** | dr, | ||
double ** | du, | ||
double ** | pl, | ||
double ** | p6, | ||
double ** | rl, | ||
double ** | r6, | ||
double ** | ul, | ||
double ** | u6, | ||
double ** | vl, | ||
double ** | v6, | ||
double ** | wl, | ||
double ** | w6, | ||
double ** | el, | ||
double ** | e6, | ||
double ** | ql, | ||
double ** | q6, | ||
double ** | deltaa, | ||
double ** | dv, | ||
double ** | dw, | ||
double ** | dq, | ||
double ** | de, | ||
double ** | scratch1, | ||
double ** | scratch2, | ||
double ** | scratch3, | ||
double ** | diffa, | ||
double ** | plft, | ||
double ** | prgh, | ||
double ** | ulft, | ||
double ** | urgh, | ||
double ** | rlft, | ||
double ** | rrgh, | ||
double ** | Cdtdx, | ||
double ** | fCdtdx, | ||
double ** | clft, | ||
double ** | crgh, | ||
double ** | plfti, | ||
double ** | prghi, | ||
double ** | pmid, | ||
double ** | pmold, | ||
double ** | wlft, | ||
double ** | wrgh, | ||
double ** | zlft, | ||
double ** | zrgh, | ||
double ** | umidl, | ||
double ** | umidr, | ||
double ** | umid, | ||
double ** | dm, | ||
double ** | dtbdm, | ||
double ** | upmid, | ||
double ** | xa1, | ||
double ** | xa2, | ||
double ** | xa3, | ||
double ** | vx_1D_old, | ||
double ** | e_int_1D, | ||
double ** | dvol, | ||
double ** | dvol0, | ||
double ** | dvol1, | ||
double ** | delta, | ||
double ** | fluxr, | ||
double ** | fluxu, | ||
double ** | fluxv, | ||
double ** | fluxw, | ||
double ** | fluxe, | ||
double ** | fluxq, | ||
double ** | dm0, | ||
double ** | steep, | ||
double ** | flat, | ||
double *** | para, | ||
double ** | rhodown, | ||
double ** | rhoup, | ||
double ** | rhofront, | ||
double ** | rhoback, | ||
double ** | vxdown, | ||
double ** | vxup, | ||
double ** | vxfront, | ||
double ** | vxback, | ||
double ** | vydown, | ||
double ** | vyup, | ||
double ** | vyfront, | ||
double ** | vyback, | ||
double ** | vzdown, | ||
double ** | vzup, | ||
double ** | vzfront, | ||
double ** | vzback, | ||
int | dimension | ||
) |
all quantities within the hydro part are freed in this function
Definition at line 174 of file init.c.
References max_array_length.
Referenced by sweep_x(), sweep_y(), and sweep_z().
int ppm_step | ( | int | i, |
int | j, | ||
int | k, | ||
int | direction, | ||
int | flag, | ||
int | nmin, | ||
int | nmax, | ||
double * | a_coef, | ||
double * | ai_coef, | ||
double * | b_coef, | ||
double * | bi_coef, | ||
double * | c_coef, | ||
double * | ci_coef, | ||
double * | d_x, | ||
double * | diffa, | ||
double * | da, | ||
double * | ar, | ||
double * | pl, | ||
double * | p6, | ||
double * | rl, | ||
double * | r6, | ||
double * | u6, | ||
double * | ul, | ||
double * | vl, | ||
double * | v6, | ||
double * | wl, | ||
double * | w6, | ||
double * | el, | ||
double * | e6, | ||
double * | ql, | ||
double * | q6, | ||
double * | dp, | ||
double * | du, | ||
double * | dr, | ||
double * | dv, | ||
double * | dw, | ||
double * | dq, | ||
double * | de, | ||
double * | scratch1, | ||
double * | scratch2, | ||
double * | scratch3, | ||
double * | plft, | ||
double * | prgh, | ||
double * | ulft, | ||
double * | urgh, | ||
double * | rlft, | ||
double * | rrgh, | ||
double * | Cdtdx, | ||
double * | fCdtdx, | ||
double * | steep, | ||
double * | flat, | ||
double ** | para, | ||
double * | clft, | ||
double * | crgh, | ||
double * | plfti, | ||
double * | prghi, | ||
double * | pmid, | ||
double * | pmold, | ||
double * | wlft, | ||
double * | wrgh, | ||
double * | zlft, | ||
double * | zrgh, | ||
double * | umidl, | ||
double * | umidr, | ||
double * | umid, | ||
double * | dm, | ||
double * | dtbdm, | ||
double * | upmid, | ||
double * | xa1, | ||
double * | xa2, | ||
double * | xa3, | ||
double * | vx_1D_old, | ||
double * | dvol, | ||
double * | dvol0, | ||
double * | dvol1, | ||
double * | delta, | ||
double * | fluxr, | ||
double * | fluxu, | ||
double * | fluxv, | ||
double * | fluxw, | ||
double * | fluxe, | ||
double * | fluxq, | ||
double * | dm0, | ||
double * | e_int_1D, | ||
double * | rho_1D, | ||
double * | pre_1D, | ||
double * | eng_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | marker_1D, | ||
double * | pressure_solid_1D, | ||
double * | dx0, | ||
double * | xa0, | ||
double * | xa, | ||
double * | dx, | ||
int | bound_checker, | ||
int | lefter, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
The inner PPM Hydro scheme. First the relevant 1D distributions from the computational domain is given to this function. The scheme is as followed: First fill the ghost cells according to the defined boundary conditions. Solve the diffusion equation if needed. Calculate the parabolic function coefficients, then the volumes and the flatten coefficients due toe the steepness of the neighboring cells. Now the parabolic interpolated quantities are determined. The quantities at the interpolated points are given. The Riemann solver solves now the shock-tube problem. Then the hydrodynamic quantities are evolved and finally re-mapped on the computational grid.
Definition at line 24 of file ppm_step.c.
References advect(), advection, evolve(), flatten(), para_coef(), parabola(), remap(), riemann(), set_boundary(), states(), viscosity(), and volume().
Referenced by sweep_x(), sweep_y(), and sweep_z().
int pre_old_copy | ( | int | x, |
int | y, | ||
int | z | ||
) |
Here the present pressure is copied into the pre_old
Definition at line 21 of file pre_old_copy.c.
References pre, pre_old, x, y, and z.
Referenced by main().
int prepare_the_dB_map | ( | int | x, |
int | y, | ||
int | z | ||
) |
Definition at line 666 of file noise_generator.c.
References dB_map, dB_Map_ready, dt_integrated, pressure_integrated, with_one_pulse, x, y, and z.
Referenced by main().
int pressure_on_solid_calc | ( | int | x, |
int | y, | ||
int | z | ||
) |
Here the total force due to the pressure on the solid is determined.
Definition at line 124 of file data_explorer.c.
References dimension, dom, output_dir, pressure_on_solid, time_sim, x, xmax, and ymax.
Referenced by main().
int read_dom | ( | int | x, |
int | y, | ||
int | z | ||
) |
Definition at line 129 of file read_ic.c.
References dom, obst_ic, x, y, and z.
Referenced by main().
int read_ic | ( | int | x, |
int | y, | ||
int | z | ||
) |
If you want to read in initial conditions from an native TYCHO file.
Definition at line 17 of file read_ic.c.
References counter, counter_restart, dens_ic, dom, gasconstant, intial_velocity_file, obst_ic, pre, rho, temp_ic, time_sim, vel_ic, vx, vy, vz, with_obstacles, x, y, and z.
Referenced by main().
int read_marker_file | ( | int | x, |
int | y, | ||
int | z | ||
) |
Definition at line 186 of file read_ic.c.
References marker, marker_ic, x, y, and z.
Referenced by main().
int read_restart | ( | int | x, |
int | y, | ||
int | z | ||
) |
Read restart from an native TYCHO file.
Definition at line 17 of file read_restart.c.
References advection, counter, counter_restart, dom, gasconstant, marker, output_dir, pre, pressure_integrated, rho, soundemitter, time_sim, vx, vy, vz, wind_marker, wind_on_off, with_obstacles, with_sound_emitter, x, y, and z.
Referenced by main().
int remap | ( | int | nmin, |
int | nmax, | ||
int | flag, | ||
double * | a_coef, | ||
double * | dx, | ||
double * | ai_coef, | ||
double * | b_coef, | ||
double * | bi_coef, | ||
double * | c_coef, | ||
double * | d_x, | ||
double ** | para, | ||
double * | ci_coef, | ||
double * | dr, | ||
double * | r6, | ||
double * | rl, | ||
double * | diffa, | ||
double * | da, | ||
double * | ar, | ||
double * | flat, | ||
double * | scratch1, | ||
double * | scratch2, | ||
double * | scratch3, | ||
double * | du, | ||
double * | u6, | ||
double * | ul, | ||
double * | dv, | ||
double * | v6, | ||
double * | vl, | ||
double * | w6, | ||
double * | wl, | ||
double * | dq, | ||
double * | q6, | ||
double * | ql, | ||
double * | de, | ||
double * | e6, | ||
double * | el, | ||
double * | xa, | ||
double * | xa0, | ||
double * | delta, | ||
double * | fluxr, | ||
double * | fluxu, | ||
double * | fluxv, | ||
double * | fluxw, | ||
double * | dw, | ||
double * | fluxe, | ||
double * | fluxq, | ||
double * | dm, | ||
double * | rho_1D, | ||
double * | dvol, | ||
double * | dvol0, | ||
double * | dm0, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | eng_1D, | ||
double * | e_int_1D, | ||
double * | pre_1D | ||
) |
Remap the new hydro-quantities on the computational domain grid.
Definition at line 30 of file remap.c.
References Gamma, max_remap(), para_coef(), parabola(), smallp, and smallr.
Referenced by ppm_step().
int reset_pressure_integrated | ( | int | x, |
int | y, | ||
int | z | ||
) |
Definition at line 634 of file noise_generator.c.
References pressure_integrated, x, y, and z.
Referenced by main().
int riemann | ( | int | nmin, |
int | nmax, | ||
double * | clft, | ||
double * | crgh, | ||
double * | rlft, | ||
double * | rrgh, | ||
double * | plfti, | ||
double * | prghi, | ||
double * | pmid, | ||
double * | pmold, | ||
double * | plft, | ||
double * | wrgh, | ||
double * | prgh, | ||
double * | wlft, | ||
double * | zlft, | ||
double * | zrgh, | ||
double * | umidl, | ||
double * | umidr, | ||
double * | umid, | ||
double * | urgh, | ||
double * | ulft | ||
) |
The Riemann solver for the shock tube problem.
Input variables are: lmin = zone number of first physical zone lmax = zone number of first ghost zone on right (lmax=nmax+1) gamma = equation of state gamma prgh = pressure state on the right side of the boundary plft = pressure state on the left side of the boundary urgh = velocity state on the right side of the boundary ulft = velocity state on the left side of the boundary vrgh = density state on the right side of the boundary rlft = density state on the left side of the boundary (rlft and vrgh are inverted to get the specific volume)
Output variables are: umid = time_sim averaged velocity at the zone interface pmid = time_sim averaged pressure at the zone interface
Definition at line 46 of file riemann.c.
References Gamma1, max_riemann(), and smallp.
Referenced by ppm_step().
int right_boundary_inflow | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
outflow boundary condition
Definition at line 1147 of file set_boundary.c.
References Gamma, gasconstant, inflow_density, inflow_temperature, and inflow_velocity.
Referenced by set_boundary().
int right_boundary_outflow | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
outflow boundary condition
Definition at line 1032 of file set_boundary.c.
Referenced by set_boundary().
int right_boundary_periodic | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
periodic boundary conditions
Definition at line 1259 of file set_boundary.c.
Referenced by set_boundary().
int right_boundary_reflecting | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Reflecting boundary condition
Definition at line 752 of file set_boundary.c.
Referenced by set_boundary().
int right_boundary_reflecting_on_obstacle | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | pressure_solid_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Reflecting boundary condition on obstacle
Definition at line 807 of file set_boundary.c.
Referenced by set_boundary().
int right_boundary_small_padding | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Small number (e.g. 1E-50) ghost cell filling
Definition at line 919 of file set_boundary.c.
References small.
Referenced by set_boundary().
int right_boundary_zero_gradient | ( | int | nmin, |
int | nmax, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
int | flag, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
Zero gradient boundary condition or if stratisfied atmosphere needed the atmosphere is extened into the ghost cells.
Definition at line 531 of file set_boundary.c.
References gasconstant, grav_acc, strat_const_atmos, y, and ymax.
Referenced by set_boundary().
int set_boundary | ( | int | nmin, |
int | nmax, | ||
int | flag, | ||
int | bound_checker, | ||
int | lefter, | ||
double * | rho_1D, | ||
double * | eng_1D, | ||
double * | pre_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | pressure_solid_1D, | ||
double * | xa0, | ||
double * | dx0, | ||
double * | xa, | ||
double * | dx, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | viscosity_on_off, | ||
int | dimension | ||
) |
The filling routine for the ghost cells according to the given boundary conditions.
Definition at line 17 of file set_boundary.c.
References boundary::back, bound, boundary::down, boundary::front, boundary::left, left_boundary_inflow(), left_boundary_outflow(), left_boundary_periodic(), left_boundary_reflecting(), left_boundary_reflecting_on_obstacle(), left_boundary_small_padding(), left_boundary_zero_gradient(), boundary::right, right_boundary_inflow(), right_boundary_outflow(), right_boundary_periodic(), right_boundary_reflecting(), right_boundary_reflecting_on_obstacle(), right_boundary_small_padding(), right_boundary_zero_gradient(), and boundary::up.
Referenced by ppm_step().
int start_file_reader | ( | char | filename[] | ) |
The read routine for the parameterfile
Definition at line 19 of file start_file_reader.c.
References advection, boundary::back, bound, courant, dens_ic, dimension, boundary::down, endtime_sim, filetype, boundary::front, Gamma1, gasconstant, gravity_on_off, inflow_density, inflow_temperature, inflow_velocity, intial_velocity_file, boundary::left, make_ics, marker_density, marker_ic, minute_in_sec, number_of_threads, obst_ic, obstacle_density, obstacle_heat_conductivity, obstacle_temperature, obstalce_absorption_coefficient, output_dir, output_frequency, restart_frequency, boundary::right, sound_frequency, sound_ic, sound_pressure_level, sound_reflexion, specific_heat_capacity_gas, specific_heat_capacity_obstacle, starting_flow, strat_const_atmos, temp_ic, boundary::up, vel_ic, viscosity_on_off, wind_ic, wind_on_off, wind_speed, with_obstacles, with_one_pulse, with_sound_emitter, x, xmax, y, ymax, z, and zmax.
Referenced by main().
int states | ( | int | nmin, |
int | nmax, | ||
int | flag, | ||
double * | pre_1D, | ||
double * | rho_1D, | ||
double * | dx, | ||
double * | Cdtdx, | ||
double * | fCdtdx, | ||
double * | plft, | ||
double * | pl, | ||
double * | dp, | ||
double * | p6, | ||
double * | ulft, | ||
double * | ul, | ||
double * | du, | ||
double * | u6, | ||
double * | rlft, | ||
double * | rl, | ||
double * | dr, | ||
double * | r6, | ||
double * | prgh, | ||
double * | urgh, | ||
double * | rrgh | ||
) |
Get hydroquantities at cell centres and faces.
Definition at line 30 of file states.c.
References dt, Gamma1, grav_acc, gravity_on_off, max_states(), smallp, and smallr.
Referenced by ppm_step().
int sweep_x | ( | int | x, |
int | y, | ||
int | z, | ||
int | flag | ||
) |
The sweeping in all three directions (if 3D). This routine does the x-sweep and is parallelized using OpenMP for shared memory machines.
Definition at line 36 of file sweep_x.c.
References advection, boundary::back, bound, dimension, dom, DOM_FLUID, DOM_SOLID, boundary::down, boundary::front, Gamma, inflow_density, inflow_velocity, init_ppm(), marker, max_sweep_x(), obstacle_density, obstacle_temperature, ppm_free(), ppm_step(), pre, pressure_on_solid, rho, rho_visc, small, smallp, boundary::up, viscosity_on_off, vx, vx_visc, vy, vy_visc, vz, vz_visc, with_obstacles, x, y, z, zdx, and zxa.
Referenced by hydro_sweeps().
int sweep_y | ( | int | x, |
int | y, | ||
int | z, | ||
int | flag | ||
) |
The sweeping in all three directions (if 3D). This routine does the y-sweep and is parallelized using OpenMP for shared memory machines.
Definition at line 35 of file sweep_y.c.
References advection, boundary::back, bound, dimension, dom, DOM_FLUID, DOM_SOLID, boundary::front, Gamma, inflow_density, inflow_velocity, init_ppm(), boundary::left, marker, max_sweep_y(), obstacle_density, obstacle_temperature, ppm_free(), ppm_step(), pre, pressure_on_solid, rho, rho_visc, boundary::right, small, smallp, viscosity_on_off, vx, vx_visc, vy, vy_visc, vz, vz_visc, with_obstacles, x, y, z, zdy, and zya.
Referenced by hydro_sweeps().
int sweep_z | ( | int | x, |
int | y, | ||
int | z, | ||
int | flag | ||
) |
The sweeping in all three directions (if 3D). This routine does the z-sweep and is parallelized using OpenMP for shared memory machines.
Definition at line 35 of file sweep_z.c.
References advection, bound, dimension, dom, DOM_FLUID, DOM_SOLID, boundary::down, Gamma, inflow_density, inflow_velocity, init_ppm(), boundary::left, marker, max_sweep_z(), obstacle_density, obstacle_temperature, ppm_free(), ppm_step(), pre, pressure_on_solid, rho, rho_visc, boundary::right, small, smallp, boundary::up, viscosity_on_off, vx, vx_visc, vy, vy_visc, vz, vz_visc, with_obstacles, x, y, z, zdz, and zza.
Referenced by hydro_sweeps().
int upper_atmos | ( | int | x, |
int | y, | ||
int | z | ||
) |
int velocity_field_analyser | ( | int | x, |
int | y, | ||
int | z | ||
) |
Quantities of the fluid are integrated here. Here you can start to implement you owns investigations of integrated quantities of the fluid.
Definition at line 25 of file data_explorer.c.
References dimension, dom, output_dir, rho, time_sim, vx, vy, vz, x, y, and z.
Referenced by main().
int viscosity | ( | int | i, |
int | j, | ||
int | k, | ||
int | flag, | ||
int | nmin, | ||
int | nmax, | ||
double * | rho_1D, | ||
double * | vx_1D, | ||
double * | vy_1D, | ||
double * | vz_1D, | ||
double * | pre_1D, | ||
double * | e_int_1D, | ||
double * | eng_1D, | ||
int | lefter, | ||
double * | rhodown, | ||
double * | rhoup, | ||
double * | rhofront, | ||
double * | rhoback, | ||
double * | vxdown, | ||
double * | vxup, | ||
double * | vxfront, | ||
double * | vxback, | ||
double * | vydown, | ||
double * | vyup, | ||
double * | vyfront, | ||
double * | vyback, | ||
double * | vzdown, | ||
double * | vzup, | ||
double * | vzfront, | ||
double * | vzback, | ||
int | dimension | ||
) |
The velocity fluxes are changed due to viscosity Velocity-Flux = -nu * grad(v_x) This routine computes the velocity fluxes from viscosity and alters the velocity fluxes from the hydro step.
Definition at line 34 of file viscosity.c.
References dt, dynamic_viscosity(), gasconstant, max_array_length, x, xmax, xmin, y, ymax, ymin, z, zmax, and zmin.
Referenced by ppm_step().
int volume | ( | int | nmin, |
int | nmax, | ||
double * | vol, | ||
double * | dx, | ||
double * | vol0, | ||
double * | dx0 | ||
) |
Definition at line 16 of file volume.c.
Referenced by diffusion(), and ppm_step().
int wind | ( | int | x, |
int | y, | ||
int | z | ||
) |
The Wind routine, which gives a constant flow with wind-speed velocity
Definition at line 39 of file wind.c.
References vx, vy, vz, wind_marker, wind_speed, x, y, and z.
Referenced by hydro_sweeps().
int write_amira | ( | int | x, |
int | y, | ||
int | z, | ||
int | counter | ||
) |
Writes the initial model in a Amira Mesh format those binaries are in BigEndian format The slow order of the nested loops is mandatory here
Definition at line 20 of file write_amira.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, ntohll(), output_dir, pre, pressure_on_solid, rho, vx, vy, vz, with_obstacles, with_sound_emitter, x, xmax, xmin, y, ymax, ymin, z, zmax, and zmin.
Referenced by main().
int write_ic_amira | ( | int | x, |
int | y, | ||
int | z | ||
) |
Writes the initial model in a Amira Mesh format those binaries are in BigEndian format The slow order of the nested loops is mandatory here
Definition at line 20 of file write_ic_amira.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, ntohll(), output_dir, pre, pressure_on_solid, rho, vx, vy, vz, with_obstacles, with_sound_emitter, x, xmax, xmin, y, ymax, ymin, z, zmax, and zmin.
Referenced by main().
int write_ic_ifrit | ( | int | x, |
int | y, | ||
int | z | ||
) |
Writes the initial model in a IFRIT uniform scale binary data file The slow order of the nested loops is mandatory here
Definition at line 20 of file write_ic_ifrit.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, output_dir, pre, pressure_on_solid, rho, vx, vy, vz, with_obstacles, with_sound_emitter, x, y, and z.
Referenced by main().
int write_ic_tyc | ( | int | x, |
int | y, | ||
int | z | ||
) |
Writes the initial model
Definition at line 17 of file write_ic_tyc.c.
References advection, counter, counter_restart, dB_map, dB_Map_ready, dom, gasconstant, marker, output_dir, pre, pressure_on_solid, rho, time_sim, vx, vy, vz, with_obstacles, with_sound_emitter, x, y, and z.
Referenced by main().
int write_ic_vtk | ( | int | x, |
int | y, | ||
int | z | ||
) |
Writes the initial model in a RECTILINEAR_GRID VTK-File VTK binaries are in BigEndian format The slow order of the nested loops is mandatory here
Definition at line 20 of file write_ic_vtk.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, minute_in_sec, ntohll(), output_dir, pre, pressure_on_solid, rho, time_sim, vx, vy, vz, with_obstacles, with_sound_emitter, x, xmax, y, ymax, z, and zmax.
Referenced by main().
int write_ifrit | ( | int | x, |
int | y, | ||
int | z, | ||
int | counter | ||
) |
Writes the initial model in a IFRIT uniform scala binary data file The slow order of the nested loops is mandatory here
Definition at line 20 of file write_ifrit.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, output_dir, pre, pressure_on_solid, rho, vx, vy, vz, with_obstacles, with_sound_emitter, x, y, and z.
Referenced by main().
int write_restart_tyc | ( | int | x, |
int | y, | ||
int | z | ||
) |
Writes the restart files
Definition at line 17 of file write_restart_tyc.c.
References advection, counter, counter_restart, dom, gasconstant, marker, output_dir, pre, pressure_integrated, rho, soundemitter, time_sim, vx, vy, vz, wind_marker, wind_on_off, with_obstacles, with_sound_emitter, x, y, and z.
Referenced by main().
int write_tyc | ( | int | x, |
int | y, | ||
int | z, | ||
int | counter | ||
) |
write_raw.c
Author: Wolfgang Kapferer
In this routine the output-files are written.
Definition at line 17 of file write_tyc.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, output_dir, pre, pressure_on_solid, rho, time_sim, vx, vy, vz, with_obstacles, with_sound_emitter, x, y, and z.
Referenced by main().
int write_vtk | ( | int | x, |
int | y, | ||
int | z, | ||
int | counter | ||
) |
Writes the initial model in a RECTILINEAR_GRID VTK-File VTK binaries are in BigEndian format The slow order of the nested loops is mandatory here
Definition at line 21 of file write_vtk.c.
References advection, dB_map, dB_Map_ready, dom, gasconstant, marker, minute_in_sec, ntohll(), output_dir, pre, pressure_on_solid, rho, time_sim, vx, vy, vz, with_obstacles, with_sound_emitter, x, xmax, y, ymax, z, and zmax.
Referenced by main().