Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical algorithm computing”

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 883 records · Page 49

Numerical pole assignment by eigenvalue Jacobian inversion

A numerical procedure for solving the linear pole placement problem is developed which operates by the inversion of an analytically determined eigenvalue Jacobian matrix. Attention is given to convergence characteristics and pathological situations. It is not concluded that the algorithm developed is suitable for computer-aided control system design with particular reference to the scan platform pointing control system for the Galileo spacecraft.

Sevaston, George E.↗

Numerical simulation of a subsonic jet in a crossflow

The aerodynamic/propulsive interaction between a subsonic jet exhausting perpendicularly through a flat plate into a crossflow is investigated numerically using an approximately factored, partially flux-split, implicit solver for the three-dimensional, thin-layer Navier-Stokes equations. This algorithm is applied to flows with a range of jet-to-crossflow velocity ratios between 4 and 8. The computations model the jet trajectory, the contrarotating vortex pair and the wake region near the plate downstream of the jet orifice. Both qualitative and quantitative agreement with the existing experimental database are demonstrated. Flow visualization is instructive for understanding the physics of this flowfield.

Roth, Karlin R.↗

Vibration of microgravity environment for rotating fluids in high and low rotating speeds

Time dependent evolutions of the profile of free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, have been studied. Numerical computations of the dynamics of bubble shapes have been carried out with the sinusoidal function vibration of gravity environment in high and low rotating cylinder speeds. The initial condition of bubble profiles was adopted from the steady-state formulations in which the computer algorithms have been developed by Hung and Leslie (1988), and Hung et al. (1988).

Hung, R. J.↗

Using deflation in the pole assignment problem with output feedback

A direct algorithm is suggested for the computation of a linear output feedback for a multi input, multi output system such that the resultant closed-loop matrix has eigenvalues that include a specified set of eigenvalues. The algorithm uses deflation based on unitary similarity transformations. Thus researchers hope the algorithm is numerically stable; however, this has not been proven as yet.

Miminis, George↗

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

Evaluation of a Multizone Impedance Eduction Method

A computational study is used to evaluate the PyCHE impedance eduction method developed at the NASA Langley Research Center. This method combines an aeroacoustic duct propagation code based on numerical solution to the convected Helmholtz equation with a global optimizer that uses the Differential Evolution algorithm. The efficacy of this method is evaluated with acoustic pressure data simulated to represent that measured with one-zone, two-zone, and three-zone liners mounted in the NASA Langley Grazing Flow Impedance Tube. The PyCHE method has a normalized impedance error of approximately 0.2 for (uniform) one-zone liners with a length of at least 5”, and produces quite reasonable results for liners as short as 2”. Whereas the impedance of the liner has an effect on eduction accuracy, the amount of attenuation is shown to be the dominant parameter. Similar results are observed for two-zone liners, for which the impedance of each zone is unique. The two-zone results also indicate it is more difficult to accurately educe resistance than reactance, and a zone length of at least 6” (slightly longer than for uniform liners) is needed to limit the normalized error to 0.2. The PyCHE method is also demonstrated to successfully educe the impedances for each zone of a three-zone liner. These results are sufficiently encouraging to warrant the continued usage of the PyCHE impedance eduction method for single and multizone liners.

Jones, M. G.↗

The change in critical technologies for computational physics

It is noted that the types of technology required for computational physics are changing as the field matures. Emphasis has shifted from computer technology to algorithm technology and, finally, to visual analysis technology as areas of critical research for this field. High-performance graphical workstations tied to a supercommunicator with high-speed communications along with the development of especially tailored visualization software has enabled analysis of highly complex fluid-dynamics simulations. Particular reference is made here to the development of visual analysis tools at NASA's Numerical Aerodynamics Simulation Facility. The next technology which this field requires is one that would eliminate visual clutter by extracting key features of simulations of physics and technology in order to create displays that clearly portray these key features. Research in the tuning of visual displays to human cognitive abilities is proposed. The immediate transfer of technology to all levels of computers, specifically the inclusion of visualization primitives in basic software developments for all work stations and PCs, is recommended.

Watson, Val↗

Sonoluminescence Bubble Measurements using Vision-Based Algorithms

Vision-based measurement methods were used to measure bubble sizes in this sonoluminescence experiment. Bubble imaging was accomplished by placing the bubble between a bright light source and a microscope-CCD camera system. A collimated light-emitting diode was operated in a pulsed model with an adjustable time delay with respect to the piezo-electric transducer drive signal. The light-emitting diode produced a bubble shadowgraph consisting of a multiple exposure made by numerous light pulses imaged onto a charge-couple device camera. Each image was transferred from the camera to a computer-controlled machine vision system via a frame grabber. The frame grabber was equipped with on-board memory to accomodate sequential image buffering while images were transferred to the host processor and analyzed. This configuration allowed the host computer to perform diameter measurements, centroid position measurements and shape estimation in "real-time" as the next image was being acquired. Bubble size measurement accuracy with an uncertainty of 3 microns was achieved using standard lenses and machine vision algorithms. Bubble centroid position accuracy was also within the 3 micron tolerance of the vision system. This uncertainty estimation accounted for the optical spatial resolution, digitization errors and the edge detection algorithm accuracy. The vision algorithms include camera calibration, thresholding, edge detection, edge position determination, distance between two edges computations and centroid position computations.

Hall, Nancy R.↗

An alternating direction implicit method for the Control Data STAR-100 vector computer

An implementation of the alternating direction implicit (ADI) method for the Control Data STAR-100 computer is presented and analyzed. Two parallel algorithms, both of which are most efficient when used to solve many independent tridiagonal systems of equations, are discussed relative to their usefulness in an ADI implementation on the STAR-100 computer. It is shown that it may be desirable to alternate between the parallel algorithms as the direction of implicitness is alternated in order to eliminate the data rearrangement which would otherwise be required. The applicability of the two parallel tridiagonal solvers to several other numerical algorithms is also discussed.

Lambiotte, J. J., Jr.↗

Numerical predictions of residence times behind a rearward facing step with transverse injection

A numerical investigation of the two-dimensional recirculation zone formed between a rear-ward facing step and a transverse hydrogen jet in a supersonic flow is presented. Calculations were made using an elliptic flow computer code with a time-relaxation finite difference integration algorithm. Cases were examined for jet static pressures between 2 and 5 times mainstream with jets located 3 to 7 step heights downstream of the step. Nominal mainstream flow conditions were Mach 2.4 air at 1 atm pressure and 1000K (1800R) static temperature. Results in the form of mean properties in the recirculation zone and residence time are presented. Calculations of the decay of hydrogen mass from the steady state are presented for selected cases. The combustion potential for each case is assessed by comparison of computed conditions in the recirculation zone with the results from a well-stirred reactor theory.

Rogers, R. C.↗

Chandrasekhar equations and computational algorithms for distributed parameter systems

The Chandrasekhar equations arising in optimal control problems for linear distributed parameter systems are considered. The equations are derived via approximation theory. This approach is used to obtain existence, uniqueness, and strong differentiability of the solutions and provides the basis for a convergent computation scheme for approximating feedback gain operators. A numerical example is presented to illustrate these ideas.

Burns, J. A.↗

Solving finite element equations on concurrent computers

This paper discusses the development of a concurrent algorithm for the solution of systems of equations arising in finite element applications. The approach is based on a hybrid of direct elimination method and preconditioned conjugate iteration. Two different preconditioners are used; diagonal scaling and a concurrent implementation of incomplete LU factorization. First, an automatic procedure is used to partition the finite element mesh into sub-structures. The particular mesh partition is chosen to minimize an estimate of the cost for evaluating the solution using this algorithm on a concurrent computer. These procedures are implemented in a finite element program on the JPL/CalTech MARK III hypercube computer. An overview of the structure of this program is presented. The performance of the solution method is demonstrated with the aid of a number of numerical test runs, and its advantages for concurrent implementations are discussed. Efficiency and speed-up factors over sequential machines for the numerical examples are highlighted.

Nour-Omid, B.↗

A forward method for optimal stochastic nonlinear and adaptive control

A computational approach is taken to solve the optimal nonlinear stochastic control problem. The approach is to systematically solve the stochastic dynamic programming equations forward in time, using a nested stochastic approximation technique. Although computationally intensive, this provides a straightforward numerical solution for this class of problems and provides an alternative to the usual dimensionality problem associated with solving the dynamic programming equations backward in time. It is shown that the cost degrades monotonically as the complexity of the algorithm is reduced. This provides a strategy for suboptimal control with clear performance/computation tradeoffs. A numerical study focusing on a generic optimal stochastic adaptive control example is included to demonstrate the feasibility of the method.

Bayard, David S.↗

Analysis of the anisotropic viscoplastic-damage response of composite laminates - Continuum basis and computational algorithms

The mathematical structure underlying the rate equations of a recently-developed constitutive model for the coupled viscoplastic-damage response of anisotropic composites is critically examined. In this regard, a number of tensor projection operators have been identified, and their properties were exploited to enable the development of a general computational framework for their numerical implementation using the Euler fully-implicit integration method. In particular, this facilitated (i) the derivation of explicit expressions of the (consistent) material tangent stiffnesses that are valid for both three-dimensional as well as subspace (e.g. plane stress) formulations, (ii) the implications of the symmetry or unsymmetry properties of these tangent operators from a thermodynamic standpoint, and (iii) the development of an effective time-step control strategy to ensure accuracy and convergence of the solution. In addition, the special limiting case of inviscid elastoplasticity is treated. The results of several numerical simulations are given to demonstrate the effectiveness of the schemes developed.

Saleeb, A. F.↗

Computer Modeling of Electrostatic Aggregation of Granular Materials in Planetary and Astrophysical Settings

Electrostatic forces strongly influence the behavior of granular materials in both dispersed (cloud) systems and semi-packed systems. These forces can cause aggregation or dispersion of particles and are important in a variety of astrophysical and planetary settings. There are also many industrial and commercial settings where granular matter and electrostatics become partners for both good and bad. This partnership is important for human exploration on Mars where dust adheres to suits, machines, and habitats. Long-range Coulombic (electrostatic) forces, as opposed to contact-induced dipoles and van der Waals attractions, are generally regarded as resulting from net charge. We have proposed that in addition to net charge interactions, randomly distributed charge carriers on grains will result in a dipole moment regardless of any net charge. If grains are unconfined, or fluidized, they will rotate so that the dipole always induces attraction between grains. Aggregates are readily formed, and Coulombic polarity resulting from the dipole produces end-to-end stacking of grains to form filamentary aggregates. This has been demonstrated in USML experiments on Space Shuttle where microgravity facilitated the unmasking of static forces. It has also been demonstrated in a computer model using grains with charge carriers of both sign. Model results very closely resembled micro-g results with actual sand grains. Further computer modeling of the aggregation process has been conducted to improve our understanding of the aggregation process, and to provide a predictive tool for microgravity experiments slated for Space Station. These experiments will attempt to prove the dipole concept as outlined above. We have considerably enhanced the original computer model: refinements to the algorithm have improved the fidelity of grain behavior during grain contact, special attention has been paid to simulation time steps to enable establishment of a meaningful, quantitative time axis, and calibration of rounding accuracies have been conducted to test cumulative numerical influences in the model. The model has been run for larger grain populations, variable initial cloud densities, and we have introduced random net charging to individual grains, as well as a net charge to the cloud as a whole. The model uses 3 positive and 3 negative charges randomly distributed on each grain, with up to 160 grains contained within various size "boxes" that define the initial number densities in the clouds. Each charge represents localized charged region on a grain, but does not necessarily imply single quantized charge carriers. The Coulomb equations are then allowed to interact for each monopole: dipoles and any higher order charge coupling is a natural product of these "free" interactions over which the modeler exerts no influence. The charges are placed on surfaces of grains at random locations. A series of runs was conducted for neutral grains that had a perfect balance of negative and positive char carriers. Runs were also conducted with grains having additional fractional charges ranging between 0 and 1. By adding fractional charges of one sign, the model created grain populations in which all grains had excess charges the same sign, giving the cloud an overall net charge. This simulates clouds subjected to ionizing radiation (e. protoplanetary debris disk around a protosun), or any other process of charge biasing in a grain population (e.g., volcanic plumes). In another run series, random fractional charges of either sign were added to the grains so th some grains had a slight net positive charge while others had a slight net negative charge. This simulates triboelectrically-charged grain populations in which acquisition of an electron by one surface is at the expense creating a hole elsewhere. This dual sign charging was applied in two ways: in one case the cloud remained neutral by ensuring that all grain excess charges added to zero; in the other case, the cloud was permitted slight net char by not imposing a charge-balance condition. Additional information is contained in the original.

Marshall, J.↗

Numerical computation of multishocked, three-dimensional supersonic flow fields with real gas effects.

A computational procedure is presented which is capable of determining the supersonic flow field surrounding three-dimensional wing-body configurations such as a delta-wing space shuttle. The governing equations in conservation-law form are solved by a finite difference method using a second-order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically, and the intersection of these shocks with the bow shock posed no difficulty. Resulting flow fields about typical blunt nose shuttle-like configurations at angle of attack are presented. The differences between perfect and real gas effects for high Mach number flows are shown.

Kutler, P.↗

Robust integration schemes for generalized viscoplasticity with internal-state variables. Part 1: Theoretical developments and applications

This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Saleeb, Atef F.↗

A PANSONIC Navier-Stokes solver

A finite-difference formulation of the full Navier-Stokes equations which demonstrates a capability to economically solve two-dimensional problems has been developed. The basic algorithm was derived from the full, Reynolds-averaged, conservative, Navier-Stokes equations expressed in curvilinear coordinates. Eddy viscosity was determined by the Baldwin and Lomax algebraic turbulence model. This non-iterative, second-order accurate, implicit, numerical algorithm is based on the approximate factorization finite-difference scheme of Beam and Warming. Results indicate a facility for solving subsonic, transonic, and supersonic (hence PANSONIC) flows about arbitrary airfoils for a wide range of Reynolds numbers, Mach numbers, and angles of attack. Current computations demonstrate that vectorized implementations of this algorithm can solve steady-state, two-dimensional problems in five to ten minutes of computer time.

Cooper, G. K.↗