Search

A numerical method is given for the determination of the impedance of an infinitely long thin-wire helix. The propagation constant of the current for zero tangential electric field is found and used in a variational expression for impedance. Asymptotic values of resistance versus pitch are compared with resistances of infinitely long straight-wire antennas.

The streamline curvature technique for the numerical solution of the equations of motion for two-dimensional flow has been analysed. By considering a simple flow case an optimum damping factor for the quasi-orthogonal method was derived as a function of grid aspect ratio, Mach number, and differentiation formula, and the number of iterations required for a given accuracy was found. The stability and accuracy properties of numerical methods for finding the second derivative of a curve were also investigated, and it was found that, in this context, finite difference methods using polynomials were the best and spline methods the worst.

In order to solve the waveguide junction problem numerically, we express the fields in the guides by truncated modal expansions and construct an error function which is a measure of the mean-square error in the matching of the boundary conditions at the junction. The minimum of this error leads to a set of linear equations for the modal amplitudes. Offset rectangular waveguides with amplitude and current excitations are studied. A weighting factor which multiplies the error contribution due to the magnetic field is studied, and a criterion for its selection given.

1. In a previous communication I proposed a formula for the wave making resistance of ships, and showed that it expressed certain general qualities of experimental results; further, notwithstanding the limitations of theory and the difficulty of interpretation of experimental data, a good numerical agreement was found in several cases with the published results of tank experiments on models when suitable numerical values were given to the coefficients in the formula. This paper records the results of a more systematic study of the numerical values of some of the coefficients, the data being taken from certain recent series of experiments; for the present the discussion is limited to those types of model whose resistance-speed curves show clearly the humps and hollows which are attributed to interference of wave-systems originating at the bow and stern. It has been remarked that although the mode of disturbance is different, the action of the bows of a ship may be roughly compared to that of a travelling pressure-point, and further, that the stern may he regarded in the same way as a negative pressure-point. This point of view originated in the well-known paper of W. Froude on the effect of the length of parallel middle body, and the theory was developed in a later paper by R. E. Froude; from an inspection of experimental results it was seen that the variations in magnitude and position of the oscillations were in directions which agreed with the above interpretation. On account of the lack of an adequate formula, the available data have not yet been examined numerically in any detail; the present investigation aims at supplying this in some measure. Section 2 is theoretical, with some necessary repetition of previous work; Sections 3 and 4 contain a numerical analysis of some available experimental curves. In Section 5 an attempt is made to estimate the effective horse-power of the "Turbinia," in order to illustrate certain points; while in Section 6 the limitations of the interference theory, in the conventional use of the term, are discussed in connection with the residuary resistance curves of finer-ended models.

A numerical-perturbation method is proposed for the determination of the nonlinear forced response of structural elements. Purely analytical techniques are capable of determining the response of structural elements having simple geometries and simple variations in thickness and properties, but they are not applicable to elements with complicated structure and boundaries. Numerical techniques are effective in determining the linear response of complicated structures, but they are not optimal for determining the nonlinear response of even simple elements when modal interactions take place due to the complicated nature of the response. Therefore, the optimum is a combined numerical and perturbation technique. The present technique is applied to beams with varying cross sections. ~ 4Y large-amplitude deflection of a beam or a plate which is restrained at its ends or along its edges results in some midplane stretching/One must account for this stretching with nonlinear strain-displacement relationships. The nonlinear equations of motion describing this situation were the basis of a number of earlier investigations and are the basis for the present paper as well. The purpose of the present paper is to present a new scheme for determining the response to a harmonic excitation. Emphasis is placed on the case when the frequency of the excitation is near a natural frequency. A convenient way to attack this nonlinear problem involves representing the deflection curve or surface with an expansion in terms of the linear, free-oscillation modes. The deflection is then determined in two steps. First, the damping, the forcing, and the nonlinear terms are deleted and the linear modes (eigenfunctions) and natural frequencies (eigenvalues) are determined. Second, the time-dependent coefficients in the expansion are obtained from a set of coupled, nonlinear, ordinary, second-order differential equations, the linear modes being used to determine the coefficients in these equations. (The procedure is described in detail in Sec. II.) Generally, one cannot obtain the linear modes analytically for structural elements having complicated boundaries and composition, and one cannot easily determine the character of the timedependent coefficients through numerical integration of the set of nonlinear equations. (The results obtained in the present numerical example are typical of the complicated manner in which the steady-state amplitudes of the various modes making up the response can vary with the amplitude and the frequency of the excitation.) Consequently, an optimal procedure involves a numerical method to determine the linear, free-oscillation modes and an analytical method to determine the time-dependent coefficients. The present procedure combines either a finiteelement or a finite-difference method with the method of multiple scales (see, for example, Ref. 1). The following brief review mentions representative examples of the work that was and is

This paper describes a fast and numerically stable method for solving the discrete Dirichlet problem for Poisson's equation in case of a rectangle (and mainly, a square). By using a special calculus for difference operators, the system of linear equations is reduced to a block-triangular system such that the diagonal blocks are heavily diagonally dominant. For a standard version of the algorithm, the number of operations and the computing time are proportional toh ?2 (h=mesh width). The method is one oftotal reduction compared with the method ofblock-cyclic reduction (odd-even reduction) [2], which we describe as a method ofpartial reduction.--Due to the developed calculus, many generalizations are possible.--In a following part II of the paper, the algorithm and numerical results will be described in detail.

Summary This paper is the first of a series that will examine the effect of earth structure on earthquake displacement, strain and tilt fields at the Earth's surface. Its purpose is to develop the numerical techniques to be applied in the papers that foliow. A general computational procedure for the evaluation of the integral expressions for the surface displacements due to an arbitrary point dislocation source in a layered medium is described. It is shown to be rapid and inexpensive to use, and its accuracy appears to be entirely adequate for practical purposes.

In recent papers, the Prufer transformation has been used to calculate the eigenvalues of second order Sturm–Liouville boundary value problems. In the present paper, a Prufer transformation for fourth order vibration problems is devised using four-dimensional spherical coordinates. Some properties of the resulting angular system are investigated. These results are then used to develop a new computational procedure for finding the eigenvalues of a vibrating beam. This procedure is used to compute the eigenvalues of several singular and nonsingular problems.

The starting point of this paper is a theorem, based on the theory of functions of a complex variable, which essentially gives a representation of the remainder in a quadrature rule as a contour integral. From this theorem it is shown how most of the well-known interpolatory and non-interpolatory quadrature rules may be derived. Some new quadrature rules are also given. The latter part of the paper is devoted to the evaluation of asymptotic estimates of the remainder, by consideration of particular examples. It is shown how this asymptotic analysis frequently gives excellent estimates of the remainder.

A separately excited d-c. motor supplied from a d-c. series generator performs an unusual type of oscillation. For a reasonably clear understanding of the electromechanical phenomenon underlying this behavior, a nonlinear differential equation was developed and its complete solution was obtained by partially graphical and numerical methods. In particular Lienard's method was adopted to examine the possibility of any periodic solutions. In this method periodic solutions are easily recognized by the presence of the so called “limit cycles,” their number being equal to the number of such loops.Of these, one or more may give stable solutions and the rest unstable.In a set of investigations with which this paper is concerned two limit cycles in general were observed of which one and only one represents stable oscillations. The existence of these limit cycles itself was found to depend on various factors such as the polar movement of inertia, the shaft load, the brush shift, the circuit resistance and excitation.The paper confines itself to the discussion of the, effects of variation of the first three factors and in the main the stability of oscillations. Any additional facts observed are also recorded.

This paper describes a numerical method for solving first-order coupled matrix differential equations. Recursive equations are used to find ‘reflection’, ‘transmission’ and ‘source’ coefficients; these coefficients are then used to construct the vector solution to the differential equations. The method can be used to solve linear and non-linear differential equations with specified initial or two-point boundary values. Numerical results for several initial value problems are given in the paper.

Consider the Cauchy transform \[T\varphi (z) = \frac{1}{{\pi i}}\int_\Gamma {\frac{{\varphi (\zeta )d\zeta }}{{\zeta - z}},} \quad z \in \Gamma ,\] in which $\Gamma $ is a simple closed curve with a continuously differentiable parameterization. In order to have direct methods for the numerical solution of equations involving $T\varphi (z)$, it is desirable to have numerical integration methods for evaluating $T\varphi (z)$. In this paper, numerical methods are investigated which are based on replacing $\varphi (z)$ by a uniformly convergent sequence $\varphi _n (z)$; if these approximations are “sufficiently smooth”, then the speed of convergence of $T_n \varphi $ to $T\varphi $ is “essentially the same” as the speed of convergence of $T\varphi _n $ to $\varphi $. Specific methods are analyzed, with $\varphi _n (z)$ defined as a piecewise linear or quadratic interpolating function to $\varphi (z)$ at a given set of nodes on $\Gamma $. Numerical examples conclude the paper.

This paper is primarily concerned with the problem of finding approximants for the solution of initial boundary value problems of parabolic type and obtaining a priori error estimates.Perhaps the most widely used method for the treatment of such problems is that of finite differences. Some difficulties arising in this method are stability of the scheme and the treatment of the boundary. Another approach to this problem is the so-called semidiscrete variational method. As in the finite difference approach, a difficulty in this method is the treatment of boundary values.The methods considered in this paper will be of the least-squares type. Roughly speaking, the method consists in finding that element of a suitable finite-dimensional subspace which best approximates the data of the initial-boundary value problem in a least-squares sense. A principal virtue of this approach is that the elements of the finite-dimensional subspace are not required to satisfy any boundary or initial conditions. A particular cho...

This paper reviews some of the more useful, current and newly developing methods for the solution of electromagnetic fields. It begins with an introduction to numerical methods in general, including specific references to the mathematical tools required for field analysis, e.g., solution of systems of simultaneous linear equations by direct and iterative means, the matrix eigenvalue problem, finite difference differentiation and integration, error estimates, and common types of boundary conditions. This is followed by a description of finite difference solution of boundary and initial value problems. The paper reviews the mathematical principles behind variational methods, from the Hilbert space point of view, for both eigenvalue and deterministic problems. The significance of natural boundary conditions is pointed out. The Rayleigh-Ritz approach for determining the minimizing sequence is explained, followed by a brief description of the finite element method. The paper concludes with an introduction to the techniques and importance of hybrid computation.

The numerical solution of linear integral equations of the types studied by Volterra has formed the subject of a recent memoir by E. T. Whittaker. Numerical methods are needed also for the solution of the linear integral equation studied by Fredholm and Hilbert and the object of this paper is to describe a method which may sometimes be useful. The linear integral equation of the second kind, f(s) = Ф(s) - λ ∫ 1 0 k(s, t) Ф(t)dt , may be solved for the unknown function Ф(t) teither directly or indirectly. In the direct methods of solution the function Ф(t) is expressed by means of infinite series of terms involving repeated integrals. In the so-called method of Neumann only one infinite series is used, while in the more complete method of Fredholm the function Ф(t) is expressed as the ratio of two infinite series which converge for all values of the parameter λ. The method of Neumann has been employed on many occasions to obtain numerical results and is generally more convenient to use than Fredholm’s method on account of the great complexity of the expressions occurring in Fredholm’s series. There are occasions, however, when the Neumann series fails to converge, or converges only slowly, and then some other method such as Fredholm’s must be used. It seems desirable if possible to devise simple approximate methods which possess some of the advantages of Fredholm’s method, because in many cases the evaluation of repeated integrals becomes very tedious and the use of the direct methods of solution becomes impracticable except for a rough approxi­mation.