Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational 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 577 records · Page 32

Aeroacoustic sensitivity analysis and optimal aeroacoustic design of turbomachinery blades

During the first year of the project, we have developed a theoretical analysis - and wrote a computer code based on this analysis - to compute the sensitivity of unsteady aerodynamic loads acting on airfoils in cascades due to small changes in airfoil geometry. The steady and unsteady flow though a cascade of airfoils is computed using the full potential equation. Once the nominal solutions have been computed, one computes the sensitivity. The analysis takes advantage of the fact that LU decomposition is used to compute the nominal steady and unsteady flow fields. If the LU factors are saved, then the computer time required to compute the sensitivity of both the steady and unsteady flows to changes in airfoil geometry is quite small. The results to date are quite encouraging, and may be summarized as follows: (1) The sensitivity procedure has been validated by comparing the results obtained by 'finite difference' techniques, that is, computing the flow using the nominal flow solver for two slightly different airfoils and differencing the results. The 'analytic' solution computed using the method developed under this grant and the finite difference results are found to be in almost perfect agreement. (2) The present sensitivity analysis is computationally much more efficient than finite difference techniques. We found that using a 129 by 33 node computational grid, the present sensitivity analysis can compute the steady flow sensitivity about ten times more efficiently that the finite difference approach. For the unsteady flow problem, the present sensitivity analysis is about two and one-half times as fast as the finite difference approach. We expect that the relative efficiencies will be even larger for the finer grids which will be used to compute high frequency aeroacoustic solutions. Computational results show that the sensitivity analysis is valid for small to moderate sized design perturbations. (3) We found that the sensitivity analysis provided important insight into how airfoils should be modified to improve their aeroelastic stability. Using this insight, we redesigned an aeroelastically unstable cascade to produce a stable cascade.

Hall, Kenneth C.↗

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

Computation of turbulent near wake for asymmetric airfoils

A numerical procedure for studying the turbulent near wake of two dimensional airfoil sections is presented. The Reynolds Navier-Stokes equations were written for flow about bodies of arbitrary geometry and solved on an arbitrary nonuniform curvilinear computational mesh. Eddy viscosity and Reynolds stress turbulence transport models are considered. Specific examples are shown for airfoil section by using an algebraic viscosity model with streamwise relaxation and the interactive Reynolds stress model.

Deiwert, G. S.↗

QUICK-GEOMETRY, a rapid response method for mathematically modeling configuration geometry

The philosophy, development, and various applications of the QUICK-GEOMETRY system were outlined. This system provides a practical method for developing the geometry models that are essential to the operation of computer-based design and manufacturing systems. Of particular interest are the various methods for modeling surface geometry that are used by aerodynamic analysis codes.

Vachris, A. F., Jr.↗

Comparison of Aero-Propulsive Performance Predictions for Distributed Propulsion Configurations

NASA's X-57 "Maxwell" flight demonstrator incorporates distributed electric propulsion technologies in a design that will achieve a significant reduction in energy used in cruise flight. A substantial portion of these energy savings come from beneficial aerodynamic-propulsion interaction. Previous research has shown the benefits of particular instantiations of distributed propulsion, such as the use of wingtip-mounted cruise propellers and leading edge high-lift propellers. However, these benefits have not been reduced to a generalized design or analysis approach suitable for large-scale design exploration. This paper discusses the rapid, "design-order" toolchains developed to investigate the large, complex tradespace of candidate geometries for the X-57. Due to the lack of an appropriate, rigorous set of validation data, the results of these tools were compared to three different computational flow solvers for selected wing and propulsion geometries. The comparisons were conducted using a common input geometry, but otherwise different input grids and, when appropriate, different flow assumptions to bound the comparisons. The results of these studies showed that the X-57 distributed propulsion wing should be able to meet the as-designed performance in cruise flight, while also meeting or exceeding targets for high-lift generation in low-speed flight.

Borer, Nicholas K.↗

The computation of mesh configurations for three-dimensional flow analysis

A geometry modeling method was developed to define three-dimensional mesh configurations, which were expressed by the set of surface intersections with a coordinate mesh. These configurations define bodies for transonic potential-flow computation in simple physical coordinates. Various configurations were prepared for 3-D potential flow analysis, including many bare inlet models in cylindrical coordinates. Trimmed surface models were used to prepare nacelle installation configurations. A variety of ducts and mixers were represented. Analytically defined mesh/surface intersections were compared with similar data obtained from the corresponding surface model. Parameters for an analytic surface definition from mesh-surface intersections obtained from a model of the surface were produced. In both studies the modeling accuracy was shown to be compatible with the flow analysis requirements. This technique fulfills the requirements of the associated three-dimensional potential-flow analysis. It preserves the flow analysis method's superior flexibility for adapting to unusual geometries. The geometry modeling technique is supported by a sufficient set of tools for effective use by general users, and sufficient accuracy was demonstrated.

Gibson, S. G.↗

Accurate In Bond Energies

InXn atomization energies are computed for n = 1-3 and X = H, Cl, and CH3. The geometries and frequencies are determined using density functional theory. The atomization energies are computed at the coupled cluster level of theory. The complete basis set limit is obtained by extrapolation. The scalar relativistic effect is computed using the Douglas-Kroll approach. While the heats of formation for InH, InCl and InCl3 are in good agreement with experiment, the current results show that the experimental value for In(CH3)3 must be wrong.

Bauschlicher, Charles W., Jr.↗

Issues and Methods Concerning the Evaluation of Hypersingular and Near-Hypersingular Integrals in BEM Formulations

It is known that higher order modeling of the sources and the geometry in Boundary Element Modeling (BEM) formulations is essential to highly efficient computational electromagnetics. However, in order to achieve the benefits of hIgher order basis and geometry modeling, the singular and near-singular terms arising in BEM formulations must be integrated accurately. In particular, the accurate integration of near-singular terms, which occur when observation points are near but not on source regions of the scattering object, has been considered one of the remaining limitations on the computational efficiency of integral equation methods. The method of singularity subtraction has been used extensively for the evaluation of singular and near-singular terms. Piecewise integration of the source terms in this manner, while manageable for bases of constant and linear orders, becomes unwieldy and prone to error for bases of higher order. Furthermore, we find that the singularity subtraction method is not conducive to object-oriented programming practices, particularly in the context of multiple operators. To extend the capabilities, accuracy, and maintainability of general-purpose codes, the subtraction method is being replaced in favor of the purely numerical quadrature schemes. These schemes employ singularity cancellation methods in which a change of variables is chosen such that the Jacobian of the transformation cancels the singularity. An example of the sin,oularity cancellation approach is the Duffy method, which has two major drawbacks: 1) In the resulting integrand, it produces an angular variation about the singular point that becomes nearly-singular for observation points close to an edge of the parent element, and 2) it appears not to work well when applied to nearly-singular integrals. Recently, the authors have introduced the transformation u(x(prime))= sinh (exp -1) x(prime)/Square root of ((y prime (exp 2))+ z(exp 2) for integrating functions of the form I = Integral of (lambda(r(prime))((e(exp -jkR))/(4 pi R) d D where A (r (prime)) is a vector or scalar basis function and R = Square root of( (x(prime)(exp2) + (y(prime)(exp2) + z(exp 2)) is the distance between source and observation points. This scheme has all of the advantages of the Duffy method while avoiding the disadvantages listed above. In this presentation we will survey similar approaches for handling singular and near-singular terms for kernels with 1/R(exp 2) type behavior, addressing potential pitfalls and offering techniques to efficiently handle special cases.

Fink, P. W.↗

Digital image processing for quantification through full flow field tracing (FFFT) in narrow geometries at low Reynolds numbers

This paper introduces a computer based image processing technique to the field of nonintrusive velocity measurements in fluid mechanics. The method is presented with two alternatives: the first is involving intense interaction between the operator and the computer system; the second is a first generation artificial intelligence based system, where a set of initially imputed rules replaces the operator. The methods are applied to flow in narrow gaps.

Braun, M. J.↗

Comparison of OVERFLOW Computational and Experimental Results for a Blunt Mars Entry Vehicle Concept During Supersonic Retropropulsion

Simulations of supersonic retropropulsion (SRP) flow over a Hypersonic Inflatable Aerodynamic Decelerator (HIAD) blunt-body vehicle were performed using the Overflow Computational Fluid Dynamics (CFD) solver. Simulation conditions and geometry were designed to match specific test runs in the Descent System Study (DSS) testing campaign. The relative accuracy of simulation predictions are assessed by direct comparison to experimental data. Computational predictions of the SRP flowfield and bow shock shape are compared to experimental schlieren images. Comparisons of the model surface pressure environment are presented for unsteady and static discrete tap data as well as time-averaged Pressure-Sensitive Paint (PSP) data.

Supersonic Retropropulsion↗

A Combined Experimental/Computational Investigation of a Rocket Based Combined Cycle Inlet

A rocket based combined cycle inlet geometry has undergone wind tunnel testing and computational analysis with Mach 4 flow at the inlet face. Performance parameters obtained from the wind tunnel tests were the mass capture, the maximum back-pressure, and the self-starting characteristics of the inlet. The CFD analysis supplied a confirmation of the mass capture, the inlet efficiency and the details of the flowfield structure. Physical parameters varied during the test program were cowl geometry, cowl position, body-side bleed magnitude and ingested boundary layer thickness. An optimum configuration was determined for the inlet as a result of this work.

Smart, Michael K.↗

Thermohydrodynamic analysis of cryogenic liquid turbulent flow fluid film bearings, phase 2

The Phase 2 (1994) Annual Progress Report presents two major report sections describing the thermal analysis of tilting- and flexure-pad hybrid bearings, and the unsteady flow and transient response of a point mass rotor supported on fluid film bearings. A literature review on the subject of two-phase flow in fluid film bearings and part of the proposed work for 1995 are also included. The programs delivered at the end of 1994 are named hydroflext and hydrotran. Both codes are fully compatible with the hydrosealt (1993) program. The new programs retain the same calculating options of hydrosealt plus the added bearing geometries, and unsteady flow and transient forced response. Refer to the hydroflext & hydrotran User's Manual and Tutorial for basic information on the analysis and instructions to run the programs. The Examples Handbook contains the test bearing cases along with comparisons with experimental data or published analytical values. The following major tasks were completed in 1994 (Phase 2): (1) extension of the thermohydrodynamic analysis and development of computer program hydroflext to model various bearing geometries, namely, tilting-pad hydrodynamic journal bearings, flexure-pad cylindrical bearings (hydrostatic and hydrodynamic), and cylindrical pad bearings with a simple elastic matrix (ideal foil bearings); (2) improved thermal model including radial heat transfer through the bearing stator; (3) calculation of the unsteady bulk-flow field in fluid film bearings and the transient response of a point mass rotor supported on bearings; and (4) a literature review on the subject of two-phase flows and homogeneous-mixture flows in thin-film geometries.

Sanandres, Luis↗

Simulation of charge exchange plasma propagation near an ion thruster propelled spacecraft

A model describing the charge exchange plasma and its propagation is discussed, along with a computer code based on the model. The geometry of an idealized spacecraft having an ion thruster is outlined, with attention given to the assumptions used in modeling the ion beam. Also presented is the distribution function describing charge exchange production. The barometric equation is used in relating the variation in plasma potential to the variation in plasma density. The numerical methods and approximations employed in the calculations are discussed, and comparisons are made between the computer simulation and experimental data. An analytical solution of a simple configuration is also used in verifying the model.

Robinson, R. S.↗

Analysis of a Channeled Centerbody Supersonic Inlet for F-15B Flight Research

The Propulsion Flight Test Fixture at the NASA Dryden Flight Research Center is a unique test platform available for use on the NASA F-15B airplane, tail number 836, as a modular host for a variety of aerodynamics and propulsion research. The first experiment that is to be flown on the test fixture is the Channeled Centerbody Inlet Experiment. The objectives of this project at Dryden are twofold: 1) flight evaluation of an innovative new approach to variable geometry for high-speed inlets, and 2) flight validation of channeled inlet performance prediction by complex computational fluid dynamics codes. The inlet itself is a fixed-geometry version of a mixed-compression, variable-geometry, supersonic in- let developed by TechLand Research, Inc. (North Olmsted, Ohio) to improve the efficiency of supersonic flight at off-nominal conditions. The concept utilizes variable channels in the centerbody section to vary the mass flow of the inlet, enabling efficient operation at a range of flight conditions. This study is particularly concerned with the starting characteristics of the inlet. Computational fluid dynamics studies were shown to align well with analytical predictions, showing the inlet to remain unstarted as designed at the primary test point of Mach 1.5 at an equivalent pressure altitude of 29,500 ft local conditions. Mass-flow-related concerns such as the inlet start problem, as well as inlet efficiency in terms of total pressure loss, are assessed using the flight test geometry.

Ratnayake, Nalin A.↗

A computer-assisted process for supersonic aircraft conceptual design

Design methodology was developed and existing major computer codes were selected to carry out the conceptual design of supersonic aircraft. A computer-assisted design process resulted from linking the codes together in a logical manner to implement the design methodology. The process does not perform the conceptual design of a supersonic aircraft but it does provide the designer with increased flexibility, especially in geometry generation and manipulation. Use of the computer-assisted process for the conceptual design of an advanced technology Mach 3.5 interceptor showed the principal benefit of the process to be the ability to use a computerized geometry generator and then directly convert the geometry between formats used in the geometry code and the aerodynamics codes. Results from the interceptor study showed that a Mach 3.5 standoff interceptor with a 1000 nautical-mile mission radius and a payload of eight Phoenix missiles appears to be feasible with the advanced technologies considered. A sensitivity study showed that technologies affecting the empty weight and propulsion system would be critical in the final configuration characteristics with aerodynamics having a lesser effect for small perturbations around the baseline.

Johnson, V. S.↗

Prediction of sonic boom at a focus

The behavior of sonic boom at a focus has been reviewed for the purpose of extending present sonic boom computational methods to include focal zones. The geometry of a focal zone - whether a smooth caustic, a cusped caustic, or a perfect focus to a point - determines the character of focused signatures. The seeming contradiction of various experimental data can be resolved by noting these differences. A ray acoustic analysis has been developed for quantitative determination of caustic geometry. The only reliable theory presently available for signatures at a focus is for a smooth caustic. There has been some controversy between theoretical and experimental values of a constant in the scaling law for this case. It has been found that this discrepancy can be resolved by accounting for the finite thickness of real sonic boom shock waves. These findings have been incorporated into an existing sonic boom computer program.

Plotkin, K. J.↗

Computation of three-dimensional viscous flows using a space-marching method

A space-marching method, developed to compute three-dimensional flows for internal geometries, has been utilized to predict viscous flows through a curved duct and over a swept wing. The Navier-Stokes equations have been posed as an initial value problem by neglecting the streamwise viscous diffusion terms and by treating the pressure gradient as a known source term. The resulting equations have been solved by a non-iterative (single pass) algorithm at each streamwise step. The results are compared with earlier computations (based on iterative methods) and the experimental data. The agreement between the present predictions, the experimental data, and the earlier predictions is good for the cases computed. The computation time is only a fraction of the iterative methods.

Murthy, K. N. S.↗

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗