Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,243 records · Page 69

Computations of three-dimensional steady and unsteady viscous incompressible flows

The INS3D family of computational fluid dynamics computer codes is presented. These codes are used to as tools in developing and assessing algorithms for solving the incompressible Navier-Stokes equations for steady-state and unsteady flow problems. This work involves applying the codes to real-world problems involving complex three-dimensional geometries. The algorithms utilized include the method of pseudocompressibility including both central and upwind differencing, several types of artificial dissipation schemes, approximate factorization, and an implicit line-relaxation scheme. These codes have been validated using a wide range of problems including flow over a backward-facing step, driven cavity flow, flow through various types of ducts, and steady and unsteady flow over a circular cylinder. Many diverse flow applications have been solved using these codes including parts of the Space Shuttle Main Engine, problems in naval hydrodynamics, low-speed aerodynamics, and biomedical fluid flows. The presentation details several of these, including the flow through a Space Shuttle Main Engine inducer, vortex shedding behind a circular cylinder, and flow through an artificial heart.

Kwak, Dochan↗

Potential flow theory and operation guide for the panel code PMARC

The theoretical basis for PMARC, a low-order potential-flow panel code for modeling complex three-dimensional geometries, is outlined. Several of the advanced features currently included in the code, such as internal flow modeling, a simple jet model, and a time-stepping wake model, are discussed in some detail. The code is written using adjustable size arrays so that it can be easily redimensioned for the size problem being solved and the computer hardware being used. An overview of the program input is presented, with a detailed description of the input available in the appendices. Finally, PMARC results for a generic wing/body configuration are compared with experimental data to demonstrate the accuracy of the code. The input file for this test case is given in the appendices.

Ashby, Dale L.↗

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗

Combined investigation of Eddy current and ultrasonic techniques for composite materials NDE

Advanced composites are not without trade-offs. Their increased designability brings an increase in the complexity of their internal geometry and, as a result, an increase in the number of failure modes associated with a defect. When two or more isotropic materials are combined in a composite, the isotropic material failure modes may also combine. In a laminate, matrix delamination, cracking and crazing, and voids and porosity, will often combine with fiber breakage, shattering, waviness, and separation to bring about ultimate structural failure. This combining of failure modes can result in defect boundaries of different sizes, corresponding to the failure of each structural component. This paper discusses a dual-technology NDE (Non Destructive Evaluation) (eddy current (EC) and ultrasonics (UT)) study of graphite/epoxy (gr/ep) laminate samples. Eddy current and ultrasonic raster (Cscan) imaging were used together to characterize the effects of mechanical impact damage, high temperature thermal damage and various types of inserts in gr/ep laminate samples of various stacking sequences.

Davis, C. W.↗

A grid generation system for multi-disciplinary design optimization

A general multi-block three-dimensional volume grid generator is presented which is suitable for Multi-Disciplinary Design Optimization. The code is timely, robust, highly automated, and written in ANSI 'C' for platform independence. Algebraic techniques are used to generate and/or modify block face and volume grids to reflect geometric changes resulting from design optimization. Volume grids are generated/modified in a batch environment and controlled via an ASCII user input deck. This allows the code to be incorporated directly into the design loop. Generated volume grids are presented for a High Speed Civil Transport (HSCT) Wing/Body geometry as well a complex HSCT configuration including horizontal and vertical tails, engine nacelles and pylons, and canard surfaces.

Jones, William T.↗

Multiwavelength study of the magnetically active T Tauri star HD 283447

We observed the luminous T Tauri star HD 283447 = V773 Tauri simultaneously at X-ray, ultraviolet, optical photometric and spectroscopic, and radio wavelengths for several hours on UT 1992 September 11. ROSAT, IUE, Very Large Array (VLA) and an intercontinental Very Long Baseline Interferometry (VLBI) network, and three optical observatories participated in the campaign. The star is known for its unusually high and variable nonthermal radio continuum emission. High levels of soft X-ray and Mg II line emission are discovered, with luminosity L(sub x) = 5.5 x 10(exp 30) ergs/s (0.2 - 2 keV) and L(sub Mg II) = 1 x 10(exp 29) ergs/s, respectively. Optically, the spectrum exhibits rather weak characteristics of `classical' T Tauri stars. A faint, broad emission line component, probably due to a collimated wind or infall, is present. During the campaign, the radio luminosity decreased by a factor of 4, while optical/UV lines and X-ray emission remained strong but constant. The large gyrosynchrotron-emitting regions are therefore decoupled from the chromospheric and coronal emission. Five models for the magnetic geometry around the star are discussed; solar-type activity, dipole magnetosphere, star-disk magnetic coupling, disk magnetic fields, and close binary interaction. The data suggest that two magnetic geometries are simultaneously present: complex multipolar fields like those on the Sun, and a large-scale field possibly associated with the circumstellar disk.

Feigelson, Eric D.↗

Microwave Scattering Model for Grass Blade Structures

In this paper, the electromagnetic scattering solution for a grass blade with complex cross-section geometry is considered. It is assumed that the blade cross section is electrically small, but its length is large compared to the incident wavelength. In a recent study it has been shown that the scattering solution for such problems, in the form of a polarizability tensor, can be obtained using the low-frequency approximation in conjunction with the method of moments. In addition, the study shows that the relationship between the polarizability tensor of a dielectric cylinder and its dielectric constant can be approximated by a simple algebraic expression. The results of this study are used to show that this algebraic approximation is valid also for cylinders with cross sections the shape of grass blades, providing that proper values am selected for each of three constants appearing in the expression. These constants are dependent on cylinder shape, and if the relationship between the constants and the three parameters describing a grass blade shape can be determined, an algebraic approximation relating polarizability tensor to blade shape, as well as dielectric constant, can be formed. Since the elements of the polarizability tensor are dependent on only these parameters, this algebraic approximation can replace the cumbersome method of moments model. A conjugate gradient method is then implemented to correctly determine the three constants of the algebraic approximation for each blade shape. A third-order polynomial fit to the data is then determined for each constant, thus providing a complete analytic replacement to the numerical (moment method) scattering model. Comparisons of this approximation to the numerical model show an average error of less than 3%.

Stiles, James M.↗

Numerical Prediction Methods (Reynolds-Averaged Navier-Stokes Simulations of Transonic Separated Flows)

During the past five years, numerous pioneering archival publications have appeared that have presented computer solutions of the mass-weighted, time-averaged Navier-Stokes equations for transonic problems pertinent to the aircraft industry. These solutions have been pathfinders of developments that could evolve into a major new technological capability, namely the computational Navier-Stokes technology, for the aircraft industry. So far these simulations have demonstrated that computational techniques, and computer capabilities have advanced to the point where it is possible to solve forms of the Navier-Stokes equations for transonic research problems. At present there are two major shortcomings of the technology: limited computer speed and memory, and difficulties in turbulence modelling and in computation of complex three-dimensional geometries. These limitations and difficulties are the pacing items of the continuing developments, although the one item that will most likely turn out to be the most crucial to the progress of this technology is turbulence modelling. The objective of this presentation is to discuss the state of the art of this technology and suggest possible future areas of research. We now discuss some of the flow conditions for which the Navier-Stokes equations appear to be required. On an airfoil there are four different types of interaction of a shock wave with a boundary layer: (1) shock-boundary-layer interaction with no separation, (2) shock-induced turbulent separation with immediate reattachment (we refer to this as a shock-induced separation bubble), (3) shock-induced turbulent separation without reattachment, and (4) shock-induced separation bubble with trailing edge separation.

Mehta, Unmeel↗

Evaluation of Perspective and Coplanar Cockpit Displays of Traffic Information to Support Hazard Awareness in Free Flight

We examined the cockpit display representation of traffic, to support the pilot in tactical planning and conflict avoidance. Such displays may support the "free flight" concept, but can also support greater situation awareness in a non-free flight environment. Two perspective views and a coplanar display were contrasted in scenarios in which pilots needed to navigate around conflicting traffic, either in the absence (low workload) or presence (high workload) of a second intruder aircraft. All three formats were configured with predictive aiding vectors that explicitly represented the predicted point of closest pass, and predicted penetration of an alert zone around ownship. Ten pilots were assigned to each of the display conditions, and each flew a series of 60 conflict maneuvers that varied in their workload and the complexity of the conflict geometry. Results indicated a tendency to choose vertical over lateral maneuvers, a tendency which was amplified with the coplanar display. Vertical maneuvers by the intruder produced an added source of workload. Importantly, the coplanar display supported performance in all measures that was equal to or greater than either of the perspective displays (i.e., fewer predicted and actual conflicts, less extreme maneuvers). Previous studies that have indicated perspective superiority have only contrasted these with UNIplanar displays rather than the coplanar display used here.

Merwin, David H.↗

Unsteady Analysis of Separated Aerodynamic Flows Using an Unstructured Multigrid Algorithm

An implicit method for the computation of unsteady flows on unstructured grids is presented. The resulting nonlinear system of equations is solved at each time step using an agglomeration multigrid procedure. The method allows for arbitrarily large time steps and is efficient in terms of computational effort and storage. Validation of the code using a one-equation turbulence model is performed for the well-known case of flow over a cylinder. A Detached Eddy Simulation model is also implemented and its performance compared to the one equation Spalart-Allmaras Reynolds Averaged Navier-Stokes (RANS) turbulence model. Validation cases using DES and RANS include flow over a sphere and flow over a NACA 0012 wing including massive stall regimes. The project was driven by the ultimate goal of computing separated flows of aerodynamic interest, such as massive stall or flows over complex non-streamlined geometries.

Pelaez, Juan↗

A Grid Generation System for Multi-Disciplinary Design Optimization

A general multi-block three-dimensional volume grid generator is presented which is suitable for Multi-Disciplinary Design Optimization. The code is fast, robust, highly automated, and written in ANSI C for platform independence. Algebraic techniques are used to generate and/or modify block face and volume grids to reflect geometric changes resulting from design optimization. Volume grids are generated/modified in a batch environment and controlled via an ASCII user input deck. This allows the code to be incorporated directly into the design loop. Generated volume grids are presented for a High Speed Civil Transport (HSCT) Wing/Body geometry as well a complex HSCT configuration including horizontal and vertical tails, engine nacelles and pylons, and canard surfaces.

Jones, William T.↗

Radiative Transfer Modeling of a Large Pool Fire by Discrete Ordinates, Discrete Transfer, Ray Tracing, Monte Carlo and Moment Methods

Five computational methods for solution of the radiative transfer equation in an absorbing-emitting and non-scattering gray medium were compared on a 2 m JP-8 pool fire. The temperature and absorption coefficient fields were taken from a synthetic fire due to the lack of a complete set of experimental data for fires of this size. These quantities were generated by a code that has been shown to agree well with the limited quantity of relevant data in the literature. Reference solutions to the governing equation were determined using the Monte Carlo method and a ray tracing scheme with high angular resolution. Solutions using the discrete transfer method, the discrete ordinate method (DOM) with both S(sub 4) and LC(sub 11) quadratures, and moment model using the M(sub 1) closure were compared to the reference solutions in both isotropic and anisotropic regions of the computational domain. DOM LC(sub 11) is shown to be the more accurate than the commonly used S(sub 4) quadrature technique, especially in anisotropic regions of the fire domain. This represents the first study where the M(sub 1) method was applied to a combustion problem occurring in a complex three-dimensional geometry. The M(sub 1) results agree well with other solution techniques, which is encouraging for future applications to similar problems since it is computationally the least expensive solution technique. Moreover, M(sub 1) results are comparable to DOM S(sub 4).

Jensen, K. A.↗

Topological Evolution of a Fast Magnetic Breakout CME in 3-Dimensions

W present the extension of the magnetic breakout model for CME initiation to a fully 3-dimensional, spherical geometry. Given the increased complexity of the dynamic magnetic field interactions in 3-dimensions, we first present a summary of the well known axisymmetric breakout scenario in terms of the topological evolution associated with the various phases of the eruptive process. In this context, we discuss the completely analogous topological evolution during the magnetic breakout CME initiation process in the simplest 3-dimensional multipolar system. We show that an extended bipolar active region embedded in an oppositely directed background dipole field has all the necessary topological features required for magnetic breakout, i.e. a fan separatrix surface between the two distinct flux systems, a pair of spine fieldlines, and a true 3-dimensional coronal null point at their intersection. We then present the results of a numerical MHD simulation of this 3-dimensional system where boundary shearing flows introduce free magnetic energy, eventually leading to a fast magnetic breakout CME. The eruptive flare reconnection facilitates the rapid conversion of this stored free magnetic energy into kinetic energy and the associated acceleration causes the erupting field and plasma structure to reach an asymptotic eruption velocity of greater than or approx. equal to 1100 km/s over an approx.15 minute time period. The simulation results are discussed using the topological insight developed to interpret the various phases of the eruption and the complex, dynamic, and interacting magnetic field structures.

Lynch, B. J.↗

Flow Simulation of Supersonic Inlet with Bypass Annular Duct

A relaxed isentropic compression supersonic inlet is a new concept that produces smaller cowl drag than a conventional inlet, but incurs lower total pressure recovery and increased flow distortion in the (radially) outer flowpath. A supersonic inlet comprising a bypass annulus to the relaxed isentropic compression inlet dumps out airflow of low quality through the bypass duct. A reliable computational fluid dynamics solution can provide considerable useful information to ascertain quantitatively relative merits of the concept, and further provide a basis for optimizing the design. For a fast and reliable performance evaluation of the inlet performance, an equivalent axisymmetric model whose area changes accounts for geometric and physical (blockage) effects resulting from the original complex three-dimensional configuration is proposed. In addition, full three-dimensional calculations are conducted for studying flow phenomena and verifying the validity of the equivalent model. The inlet-engine coupling is carried out by embedding numerical propulsion system simulation engine data into the flow solver for interactive boundary conditions at the engine fan face and exhaust plane. It was found that the blockage resulting from complex three-dimensional geometries in the bypass duct causes significant degradation of inlet performance by pushing the terminal normal shock upstream.

Kim, HyoungJin↗

Annular MHD Physics for Turbojet Energy Bypass

The use of annular Hall type MHD generator/accelerator ducts for turbojet energy bypass is evaluated assuming weakly ionized flows obtained from pulsed nanosecond discharges. The equations for a 1-D, axisymmetric MHD generator/accelerator are derived and numerically integrated to determine the generator/accelerator performance characteristics. The concept offers a shockless means of interacting with high speed inlet flows and potentially offers variable inlet geometry performance without the complexity of moving parts simply by varying the generator loading parameter. The cycle analysis conducted iteratively with a spike inlet and turbojet flying at M = 7 at 30 km altitude is estimated to have a positive thrust per unit mass flow of 185 N-s/kg. The turbojet allowable combustor temperature is set at an aggressive 2200 deg K. The annular MHD Hall generator/accelerator is L = 3 m in length with a B(sub r) = 5 Tesla magnetic field and a conductivity of sigma = 5 mho/m for the generator and sigma= 1.0 mho/m for the accelerator. The calculated isentropic efficiency for the generator is eta(sub sg) = 84 percent at an enthalpy extraction ratio, eta(sub Ng) = 0.63. The calculated isentropic efficiency for the accelerator is eta(sub sa) = 81 percent at an enthalpy addition ratio, eta(sub Na) = 0.62. An assessment of the ionization fraction necessary to achieve a conductivity of sigma = 1.0 mho/m is n(sub e)/n = 1.90 X 10(exp -6), and for sigma = 5.0 mho/m is n(sub e)/n = 9.52 X 10(exp -6).

Schneider, Steven J.↗

Asteroid Bennu Temperature Maps for OSIRIS-REx Spacecraft and Instrument Thermal Analyses

A thermophysical model has been developed to generate asteroid Bennu surface temperature maps for OSIRIS-REx spacecraft and instrument thermal design and analyses at the Critical Design Review (CDR). Two-dimensional temperature maps for worst hot and worst cold cases are used in Thermal Desktop to assure adequate thermal design margins. To minimize the complexity of the Bennu geometry in Thermal Desktop, it is modeled as a sphere instead of the radar shape. The post-CDR updated thermal inertia and a modified approach show that the new surface temperature predictions are more benign. Therefore the CDR Bennu surface temperature predictions are conservative.

thermal analysis↗

Neutralizer Hollow Cathode Simulations and Comparisons with Ground Test Data

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially ionized gas in these devices. A summary of the comparisons between simulation results and Langmuir probe measurements is provided. The model has also been employed to provide insight into recent ground test observations of the neutralizer cathode in NEXT. It is found that a likely cause of the observed keeper voltage drop is cathode orifice erosion. However, due to the small magnitude of this change, is approx. 0.5 V (less than 5% of the beginning-of-life value) over 10 khrs, and in light of the large uncertainties of the cathode material sputtering yield at low ion energies, other causes cannot be excluded. Preliminary simulations to understand transition to plume mode suggest that in the range of 3-5 sccm the existing 2-D model reproduces fairly well the rise of the keeper voltage in the NEXT neutralizer as observed in the laboratory. At lower flow rates the simulation produces oscillations in the keeper current and voltage that require prohibitively small time-steps to resolve with the existing algorithms.

electric propulsion↗

Neutralizer Hollow Cathode Simulations and Validation with Experiments

The fidelity of electric propulsion physics-based models depends largely on the validity of their predictions over a range of operating conditions and geometries. In general, increased complexity of the physics requires more extensive comparisons with laboratory data to identify the region(s) that lie outside the validity of the model assumptions and to quantify the uncertainties within its range of application. This paper presents numerical simulations of neutralizer hollow cathodes at various operating conditions and orifice sizes. The simulations were performed using a two-dimensional axisymmetric model that solves numerically a relatively extensive system of conservation laws for the partially-ionized gas in these devices. The results for the plasma are compared directly with Langmuir probe measurements. The computed keeper voltages are also compared with the observed values. Wherever model inputs and/or specific physics of the cathode discharge are uncertain, additional sensitivity calculations have been performed to quantify the uncertainties.

partially ionized gases↗