Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hyperbolic”

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 523 records · Page 29

Study of multi-kilowatt solar arrays for Earth orbit applications

A miniaturized Cassegrainian concentrator (MCC) solar array concept is being developed with the objective of significantly reducing the recurring cost of multikilowatt solar arrays. The desired cost reduction is obtained as a result of using very small high efficiency solar cells in conjuction with low cost optics. The MCC single element concept incident slar radiation is reflected rom a primary parabolic reflector to a secondary hyperbolic reflector and finally to a 4 millimeter diameter solar cell. A light catcher cone is used to improve off axis performance. The solar cell is mounted to a heat fin. An element is approximately 13 millimeters thick which permits efficient launch stowage of the concentrator system panels without complex optical component deployments or retractions. The MCC elements are packed in bays within graphite epoxy frames and are electrically connected into appropriate series-parallel circuits. A MCC sngle element with a 21 sq cm entrance aperture and a 20 efficient, 0.25 sq cm gallium arsenide solar cell has the same power output as 30 sq cm of 11-percent efficiency (at 68 C) silicon solar cells.

Patterson, R. E.↗

Computational methods for ideal compressible flow

Conservative dissipative difference schemes for computing one dimensional flow are introduced, and the recognition and representation of flow discontinuities are discussed. Multidimensional methods are outlined. Second order finite volume schemes are introduced. Conversion of difference schemes for a single linear convection equation into schemes for the hyperbolic system of the nonlinear conservation laws of ideal compressible flow is explained. Approximate Riemann solvers are presented. Monotone initial value interpolation; and limiters, switches, and artificial dissipation are considered.

Vanleer, B.↗

Universal single level implicit algorithm for gasdynamics

A single level effectively explicit implicit algorithm for gasdynamics is presented. The method meets all the requirements for unconditionally stable global iteration over flows with mixed supersonic and supersonic zones including blunt body flow and boundary layer flows with strong interaction and streamwise separation. For hyperbolic (supersonic flow) regions the method is automatically equivalent to contemporary space marching methods. For elliptic (subsonic flow) regions, rapid convergence is facilitated by alternating direction solution sweeps which bring both sets of eigenvectors and the influence of both boundaries of a coordinate line equally into play. Point by point updating of the data with local iteration on the solution procedure at each spatial step as the sweeps progress not only renders the method single level in storage but, also, improves nonlinear accuracy to accelerate convergence by an order of magnitude over related two level linearized implicit methods. The method derives robust stability from the combination of an eigenvector split upwind difference method (CSCM) with diagonally dominant ADI(DDADI) approximate factorization and computed characteristic boundary approximations.

Lombard, C. K.↗

Stability of finite difference models containing two boundaries or interfaces

The stability of finite difference models of hyperbolic initial boundary value problems is connected with the propagation and reflection of parasitic waves. Wave propagation ideas are applied to models containing two boundaires or interfaces, where repeated reflection of trapped wave packets is a potential new source of instability. Various known instability phenomena are accounted for in a unified way. Results show: (1) dissipativity does not ensure stability when three or more formulas are concatenated at a boundary or internal interface; (2) algebraic GKS instabilities can be converted by a second boundary to exponential instabilities only when an infinite numerical reflection coefficient is present; and (3) GKS-stability and P-stability can be established in certain problems by showing that all numerical reflection coefficients have modulus less than 1.

Trefethen, L. N.↗

Spectral methods for partial differential equations

Origins of spectral methods, especially their relation to the Method of Weighted Residuals, are surveyed. Basic Fourier, Chebyshev, and Legendre spectral concepts are reviewed, and demonstrated through application to simple model problems. Both collocation and tau methods are considered. These techniques are then applied to a number of difficult, nonlinear problems of hyperbolic, parabolic, elliptic, and mixed type. Fluid-dynamical applications are emphasized.

Hussaini, M. Y.↗

A computational method for viscous incompressible flows

An implicit, finite-difference procedure for numerically solving viscous incompressible flows is presented. The pressure-field solution is based on the pseudocompressibility method in which a time-derivative pressure term is introduced into the mass-conservation equation to form a set of hyperbolic equations. The pressure-wave propagation and the spreading of the viscous effect is investigated using simple test problems. Computed results for external and internal flows are presented to verify the present method which has proved to be very robust in simulating incompressible flows.

Kwak, D.↗

Analysis of lateral and torsional vibration characteristics of beams and shafts with end located rotational masses

Partial differential equations are derived for free lateral and torsional vibration of a uniform free-free beam with a rotational mass attached to each extremity. For appropriate boundary conditions, nonlinear algebraic equations are obtained using a symbolic manipulation computer program, the solutions of which enable the computation of the neutral frequencies and mode-shapes. The mode-shapes are linear combinations of trigonometric and hyperbolic sine and cosine functions. A computer program is written for the numerical solution of the algebraic equations mentioned above, which can compute the natural frequencies, mode-shapes, and node points for any given set of parameters, for any given number of modes.

Robertson, D. K.↗

Deployment of a Secondary Concentrator to Increase the Intercept Factor of a Dish with Large Slope Errors

The testing of a hyperbolic trumpet non-imaging secondary concentrator with a parabolic dish having slope errors of about 10 mrad is reported. The trumpet, which has a concentration ratio of 2.1, increased the flux through a 141-mm focal aperture by 72%, with an efficiency of 96%, thus demonstrating its potential for use in tandem with cheap dishes having relatively large slope errors.

Ortabasi, U.↗

Enthalpy damping for the steady Euler equations

For inviscid steady flow problems where the enthalpy is constant at steady state, it was previously proposed to use the difference between the local enthalpy and the steady state enthalpy as a driving term to accelerate convergence of iterative schemes. This idea is analyzed, both on the level of the partial differential equation and on the level of a particular finite difference scheme. It is shown that for the two-dimensional unsteady Euler equations, a hyperbolic system with eigenvalues on the imaginary axis, there is no enthalpy damping strategy which moves all the eigenvalues into the open left half plane. For the numerical scheme, however, the analysis shows and examples verify that enthalpy damping is potentially effective in accelerating convergence to steady state.

Jespersen, D. C.↗

Second-order explicit finite-difference methods for transient-flow analysis

Three second-order accurate numerical methods - MacCormack's method, Lambda scheme and Gabutti scheme - are introduced to solve the quasi-linear, hyperbolic partial differential equations describing transient flows in closed conduits. The details of these methods and the treatment of boundary conditions are presented and the results computed by using these methods for a typical piping system are compared. It is shown that for the same accuracy, second-order methods require considerably lesser number of computational nodes and computer time as compared to those required by the first-order methods.

Chaudhry, M. H.↗

Numerical approximations to nonlinear conservation laws with locally varying time and space grids

Numerical approximations to the initial value problem for nonlinear systems of conservation laws are considered. The considered system is said to be hyperbolic when all eigenvalues of every real linear combination of the Jacobian matrices are real. Solutions may develop discontinuities in finite time, even when the initial data are smooth. In the investigation, explicit finite difference methods which use locally varying time grids are considered. The global CFL restriction is replaced by a local restriction. The numerical flux function is studied from a finite volume viewpoint, and a differencing technique is developed at interface points between regions of distinct time increments.

Osher, S.↗

An incompressible Navier-Stokes flow solver in three-dimensional curvilinear coordinate systems using primitive variables

An implicit, finite-difference computer code has been developed to solve the incompressible Navier-Stokes equations in a three-dimensional, curvilinear coordinate system. The pressure-field solution is based on the pseudo compressibility approach in which the time derivative pressure term is introduced into the mass conservation equation to form a set of hyperbolic equations. The solution procedure employs an implicit, approximate factorization scheme. The Reynolds stresses, that are uncoupled from the implicit scheme, are lagged by one time-step to facilitate implementing various levels of the turbulence model. Test problems for external and internal flows are computed, and the results are compared with existing experimental data. The application of this technique for general three-dimensional problems is then demonstrated.

Kwak, D.↗

Relaxation methods for unfactored implicit upwind schemes

Relaxation methods are presented for unfactored implicit upwind schemes for hyperbolic equations. The theoretical bases are explained using linear and nonlinear scalar equations; construction of the method for the unsteady Euler equations (nonlinear system) is but a natural extension. One of the important advantages of the above methods vis a vis factored implicit schemes is the possibility of faster convergence to steady state, as illustrated by the results. Several classes of relaxation schemes such as pointwise, linewise, Gauss-Seidel, and non-Gauss-Seidel methods are discussed, along with various strategies for convergence.

Chakravarthy, S. R.↗

On the method of pseudo compressibility for numerically solving incompressible flows

Pseudo compressibility is used for numerically solving incompressible flows to achieve computational efficiency. The use of pseudo compressibility results in a system of hyperbolic-type equations of motion that introduce waves of finite speed. The interactions of the wave propagation and the vorticity spreading are analyzed. A criterion governing the dependence of the pseudo compressiblity on the Reynolds number and the characteristic length of the flow geometry is obtained that allows for a proper convergence. It is demonstrated that the solution does tend to the incompressible limit. External and internal viscous flow test problems are presented to verify the theory.

Chang, J. L. C.↗

Wrinkling phenomenon in structures. II - Wrinkling criteria

The wrinkling phenomenon in elasticity is mathematically defined and described as a special type of instability associated with the loss of hyperbolicity of the governing equations when some of the characteristic speeds of elastic wave propagation become imaginary. The developed criteria of wrinkling are applied to such structural elements as membranes, strings, soft shells, soft rods, and laminated material. The cause of wrinkling in these elements is due to compression in longitudinal directions that exceeds the through-the-thickness shear modulus. For elements with internal fluid flow the critical compression is reduced due to the destabilizing effect of the flow. The most effective mathematical analysis technique is a combination of tensor analysis and index notations.

Zak, M.↗

A comparative study of the parabolized Navier-Stokes code using various grid-generation techniques

The parabolized Navier-Stokes (PNS) equations are used to calculate the flow-field characteristics about the hypersonic research aircraft X-24C. A comparison of the results obtained using elliptic, hyperbolic, and algebraic grid generators is presented. The outer bow shock is treated as a sharp discontinuity, and the discontinuities within the shock layer are captured. Surface pressures and heat-transfer results at angles of attack of 6 deg and 20 deg, obtained using the three grid generators, are compared. The PNS equations are marched downstream over the body in both Cartesian and cylindrical base coordinate systems, and the results are compared. A robust marching procedure is demonstrated by successfully using large marching step sizes with the implicit shock fitting procedure. A correlation is found between the marching step size, Reynolds number, and the angle of attack at fixed values of smoothing and stability coefficients for the marching scheme.

Kaul, U. K.↗

Dynamical evolution of the Oort cometary cloud

The dynamical evolution of comets in the Oort cloud under the influence of stellar perturbations has been modeled using Monte Carlo techniques. It is shown that the cloud has been depleted over the history of the solar system. Comets are lost from the cloud by direct ejection due to close stellar encounters, diffusion of aphelia to distances beyond the sun's sphere of influence, or diffusion of perihelia into the planetary region where Jupiter and Saturn perturbations either eject them on hyperbolic trajectories or capture them to short-period orbits. The population of the cloud is estimated to be 1.0 - 1.5 x 10 to the 12th comets and the total mass is on the order of 1.9 earth masses. In addition to random passing stars, less frequent encounters with giant molecular clouds may play a significant role in randomizing the orbits of comets in the cloud and reducing the effective radius of the sun's sphere of influence.

Weissman, P. R.↗

A minimal residual method for transonic potential flows

For transonic flow calculations, a combination of the successive line over-relaxation (SLOR) and the preconditioned conjugate gradient (CG) method has been suggested by Wong and Hafez (1981). This paper studies the method of minimal residual (MR) which avoids a combined iteration. This method is closely related to the CG method, may be regarded as a first-order gradient method, and is applicable to symmetric and nonsymmetric matrices. The problem is formulated as a nonlinear mixed elliptic-hyperbolic partial differential equation which includes an artificial viscosity and a switching function which is zero in subsonic regions and nonzero in supersonic regions. Alternatives to the SLOR method which provide faster convergence rates are introduced. The preconditioned MR algorithm is developed, and transonic potential flows around NACA 0012 airfoil are calculated for different Mach numbers and angles of attack. Preliminary results are presented, demonstrating that the MR algorithm requires no parameter estimation and rapidly converges for subsonic flows.

Wong, Y. S.↗