An algorithm for reconstructing transient rocket chamber pressures from attenuated time histories is presented. This method is presented as an alternative to traditional methods, in which motor pressures are sensed by transducers mounted close-coupled to the motor casing. This close-coupled installation presents several significant measurement issues, including potential structural weakening of chamber walls, overheating of the pressure sensing element, motor-case strain biasing of the sensing element, and resonance within the measurement port. A less complex, and consequently less risky, installation is to tap the case at desired locations using very small pressure ports and then transmit pressure from the port to a pressure transducer using a significant length of pneumatic transmission tubing. This installation allows the transducer to be mounted in a controlled environment and virtually eliminates any sensing errors due to motor-case strain. Very small pressure taps are far easier to integrate with motor-case walls and allow multiple longitudinal measurements to be obtained without significantly reducing the structural integrity of the motor case. The method, based on optimal deconvolution theory, amplifies attenuated pressure signals while rejecting additive noise. The method is validated using laboratory-derived data and then applied to reconstructing transient chamber pressure for a small-scale solid rocket motor.
Nomenclature= characteristic velocity, m=s c = sonic velocity, m=s D = tube diameter, mm E = expectation operator G = optimal deconvolution function i = frequency index J = cost functional j = complex constant, 1 1=2 k = time index L = tube length, cm M = number of data points in record N = frequency-domain noise, kPa=s kNk 2 = measurement variance, kPa 2 n = discrete-time advance index P c = chamber pressure, kPa P L = frequency-domain response (sensed) pressure, kPa kP LL k 2 = measured pressure variance, kPa 2 . P 0 = frequency-domain input pressure to pneumatic tubing, kPa kP 0L k = covariance, kPa 2 . kP 00 k 2 = measured pressure variance, kPa 2 . p L = time-domain response pressure, kPa p 0 = time-domain input pressure to pneumatic tubing, kPatransducer volume, cm 3 x = longitudinal coordinate, cm z = discrete-frequency response variable, 1=s t = sample interval, s p = wave propagation factor = ratio of specific heats = damping ratio = time-domain noise function, kPa = dynamic viscosity, N s=m 2 = time-domain convolution function = polytropic exponent = fluid density, kg=m 3 0 = mean fluid density in tube, kg=m 3 = time lag or dummy variable, s, = frequency-domain convolution function ! = angular frequency, rad=s ! n = radian natural frequency, rad=ŝ = optimal estimate Subscripts I = imaginary component of complex number R = real component of complex number Superscript = complex conjugate