research/conditional-analysis-dns
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

implemented binning 2nd order dervative

Browse source
This commit is contained in:
ignis 2021-08-12 03:50:07 +09:00
parent 44ef8da3b9
commit dab488759b
1 changed files with 5 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

17
binning_ddxc_d2dxc_given_c_from_raw.py
View file

@ -121,23 +121,21 @@ cs = CompactScheme(nx, ny, nz, False, True, True, 4, 2, 2)
# c mean at sampling x coordinates (verification purpose)
cmean_verify = np.zeros(len(x_indices))
avg_objs_u = [ConditionalMeanBinning(name=str(xi)) for xi in x_indices]
avg_objs_ddxc = [ConditionalMeanBinning(name=str(xi)) for xi in x_indices]
avg_objs_uddxc = [ConditionalMeanBinning(name=str(xi)) for xi in x_indices]
avg_objs_d2dxc = [ConditionalMeanBinning(name=str(xi)) for xi in x_indices]
for i, fname in enumerate(case.data_files[:n_timestep]):
print(datetime.datetime.now(), '{:5d}/{}'.format(i+1, n_timestep), fname)
fpath = os.path.join(case.case_root, fname)
time, u = case.read_single_field(fpath, u_idx)
time, y = case.read_single_field(fpath, y_idx)
c = (1. - y)
ddxc = cs.ddx(c)
d2dxc = cs.d2dx(c)
for j, xi in enumerate(x_indices):
avg_objs_u[j].feed_samples(c[:,:,xi], u[:,:,xi])
avg_objs_ddxc[j].feed_samples(c[:,:,xi], ddxc[:,:,xi])
avg_objs_uddxc[j].feed_samples(c[:,:,xi], u[:,:,xi]*ddxc[:,:,xi])
avg_objs_d2dxc[j].feed_samples(c[:,:,xi], d2dxc[:,:,xi])
cmean_verify[j] += np.sum(c[:,:,xi])
print(datetime.datetime.now(), 'Finished')
@ -149,17 +147,12 @@ for xi_cmean_pair in zip(x_indices, cmean_verify):
print()
print("Save the result")
for j, xi in enumerate(x_indices):
filename = "cavg_u_c_{}_{}".format(casename, avg_objs_u[j].name)
np.savez(filename, **avg_objs_u[j].average())
print("saved ", filename)
for j, xi in enumerate(x_indices):
filename = "cavg_ddxc_c_{}_{}".format(casename, avg_objs_ddxc[j].name)
np.savez(filename, **avg_objs_ddxc[j].average())
print("saved ", filename)
for j, xi in enumerate(x_indices):
filename = "cavg_uddxc_c_{}_{}".format(casename, avg_objs_uddxc[j].name)
np.savez(filename, **avg_objs_uddxc[j].average())
filename = "cavg_d2dxc_c_{}_{}".format(casename, avg_objs_d2dxc[j].name)
np.savez(filename, **avg_objs_d2dxc[j].average())
print("saved ", filename)