1 pv . vdV, where p (x, y, z) and v (x, y, z) are the pointwise fluid density and velocity, respectively. Fluid with uniform density 1 Kinetic energy of a fluid flow can be computed by flows in the domain bounded by x? + z? = 9 and 0

(a) For a circular piston, plot the on-axis pressure amplitude as a function of scaled distance r/a for several values of ka between 3 and 12. (b) Plot the range beyond which the pressure amplitude is within 10% of the asymptotic form (7.4.7). (c) For a piston of 20 cm radius operating at 4 kHz in water, find the distance corresponding to (b).

% ParametersD = 0.1; % Diameter of the tube (m)L = 1.0; % Length of the tube bundle (m)N = 8; % Number of tubes in the bundleU = 1.0; % Inlet velocity (m/s)rho = 1.2; % Density of the fluid (kg/m^3)mu = 0.01; % Dynamic viscosity of the fluid (Pa.s) % Define the grid size and time stepdx = D/10; % Spatial step size (m)dy = L/10; % Spatial step size (m)dt = 0.01; % Time step size (s) % Calculate the number of grid points in each directionnx = ceil(D/dx) + 1;ny = ceil(L/dy) + 1; % Create the velocity matrixU_matrix = U * ones(nx, ny); % Perform the iterationsfor iter = 1:100    % Calculate the velocity gradients    dUdx = (U_matrix(:, 2:end) - U_matrix(:, 1:end-1)) / dx;    dUdy = (U_matrix(2:end, :) - U_matrix(1:end-1, :)) / dy;        % Calculate the pressure gradients    dpdx = -mu * dUdx;    dpdy = -mu * dUdy;        % Calculate the change in velocity    dU = dt * (dpdx / rho);        % Update the velocity matrix    U_matrix(:, 2:end-1) = U_matrix(:, 2:end-1) + dU;        % Apply…

Need help with the calculations i keep getting lost in the conversions