def soln_specific_internal_energy(face_mesh_velocity,
outward_directed_normal_covector, coords,
time, dim):
return hydro.specific_internal_energy(coords, time,
_soln_mean_velocity(dim),
_soln_wave_vector(dim),
_soln_pressure,
_soln_adiabatic_index,
_soln_perturbation_size)
def soln_energy_density(face_mesh_velocity, outward_directed_normal_covector,
coords, time, dim):
velocity = hydro.spatial_velocity(coords, time, _soln_mean_velocity(dim),
_soln_wave_vector(dim), _soln_pressure,
_soln_adiabatic_index,
_soln_perturbation_size)
int_energy = hydro.specific_internal_energy(coords, time,
_soln_mean_velocity(dim),
_soln_wave_vector(dim),
_soln_pressure,
_soln_adiabatic_index,
_soln_perturbation_size)
return hydro.rest_mass_density(
coords, time, _soln_mean_velocity(dim), _soln_wave_vector(dim),
_soln_pressure, _soln_adiabatic_index, _soln_perturbation_size) * (
0.5 * np.dot(velocity, velocity) + int_energy)
def soln_specific_internal_energy(coords, time):
return hydro.specific_internal_energy(coords, time, _soln_mean_velocity(),
_soln_wave_vector(), _soln_pressure,
_soln_adiabatic_index,
_soln_perturbation_size)