THE UNIVERSITY OF MICHIGAN 2871-4-T NON-LINEAR MODELING OF MAXWELL'S EQUATIONS by J. E.\elyea, R. D. Low, and K. M. Siegel December 1959 Report No. 2871-4-T on Contract AF 19(604)-4993 Prepared for ELECTRONICS RESEARCH DIRECTORATE AIR FORCE CAMBRIDGE RESEARCH CENTER AIR RESEARCH AND DEVELOPMENT COMMAND UNITED STATES AIR FORCE BEDFORD, MASSACHUSETTS

THE UNIVERSITY OF MICHIGAN 2871-4-T Requests for additional copies by Agencies of the Department of Defense, their contractors, and other Government agencies should be directed to the: ARMED SERVICES TECHNICAL INFORMATION AGENCY ARLINGTON HALL STATION ARLINGTON 12, VIRGINIA Department of Defense contractors must be established for ASTIA services or have their "need-to-know" certified by the cognizant military agency of their project or contract. All other persons and organizations should apply to the: U. S. DEPARTMENT OF COMMERCE OFFICE OF TECHNICAL SERVICES WASHINGTON 25, D. C. ii

THE UNIVERSITY OF MICHIGAN 2871-4-T STUDIES IN RADAR CROSS SECTIONS I "Scattering by a Prolate Spheroid", F.V. Schultz (UMM-42, March 1950), W-33(038)-ac-14222. UNCLASSIFIED. 65 pgs. m II "The Zeros of the Associated Legendre Functions P (p') of Non-Integral Degree", K.M. Siegel, D.M. Brown, H.E. Hunter, H.A. Alperin and C.W. Quillen (UMM-82, April 1951), W-33(038)-ac-14222. UNCLASSIFIED. 20 pgs. III "Scattering by a Cone", K. M. Siegel and H.A. Alperin (UMM-87, January 1952), AF-30(602)-9. UNCLASSIFIED. 56 pgs. IV "Comparison between Theory and Experiment of the Cross Section of a Cone", K.M. Siegel, H.A. Alperin, J.W. Crispin, Jr., H.E. Hunter, R. E. Kleinman, W.C. Orthwein and C. E. Schensted (UMM-92, February 1953), AF-30(602)-9. UNCLASSIFIED. 70 pgs. V "An Examination of Bistatic Early Warning Radars", K. M. Siegel (UMM-98, August 1952), W-33(038)-ac-14222. SECRET. 25 pgs. VI "Cross Sections of Corner Reflectors and Other Multiple Scatterers at Microwave Frequencies", R.R. Bonkowski, C. R. Lubitz and C. E. Schensted (UMM-106, October 1953), AF-30(602)-9. SECRET - Unclassified when appendix is removed. 63 pgs. VII "Summary of Radar Cross Section Studies under Project Wizard", K.M. Siegel, J.W. Crispin,Jr. and R.E. Kleinman (UMM-108, November 1952), W-33(038)-ac-14222. SECRET. 75 pgs. VIII "Theoretical Cross Section as a Function of Separation Angle between Transmitter and Receiver at Small Wavelengths", K. M. Siegel, H.A. Alperin, R.R. Bonkowski, J.W. Crispin, Jr., A.L. Maffett, C.E. Schensted and I. V. Schensted (UMM-115, October 1953), W-33(038)-ac14222. UNCLASSIFIED. 84 pgs. IX "Electromagnetic Scattering by an Oblate Spheroid", L. M. Rauch UMM-116, October 1953), AF-30(602)-9. UNCLASSIFIED. 38 pgs. iii

THE UNIVERSITY OF MICHIGAN 2871-4-T X "Scattering of Electromagnetic Waves by Spheres", H. Weil, M.L. Barasch and T. A. Kaplan (2255-20-T, July 1956), AF-30(602)-1070. UNCLASSIFIED. 104 pgs. XI "The Numerical Determination of the Radar Cross Section of a Prolate Spheroid", K. M. Siegel, B. H. Gere, I. Marx and F. B. Sleator (UMM-126, December 1953), AF-30(602)-9. UNCLASSIFIED. 75 pgs. XII "Summary of Radar Cross Section Studies under Project MIRO", K. M. Siegel, M.E. Anderson, R.R. Bonkowski and W.C. Orthwein (UMM127, December 1953), AF-30(602)-9. SECRET. 90 pgs. XIII "Description of a Dynamic Measurement Program", K. M. Siegel and J. M. Wolf (UMM-128, May 1954), W-33(038)-ac-14222. CONFIDENTIAL. 152 pgs. XIV "Radar Cross Section of a Ballistic Missile "' K. M. Siegel, M. L. Barasch, J.W. Crispin, Jr., W.C. Orthwein, I.V. Schensted and H. Weil (UMM134, September 1954), W-33(038)-ac-14222. SECRET. 270 pgs. XV "Radar Cross Sections of B-47 and B-52 Aircraft", C. E. Schensted, J.W. Crispin, Jr. and K.M. Siegel (2260-1-T, August 1954), AF-33(616)2531. CONFIDENTIAL. 155 pgs. XVI "Microwave Reflection Characteristics of Buildings", H. Weil, R. R. Bonkowski, T.A. Kaplan and M. Leichter (2255-12-T, May 1955), AF-30(602)-1070. SECRET. 148 pgs. XVII "Complete Scattering Matrices and Circular Polarization Cross Sections for the B-47 Aircraft at S-band", A. L. Maffett, M. L. Barasch, W. E. Burdick, R. F. Goodrich, W.C. Orthwein, C. E. Schensted and K. M. Siegel (2260-6-T, June 1955), AF-33(616)-2531. CONFIDENTIAL. 157 pgs. XVIII "Airborne Passive Measures and Countermeasures", K.M. Siegel, M. L. Barasch, J. W. Crispin, Jr., R. F. Goodrich, A. H. Halpin, A. L. Maffett, W.C. Orthwein, C. E. Schensted and C.J. Titus (226029-F, January 1956), AF-33(616)-2531. SECRET. 177 pgs. XIX "Radar Cross Section of a Ballistic Missile II", K. M. Siegel, M. L. Barasch, H. Brysk, J.W. Crispin, Jr., T.B. Curtz and T.A. Kaplan (2428-3-T, January 1956), AF-04(645)-33. SECRET. 189 pgs. iv

THE UNIVERSITY OF MICHIGAN 2871-4-T XX "Radar Cross Section of Aircraft and Missiles", K.M. Siegel, W.E. Burdick, J.W. Crispin, Jr. and S. Chapman (WR-31-J, March 1956), SECRET. 151 pgs. XXI "Radar Cross Section of a Ballistic Missile III", K.M. Siegel, H. Brysk, J.W. Crispin, Jr. and R.E. Kleinman (2428-19-T, October 1956), AF-04(645)-33. SECRET. 125 pgs. XXII "Elementary Slot Radiators", R. F. Goodrich, A.L. Maffett, N.E. Reitlinger, C. E. Schensted and K. M. Siegel (2472-13-T, November 1956), AF-33(038)-28634, HAC-PO L-265165-F31. UNCLASSIFIED. 100 pgs. XXIII "A Variational Solution to the Problem of Scalar Scattering by a Prolate Spheroid", F.B. Sleator (2591-1-T, March 1957), AF-19(604)-1949, AFCRC-TN-57-586, AD 133631. UNCLASSIFIED. 67 pgs. XXTV "Radar Cross Section of a Ballistic Missile - IV The Problem of Defense", M.L. Barasch, H. Brysk, J.W. Crispin, Jr., B.A. Harrison, T.B.A. Senior, K.M. Siegel, H. Weil and V.H. Weston (2778-1-F, April 1959), AF-30(602)-1853. SECRET. 362 pgs. XXV "Diffraction by an Imperfectly Conducting Wedge", T. B.A. Senior (2591-2-T, October 1957), AF-19(604)-1949, AFCRC-TN-57-591, AD 133746. UNCLASSIFIED. 71 pgs. XXVI "Fock Theory", R. F. Goodrich (2591-3-T, July 1958), AF-19(604)-1949, AFCRC-TN-58-350, AD 160790. UNCLASSIFIED. 73 pgs. KXVII "Calculated Far Field Patterns from Slot Arrays on Conical Shapes", R.E. Doll, R. F. Goodrich, R.E. Kleinman, A. L. Maffett, C. E. Schensted and K.M. Siegel (2713-1-F, February 1958), AF-33(038)28634 and 33(600)-36192; HAC-POs L-265165-F47, 4-500469-FC-47-D and 4-526406-FC-89-3. UNCLASSIFIED. 115 pgs. XVIII "The Physics of Radio Communication via the Moon", M.L. Barasch, H. Brysk, B.A. Harrison, T.B.A. Senior, K.M. Siegel and H. Weil (2673-1-F, March 1958), AF-30(602)-1725. UNCLASSIFIED. 86 pgs. XXIJX "The Determination of Spin, Tumbling Rates and Sizes of Satellites and Missiles", M.L. Barasch, W.E. Burdick, J.W. Crispin, Jr., B.A. Harrison, R.E. Kleinman, R.J. Leite, D.M. Raybin, T.B.A. Senior, K.M. Siegel and H. Weil (2758-1-T, April 1959), AF-33(600)-36793. SECRET. 180 pgs. v

THE UNIVERSITY OF MICHIGAN 2871-4-T XXX "The Theory of Scalar Diffraction with Application to the Prolate Spheroid", R.K. Ritt (with Appendix by N.D. Kazarinoff), (2591-4-T, August 1958), AF-19(604)-1949, AFCRC-TN-58-531, AD 160791. UNCLASSIFIED. 66 pgs. XXXI "Diffraction by an Imperfectly Conducting Half-Plane at Oblique Incidence", T.B.A. Senior (2778-2-T, February 1959), AF-30(602)-1853. UNCLASSIFIED. 35 pgs. XXXII "On the Theory of the Diffraction of a Plane Wave by a Large Perfectly Conducting Circular Cylinder", P.C. Clemmow (2778-3-T, February 1959), AF-30(602)-1853. UNCLASSIFIED. 29 pgs. XXXIII "Exact Near-Field and Far-Field Solution for the Back Scattering of a Pulse from a Perfectly Conducting Sphere", V. H. Weston (2778-4-T, April 1959), AF-30(602)-1853. UNCLASSIFIED. 61 pgs. XXXIV "An Infinite Legendre Integral Transform and Its Inverse", P.C. Clemmow (2778-5-T, March 1959), AF-30(602)-1853. UNCLASSIFIED. 35 pgs. XXXV "On the Scalar Theory of the Diffraction of a Plane Wave by a Large Sphere", P.C. Clemmow (2778-6-T, April 1959), AF-30(602)-1853. UNCLASSIFIED. 39 pgs. XXXVI "Diffraction of a Plane Wave by an Almost Circular Cylinder", P.C. Clemmow and V.H. Weston (2871-3-T, September 1959), AF19(604)-4993. AFCRC-TN-59-955. UNCLASSIFIED. 47 pgs. XXXVII "Enhancement of Radar Cross Sections of Warheads and Satellites by the Plasma Sheath", C. L. Dolph and H. Weil (2778-2-F, December 1959), AF-30(602)-1853. SECRET. XXXVIII "Non-Linear Modeling of Maxwell's Equations", J. E. Belyea, R. D. Low and K.M. Siegel (2871-4-T, December 1959), AF-19(604)-4993. UNCLASSIFIED. vi

THE UNIVERSITY OF MICHIGAN 2871-4-T Preface This is the thirty-eighth in a series of reports growing out of the study of radar cross sections at The Radiation Laboratory of The University of Michigan. Titles of the reports already published or presently in process of publication are listed on the preceding pages. When the study was first begun, the primary aim was to show that radar cross sections can be determined theoretically, the results being in good agreement with experiment. It is believed that by and large this aim has been achieved. In continuing this study, the objective is to determine means for computing the radar cross section of objects in a variety of different environments. This has led to an extension of the investigation to include not only the standard boundary-value problems, but also such topics as the emission and propagation of electromagnetic and acoustic waves, and phenomena connected with ionized media. Associated with the theoretical work is an experimental program which embraces (a) measurement of antennas and radar scatterers in order to verify data determined theoretically; (b) investigation of antenna behavior and cross section problems not amenable to theoretical solution; (c) problems associated with the design and development of microwave absorbers; and (d) low and high density ionization phenomena. K. M. Siegel vii

THE UNIVERSITY OF MICHIGAN 2871-4-T Table of Contents Page Abstract 1 Introduction 3 Non-linear Modeling of the Non-linear Spring Equation 11 General Non-linear Oscillator 16 Euler' s Equation 21 Boltzmann Equation 23 Non-linear Modeling in Acoustics 25 Maxwell's Equations 31 Conclusions 35 Appendix 37 References 39 ix

THE UNIVERSITY OF MICHIGAN 2871-4-T Non-Linear Modeling of Maxwell's Equations Abstract A basic stumbling block to carrying out many experiments in the laboratory is the reliance which is placed on the theory of linear modeling. In aerodynamics and electromagnetics one is seldom able to construct a precise linear model; one must appeal to physical reasoning to show that the linear model one is able to construct gives results in close agreement with those which would have been obtained with a precise one. For example, as long as good conductors are used, conductivity is not modeled. Until recent years the limitation on indoor radar cross section ranges involved the size of the model that could be handled, and the magnitude of the frequency of the coherent source. Even now no precise linear modeling experiments exist for bomber cross sections at very high microwave frequencies. With the growth of interest in I. C. B. M. programs it has become desirable (as measured in dollars) to use modeling to determine the radar 1

THE UNIVERSITY OF MICHIGAN 2871-4-T cross sections of ablating warheads, and to determine attenuation effects and scattering effects in or by plasmas. Linear modeling is of very limited utility for such experiments. One must model a lossy dielectric at high frequencies or not do the experiment in the laboratory. One must model some plasma parameters, or risk missing crucial first order effects. Previous analysis predicted the feasibility of non-linear modeling. The first such modeling results are given in Reference 1. In this paper several non-linear models of equations of mathematical physics are presented. Maxwell's equations are non-linearly modeled. Experiments making use of non-linear models are discussed. 2

THE UNIVERSITY OF MICHIGAN 2871-4-T Introduction Just exactly what is meant by modeling is not at all simple to define. By modeling we mean more than what the theorist implies when he reasons by analogy with some known and believed "similar" physical problem. If we are faced with finding a solution to a physical problem which (1) cannot be solved analytically by known techniques, (2) could be solved analytically but at too great a cost, or (3) would involve measurements, in an experimental approach, which are either too costly or too difficult to do on a full scale basis (at least initially), we must turn to a model-experiment approach. Thus we see the economic factor usually comes into the question of modeling. An aerodynamicist would not consider building a new aircraft until he had made wind tunnel tests,or for that matter would he obtain financing without such tests. We should also keep in mind that in this country it might be, and probably would be, easier to "'sell" results based on "modeling" experiments than it would be to "sell" results based on reasonable approximations to exact analytic results. In the latter case most people would require an experimental "check". When we discuss modeling we state one measurement system on a small model will allow us to obtain results which we can relate to the difficult or expensive full scale experimental results desired. The way this relationship is usually made is via linear transformations. To 3

THE UNIVERSITY OF MICHIGAN 2871-4-T clarify the meaning of linear modeling we will now go through a standard textbook example (Ref. 2). A fundamental example in fluid dynamics would be based on the Navier-Stokes equation. Consider the equation of motion aq 2 -q+ (q-7) q -v(')+ 2V2 q (1) at P where q is the velocity, p the pressure, p the density, t the time and V the kinematic viscosity. Now consider a scaled motion which differs only in the scales of length and time. Let us denote the quantities describing the second motion with dashes. at + (q'.V') q' =-V' )+ w' V'2 q' (2) ati' - p The motion is said to be "similar" if we can proceed from equation (1) to equation (2) by multiplying every term in equation (1) by the same constant factor a. By hypothesis we consider only scaling length and time r =Xr', t =Kt'; (3) for similarity we must have -q q q, 2 2 p, P t' t r' r p'r' p r v' q' avq r,2 r2

THE UNIVERSITY OF MICHIGAN 2871-4-T We note by division q' r'/t' p' p q2 q' r' qr q r/t P p q2' where R is the Reynolds number. This must be the same for the two motions. Since the equation of continuity as well as the first two of the above conditions is satisfied by (3), it follows that the equality of R is both necessary and sufficient for the similarity of motions. In experiments made with models in the wind tunnel the quantities q' and r' for the tunnel and model are less than for the full scale case, while W is the same for both. This has led to the use of compressed air wind tunnels where v = -- [, =coefficient of viscosity ] is decreased by p taking an increased value of p. Now that we understand what similarity implies, and the methods which are used in linear modeling, let us go on to non-linear modeling. Let us first consider an elementary example, the simple harmonic oscillator, to clarify the method of approach. Consider the initial dynamic system as m' = - k2 x (4) The modeling transformations are m = am (5) k1 = ~'k t1 = At 5

THE UNIVERSITY OF MICHIGAN 2871-4-T An unnecessary but nevertheless perhaps desirable requirement would be that the experiment in the laboratory be governed by the same equation of the motion. Thus we require d x1 2 m =- kl x1. (6) dt1 Now it is clear that if (x) were equal to a constant times x we would have linear modeling. Since there is nothing which a priori precludes non-linear modeling we shall attempt it here, but first we will present some ground rules. We pretend we cannot solve (4), and we perform transformations (5) and attempt to determine the values of ai such that we arrive at (6). We again pretend we cannot solve (6), and that the purpose of all this is that the experiment given by (6) is more economical than the experiment given by (4). Since we know the solutions to (4) and (6) we are in a position to check the concepts of non-linear modeling but we must not use this knowledge of the solutions to (4) and (6) in deriving the a.. Substituting relationships (5) into (6) we obtain 1 am _k2 (2k2 If we let 22 0t2 = 2, then m x =- 2k2. 6

THE UNIVERSITY OF MICHIGAN 2871-4-T By use of the chain rule, and of (4), this becomes m d2 + (-k) 2 (7) dx2 m dx Differentiating (7) with respect to time we obtain d3 + 2k d2 -k2 2. m - 3 r ___ (kx )+() _..k +( X)__ [ dx3 dx2 m m dx m dx2 (8) -2 2dp. dx Dividing through by -k2 x and collecting terms we obtain M d% 22 +x2 m2 ___2+3 x + (9) k2 d dx dx J dx.2 We can eliminate x between (7) and (9) to obtain 3 d x+(1X 2) d3 dx3 dx =0. (lOa) + =0. (1 Oa) d2o 2 k _ do x dx2 dx Now (10 Oa) may be written

THE UNIVERSITY OF MICHIGAN 2871-4-T d3 d dx2 d. =_.. (lOb) dx dA x_2p dx so that if we let d x2 v; 2~2~ _ dp x dx (10b) may be written 1 dv = -2 vx v dx A first integral of this is 1 =x2 12 v Substituting for v: X2 -c2) d p x 2 0. dx dx This equation possesses a general solution

THE UNIVERSITY OF MICHIGAN 2871-4-T -1 X-x =C2 cos(Xcos- ) + c sin( sin ). (11) We note if A= 1 we have linear modeling, since (C2 + c3) f =___ x. (Cl) The condition X = 1 forces 2/ 2 =aa. This is the same type similarity condition as we discussed for the Navier-Stokes equation. Note that (11) defines an infinite set of non-linear modeling functions. In the event that Xis an integer and c3 = 0 we see that 0 is a Tschebyscheff polynomial of the first kind of order X. The simplest example of this, that for A= 2, is 1 =0(x) =c 2 o (12a) Let us use this to test the validity of our model. A solution of (4) is x =A sin kt + B cos -k From(12a) we have[ with c2 = 1 r2(Asin k +Bcos -- ) x1 =L. (12b) 9

THE UNIVERSITY OF MICHIGAN 2871-4-T Now we use (5) to see if this is really a solution of (6). Since we let X= 2 to obtain (12a), we must use this restriction on the modeling constants: 2 2 _2-2- = 4. (13) Equation (6) has the form d2l dx =- k12 x, (6) dtl2 which becomes upon the use of (12b) 2 kt kt 2 d (Asin kt +B cos 2 am Ji_ __ i c12 /2 dt2 kt kt 2 2 2 2(A sin + Bcos k 1 414) C12 L2 and by use of (13) kt kt 2 2(A sin -z+ B cO- ) m d2 kt kt 2 2 sm+ Jm — (Asin-+Bcos =-2k cl2 dt2 c2 (15) Performing the indicated differentiation, and dividing through by 2 k we obtain 10

THE UNIVE RSITY O r MICHIGAN 2871-4-T 2 2 A + = 1. (16) c12 c12 2Ac 2 c 2 Thus if c12 =A + B we see that (12b) is a solution to (6). We have presented what we have called the simplest type of non-linear modeling. It will be recalled that the equation of motion was the same in both systems in this example. In the next section we shall consider a case of non-linear modeling in which the two systems are governed by different equations of motion. Non-linear Modeling of the Non-linear Spring Equation. The non-linear spring equation is'x + h(x2 -1) x =0. (17) In an attempt to model equation (17) non-linearly with respect to displacement, we use the modeling transformations x1 = (x) p bh t1 = ct. Let us assume that the model system is governed by dxl2 -d - p X1 = o0. (18) dt12 1 11

THE UNIVERSITY OF MICHIGAN 2871-4-T Substitution of the transformations into (18) yields -1 L-(- -k) - bh =1 d.2 + p - bh = c2 dt dx c2 d dx (19) Let ~ =bc2; then (19) becomes d_ _ 2+ d_ _A _ - ho =0. (20) dx2 dx Now by (17) we obtain d2 d x 2 d ) x2 _ d h(x -1) x -h =0. (21) dx2 dx We differentiate (21) with respect to time obtaining d3~ 2 3 d2 2 x3 + d 2x x[- h(x - l)x dx3 dx2 dx2 _ d h(x2-)1)_ d- h(2x )x - h d =0. dx dx dx (22) Dividing through by x, we obtain d3P.2 d2p 2x d2d h(x2 + (-1)+ 2d h2x2+hd dx3 dx2 dx2 dx dx dx (23) 12

THE UNIVERSITY OF MICHIGAN 2871-4-T Eliminating x2 between (21) and (23), we obtain 3 d2 -2x + h(x + -1)+h2x2 +h] _dx dx d22p dofd h(x2 -1) x + rhA dx dx (24) Since dd [- h(x2 -l)x+ hP d2 h(2-x)x+ h+ h(X-1)+2hx, dx Ldx iJdx2 dx L (25) (24) becomes, d~ 2~ d din I d lFg h(x21-)x+1h0l dx2 d x dx (:l d~in ln-x2 h -1)X+' dx (26) Using the temporary change of variables d2 v= (27) dJ h(x2-l1)x+' hp dx equation (26) takes on the form dd [lnv =-v.21. (28a)

THE UNIVERSITY OF MICHIGAN 2871-4-T We can rewrite (28a) in the form dv =2 [hx3 - hx dx (28b) v2 which upon integration becomes 1 2 hx4 2 1 = 2h x 2 + c. (28c) v 4 Replacing v by the expression (27), (28c) becomes 4 2 hx hX d + + (hx3 -hx) d +'hp =0. (29) 2 d2 dx Dividing through by h and letting c1 = we obtain (x4 2 l) d2 + (x - x) + = (30) ~~~ - x + C ~1) _ x) + ep =(30) 2 dx dx an ordinary linear equation with non-constant coefficients. In order to solve (30), we can take advantage of the fact that it may easily be thrown into self-adjoint form: dividing it by 2 -x +cl we obtain X4 x2+ e d _ - d_ + =0, 2 2 d2 xx 2 d 2 x +X -C1 - -x +c 42 12 (31) 14

THE UNIVERSITY OF MICHIGAN 2871-4-T which can be put in the form x2-x +eI —- d-x 2+ d + C =o. (32) 2 1 dx 2 1 dx If we define: x 2 d d -— x + c __ = 2 1dx dy then (32) becomes d20 + B =0 pc e (33a) 2 2 dy The above definition of y implies that dx dy= 4 X -x +c 2 1 an expression which may be integrated by quadrature to yield x x/f/t dt - IS ds 0 - 24 + c' ~ \(1- 02) 15 ~ 1s)1~

THE UNIVERSITY OF MICHIGAN 2871-4-T t where c' = 2 c1 ab, b + a = 2, and s. Since the integral appearing \Fa above is in the standard form of the elliptic integral of the first kind with modulus JI, expressed F ( ), the'modeling function we seek b 2 W c LC exp F (33b),( 2 b,I 33a If c' =1 we have the following simple special case. Y25i 2 =IZ __ __ __Q-<d t + 1+t 1 1+ 1-t 0 i 1+X W c2=c ( 1 -x Now we are ready to discuss more general systems. General Non-linear Oscillator Consider as a prototype a system governed by the equation 2 dx + f(x). = 0, (34) which describes the free, undamped, non-linear ~vibrations of quite general systems having one degree of freedom. If the model is required

THE UNIVERSITY OF MICHIGAN 2871-4-T to satisfy the equation 2 dy +f (y)=0, (35) ds in which f* differs from f only by a linear change in parameters, then it is desired to obtain the functional relation y = F(x) when the independent variables are related by s = at. This last relation permits one to write (35) in the form d2y 2 +a f* (y) =0, (36) dt and by use of the chain rule this can be written as 2 dY dx 2 dy d +y 2 d dx __) + ___ x +a f*(y) = 0. (37) d dx t dt2 2 d x Substitution for 2 from (34) and use of the identity dt d2x 1d dx 2 dt 2 2 dx dt permits one to write (37) in the form d2 f(x) dy a2 f* (y) (38) dY + dx 2 2Sf(x) dx 2 f(x) dx 17

THE UNIVERSITY OF MICHIGAN 2871-4-T Equation (38) is the differential equation whose solution gives the desired functional relationship between x and y. It does not seem possible to obtain its solution for arbitrary f; however, the equation can be written in a form more suitable for solution once f is prescribed. Setting g(x) = 5 f(x) dx (38) can be put in the form d 2 g(x) (g(x) ) _ a f (y) = 0, (39) dx dx 2 and by the change of variable dx dmg(x) (39) is reduced to d2y 2 a2 f (y) 0. (40) dm2 Equation (40) can be written in the form d (dy)2 2 d ( dy) - a f*(y) =0, dy dm which is solvable for m. Hence m i S dy (41) a g*(y) 18

THE UNIVERSITY OF MICHIGAN 2871-4-T where g* (y)= f *(y) dy. The definition of m together with (41) imply that Idx g S(dy) (42) g(x) a g* (y) which is the desired functional relationship between x and y. The non-linear spring equation (17) is of course a special case of the results just presented. For if in (34) we let 3 3 f(x) =,u(x3 -x), then f* (y) =,' (y -y), whereu' =b,u and (42) becomes 9 ~dx =1 *4dy.(43) 5 (x3 x) dx (y y) dy If the initial conditions are dx 0 when x =x dt 0 dy -O when y = y, dt 0 and if y = y1 when x = O, then (43) becomes 19

THE UNIVERSITY OF MICHIGAN 2871-4-T x Y S du 1 du x Y 0 v3 _ v) dv Y (v3 -v) dv Xo Yo (44) By performing two of the integrals in (44), and defining A, B=1+ 1+x4- 2x 2 C, D =1 \ + y4 - 2yo we obtain the expression F QC4 y) =F (\, Y)+ a b F Q) where F (, ) is the elliptic integral of the first kind as defined above in equation (33b). It should be noted that the above treatment of the nonlinear spring differs from that in the preceding section in that both model and prototype satisfy equations of motion of the same form. One reason for its inclusion was to show how easily boundary conditions may be handled. 20

THE UNIVERSITY OF MICHIGAN 2871-4-T Euler's Equation Consider a hydrodynamic system governed by Euler's equation d; 1 - Vp + G. (45) dt In the event that this system is either an adiabatic or an isothermal gas, p may be expressed as a function of p by an equation of state, p = f(p), and (45) may be written dv = 1 Vp+G(45) dt f(p) If for this prototype we specify a model system obeying d1' q+ G', (46) ds f *(q) it is desired to obtain the functional relation p = F(q) which arises when the other variables of the two systems are related by w= a s =bt, G = eG, c and when f * differs from f only by some linear change of parameters. These relations permit us to rewrite (46) as a dv 1 dq 7p+eG. (47) b dt c f*(q) dp 21

THE UNIVERSITY OF MICHIGAN 2871-4-T Solution of (47) for G yields G a dv 1 dq > + 1= sL p, be;t ce f* (q) dp and substitution of this result in (45) gives us di; a 1 dq 1 (1 a )+ Vp ( 1 dq 0. dt be f(p) dp ce f (q) If the stipulation is made that a =be, the above reduces to dq ce f* (q) dq - O, (48) dp f(p) since we may assume that Vp is not identically zero. The variables in this equation are separable, so that if f(p) is specified then p = F(q) may be found. For example, in the event that the system is an ideal adiabatic gas, then p =Kp and (48) becomes K' q dq ac 0. dp b IKpl/) This is immediately integrable to yield 22

THE UNIVERSITY OF MICHIGAN 2871-4-T P'('-1) bK q + C. _ ('-1) acK' Note that the modeling is linear if both model and prototype have the same specific heat ratios (' = ='). Boltzmann Equation Consider a statistical distribution of particles whose distribution function, F, satisfies the Boltzmann equation + V F + A. V F + I(F) = 0, (49) at v where I(F) is the collision integral of the system. Now the first three terms in (49) are simply the expanded form of the material DF derivative. Hence (49) can be written in the form Dt DF + I(F) 0. (50) Dt This form lends itself to modeling more rapidly than does (49). Hence if the model is required to satisfy an equation of the form - + I* (G)=O, (51) Ds 23

THE UNIVERSITY OF MICHIGAN 2871-4-T where I* differs from I only by a linear change in parameters, and if S = at, then a simple manipulation permits (50) to be written in the form a I* (G) F _ I(F) =0. (52) dG The variables are separable in (52); hence it may presumably be solved once I is prescribed. For example, in the event that both systems are neutral or weakly ionized gases, then I is given approximately by F-F I(F) where T is the characteristic time of the system, and F is its Maxwellian 0 distribution. Hence G-G I*(G) = 0 and the solution of (52.) is'l/aT F =F + C(G-G ), (53) O o where C is a constant of integration. Inspection of (53) shows that the modeling is non-linear except when'- = a Z. 24

THE UNIVERSITY OF MICHIGAN 2871-4-T Non-linear Modeling in Acoustics Let us next consider the scalar wave equation V2 p +k2 P =0, k =2wr/X (54) where X = wavelength. This equation we plan to non-linearly model by means of an equation of the same form v 2 + g= 0 (55) wherein we take xi' = a xi and Ag = b X. If we find we cannot do an economical experiment when b = a, then linear modeling is out of the question, since two distances X and the xi are modeled differently. In all cases known to the authors the experimentalists have used the laws of linear modeling and have forced the model size to be in the same ratio to the true size as the modeled wavelength is to the true wavelength. However, since in our case we cannot do so, let us find the non-linear modeling function P =p (e/) (56) which allows a andX to be arbitrary. From (56) we obtain vat =t'v yf, V2 =" (V7 )2 + p 2(57) = B, ( 7 t)2 5 f2 t, 25

THE UNIVERSITY OF MICHIGAN 2871-4-T where (V f )2 means 7 1 ~ 7 and use has been made of the fact that V2' =- Vj 2. Substitution of (57) into (54) gives (Vq)2 ",, -1_2 V'p + k2 =0. (58) To solve (58) one may proceed as follows: First, take the gradient of (58) to obtain "'VI(VI/2 V/ /+20,1(V V # ) WV -L2 V,,7+(k2-_12),'V:= 0. (59) Next, divide (59) by I| I, on the assumption that V O' 0, to get p' (V ) + 2 t+" (VV) -' 2 n, n+(k2 ~2)x, n =0, where n = is the unit normal to any solution surface I- = constant. Introducing the unit dyadic II, one can write this last equation in the form L0 )II+2P H+ V2 H+(k21. (6 0) Now since the surface ~ = const,, to which n is the normal, is quite arbitrary it may be concluded from (60) that b1"(V 4, )2 11+ 2p" V Vq-e2 42 c,, II+(k2_~2), II =0. (61) Next, take the divergence of (61) to get 0iv(Vt) 2Vp+ 40"?' (7V+q) V* +(k2 - 4L2)'V-2 2'V = 0. (62) In obtaining (62) the following results have been used:

THE UNIVERSITY OF MICHIGAN 2871-4-T (i) V ~ i=o, (ii) A * = A, for all A, (iii) V ( V V) =V(V2 ) = _e2 T (iv);V V? is self-conjugate. Of these results, (, (ii), and (iv) follow directly from the properties and definitions of the dyadics involved. The first part of (iii) follows directly upon expansion of 7 x ( V x 7 3) = 0, and the second part is a consequence of (55). Next, from (58) one finds that (V'),2 =e2(' -'X ),2 (63) where 2X = k/ 2, and from (59) and (63) it follows that 2 ( V...) / t= 2 t - 2 (64) Substitution of (63) and (64) into (62) gives (after division by v 2) 2+ ( { 2 4),, _t v,2= O 27

THE UNIVERSITY OF MICHIGAN 2871-4-T Now it has been assumed that 7V7/P 0, hence this last result implies, after some rearrangement, that C(VP'- X2 )+2(1-_~2) PIP"' -2 — tt2 (Vt, (2)+( -4)9'+3 t0. (65) A first integral of (65) is 2" (Y0, -_ 2 ) =(1 — 2). + 3?+ 2 C1 (66) ~,,2 0", where 2 C1 is a constant of integration. Multiplication of (66) by t" and a rearrangement of the result gives b,, Wb"+(1 _ %2)j, P" +(2 C + 2Y) P"l u'-rX22 1 p, 2 and if the substitution ~,, u = (67) is made, the previous equation becomes d (u)=(2C1 +2P) u. dY Integration of this equation gives 1 -C2 2 C 2 2C1, (68) See Appendix 28

THE UNIVERSITY OF MICHIGAN 2871-4-T where the constant of integration has been taken as -C2. Equations (67) and (68) now give (C22 _'2 -2C1T) _P +12 =0. (69) Setting S ='+ C1 reduces (69) to (m - d2_ + (C2 S = 0, (70) dSm dS 2 2 2 where m2 =C12 + C2 Equation (70) is similar in form to the standard hypergeometric equation and under the change of independent variable S = m (1 - 2x), it becomes x(l-x) d2 +( C1 - x) d- +72 =0. (71) dx2 2 2m dx The hypergeometric equation referred to is of the form 2 d( x(l -x) d + IY-(a+3+1) x] Y- -I3y=0, (72) dx2 dx and comparison of (72) and (71) shows that a =+, 3 =X, and ~ = 1/ 2 - C1/ 2m. The analysis from this point depends quite heavily on the value of' which in turn depends on the constants C1 and C2. It is known that the hypergeometric equation has two linearly independent 29

THE UNIVERSITY OF MICHIGAN 2871-4-T solutions regardless of the value of )/, but when /= O, 0 1, + 2, one or both of these solutions possess singularities of various types. Since the experimental parameters which make up 2 are under our control, we may avoid these troublesome integral values, however. The general solution of (72), valid for Ixl <1, is y =A 2F1(, /; x; x) + Bx -F (a-+ 1, 2-3(+ 1; 2 -4; x), m-S hence the general solution of (71), valid for < 1, is 2m 1 1 C1 m-S 3" F(22 2m 2m 1 C1 i-S 2 m 1 C1 3 C1 m-S +C (_ 2m F (+1+, ++ 1+; ) 24 m 2 1 2 2m 2 2m 2 2m 2m mI-C1 - y Thus, if < 1 then the desired functional relationship 2m between A and y becomes 6 =C2 F1 (XY, _; p C1 C m -C1 - 2m 1 C1 _ 2m 1 C1 1 C1 3 C1 m-C1-W 2mJl22m 22m2 2m 2m (73) 30

THE UNIVERSITY OF MICHIGAN 2871-4-T m - C1 -, The restriction < 1 may seem rather severe, 2m but since C1 is subject to our choice, we can produce a value such that the inequality is obeyed if we so desire. As a further remark, attention is called to the fact that if C1 = 0, equation (69) becomes (C22-_ Y2) o,, _ Vo + X 2 = 0, the general solution of which is = C3 cos (X Cos 1 L )+ C4 sin(7 sin (74) 2 C2 This is the same result obtained by us in the Introduction, equation (11) which is the one-dimensional analogue of equation (54). This suggests that for the usual three-dimensional cases C1 / 0. Maxwell' s Equations Following Stratton (Ref. 3), we will first go through the linear modeling procedure for electromagnetic phenomena. Maxwell's equations, which describe these phenomena, may be written Vx E + E) =0 (75) xH-E- - C a E =0. (76) @t 31

THE UNIVERSITY OF MICHIGAN 2871-4-T Suppose we construct a model system governed by D HI V' x E+u' =0 (77) at' and V' xlHI - -t -a' E' =0. (78);qt' whose variables and parameters are related to those of the system of (75) - (76) by E =a E' H =f H';> = ~/~' -(79) = al' x. =a x. 1 2 t = a t' 3 Substitution of (79) into (75) gives V' xE' +' H' 0; ( 80) a a Also, substitution into (76) gives 32

THE UNIVERSITY OF MICHIGAN 2871-4-T V3 xH' _- _ _ _ - a ac' E' =0 (81) a2 a3 Multiplying by a2//3: 2 fI aa1 a2 7'x H' HI - I' -'?E' =0 (82) a38 at' / II III I, II, III are the invariants of the linear modeling. We can eliminate the common ratio a// among them to obtain - 2 = constant(83) (83) a22 s a = constant 1 a3 If the two quantities of (83) cannot be made invariant for a particular modeling application, exact linear modeling is impossible. One such case arises when one attempts to obtain the X-band radar cross section of a B-70 aircraft, where object dimensions and wavelength cannot be scaled in the same manner for a reasonable sized experiment. The method of "approximately linear" models, which has been used in the past, is undesirable. 33

THE UNIVERSITY OF MICHIGAN 2871-4-T However, non-linear modeling of Maxwell's equations is now possible because of our results on the scalar wave equation. It has been found by Schelkunoff (Ref. 4) among others, that the most general electromagnetic field in a source-free region can be described in terms of two scalar wave functions. By this method, if 01 and 02 are two scalar functions satisfying V2 i+kk2 =O i=1, 2 (84) then the electric and magnetic vectors are given by 1 E = 1E = (_2 +k2_ ) (85) y,E ( +k Ex E' z 2 1 +iw6 axDz Y +i(w eyaz +iwz D (z ~2~2,H H 2)2 H 2 02 H ay ( - a2 +k k2 ) (86) il ai xz Y i0 ay yaz ic0p a z Since both 01 and 0Z are scalar wave functions they can be modeled as was shown in the previous section. This leads us to believe that we may non-linearly model electromagnetic systems in source-free regions, since 01 and 02 uniquely represent such systems. The actual techniques which would have to be used to perform this modeling may be summarized as follows: The electric and magnetic fields for the model system are measured. The associated potentials 01 and 02 are then found using equations of the form (85), (86). The potentials of the full-scale system may be determined from these by the scalar wave function modeling transformations'1= Y/l(01) and UY2= ~2(3) given on page 30. The electric and magnetic 34

THE UNIVERSITY OF MICHIGAN 2871-4-T fields of the full-scale system are found from these by using another set of equations of the form (85), (86). It should be noted that it is possible to represent an electromagnetic field in terms of two scalar wave functions by other methods which are presumably different from the one described above. Readers are referred to References 5 and 6. The relation of these different scalar wave functions to each other, and to the Debye potentials will be discussed in a forthcoming report. This paper represents the sum total of our results in non-linear modeling to date. During the next few months we plan to verify the results of this section as was done with the simple harmonic oscillator result in the introduction. If this is successful we plan to perform some experiments using non-linear modeling. This experimental program will be under the direction of Ralph E. Hiatt. Conclusions In this paper we have non-linearly modeled the simple harmonic oscillator, the non-linear spring equation, Euler's equation, the Boltzmann equation, a general non-linear oscillator equation, and the scalar wave equation. We have exhibited a method which shows how to nonlinearly model Maxwell's equations by utilizing the non-linear modeling 35

THE UNIVERSITY OF MICHIGAN 2871-4-T of the scalar wave equation. Despite the a priori expectation that non-linear modeling was at best mathematically very complex, if not impossible, it now appears that there is a possibility of non-linearly modeling more difficult equations, such as the coupled set consisting of Maxwell's equations and the Boltzmann transport equation. Of course, if these non-linear models can be made useful and if results remain general, or at least if the physically interesting cases can be solved, then it appears that one should be able to eliminate a great many expensive facilities which are presently required. It also appears that certain experiments previously considered beyond our present capabilities will now become feasible. The authors would like to gratefully acknowledge the many helpful discussions which they had with Professors C. L. Dolph and R. K. Ritt. 36

THE UNIVERSITY OF MICHIGAN 2871-4-T Appendix The analysis leading to the first integral, (66) of equation (65), although lengthy, adds little to the main body of the section and is thus included here for the interested reader. Recall that equation (65) is iv (p, 12)+2( -12) -2 (' )+((2 -4)0 + ~,,2 If this equation is multiplied by and the result slightly _;t2p rearranged, it becomes P'8t Pt_2 3....3, +(1_7(2)P. IV +(1-) -- - 0- -- 2f 1 0 2132=o. Noting that the coefficient of ~" p"' in the fourth term of this equation is the derivative, with respect to 3//, of A ( i/%t -_( 2 ), one obtains, after solving the equation for this derivative, the result ad, 32 2pp' 0,,i 2 _p,2 A-~ (8I pi 0 )-+' - (1 _-(2)'-e do p? 1 6,, 6"' p"'(Pf'- 2 p) Now the first two terms on the right side of this last equation are the derivatives of An(,2) and -a2 " "' respectively; hence combining them with the logarithm term on the left, one obtains

THE UNIVERSITY OF MICHIGAN 2871-4-T d 2~~2P ~,,, t(~,_20) 3 ~,,2,2 0 _p p,,2 1 1,,, ((A,_2)1,,, The last term on the right can be written as _(1 2) _ 2 d ( ), hence, introduction of this result into (A-1) followed by multiplication by 0"'(~p't -x2b) gives ~,,2 d ptp 2x p) d P 0 d 11,,2 d (p102 d j d d' (3 + 2 Cl) + (1 _2) ( ), d3,, where 2C1 is a constant of integration. Equation(66) readily follows. 38

THE UNIVERSITY OF MICHIGAN 2871-4-T References 1. Ritt, R. K., "The Modeling of Physical Systems", IRE Trans. on Antennas and Propagation, Vol. AP-4, No. 3, July 1958. 2. Milne-Thomson, L. M., Theoretical Hydrodynamics, 2nd Edition, The MacMillan Company, New York (1950). 3. Stratton, J. A., Electromagnetic Theory, First Edition, McGrawHill Book Company, Inc., New York and London, (1941). 4. Schelkunoff, S. A., Electromagnetic Waves, D. VanNostrand Company, Inc., New York, (1943). 5. Green, H. S. and Wolf, E., Proc. Phys. Soc., A, 66, (1953), 1129. 6. Nisbet, A., Proc. Roy. Soc., A, 231 (1955), 251. 39

UNIVERSITY OF MICHIGAN 3 901 5 02229 3024 3 0502229 3024