Search NASASearch

SEARCH · Search NASA

Results for “Complex Geometry”

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 91 records · Page 5

Annual Research Briefs

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Spinks, Debra

A Near-Wall Reynolds-Stress Closure Without Wall Normals

Turbulent wall-bounded complex flows are commonly encountered in engineering practice and are of considerable interest in a variety of industrial applications. The presence of a wall significantly affects turbulence characteristics. In addition to the wall effects, turbulent wall-bounded flows become more complicated by the presence of additional body forces (e.g. centrifugal force and Coriolis force) and complex geometry. Most near-wall Reynolds stress models are developed from a high-Reynolds-number model which assumes turbulence is homogenous (or quasi-homogenous). Near-wall modifications are proposed to include wall effects in near-wall regions. In this process, wall normals are introduced. Good predictions could be obtained by Reynolds stress models with wall normals. However, ambiguity arises when the models are applied in flows with multiple walls. Many models have been proposed to model turbulent flows. Among them, Reynolds stress models, in which turbulent stresses are obtained by solving the Reynolds stress transport equations, have been proved to be the most successful ones. To apply the Reynolds stress models to wall-bounded flows, near-wall corrections accounting for the wall effects are needed, and the resulting models are called near-wall Reynolds stress models. In most of the existing near-wall models, the near-wall corrections invoke wall normals. These wall-dependent near-wall models are difficult to implement for turbulent flows with complex geometry and may give inaccurate predictions due to the ambiguity of wall normals at corners connecting multiple walls. The objective of this study is to develop a more general and flexible near-wall Reynolds stress model without using any wall-dependent variable for wall-bounded turbulent flows. With the aid of near-wall asymptotic analysis and results of direct numerical simulation, a new near-wall Reynolds stress model (NNWRS) is formulated based on Speziale et al.'s high-Reynolds-stress model with wall-independent near-wall corrections. Moreover, only one damping function is used for flows with a wide range of Reynolds numbers to ensure that the near-wall modifications diminish away from the walls.

Yuan, S. P.

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.

Additive Manufacture of Novel Tungsten Wicking and Heat Transfer Structures

Additive manufacturing (AM) is utilized for fabricating complex geometries from refractory metals, such as tungsten and molybdenum. In this study, plates with ultra-fine and fine lattice structures were manufactured from tungsten to demonstrate printability and the component fabrication advantages of AM. Previous work was completed on individual tungsten lattice structures; this work demonstrates the advantages of printing these structures on tubing and first-wall materials. Initial work using stainless steel demonstrated increased heat transfer properties promoted by these lattice structures. Applications for these structures include wicking structures, heat pipes, and first-wall materials in fusion devices. A major advantage of these designs is the capability to produce a monolithic component with tailored heat transfer and wicking properties. This research helps promote the ability to use liquid lithium in fusion devices since these lattice structures manufactured from tungsten could act as a high temperature first-wall or divertor material and as a wicking structure for liquid lithium. Different lattice structures with varying unit cell sizes and strut thicknesses were printed and the effects of these structures on heat transfer properties were measured. Abrasion resistance of these structures was also investigated. This presentation will discuss the results of these experiments and their implications related to heat transfer applications. This work will help inform the development and implementation of complex geometries made from tungsten for future heat pipe and fusion first-wall or divertor component designs.

Carly J. Romnes

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar

Application of a full potential method for computation of three-dimensional supersonic flows

A nonlinear aerodynamic analysis technique based on the full potential equation in conservative form has been modified to permit treatment of supersonic flows with embedded subsonic regions (typically near the fuselage-canopy juncture and the wing leading edge). Solution procedures for the equations do not require any specific form of geometry or physical grid system. This results in the capability to analyze easily very complex geometries provided the posed problem lies within the isentropic restrictions of the full potential theory. Characteristic signal propagation theory is used to monitor the type dependent flow and a conservative switching scheme is employed to transition from the supersonic marching algorithm to a subsonic relaxation procedure and vice versa. An implicit approximate factorization scheme is used to solve the finite-difference equations. These modifications now permit analysis of fully three-dimensional flowfields including the interference effects due to lifting surface wakes. Improved grid generation capability allows analysis of complete complex aircraft geometries (fuselage, wing, tail, wing wake, and tail wake). Results are presented showing very good correlations with experimental surface pressure data and aerodynamic force data at both design and off-design operating points. Configurations examined include several waverider concepts, an arrow wing-body with wake, an advanced tactical fighter concept, and a fighter forebody-canard configuration.

Jones, K. M.

An interactive graphics program for manipulation and display of panel method geometry

Modern aerodynamic panel methods that handle large, complex geometries have made evident the need to interactively manipulate, modify, and view such configurations. With this purpose in mind, the GEOM program was developed. It is a menu driven, interactive program that uses the Tektronix PLOT 10 graphics software to display geometry configurations which are characterized by an abutting set of networks. These networks are composed of quadrilateral panels which are described by the coordinates of their corners. GEOM is divided into fourteen executive controlled functions. These functions are used to build configurations, scale and rotate networks, transpose networks defining M and N lines, graphically display selected networks, join and split networks, create wake networks, produce symmetric images of networks, repanel and rename networks, display configuration cross sections, and output network geometry in two formats. A data base management system is used to facilitate data transfers in this program. A sample session illustrating various capabilities of the code is included as a guide to program operation.

Hall, J. F.

Development of a Robust and Efficient Parallel Solver for Unsteady Turbomachinery Flows

The traditional design and analysis practice for advanced propulsion systems relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. In the design of advanced propulsion systems, CFD plays a major role in defining the required performance over the entire flight regime, as well as in testing the sensitivity of the design to the different modes of operation. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment. The turbomachinery simulation capability presented here is being developed in a computational tool called Loci-STREAM [1]. It integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci [2] which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective is to be able to routinely simulate problems involving complex geometries requiring large unstructured grids and complex multidisciplinary physics. An immediate application of interest is simulation of unsteady flows in rocket turbopumps, particularly in cryogenic liquid rocket engines. The key components of the overall methodology presented in this paper are the following: (a) high fidelity unsteady simulation capability based on Detached Eddy Simulation (DES) in conjunction with second-order temporal discretization, (b) compliance with Geometric Conservation Law (GCL) in order to maintain conservative property on moving meshes for second-order time-stepping scheme, (c) a novel cloud-of-points interpolation method (based on a fast parallel kd-tree search algorithm) for interfaces between turbomachinery components in relative motion which is demonstrated to be highly scalable, and (d) demonstrated accuracy and parallel scalability on large grids (approx 250 million cells) in full turbomachinery geometries.

West, Jeff

Development of a gas pressure bonded four-pole alternator rotor

Methods were developed for fabrication of a solid four pole alternator rotor by hot isostatic pressure welding. The rotor blanks welded in this program had complex geometrical mating interfaces and were of considerable bulk, being approximately 3-1/2 inches (0.089 meters) in diameter and 14 inches (0.356 meters) long. Magnetic end pieces were machined from AlSl 4340 steel, while the non-magnetic central section was of Inconel 718. Excellent welds were produced which were shown to be responsive to post weld heat treatments which substantially improved joint strength. Prior to welding the rotors, test specimens of complex geometry were welded to demonstrate that complex surfaces with intentional mechanical misfit could be readily joined using HIP welding. This preliminary work demonstrated not only that interface compliance is achieved during welding but that welding pressure is developed in these thick sections sufficient to produce sound joints. Integral weld-heat treatment cycles were developed that permitted the attainment of magnetic properties while minimizing residual stress associated with the allotropic transformation of 4340 steel.

Lessmann, G. G.

Atmospheric and Science Complexity Effects on Surface Bidirectional Reflectance

Among the tools used in passive remote sensing of Earth resources in the visible and near-infrared spectral regions are measurements of spectral signature and bidirectional reflectance functions (BDRFs). Determination of surface properties using these observables is complicated by a number of factors, including: (1) mixing of surface components, such as soil and vegetation, (2) multiple reflections of radiation due to complex geometry, such as in crop canopies, and (3) atmospheric effects. In order to bridge the diversity in these different approaches, there is a need for a fundamental physical understanding of the influence of the various effects and a quantiative measure of their relative importance. In particular, we consider scene complexity effects using the example of reflection by vegetative surfaces. The interaction of sunlight with a crop canopy and interpretation of the spectral and angular dependence of the emergent radiation is basically a multidimensional radiative transfer problem. The complex canopy geometry, underlying soil cover, and presence of diffuse as well as collimated illumination will modify the reflectance characteristics of the canopy relative to those of the individual elements.

Diner, D. J.

Numerical Study of Base Pressure Characteristic Curve for a Four-Engine Clustered Nozzle Configuration

Excessive base heating has been a problem for many launch vehicles. For certain designs such as the direct dump of turbine exhaust in the nozzle section and at the nozzle lip of the Space Transportation Systems Engine (STME), the potential burning of the turbine exhaust in the base region has caused tremendous concern. Two conventional approaches have been considered for predicting the base environment: (1) empirical approach, and (2) experimental approach. The empirical approach uses a combination of data correlations and semi-theoretical calculations. It works best for linear problems, simple physics and geometry. However, it is highly suspicious when complex geometry and flow physics are involved, especially when the subject is out of historical database. The experimental approach is often used to establish database for engineering analysis. However, it is qualitative at best for base flow problems. Other criticisms include the inability to simulate forebody boundary layer correctly, the interference effect from tunnel walls, and the inability to scale all pertinent parameters. Furthermore, there is a contention that the information extrapolated from subscale tests with combustion is not conservative. One potential alternative to the conventional methods is computational fluid dynamics (CFD), which has none of the above restrictions and is becoming more feasible due to maturing algorithms and advancing computer technology. It provides more details of the flowfield and is only limited by computer resources. However, it has its share of criticisms as a predictive tool for base environment. One major concern is that CFD has not been extensively tested for base flow problems. It is therefore imperative that CFD be assessed and benchmarked satisfactorily for base flows. In this study, the turbulent base flowfield of a experimental investigation for a four-engine clustered nozzle is numerically benchmarked using a pressure based CFD method. Since the cold air was the medium, accurate prediction of the base pressure distributions at high altitudes is the primary goal. Other factors which may influence the numerical results such as the effects of grid density, turbulence model, differencing scheme, and boundary conditions are also being addressed. Preliminary results of the computed base pressure agreed reasonably well with that of the measurement. Basic base flow features such as the reverse jet, wall jet, recompression shock, and static pressure field in plane of impingement have been captured.

Wang, Ten-See

Tensoral Present and Future

The coding of high-performance fluids simulations requires significant knowledge of both numerical and computational details. The magnitude and complexity of low-level details is often enough to discourage many users of turbulence data wishing to study more important, higher-level fluid dynamical questions. These same complexities are often a practical barrier to simulation experts who develop, verify and maintain the codes which generate this data. Future fluids codes, with high resolution and complex geometries, are likely to involve far more coding complexity. The research--the design and implementation of the Tensoral computer language--aims to greatly ease the coding of today's simulation and post-processing codes and at the same time provide a general computational tool for future simulations.

Dresselhaus, Eliot

Fiber Optic Strain Sensor for Planetary Gear Diagnostics

This paper presents a new sensing approach for helicopter damage detection in the planetary stage of a helicopter transmission based on a fiber optic strain sensor array. Complete helicopter transmission damage detection has proven itself a difficult task due to the complex geometry of the planetary reduction stage. The crowded and complex nature of the gearbox interior does not allow for attachment of sensors within the rotating frame. Hence, traditional vibration-based diagnostics are instead based on measurements from externally mounted sensors, typically accelerometers, fixed to the gearbox exterior. However, this type of sensor is susceptible to a number of external disturbances that can corrupt the data, leading to false positives or missed detection of potentially catastrophic faults. Fiber optic strain sensors represent an appealing alternative to the accelerometer. Their small size and multiplexibility allows for potentially greater sensing resolution and accuracy, as well as redundancy, when employed as an array of sensors. The work presented in this paper is focused on the detection of gear damage in the planetary stage of a helicopter transmission using a fiber optic strain sensor band. The sensor band includes an array of 13 strain sensors, and is mounted on the ring gear of a Bell Helicopter OH-58C transmission. Data collected from the sensor array is compared to accelerometer data, and the damage detection results are presented

Kiddy, Jason S.

Quantitative Electron Probe Microanalysis: State of the Art

Quantitative electron-probe microanalysis (EPMA) has improved due to better instrument design and X-ray correction methods. Design improvement of the electron column and X-ray spectrometer has resulted in measurement precision that exceeds analytical accuracy. Wavelength-dispersive spectrometer (WDS) have layered-dispersive diffraction crystals with improved light-element sensitivity. Newer energy-dispersive spectrometers (EDS) have Si-drift detector elements, thin window designs, and digital processing electronics with X-ray throughput approaching that of WDS Systems. Using these systems, digital X-ray mapping coupled with spectrum imaging is a powerful compositional mapping tool. Improvements in analytical accuracy are due to better X-ray correction algorithms, mass absorption coefficient data sets,and analysis method for complex geometries. ZAF algorithms have ban superceded by Phi(pz) algorithms that better model the depth distribution of primary X-ray production. Complex thin film and particle geometries are treated using Phi(pz) algorithms, end results agree well with Monte Carlo simulations. For geological materials, X-ray absorption dominates the corretions end depends on the accuracy of mass absorption coefficient (MAC) data sets. However, few MACs have been experimentally measured, and the use of fitted coefficients continues due to general success of the analytical technique. A polynomial formulation of the Bence-Albec alpha-factor technique, calibrated using Phi(pz) algorithms, is used to critically evaluate accuracy issues and can be also be used for high 2% relative and is limited by measurement precision for ideal cases, but for many elements the analytical accuracy is unproven. The EPMA technique has improved to the point where it is frequently used instead of the petrogaphic microscope for reconnaissance work. Examples of stagnant research areas are: WDS detector design characterization of calibration standards, and the need for more complete treatment of the continuum X-ray fluorescence correction.

Carpernter, P. K.

Conservative zonal schemes for patched grids in 2 and 3 dimensions

The computation of flow over complex geometries, such as realistic aircraft configurations, poses difficult grid generation problems for computational aerodynamicists. The creation of a traditional, single-module grid of acceptable quality about an entire configuration may be impossible even with the most sophisticated of grid generation techniques. A zonal approach, wherein the flow field is partitioned into several regions within which grids are independently generated, is a practical alternative for treating complicated geometries. This technique not only alleviates the problems of discretizing a complex region, but also facilitates a block processing approach to computation thereby circumventing computer memory limitations. The use of such a zonal scheme, however, requires the development of an interfacing procedure that ensures a stable, accurate, and conservative calculation for the transfer of information across the zonal borders.

Hessenius, Kristin A.

Progress in Grid Generation: From Chimera to DRAGON Grids

Hybrid grids, composed of structured and unstructured grids, combines the best features of both. The chimera method is a major stepstone toward a hybrid grid from which the present approach is evolved. The chimera grid composes a set of overlapped structured grids which are independently generated and body-fitted, yielding a high quality grid readily accessible for efficient solution schemes. The chimera method has been shown to be efficient to generate a grid about complex geometries and has been demonstrated to deliver accurate aerodynamic prediction of complex flows. While its geometrical flexibility is attractive, interpolation of data in the overlapped regions - which in today's practice in 3D is done in a nonconservative fashion, is not. In the present paper we propose a hybrid grid scheme that maximizes the advantages of the chimera scheme and adapts the strengths of the unstructured grid while at the same time keeps its weaknesses minimal. Like the chimera method, we first divide up the physical domain by a set of structured body-fitted grids which are separately generated and overlaid throughout a complex configuration. To eliminate any pure data manipulation which does not necessarily follow governing equations, we use non-structured grids only to directly replace the region of the arbitrarily overlapped grids. This new adaptation to the chimera thinking is coined the DRAGON grid. The nonstructured grid region sandwiched between the structured grids is limited in size, resulting in only a small increase in memory and computational effort. The DRAGON method has three important advantages: (1) preserving strengths of the chimera grid; (2) eliminating difficulties sometimes encountered in the chimera scheme, such as the orphan points and bad quality of interpolation stencils; and (3) making grid communication in a fully conservative and consistent manner insofar as the governing equations are concerned. To demonstrate its use, the governing equations are discretized using the newly proposed flux scheme, AUSM+, which will be briefly described herein. Numerical tests on representative 2D inviscid flows are given for demonstration. Finally, extension to 3D is underway, only paced by the availability of the 3D unstructured grid generator.

Liou, Meng-Sing

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Britcher, Colin P.