Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

An instability associated with a magnetosphere-disk interaction

The evolution of a thin accretion disk surrounding a rapidly rotating magnetosphere is considered. By taking account of the variations of the magnetospheric boundary in response to the conditions at the inner edge of the disk, we find from linear analysis and numerical computation that the accretion disk can become unstable. Mass can be accreted by the central object in a cyclic fashion, with the cycle involving the storage and release of mass in the inner parts of the disk. The physical origin of the instability is associated with the variations of the magnetospheric boundary about corotation. The recurrence time scale of the cycle can vary by several orders of magnitude depending on the details of the conditions at the magnotosphere. The possible applicability of this instability process to the 'rapid burster' MXB 1730-335 is briefly discussed.

Spruit, H. C.↗
NASA NTRS1993

Further development towards a new virtual-pulse time integral methodology for general nonlinear transient thermal analysis

New developments describing the theoretical basis towards an effective virtual-pulse time integral procedure are presented for general nonlinear transient heat transfer problems. To validate the proposed methodology of computation, numerical test examples are given and comparisons are drawn with the implicit second-order accurate Crank-Nicolson technique. The proposed method shows superior or improved accuracy and stability characteristics for the models tested.

Tamma, K. K.↗
NASA NTRS1993

Centrifuge in space fluid flow visualization experiment

A prototype flow visualization system is constructed to examine buoyancy driven flows during centrifugation in space. An axial density gradient is formed by imposing a thermal gradient between the two ends of the test cell. Numerical computations for this geometry showed that the Prandtl number plays a limited part in determining the flow.

Arnold, William A.↗
NASA NTRS1992

Prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations

Steady and unsteady asymmetric vortical flows around slender bodies at high angles of attack are solved using the unsteady, compressible, thin-layer Navier-Stokes equations. An implicit, upwind-biased, flux-difference splitting, finite-volume scheme is used for the numerical computations. For supersonic flows past point cones, the locally conical flow assumption have been used for efficient computational studies of this phenomenon. Asymmetric flows past a 5-deg semiapex-angle circular cone at different angles of attack, free-stream Mach numbers, and Reynolds numbers have been studied in responses to different sources of disturbances. The effects of grid fineness and computational domain size have also been investigated. Next, the responses of three-dimensional supersonic asymmetric flow around a 5-deg circular cone at different angles of attack and Reynolds numbers to short-duration sideslip disturbances are presented. The results show that flow asymmetry becomes stronger as the Reynolds number and angles of attack are increased. One of the cases of flow over a cone-cylinder configuration is validated fairly well by experimental data.

Liu, C. H.↗
NASA NTRS1993

Neptune's magnetic field - Calculation of field geometric invariants derived from the I8E1 GSFC/BRI model

L shell values along the Voyager 2 encounter trajectory and those associated with the N1 through N6 moons and N1R through N6R rings of Neptune are computed numerically on the basis of a simplified description of the Neptunian magnetic field derived from the Goddard Space Flight Center/Bartol Research Institute I8E1 model, which includes internal terms up to and including the octupole (but no external terms). Like Uranus, the large tilt between the dipole term and the rotation axis causes the moons and rings to sweep a very large range of L shells. Their orbital motion introduces additional periodicities, causing the maxima and minima in L space to vary systematically with time.

Acuna, M. H.↗
NASA NTRS1993

Experimental and numerical investigation of supersonic turbulent flow in an annular duct

Experimental and numerical results are presented for developing supersonic turbulent flow in an annular duct formed by a circular centerbody and outer shroud. The experimental results are based on data taken in a new flow facility that was designed to generate a shock-free, supersonic annular flow. Numerical computations were performed using the Baldwin-Lomax turbulence model for comparison with experimentally measured profiles. The results demonstrate that computed and measured profiles are in excellent agreement, so that studies can now be conducted of shock wave/boundary layer interaction phenomena within the duct, such as those induced by changes in downstream duct geometry or by the placement of struts between the duct walls.

Williams, K. E.↗
NASA NTRS1993

Coronal heating by the resonant absorption of Alfven waves - Importance of the global mode and scaling laws

Numerical simulations of the MHD equations for a fully compressible, low-beta, resistive plasma are used to study the resonance absorption process for the heating of coronal active region loops. Comparisons with more approximate analytic models show that the major predictions of the analytic theories are, to a large extent, confirmed by the numerical computations. The simulations demonstrate that the dissipation occurs primarily in a thin resonance layer. Some of the analytically predicted features verified by the simulations are (a) the position of the resonance layer within the initial inhomogeneity; (b) the importance of the global mode for a large range of loop densities; (c) the dependence of the resonance layer thickness and the steady-state heating rate on the dissipation coefficient; and (d) the time required for the resonance layer to form. In contrast with some previous analytic and simulation results, the time for the loop to reach a steady state is found to be the phase-mixing time rather than a dissipation time. This disagreement is shown to result from neglect of the existence of the global mode in some of the earlier analyses. The resonant absorption process is also shown to behave similar to a classical driven harmonic oscillator.

Steinolfson, Richard S.↗
NASA NTRS1924

Some tables of the factor of apparent additional mass

This note, prepared for publication by the National Advisory Committee for Aeronautics, is a collection of the tables of the factor of apparent mass that have been published up to now. The theory of the motion of solids in a perfect fluid is of the greatest value for the study of most aerodynamic problems, and the additional apparent mass of an immersed solid is the most important characteristic for such theoretical numerical computations. It will therefore be helpful to have the most important values of the apparent mass - for some elementary cases - collected in a convenient form.

Munk, Max M↗
NASA NTRS1938

A study of the torque equilibrium of an autogiro rotor

Two improvements have been made in the method developed in NACA Reports nos. 487 and 591 for the estimation of the inflow velocity required to overcome a given decelerating torque in an autogiro rotor. At low tip-speed ratios, where the assumptions necessary for the analytical integrations of the earlier papers are valid, the expressions therein derived are greatly simplified by combining and eliminating terms with a view of minimizing the numerical computations required. At high tip-speed ratios, by means of charts based on graphical integrations, errors inherent in the assumptions associated with the analytical method are largely eliminated. The suggested method of estimating the inflow velocity presupposes a knowledge of the decelerating torque acting on the rotor; all available full-scale experimental information on this subject is included.

Bailey, F J , Jr↗
NASA NTRS1947

Isolated and cascade airfoils with prescribed velocity distribution

An exact solution of the problem of designing an airfoil with a prescribed velocity distribution on the suction surface in a given uniform flow of an incompressible perfect fluid is obtained by replacing the boundary of the airfoil by vortices. By this device, a method of solution is developed that is applicable both to isolated airfoils and to airfoils in cascade. The conformal transformation of the designed airfoil into a circle can then be obtained and the velocity distribution at any angle of attack computed. Numerical illustrations of the method are given for the airfoil in cascade.

Goldstein, Arthur W↗
NASA NTRS1952

An empirically derived basis for calculating the area, rate, and distribution of water-drop impingement on airfoils

An empirically derived basis for predicting the area, rate, and distribution of water-drop impingement on airfoils of arbitrary section is presented. The concepts involved represent an initial step toward the development of a calculation technique which is generally applicable to the design of thermal ice-prevention equipment for airplane wing and tail surfaces. It is shown that sufficiently accurate estimates, for the purpose of heated-wing design, can be obtained by a few numerical computations once the velocity distribution over the airfoil has been determined. The calculation technique presented is based on results of extensive water-drop trajectory computations for five airfoil cases which consisted of 15-percent-thick airfoils encompassing a moderate lift-coefficient range. The differential equations pertaining to the paths of the drops were solved by a differential analyzer.

Bergrun, Norman R↗
NASA NTRS1990

An alternative formulation of the global sensitivity equations

To optimize the performance of any system, the sensitivity derivatives of the system's output variables with respect to its input variables must be readily available. It is also desirable that these derivatives be inexpensive to calculate as the optimization process requires many evaluations of the output variables and their derivatives. Optimization methods that have been developed for use in automated structural design programs may not be extended for use in integrated multidisciplinary design programs until adequate means of calculating accurate sensitivity derivatives of complex, internally coupled systems have been developed. Until the development of the Global Sensitivity Equations (GSE), the only method of determining the sensitivity derivatives of coupled systems has been by using finite differences. Analytical or semi-analytical derivatives do not exist as there is no analytical solution to the coupled problem. Also, difficulties arise because the finite difference method is expensive as the system has to iterate to a converged solution for each incremental input variable. The method may not be accurate, and the choice of the input variable increment may cause the difference in the output variable to be insignificant compared to computer numerical error if the choice is too small, or the process may not predict the true value of the output variable if the increment is too large. The GSE allow the system's sensitivity derivatives to be calculated as functions of the component subsystem's (local) sensitivity derivatives. These local sensitivity derivatives are calculated from specifically decoupled subsystems, whereas the GSE account for total system coupling. Since the subsystems are decoupled, it may be possible for the local derivatives to be calculated by analytical or semi-analytical methods, which generally reduce cost and improve accuracy. Several academic problems have been solved using GSE and have demonstrated encouraging results.

James, Benjamin B.↗
NASA NTRS1993

Sloshing dynamics on rotating helium dewar tank

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics were based on the non-inertia frame spacecraft bound coordinate, and solve time dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers were derived. Results were widely published in the open journals.

Hung, R. J.↗
NASA NTRS1993

Size distribution of interplanetary iron and stony particles related with deep-sea spherules

To study origin and evolution of the interplanetary dust, it is very important to investigate the size distribution. Here the changes of the size distributions of meteoroid particles due to the ablative effects during atmospheric entry were investigated by numerical computer simulation. Using the results, the pre-atmospheric size distributions of the interplanetary dust particles could be estimated from that of ablated spherules taken from deep-sea sediments. We are now analyzing deep-sea spherules from some aspects and examining if we could get any information about the interplanetary dust.

Matsuzaki, H.↗
NASA NTRS1993

An efficient formulation of the coupled finite element-integral equation technique for solving large 3D scattering problems

It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.

Cwik, T.↗
NASA NTRS1993

Multigrid schemes for viscous hypersonic flows

Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm employs upwind spatial discretization with explicit multistage time stepping. Two-level versions of the various multigrid algorithms are applied to the two-dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving two different hypersonic flow problems. Some new multigrid schemes, based on semicoarsening strategies, are shown to be quite effective in relieving the stiffness caused by the high-aspect-ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6).

Swanson, R. C.↗
NASA NTRS1994

Superfluid helium sloshing dynamics induced oscillations and fluctuations of angular momentum, force and moment actuated on spacecraft driven by gravity gradient or jitter acceleration associated with slew motion

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by the gravity gradient and jitter accelerations associated with slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with the Advanced X-Ray Astrophysics Facility-Spectroscopy (AXAF-S) for slew motion which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics is based on the non-inertia frame spacecraft bound coordinate, and solve time-dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid-vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers have also been derived. Examples are also given for cases applicable to the AXAF-S spacecraft sloshing dynamics associated with slew motion.

Hung, R. J.↗
NASA NTRS1994

Generation of gear tooth surfaces by application of CNC machines

This study will demonstrate the importance of application of computer numerically controlled (CNC) machines in generation of gear tooth surfaces with new topology. This topology decreases gear vibration and will extend the gear capacity and service life. A preliminary investigation by a tooth contact analysis (TCA) program has shown that gear tooth surfaces in line contact (for instance, involute helical gears with parallel axes, worm gear drives with cylindrical worms, etc.) are very sensitive to angular errors of misalignment that cause edge contact and an unfavorable shape of transmission errors and vibration. The new topology of gear tooth surfaces is based on the localization of bearing contact, and the synthesis of a predesigned parabolic function of transmission errors that is able to absorb a piecewise linear function of transmission errors caused by gear misalignment. The report will describe the following topics: description of kinematics of CNC machines with six degrees of freedom that can be applied for generation of gear tooth surfaces with new topology. A new method for grinding of gear tooth surfaces by a cone surface or surface of revolution based on application of CNC machines is described. This method provides an optimal approximation of the ground surface to the given one. This method is especially beneficial when undeveloped ruled surfaces are to be ground. Execution of motions of the CNC machine is also described. The solution to this problem can be applied as well for the transfer of machine tool settings from a conventional generator to the CNC machine. The developed theory required the derivation of a modified equation of meshing based on application of the concept of space curves, space curves represented on surfaces, geodesic curvature, surface torsion, etc. Condensed information on these topics of differential geometry is provided as well.

Litvin, F. L.↗