SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Lighthill method with Lagrange expansion for approximative solutions of partial differential equations
It is shown that a general class of nonlinear integral equations may be transformed into a Cauchy system. That this leads to an effective numerical scheme is demonstrated by solving the Ambarzumian integral equation. The new method does not involve successive approximations or series expansions.
A scheme is introduced which makes it feasible to make completely self-consistent Brueckner-Hartree-Fock (BHF) and renormalized BHF calculations for spherical, closed-shell and axially-symmetric deformed nuclei. The usual requirement or orbital self-consistency has been imposed, as well as self-consistency in the starting energies and occupation probabilities. Previously, only approximate forms were used for the Pauli operator. This approximation is removed and a method for making the necessary Pauli corrections to the reaction matrix during the approach to self-consistency is presented. A discussion of the symmetries which reduce the problem to one of manageable proportions is included.
A solution procedure was developed using linear small deflection theory for the flutter of simply supported laminated plates. For such plates, the bending and extensional governing equations are coupled and have cross-stiffness terms which do not appear in classical plate theory. An extended Galerkin method is used to obtain approximate solutions to the governing equations, and the aerodynamic pressure loading used in the analysis is that given by linear piston theory with flow at arbitrary cross-flow angle. A limited parametric study was conducted for typical laminated composite plates. The calculations show that both the bending-extensional coupling and the cross-stiffness terms have a large destabilizing effect on flutter. Since classical plate theory does not consider bending-extensional coupling and cross stiffness terms, it usually gives inaccurate and nonconservative flutter boundaries for laminated plates.
A model for the motion of flare-generated shocks through interplanetary streams is presented which illustrates the effects of a stream-shock interaction on the shock strength and geometry. It is a gasdynamic calculation based on Whitham's (1958) method and on an empirical approximation for the relevant characteristics of streams. The results show that the Mach number of a shock can decrease appreciably to near unity in the interaction region ahead of streams and that the interaction of a spherically symmetric shock with a spiral-shaped corotating stream can cause significant distortions of the initial shock-front geometry. The geometry of a shock discussed by Lepping and Chao (1972) is qualitatively explained by this model.
An anisotropic buckling and flutter analysis is developed with allowance for both bending-extensional coupling and bending-twisting coupling within the framework of linear small deflection theory for simply supported general laminated plates. The extended Galerkin method is used to obtain approximate solutions to the coupled governing equations. The effects of various anisotropic stiffness parameters on the static and dynamic stability of laminated plates are evaluated, with particular emphasis on assessing the range of applicability of classical orthotropic plate theory. It is shown that bending-extensional coupling and bending-twisting stiffness terms have a destabilizing effect on buckling and flutter, the effect being more pronounced for a small number of layers. For symmetric plates, the number of layers required for orthotropic plate theory to be applicable is generally less for the buckling problem than for flutter. For square plates, aligning the fibers with the direction of airflow over the plate surface results in the highest flutter dynamic pressure.
Finite-element techniques, mathematical programming methods, and innovative collection of approximation techniques are combined in program for finding efficient minimum-weight optimum designs for significant class of structural synthesis problems.
The problem of calculating the steady-state free-electron energy distribution in a hydrogen gas is considered in order to study departures of that distribution from a Maxwellian at sufficiently low degrees of ionization. A model kinetic equation is formulated and solved analytically for the one-particle electron distribution function in a steady-state partially ionized hydrogen gas, and it is shown that the formal solution can be accurately approximated by using the WKB method. The solutions obtained indicate that the high-energy tail of the distribution is susceptible to distortion by imbalanced inelastic collisions for ionization fractions not exceeding about 0.1 and that such departures from a Maxwellian can lead to significant changes in the collisional excitation and ionization rates of ground-state hydrogen atoms. Expressions for the electron-hydrogen collision rates are derived which explicitly display their dependence on the hydrogen departure coefficients. The results are applied in order to compare self-consistent predictions with those based on the a priori assumption of a Maxwellian distribution for models of the thermal ionization equilibrium of hydrogen in the optically thin limit, spectral-line formation by a gas consisting of two-level atoms, and radiative transfer in finite slabs by a gas of four-level hydrogen atoms.
A review of recent theoretical modeling of the effects of turbulence on acoustic wave propagation is presented. The Born approximation limit of the smoothing method of Keller and of the kinetic theory method of Howe are compared with the classical perturbation results of Chernov. Secular growth is observed in the higher order terms of the kinetic theory method of Howe when multiple scaling is ignored. The forward propagation - beam broadening model of Brown and Clifford is discussed as well as the single scatter model of DeLoach and of Rudd. Each model is examined with respect to its prediction of conservation of energy.
A method of finding a first approximation to a crustal magnetization distribution from inversion of satellite magnetic anomaly data is described. Magnetization is expressed as a Fourier Series in a segment of spherical shell. Input to this procedure is an equivalent source representation of the observed anomaly field. Instability of the inversion occurs when high frequency noise is present in the input data, or when the series is carried to an excessively high wave number. Preliminary results are given for the United States and adjacent areas.
The Galerkin finite-element method is used to obtain approximate solutions for the three-dimensional induction problem. A rectangular conductive prism is considered as an example, and solutions are obtained for linear and circularly polarized incident plane-wave fields. Magnetotelluric tensor impedances and magnetic transfer functions are computed. Polar diagrams of the tensor impedances and magnetic transfer functions along with their amplitude contour maps are presented. The dimensionality parameter, skew, is contoured at the surface of the earth. It is shown that the relative amplitudes and shapes of the additional and principal impedance polar diagrams can be used to determine the dimensionality of geoelectrical structures. Stations with skew values greater than 0.2 are significantly influenced by the three-dimensionality of the geoelectric structure. The amplitudes of the magnetic transfer function and the orientations of its polar diagrams exhibit large anomalies in the vicinity of the intersection of the lateral contacts.
A method of finding a first approximation to a crustal magnetization distribution from inversion of satellite magnetic anomaly data is described. Magnetization is expressed as a Fourier series in a segment of spherical shell. Input to this procedure is an equivalent source representation of the observed anomaly field. Instability of the inversion occurs when high frequency noise is present in the input data, or when the series is carried to an excessively high wave number. Preliminary results are given for the United States and adjacent areas.
An extrapolation method based on a screening approximation, applied to available initial values of polarizability for low stages of ionization, is used to obtain dipole and quadrupole polarizabilities for more highly ionized members of many isoelectronic sequences. It is suggested that the derived screening constants x sub L and limiting ratios F sub L may have significant physical meaning, especially the latter which may have an interpretation in terms of hydrogenic polarizabilities.
A study of the response of a clamped panel to a supersonic turbulent boundary layer, based on a recently developed theoretical model, is presented. It is demonstrated that the model incorporates the effect of coupling between the panel motion and the flow of the surrounding fluid. Further, a Ritz-Galerkin method is used to obtain approximate solutions for the statistics of the panel response to the turbulence. The numerical results are compared with previous experimental data and the theoretical model is assessed.
The objective of this paper is to develop a new and more accurate local linearization algorithm for numerically solving sets of linear time-varying differential equations. Of special interest is the application of this algorithm to the quaternion rate equations. The results are compared, both analytically and experimentally, with previous results using local linearization methods. The new algorithm requires approximately one-third more calculations per step than the previously developed local linearization algorithm; however, this disadvantage could be reduced by using parallel implementation. For some cases the new algorithm yields significant improvement in accuracy, even with an enlarged sampling interval. The reverse is true in other cases. The errors depend on the values of angular velocity, angular acceleration, and integration step size. One important result is that for the worst case the new algorithm can guarantee eigenvalues nearer the region of stability than can the previously developed algorithm.
Differential, integral, momentum transfer, and partial cross sections have been calculated for elastic scattering and rotational excitation of C2H2 by 10-eV electrons. The effective potential includes static, exchange, and polarization interactions calculated by the INDOX/1s method and the semiclassical exchange approximation with adiabatic polarization at large electron-molecule distances. The scattering is treated by well converged rotational close coupling using the centrifugal dominant scheme to select the channels included and including up to 32 coupled channels for a given total angular momentum. The calculated integral cross sections for pure elastic scattering and rotation excitation are 54.5 and 41.4 a(0)squared, respectively. These are much larger than the values (34.4 and 18.6 a(0)squared) previously (Onda and Truhlar, 1979) calculated for the isoelectronic molecule N2, at this energy. This illustrates how the greater spatial extent of C2H2 greatly increases the cross sections for pure elastic and rotationally inelastic scattering.
The linear instability of a non-zonal flow can be reduced to an eigenvalue-eigenfunction problem, governed by a nonseparable partial differential equation (Niehaus, 1980). Approximate solutions, found by the method of multiple scales, are derived here and compared with earlier results found using a spectral method. The amplitude maxima are correctly located. The zonal variations of local wavenumber and of amplitude are qualitatively correct, but not sufficiently extreme. Because the method is oversensitive to local conditions, and less sensitive to global constraints, this comparison provides theoretical limits to the possibility of parameterizing transient eddies in terms of the local time mean state of the atmosphere. The method can be extended easily to flows with more realistic vertical structure.
Improved calculation method promises to reduce computer time by factor of 10 for vibration analysis of complex structures. Method originally developed to predict mechanical response of spacecraft to maneuvers in various cases in which same booster is used with different payloads. Method is more exact than conventional normal-mode methods; it involves no initial approximations of assumptions. Calculation cost is reduced without loss of accuracy.