Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Computational Number Theory”

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 55 records · Page 3

Electron Stark Broadening Database for Atomic N, O, and C Lines

A database for efficiently computing the electron Stark broadening line widths for atomic N, O, and C lines is constructed. The line width is expressed in terms of the electron number density and electronatom scattering cross sections based on the Baranger impact theory. The state-to-state cross sections are computed using the semiclassical approximation, in which the atom is treated quantum mechanically whereas the motion of the free electron follows a classical trajectory. These state-to-state cross sections are calculated based on newly compiled line lists. Each atomic line list consists of a careful merger of NIST, Vanderbilt, and TOPbase line datasets from wavelength 50 nm to 50 micrometers covering the VUV to IR spectral regions. There are over 10,000 lines in each atomic line list. The widths for each line are computed at 13 electron temperatures between 1,000 K 50,000 K. A linear least squares method using a four-term fractional power series is then employed to obtain an analytical fit for each line-width variation as a function of the electron temperature. The maximum L2 error of the analytic fits for all lines in our line lists is about 5%.

Liu, Yen

A numerical simulation of hypersonic viscous flow over arbitrary geometries at angle of attack

An implicit conservative, noniterative, finite-difference algorithm that predicts the supersonic, laminar or turbulent viscous flow about arbitrary geometries at large angles of attack is presented. The three-dimensional parabolized form of the thin-layer Navier-Stokes equations are written in generalized coordinates. These equations are solved using the delta form of the Beam-Warming implicit algorithm. Flow field simulations have been obtained for a blunt biconic with windward and leeward cuts and an X-24C lifting body for both laminar and turbulent flow at various Mach numbers and angles of attack. When compared with experiment or with previous theories, these computational predictions show good agreement.

Chaussee, D. S.

Wall-temperature effects on the aerodynamics of a hydrogen-fueled transport concept in Mach 8 blowdown and shock tunnels

Results are presented from two separate tests on the same blended wing-body hydrogen fueled transport model at a Mach number of about 8 and a range of Reynolds numbers (based on theoretical body length) of 0.597 x 10 to the 6th power to about 156.22 x 10 to the 6th power. Tests were made in conventional hypersonic blowdown tunnel and a hypersonic shock tunnel at angles of attack of -2 deg to about 8 deg, with an extensive study made at a constant angle of attack of 3 deg. The model boundary-layer flow varied from laminar at the lower Reynolds numbers to predominantly turbulent at the higher Reynolds numbers. Model wall temperatures and stream static temperatures varied widely between the two tests, particularly at the lower Reynolds numbers. These temperature differences resulted in marked variations of the axial-force coefficients between the two tests, due in part to the effects of induced pressure and viscous interaction variations. The normal-force coefficient was essentially independent of Reynolds number. Analysis of results utilized current theoretical computer programs and basic boundary-layer theory.

Penland, J. A.

Absolute and convective instability of a cylinder wake

The stability behavior of a circular cross-section cylinder's wake at Reynolds number values of up to 45 is presently investigated by means of local linear stability theory. The steady-wake profiles computed are Navier-Stokes solutions of a uniform, incompressible viscous flow around a cylinder obtained by a spectral method. An absolutely unstable region is found to begin to form at a Reynolds number of about 20, and grows with incresing Reynolds number. The onset of global instability response must be characterized by a critical length of an absolutely unstable region; a critical Reynolds number criterion and preferred frequency are proposed based on linear stability analysis.

Yang, X.

Predicting chaos for infinite dimensional dynamical systems: The Kuramoto-Sivashinsky equation, a case study

The results of extensive computations are presented in order to accurately characterize transitions to chaos for the Kuramoto-Sivashinsky equation. In particular, the oscillatory dynamics in a window that supports a complete sequence of period doubling bifurcations preceding chaos is followed. As many as thirteen period doublings are followed and used to compute the Feigenbaum number for the cascade and so enable, for the first time, an accurate numerical evaluation of the theory of universal behavior of nonlinear systems, for an infinite dimensional dynamical system. Furthermore, the dynamics at the threshold of chaos exhibit a fractal behavior which is demonstrated and used to compute a universal scaling factor that enables the self-similar continuation of the solution into a chaotic regime.

Smyrlis, Yiorgos S.

Lifting surface theory for rectangular wings

A new incompressible lifting-surface theory is developed for thin rectangular wings. The solution requires the downwash equation to be in the form of Cauchy-type integrals. Lan's method is employed for the chordwise integrals since it properly accounts for the leading-edge singularity, Cauchy singularity and Kutta condition. The Cauchy singularity in the spanwise integral is also accounted for by using the midpoint trapezoidal rule and theory of Chebychev polynomials. The resulting matrix equation, formed by satisfying the boundary condition at control points, is simpler and quicker to compute than other lifting surface theories. Solutions were found to converge with only a small number of control points and to compare favorably with results from other methods.

Dejarnette, F. R.

The economics of time shared computing: Congestion, user costs and capacity

Time shared systems permit the fixed costs of computing resources to be spread over large numbers of users. However, bottleneck results in the theory of closed queueing networks can be used to show that this economy of scale will be offset by the increased congestion that results as more users are added to the system. If one considers the total costs, including the congestion cost, there is an optimal number of users for a system which equals the saturation value usually used to define system capacity.

Agnew, C. E.

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: a gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.

Solution-Adaptive Cartesian Cell Approach for Viscous and Inviscid Flows

A Cartesian cell-based approach for adaptively refined solutions of the Euler and Navier-Stokes equations in two dimensions is presented. Grids about geometrically complicated bodies are generated automatically, by the recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, polygonal cut cells are created using modified polygon-clipping algorithms. The grid is stored in a binary tree data structure that provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite volume formulation. The convective terms are upwinded: A linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The results of a study comparing the accuracy and positivity of two classes of cell-centered, viscous gradient reconstruction procedures is briefly summarized. Adaptively refined solutions of the Navier-Stokes equations are shown using the more robust of these gradient reconstruction procedures, where the results computed by the Cartesian approach are compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.

An Enabling Platform for Achieving Multiscale Multiphysics Analysis of Multiphase Materials

This paper will address current NASA GRC research efforts utilizing the NASA Multiscale Analysis Tool (NASMAT) which address technical gaps in three of the Vision 2040 key discipline areas. NASMAT is a versatile platform for performing computationally efficient multiscale analyses of heterogeneous materials and is available free through the NASA software catalog. It offers the user flexibility to define an arbitrary number of length scales (levels) where a variety of micromechanics theories can be implemented at each level. Micromechanics theories can be selected to balance accuracy and computational efficiency and range from analytical (Mori-Tanaka) to several semi-analytical (method of cells) formulations. The resulting anisotropic, evolving nonlinear, thermomechanical constitutive model can also be coupled with external software and used to perform multiscale analyses of more complex structures. A recent application to model a complex, three-dimensional woven composite, with a particular emphasis placed on multiscale measurements utilized to enhance the quality of the multiscale analysis will be discussed. Since typical NASMAT analyses can be completed on the order of seconds to minutes, a second example will demonstrate the ability to generate large quantities of data useful for sensitivity analysis, uncertainty quantification, or machine learning applications. Current progress on implementing NASMAT within a multiscale digital thread/digital twin framework will also be addressed. These examples will demonstrate that NASMAT is an enabling platform for achieving the goals in Vision 2040.

Multiscale

Theory and operation of the real-time data acquisition system for the NASA-LaRC differential absorption lidar (DIAL)

The improvement of computer hardware and software of the NASA Multipurpose Differential Absorption Lidar (DIAL) system is documented. The NASA DIAL system is undergoing development and experimental deployment at NASA Langley Research Center for the remote measurement of atmospheric trace gas concentrations from ground and aircraft platforms. A viable DIAL system was developed capable of remotely measuring O3 and H2O concentrations from an aircraft platform. Test flights of the DIAL system were successfully performed onboard the NASA Goddard Flight Center Electra aircraft from 1980 to 1985. The DIAL Data Acquisition System has undergone a number of improvements over the past few years. These improvements have now been field tested. The theory behind a real time computer system as it applies to the needs of the DIAL system is discussed. This report is designed to be used as an operational manual for the DIAL DAS.

Butler, C.

Kinematics and dynamics of vortex rings in a tube

Kinematic theory and flow visualization experiments were combined to examine the dynamic processes which control the evolution of vortex rings from very low to very high Reynolds numbers, and to assess the effects of the wall as a vortex ring travels up a tube. The kinematic relationships among the size, shape, speed, and strength of vortex rings in a tube were computed from the theory. Relatively simple flow visualization measurements were used to calculate the total circulation of a vortex rings at a given time. Using this method, the strength was computated and plotted as a function of time for experimentally produced vortex rings. Reynolds number relationships are established and quantitative differences among the three Reynolds number groups are discussed.

Brasseur, J. G.

Nonlinear Theory and Breakdown

The main points of recent theoretical and computational studies on boundary-layer transition and turbulence are to be highlighted. The work is based on high Reynolds numbers and attention is drawn to nonlinear interactions, breakdowns and scales. The research focuses in particular on truly nonlinear theories, i.e. those for which the mean-flow profile is completely altered from its original state. There appear to be three such theories dealing with unsteady nonlinear pressure-displacement interactions (I), with vortex/wave interactions (II), and with Euler-scale flows (III). Specific recent findings noted for these three, and in quantitative agreement with experiments, are the following. Nonlinear finite-time break-ups occur in I, leading to sublayer eruption and vortex formation; here the theory agrees with experiments (Nishioka) regarding the first spike. II gives rise to finite-distance blowup of displacement thickness, then interaction and break-up as above; this theory agrees with experiments (Klebanoff, Nishioka) on the formation of three-dimensional streets. III leads to the prediction of turbulent boundary-layer micro-scale, displacement-and stress-sublayer-thicknesses.

Smith, Frank

Effects of viscosity on transonic-aerodynamic and aeroelastic characteristics of oscillating airfoils

Studies were made to investigate the effects of viscosity on aerodynamic and aeroelastic characteristics of oscillating airfoils. The computer code LTRAN2 (viscous), which is based on a small disturbance aerodynamic theory, was used to make aerodynamic computations. Two viscous models, the viscous-ramp model and the lag-entrainment model were considered. The unsteady viscous effects were obtained by use of the quasi-steady assumptions. Two cases, a conventional airfoil, NACA 64A010, and a supercritical airfoil, MBB-A3, were studied at Mach numbers 0.796 and 0.7557, respectively. For both the airfoils, steady and unsteady aerodynamic computations were made by using inviscid and viscous theories. The steady and the unsteady results for the NACA 64A010 airfoil and the steady results for the MBB-A3 airfoil were compared with the available wind-tunnel results. Flutter speeds were computed for both airfoils by using the U-g method and the effects of viscosity on the airfoils were studied. Results from the viscous methods show improvements over the inviscid method.

Guruswamy, P.

Applications of potential theory computations to transonic aeroelasticity

Unsteady aerodynamic and aeroelastic stability calculations based upon transonic small disturbance (TSD) potential theory are presented. Results from the two-dimensional XTRAN2L code and the three-dimensional XTRAN3S code are compared with experiment to demonstrate the ability of TSD codes to treat transonic effects. The necessity of nonisentropic corrections to transonic potential theory is demonstrated. Dynamic computational effects resulting from the choice of grid and boundary conditions are illustrated. Unsteady airloads for a number of parameter variations including airfoil shape and thickness, Mach number, frequency, and amplitude are given. Finally, samples of transonic aeroelastic calculations are given. A key observation is the extent to which unsteady transonic airloads calculated by inviscid potential theory may be treated in a locally linear manner.

Edwards, J. W.

Large amplitude flutter of a low aspect ratio panel at low supersonic speeds comparison of theory and experiment

Flutter boundaries, as well as flutter limit cycle amplitudes, frequencies and stresses were computed for a panel of length-width ratio 4.48 exposed to applied in-plane and transverse loads. The Mach number range was 1.1 to 1.4. The method used involved direct numerical integration of modal equations of motion derived from the nonlinear plate equations of von Karman, coupled with linearized potential flow aerodynamic theory. The flutter boundaries agreed reasonably well with experiment, except when the in-plane loading approached the buckling load. Structural damping had to be introduced, to produce frequencies comparable to the experimental values. Attempts to compute panel deflections or stress at a given point met with limited success. There is some evidence, however, that deflection and stress maxima can be estimated with somewhat greater accuracy.

Ventres, C. S.

Direct simulation of a turbulent oscillating boundary layer

The turbulent boundary layer driven by a freestream velocity that varies sinusoidally in time around a zero mean is considered. The flow has a rich behavior including strong pressure gradients, inflection points, and reversal. A theory for the velocity and stress profiles at high Reynolds number is formulated. Well-resolved direct Navier-Stokes simulations are conducted over a narrow range of Reynolds numbers, and the results are compared with the theoretical predictions. The flow is also computed over a wide range of Reynolds numbers using a new algebraic turbulence model; the results are compared with the direct simulations and the theory.

Spalart, Philippe R.