Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer Based Procedures”

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 145 records · Page 8

The effects of Reynolds number, rotor incidence angle, and surface roughness on the heat transfer distribution in a large-scale turbine rotor passage

A combined experimental and computational program was conducted to examine the heat transfer distribution in a turbine rotor passage geometrically similiar to the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump (HPFTP). Heat transfer was measured and computed for both the full-span suction and pressure surfaces of the rotor airfoil as well as for the hub endwall surface. The primary objective of the program was to provide a benchmark-quality data base for the assessment of rotor passage heat transfer computational procedures. The experimental portion of the study was conducted in a large-scale, ambient temperature, rotating turbine model. Heat transfer data were obtained using thermocouple and liquid-crystal techniques to measure temperature distributions on the thin, electrically-heated skin of the rotor passage model. Test data were obtained for various combinations of Reynolds number, rotor incidence angle and model surface roughness. The data are reported in the form of contour maps of Stanton number. These heat distribution maps revealed numerous local effects produced by the three-dimensional flows within the rotor passage. Of particular importance were regions of local enhancement produced on the airfoil suction surface by the main-passage and tip-leakage vortices and on the hub endwall by the leading-edge horseshoe vortex system. The computational portion consisted of the application of a well-posed parabolized Navier-Stokes analysis to the calculation of the three-dimensional viscous flow through ducts simulating the a gas turbine passage. These cases include a 90 deg turning duct, a gas turbine cascade simulating a stator passage, and a gas turbine rotor passage including Coriolis forces. The calculated results were evaluated using experimental data of the three-dimensional velocity fields, wall static pressures, and wall heat transfer on the suction surface of the turbine airfoil and on the end wall. Particular attention was paid to an accurate modeling of the passage vortex and to the development of the wall boundary layers including crossflow.

Blair, Michael F.↗

Direct computation of turbulence and noise

Jet exhaust turbulence noise is computed using a time dependent solution of the three dimensional Navier-Stokes equations to supply the source terms for an acoustic computation based on the Phillips convected wave equation. An extrapolation procedure is then used to determine the far field noise spectrum in terms of the near field sound. This will lay the groundwork for studies of more complex flows typical of noise suppression nozzles.

Berman, C.↗

Numerical Determination of Critical Conditions for Thermal Ignition

The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.

Luo, W.↗

Numerical Investigation of the Interaction of Counterflowing Jets and Supersonic Capsule Flows

Use of counterflowing jets ejected into supersonic freestreams as a flow control concept to modify the external flowfield has gained renewed interest with regards to potential retropropulsion applications pertinent to entry, descent, and landing investigations. This study describes numerical computations of such a concept for a scaled wind-tunnel capsule model by employing the space-time conservation element solution element viscous flow solver with unstructured meshes. Both steady-state and time-accurate computations are performed for several configurations with different counterflowing jet Mach numbers. Axisymmetric computations exploring the effect of the jet flow rate and jet Mach number on the flow stability, jet interaction with the bow shock and its subsequent impact on the aerodynamic and aerothermal loads on the capsule body are carried out. Similar to previous experimental findings, both long and short penetration modes exist at a windtunnel Mach number of 3.48. It was found that both modes exhibit non-stationary behavior and the former is much more unstable than the latter. It was also found that the unstable long penetration mode only exists in a relatively small range of the jet mass flow rate. Solution-based mesh refinement procedures are used to improve solution accuracy and provide guidelines for a more effective mesh generation procedure for parametric studies. Details of the computed flowfields also serve as a means to broaden the knowledge base for future retropropulsion design studies.

Venkatachari, Balaji Shankar↗

Minimization of the vibration energy of thin-plate structure

An optimization method is proposed to reduce the vibration of thin plate structures. The method is based on a finite element shell analysis, a modal analysis, and a structural optimization method. In the finite element analysis, a triangular shell element with 18 dof is used. In the optimization, the overall vibration energy of the structure is adopted as the objective function, and it is minimized at the given exciting frequency by varying the thickness of the elements. The technique of modal analysis is used to derive the sensitivity of the vibration energy with respect to the design variables. The sensitivity is represented by the sensitivities of both eigenvalues and eigenvectors. The optimum value is computed by the gradient projection method and a unidimensional search procedure under the constraint condition of constant weight. A computer code, based on the proposed method, is developed and is applied to design problems using a beam and a plate as test cases. It is confirmed that the vibration energy is reduced at the given exciting frequency. For the beam excited by a frequency slightly less than the fundamental natural frequency, the optimized shape is close to the beam of uniform strength.

Inoue, Katsumi↗

An experimental modal testing/identification technique for personal computers

A PC-based system for mode shape identification is evaluated. A time-domain modal identification procedure is utilized to identify the mode shapes of a beam apparatus from discrete time-domain measurements. The apparatus includes a cantilevered aluminum beam, four accelerometers, four low-pass filters, and the computer. The method's algorithm is comprised of an identification algorithm: the Eigensystem Realization Algorithm (ERA) and an estimation algorithm called Minimum Model Error (MME). The identification ability of this algorithm is compared with ERA alone, a frequency-response-function technique, and an Euler-Bernoulli beam model. Detection of modal parameters and mode shapes by the PC-based time-domain system is shown to be accurate in an application with an aluminum beam, while mode shapes identified by the frequency-domain technique are not as accurate as predicted. The new method is shown to be significantly less sensitive to noise and poorly excited modes than other leading methods. The results support the use of time-domain identification systems for mode shape prediction.

Roemer, Michael J.↗

Simulation System for Training in Laparoscopic Surgery

A computer-based simulation system creates a visual and haptic virtual environment for training a medical practitioner in laparoscopic surgery. Heretofore, it has been common practice to perform training in partial laparoscopic surgical procedures by use of a laparoscopic training box that encloses a pair of laparoscopic tools, objects to be manipulated by the tools, and an endoscopic video camera. However, the surgical procedures simulated by use of a training box are usually poor imitations of the actual ones. The present computer-based system improves training by presenting a more realistic simulated environment to the trainee. The system includes a computer monitor that displays a real-time image of the affected interior region of the patient, showing laparoscopic instruments interacting with organs and tissues, as would be viewed by use of an endoscopic video camera and displayed to a surgeon during a laparoscopic operation. The system also includes laparoscopic tools that the trainee manipulates while observing the image on the computer monitor (see figure). The instrumentation on the tools consists of (1) position and orientation sensors that provide input data for the simulation and (2) actuators that provide force feedback to simulate the contact forces between the tools and tissues. The simulation software includes components that model the geometries of surgical tools, components that model the geometries and physical behaviors of soft tissues, and components that detect collisions between them. Using the measured positions and orientations of the tools, the software detects whether they are in contact with tissues. In the event of contact, the deformations of the tissues and contact forces are computed by use of the geometric and physical models. The image on the computer screen shows tissues deformed accordingly, while the actuators apply the corresponding forces to the distal ends of the tools. For the purpose of demonstration, the system has been set up to simulate the insertion of a flexible catheter in a bile duct. [As thus configured, the system can also be used to simulate other endoscopic procedures (e.g., bronchoscopy and colonoscopy) that include the insertion of flexible tubes into flexible ducts.] A hybrid approach has been followed in developing the software for real-time simulation of the visual and haptic interactions (1) between forceps and the catheter, (2) between the forceps and the duct, and (3) between the catheter and the duct. The deformations of the duct are simulated by finite-element and modalanalysis procedures, using only the most significant vibration modes of the duct for computing deformations and interaction forces. The catheter is modeled as a set of virtual particles uniformly distributed along the center line of the catheter and connected to each other via linear and torsional springs and damping elements. The interactions between the forceps and the duct as well as the catheter are simulated by use of a ray-based haptic-interaction- simulating technique in which the forceps are modeled as connected line segments.

Basdogan, Cagatay↗

Progress in Parallel Schur Complement Preconditioning for Computational Fluid Dynamics

We consider preconditioning methods for nonself-adjoint advective-diffusive systems based on a non-overlapping Schur complement procedure for arbitrary triangulated domains. The ultimate goal of this research is to develop scalable preconditioning algorithms for fluid flow discretizations on parallel computing architectures. In our implementation of the Schur complement preconditioning technique, the triangulation is first partitioned into a number of subdomains using the METIS multi-level k-way partitioning code. This partitioning induces a natural 2X2 partitioning of the p.d.e. discretization matrix. By considering various inverse approximations of the 2X2 system, we have developed a family of robust preconditioning techniques. A computer code based on these ideas has been developed and tested on the IBM SP2 and the SGI Power Challenge array using MPI message passing protocol. A number of example CFD calculations will be presented to illustrate and assess various Schur complement approximations.

Barth, Timothy J.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Finite dynamic element formulation for a plane triangular element

Higher order correction terms for the stiffness and inertia matrices associated with a triangular plane stress-strain finite dynamic element are developed in detail. Numerical results are presented which indicate that the adoption of these matrices along with a suitable quadratic matrix eigenproblem solver effects a significant economy in the free vibration solution of structure when compared with the analysis based on the usual finite element procedure. Finally, a FORTRAN IV computer program listing of the various relevant matrices is given.

Gupta, K. K.↗

Computer-aided design of bevel gear tooth surfaces

This paper presents a computer-aided design procedure for generating bevel gears. The development is based on examining a perfectly plastic, cone-shaped gear blank rolling over a cutting tooth on a plane crown rack. The resulting impression on the plastic gear blank is the envelope of the cutting tooth. This impression and envelope thus form a conjugate tooth surface. Equations are presented for the locus of points on the tooth surface. The same procedures are then extended to simulate the generation of a spiral bevel gear. The corresponding governing equations are presented.

Shuo, Hung Chang↗

Computer aided design of bevel gear tooth surfaces

This paper presents a computer-aided design procedure for generating bevel gears. The development is based on examining a perfectly plastic, cone-shaped gear blank rolling over a cutting tooth on a plane crown rack. The resulting impression on the plastic gear blank is the envelope of the cutting tooth. This impression and envelope thus form a conjugate tooth surface. Equations are presented for the locus of points on the tooth surface. The same procedures are then extended to simulate the generation of a spiral bevel gear. The corresponding governing equations are presented.

Chang, S. H.↗

Interaction dynamics of on-orbit construction

Deployment and assembly of large structures in orbit is a critical technology to the overall problem of orbital construction. The attendant large configuration changes of structures will cause significant changes in the dynamic characteristics of the entire system, and perturbation to the orbital dynamics of the spacecraft from which the structures are deployed and/or assembled. To better design structures for deployment and assembly, and to better design controlled deployment/assembly processes, accurate modeling techniques are absolutely essential. The problem of modeling the dynamics of deploying and retrieving beam-like structures from a rotating base was addressed. A methodology for discrete modeling, and a computational procedure were developed. These results give us the capability of understanding and predicting the effects on the overall satellite motion of deploying flexible appendages. This is an initial step towards a general capability of treating axially moving three-dimensional beams. The interaction dynamics of the orbiter, its flexible manipulator, and the structures to be assembled/deployed, as a prerequisite in order to simulate incremental in-space structural construction processes are investigated. Preliminary results so obtained indicate that, as the inertia properties of the flexible large space structure under construction change during the space assembly/construction process, the interaction dynamics undergo significant changes in their characteristics, thus revealing the need for a variety of control strategies throughout construction.

Park, K. C.↗

On some recent advances in multidisciplinary analysis of hypersonic vehicles

This paper presents pertinent details of the development and application of an integrated, multidisciplinary finite element analysis tool for the effective modeling and simulation of aerospace vehicles including hypersonic ones. Recent advances in this connection include more efficient CFD solution techniques and also aeroservoelastic stability analysis methods that employ unsteady aerodynamic forces for computation of flutter and divergence speeds. An accelerated Euler solution procedure based on the Aitken acceleration technique has recently been implemented that effects considerable improvement in solution efficiency for the same level of solution accuracy, as evidenced by example problems presented in the paper. A summary of finite element numerical formulations for the various individual disciplines as well as for the unified aeroservoelastic analysis and some relevant numerical examples are also presented in the paper.

Gupta, K. K.↗

User's guide for ENSAERO: A multidisciplinary program for fluid/structural/control interaction studies of aircraft (release 1)

Strong interactions can occur between the flow about an aerospace vehicle and its structural components resulting in several important aeroelastic phenomena. These aeroelastic phenomena can significantly influence the performance of the vehicle. At present, closed-form solutions are available for aeroelastic computations when flows are in either the linear subsonic or supersonic range. However, for aeroelasticity involving complex nonlinear flows with shock waves, vortices, flow separations, and aerodynamic heating, computational methods are still under development. These complex aeroelastic interactions can be dangerous and limit the performance of aircraft. Examples of these detrimental effects are aircraft with highly swept wings experiencing vortex-induced aeroelastic oscillations, transonic regime at which the flutter speed is low, aerothermoelastic loads that play a critical role in the design of high-speed vehicles, and flow separations that often lead to buffeting with undesirable structural oscillations. The simulation of these complex aeroelastic phenomena requires an integrated analysis of fluids and structures. This report presents a summary of the development, applications, and procedures to use the multidisciplinary computer code ENSAERO. This code is based on the Euler/Navier-Stokes flow equations and modal/finite-element structural equations.

Guruswamy, Guru P.↗

Theoretical Evaluation of the Pressures, Forces, and Moments at Hypersonic Speeds Acting on Arbitrary Bodies of Revolution Undergoing Separate and Combined Angle-of-attack and Pitching Motions

Equations based on Newtonian impact theory have been derived and a computational procedure developed with the aid of several design-type charts which enable the determination of the aerodynamic forces and moments acting on arbitrary bodies of revolution undergoing either separate or combined angle-of-attack and pitching motions. Bodies with axially increasing and decreasing cross-sectional area distributions are considered; nose shapes may be sharp, blunt, or flat faced. The analysis considers variations in angle of attack from -90 degrees to 90 degrees and allows for both positive and negative pitching rates of arbitrary magnitude. The results are also directly applicable to bodies in either separate or combined sideslip and yawing maneuvers.

Margolis, Kenneth↗

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗

Computer aided design and analysis of gear tooth geometry

A simulation method for gear hobbing and shaping of straight and spiral bevel gears is presented. The method is based upon an enveloping theory for gear tooth profile generation. The procedure is applicable in the computer aided design of standard and nonstandard tooth forms. An inverse procedure for finding a conjugate gear tooth profile is presented for arbitrary cutter geometry. The kinematic relations for the tooth surfaces of straight and spiral bevel gears are proposed. The tooth surface equations for these gears are formulated in a manner suitable for their automated numerical development and solution.

Chang, S. H.↗