CapacitorSkinEffect

Determines the skin effect current density for a capacitor discharge

Input values that describe a capacitor discharge function. This code will output the heterogeneous current density at specific query points of time and will also create a video to demonstrate how the current flow behaves during the pulse.

This code uses an equation derived from literature:

Haines, Malcolm G. "The inverse skin effect." Proceedings of the Physical Society (1958-1967) 74.5 (1959): 576.

[SkinEffect] = CapacitorSkinEffect(t_Peak_Desired, PulseWidth_Desired, PulseAmplitude_Desired, Pulse_Off_Desired, use_Parallel, r, QueryPoints, num_Frames, Video_on_off, VideoName)

% Finds the variables a, b, and c for the following equation based on desired parameters:

% y(x) = c*(exp(-a*x)-exp(-b*x))

% t_Peak_Desired: time at which Y is max

% PulseWidth_Desired: time at which Y decays to the Pulse Off limit

% PulseAmplitude_Desired: Amplitude of maximum Y value

% Pulse_Off_Desired: Pulse off limit set as percentage of maximum

% amplitude - possible range: 0 - 1

% use_Parallel: 1=use parfor; 0=do not use parfor

% If nothing is plotted, t_peak is likely too close to pulse width. Reduce

% t_peak or increase pulse width

% QueryPoints = matrix of time points to determine skin effect for (should

% be between 0-PulseWidth_Desired)

% r = radius of conductor

% Video_on_off = 1 - save video

% Video_on_off = 0 - do not save video

% VideoName = file name for video file