Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “domain structure”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Derivation of Effective Properties Based on Porous Scale Simulations Using Filtering Techniques

This study presents a method for derivation of effective properties at the interface and in-depth of porous materials. The method defines a Representative Elementary Volume (REV) and applies filtering techniques to computer effective properties such as porosity and flow quantities, such as velocity and pressure. The script, developed to process the data was tested on the VTK type files that contain the mesh information and the flow solution. The method allows to choose between two types of filters, such as cellular and top-hat and define the size of the REV and number of samples along the domain. Extraction of the REV from the domain is performed to exact boundaries requested for the user. This is done using a triangulation technique and cutting through the cells to comply to the requested boundaries of the volume. The method can be applied to both structured and unstructured meshes. Filtering the material porosity and flow quantities involves integration of the numerical data. The algorithm provides three integration methods, such as Riemann sum, Monte Carlo and Quadrature rule to perform the integration. The Monte-Carlo technique permits the use of either uniform or linearly spaced distribution of points. The Quadrature rule is currently applicable to tetrahedral element types. The Monte Carlo and Quadrature rule methods require interpolation of the flow quantities at the sample points. For interpolation, two methods were tested and are readily available, Gaussian interpolation and re-sampling. It has been shown that re-sampling method has better consistency and acceptable accuracy in interpolation of the data. The algorithm was written in Python language and uses a number of modules. The major module besides numpy is PyVista. It is used to process the computational domain, clip the REV and interpolate the data. Quadrature rule integration was performed using a quadpy module. ParaView software was used externally to convert the flow solution to the VTK (or more specifically VTU) format. Integration of ParaView in the same environment with PyVista encountered problems and could not be implemented in this work. The developed algorithm is expected to be applicable to unstructured meshes and more complex porous structures as soon as the data can be passed in VTK type format. With the report is provided Python script for filtering the solution and a Matlab script for simple generation and processing of 2-D and 3-D porous channel geometries. The two scripts don't communicate.

Alexsander Zibitsker↗
NASA NTRS1991

Kinematic domain decomposition to simulate flows past moving objects

A new method is developed to solve the unsteady Navier-Stokes equations on a composite grid which consists of subdomain grids moving with respect to each other. These subdomains are structured grids with different topologies. However, their interfaces are unstructured in the sense that they may or may not arbitrarily overlap with each other. This method captures the boundary-motion-induced flow component. The method is demonstrated for a transonic flow past an airfoil which experiences a combined motion of pitching and plunging.

Baysal, Oktay↗
NASA NTRS1992

Sub-domain decomposition methods and computational controls for multibody dynamical systems

This paper presents a concurrent methodology to simulate the dynamics of flexible multibody systems with a large number of degrees of freedom. A general class of open-loop structures is treated and a redundant coordinate formulation is adopted. A range space method is used in which the constraint forces are calculated using a preconditioned conjugate gradient method. By using a preconditioner motivated by the regular ordering of the directed graph of the structures, it is shown that the method is order N in the total number of coordinates of the system. The overall formulation has the advantage that it permits fine parallelization and does not rely on system topology to induce concurrency. It can be efficiently implemented on the present generation of parallel computers with a large number of processors. Validation of the method is presented via numerical simulations of space structures incorporating large number of flexible degrees of freedom.

Menon, R. G.↗
NASA NTRS1991

Kinematic domain decomposition for boundary-motion-induced flow simulations

A method is developed to solve the unsteady Navier-Stokes equations on a composite grid, which consists of subdomain grids moving with respect to each other. These subdomains are structured grids with different topologies. This method eliminates assuming the moving components to be instantaneously stationary, where deciding on the particular frozen instants is difficult and affects the solution adversely. Moreover, this method captures the boundary-motion-induced flow component. The method is demonstrated through a transonic flow past an airfoil, which experiences a combined motion of pitching and plunging. An O-grid around the airfoil is overlapped on a fine Cartesian grid, which is zonally embedded in a coarse Cartesian grid. The coarse grid is stationary but the other two grids are plunging. Only the O-grid is also sinusoidally pitching. The results are compared successfully with the experimental data.

Baysal, Oktay↗
NASA NTRS1998

Fuzzy Versions of Epistemic and Deontic Logic

Epistemic and deontic logics are modal logics, respectively, of knowledge and of the normative concepts of obligation, permission, and prohibition. Epistemic logic is useful in formalizing systems of communicating processes and knowledge and belief in AI (Artificial Intelligence). Deontic logic is useful in computer science wherever we must distinguish between actual and ideal behavior, as in fault tolerance and database integrity constraints. We here discuss fuzzy versions of these logics. In the crisp versions, various axioms correspond to various properties of the structures used in defining the semantics of the logics. Thus, any axiomatic theory will be characterized not only by its axioms but also by the set of properties holding of the corresponding semantic structures. Fuzzy logic does not proceed with axiomatic systems, but fuzzy versions of the semantic properties exist and can be shown to correspond to some of the axioms for the crisp systems in special ways that support dependency networks among assertions in a modal domain. This in turn allows one to implement truth maintenance systems. For the technical development of epistemic logic, and for that of deontic logic. To our knowledge, we are the first to address fuzzy epistemic and fuzzy deontic logic explicitly and to consider the different systems and semantic properties available. We give the syntax and semantics of epistemic logic and discuss the correspondence between axioms of epistemic logic and properties of semantic structures. The same topics are covered for deontic logic. Fuzzy epistemic and fuzzy deontic logic discusses the relationship between axioms and semantic properties for these logics. Our results can be exploited in truth maintenance systems.

Ramasamy S Gounder↗
NASA NTRS2007

The Spatiotemporal Structure of 20th Century Climate Variations in Observations and Reanalyses: Long-term Trend - Part 1

The dominant interannual El Nino-Southern Oscillation phenomenon (ENSO) and the short length of climate observation records make it difficult to study long-term climate variations in the spatiotemporal domain. Based on the fact that the ENS0 signal spreads to remote regions and induces delayed climate variation through atmospheric teleconnections, we develop an ENSO-removal method through which the ENS0 signal can be approximately removed at the grid box level from the spatiotemporal field of a climate parameter. After this signal is removed, long-term climate variations, namely, the global warming trend (GW) and the Pacific pan-decadal variability (PDV), are isolated at middle and low latitudes in the climate parameter fields from observed and reanalyses datasets. Except for known GW characteristics, the warming that occurs in the Pacific basin (approximately 0.4K in the 2oth century) is much weaker than in surrounding regions and the other two ocean basins (approximately 0.8K). The modest warming in the Pacific basin is likely due to its dynamic nature on the interannual and decadal time scales and/or the leakage of upper ocean water through the Indonesian Throughflow. Based on NCEP/NCAR and ERA-40 reanalyses, a comprehensive atmospheric structure associated with GW is given. Significant discrepancies exist between the two datasets, especially in the tightly coupled dynamic and water vapor fields. The dynamic field based on NCEP/NCAR reanalysis, which shows a change in the Walker Circulation, is consistent with the GW change in the surface temperature field. However, intensification in the Hadley Circulation is associated with GW trend in the ERA-40 reanalysis.

Chen, Junye↗
NASA NTRS2004

Realtime Knowledge Management (RKM): From an International Space Station (ISS) Point of View

We are developing automated methods to provide realtime access to spacecraft domain knowledge relevant a spacecraft's current operational state. The method is based upon analyzing state-transition signatures in the telemetry stream. A key insight is that documentation relevant to a specific failure mode or operational state is related to the structure and function of spacecraft systems. This means that diagnostic dependency and state models can provide a roadmap for effective documentation navigation and presentation. Diagnostic models consume the telemetry and derive a high-level state description of the spacecraft. Each potential spacecraft state description is matched against the predictions of models that were developed from information found in the pages and sections in the relevant International Space Station (ISS) documentation and reference materials. By annotating each model fragment with the domain knowledge sources from which it was derived we can develop a system that automatically selects those documents representing the domain knowledge encapsulated by the models that compute the current spacecraft state. In this manner, when the spacecraft state changes, the relevant documentation context and presentation will also change.

Robinson, Peter I.↗
NASA NTRS2003

Fast-Running Aeroelastic Code Based on Unsteady Linearized Aerodynamic Solver Developed

The NASA Glenn Research Center has been developing aeroelastic analyses for turbomachines for use by NASA and industry. An aeroelastic analysis consists of a structural dynamic model, an unsteady aerodynamic model, and a procedure to couple the two models. The structural models are well developed. Hence, most of the development for the aeroelastic analysis of turbomachines has involved adapting and using unsteady aerodynamic models. Two methods are used in developing unsteady aerodynamic analysis procedures for the flutter and forced response of turbomachines: (1) the time domain method and (2) the frequency domain method. Codes based on time domain methods require considerable computational time and, hence, cannot be used during the design process. Frequency domain methods eliminate the time dependence by assuming harmonic motion and, hence, require less computational time. Early frequency domain analyses methods neglected the important physics of steady loading on the analyses for simplicity. A fast-running unsteady aerodynamic code, LINFLUX, which includes steady loading and is based on the frequency domain method, has been modified for flutter and response calculations. LINFLUX, solves unsteady linearized Euler equations for calculating the unsteady aerodynamic forces on the blades, starting from a steady nonlinear aerodynamic solution. First, we obtained a steady aerodynamic solution for a given flow condition using the nonlinear unsteady aerodynamic code TURBO. A blade vibration analysis was done to determine the frequencies and mode shapes of the vibrating blades, and an interface code was used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor was used to interpolate the mode shapes from the structural dynamic mesh onto the computational dynamics mesh. Then, we used LINFLUX to calculate the unsteady aerodynamic forces for a given mode, frequency, and phase angle. A postprocessor read these unsteady pressures and calculated the generalized aerodynamic forces, eigenvalues, and response amplitudes. The eigenvalues determine the flutter frequency and damping. As a test case, the flutter of a helical fan was calculated with LINFLUX and compared with calculations from TURBO-AE, a nonlinear time domain code, and from ASTROP2, a code based on linear unsteady aerodynamics.

Reddy, T. S. R.↗
NASA NTRS1981

Application of a flight test and data analysis technique to flutter of a drone aircraft

Modal identification results presented were obtained from recent flight flutter tests of a drone vehicle with a research wing (DAST ARW-1 for Drones for Aerodynamic and Structural Testing, Aeroelastic Research Wing-1). This vehicle is equipped with an active flutter suppression system (FSS). Frequency and damping of several modes are determined by a time domain modal analysis of the impulse response function obtained by Fourier transformations of data from fast swept sine wave excitation by the FSS control surface on the wing. Flutter points are determined for two different altitudes with the FSS off. Data are given for near the flutter boundary with the FSS on.

Bennett, R. M.↗
NASA NTRS1990

Evolution of an Intelligent Information Fusion System

Consideration is given to the hardware and software needed to manage the enormous amount and complexity of data that the next generation of space-borne sensors will provide. An anthology is presented illustrating the evolution of artificial intelligence, science data processing, and management from the 1960s to the near future. Problems and limitations of technologies, data structures, data standards, and conceptual thinking are addressed. The development of an end-to-end Intelligent Information Fusion System that embodies knowledge of the user's domain-specific goals is proposed.

Campbell, William J.↗
NASA NTRS1994

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods for the analysis of complete jet engines. We have entitled this program by applying the two dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition, and solution capabilities were successfully tested. We then focused attention on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion that results from these structural displacements. This is treated by a new arbitrary Lagrangian-Eulerian (ALE) technique that models the fluid mesh motion as that of a fictitious mass-spring network. New partitioned analysis procedures to treat this coupled three-component problem are developed. These procedures involved delayed corrections and subcycling. Preliminary results on the stability, accuracy, and MPP computational efficiency are reported.

Felippa, C. A.↗
NASA NTRS1994

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations

A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: a gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment, and other accepted computational results for a series of low and moderate Reynolds number flows.

Coirier, William J.↗
NASA NTRS1994

Simulations of bypass transition for spatially evolving disturbances

The spatial evolution of disturbances in plane Poiseuille flow and zero pressure gradient boundary layer flow is considered. For disturbances governed by the linearized equations, potential for significant transient growth of the amplitude is demonstrated. The maximum amplification occurs for disturbances with zero or near zero frequencies. Spatial numerical simulations of the transition scenario involving a pair of oblique waves has been conducted for both flows. A fully spectral solver using a simple but efficient fringe region technique allowed the flows to be computed with high resolution into the fully turbulent domain. A modal decomposition of the simulation results indicates that non-linear excitation of the transient growth is responsible for the rapid emergence of low-frequency structures. Physically, this corresponds to streaky flow structures, as seen from the results of a numerical amplitude expansion. Thus, this spatial transition scenario has been found to be similar to the corresponding temporal one. In the boundary layer simulations the streaks are seen to break down from what appears to be a secondary instability.

Lundbladh, A.↗
NASA NTRS1995

High-performance parallel analysis of coupled problems for aircraft propulsion

This research program deals with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a by-pass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by an ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled 3-component problem were developed in 1994. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 has been developed. It is planned to use the steady-state global solution provided by ENG10 as input to a localized three-dimensional FSI analysis for engine regions where aeroelastic effects may be important.

Felippa, C. A.↗
NASA NTRS1994

The deep structure of Venusian plateau highlands

Magellan gravity data confirm that several of the large, tectonically deformed, plateau-like highlands on venus are shallowly compensated, most likely by crustal thickness variations. Apparent depths of isostatic compensation, computed in the spatial domain, range from 30 to 50 km for Alpha, Tellus, Ovda, and Thetis Regiones. Using a two-layer model for isostatic compensation, Alpha, Tellus, and Ovda are best represented as nearly completely compensated in crust that is regionally 20-40 km thick around these highlands, with little contribution from deeper mantle sources. In contrast to these three areas, a stronger regional gravity high associated with Thetis requires a significant upper mantle component to compensation. This is evident in the spectral admittance as a pronounced deep, long-wavelength anomaly. In the two-layer isostatic model, a broad, deeply compensated upland underlies a shallowly compensated central block of Thetis. If this deep component is interpreted as a thermal anomaly, the loci of maximum upwelling agree well with sites of recent extension. The plateau highlands are thus physiographically and isostatically equivalent to terrestrial continents, though probably not compositionally. They also share the record of a long tectonic history. The large regional gravity anomaly of Thetis indicates that active mantle proceses continue even beneath some areas (tessera) thought to be a relic of a former geological regime. The excellent agreement of modeled crustal thicknesses around Alpha, Tellus, and Ovda Regiones suggests that 20-40 km is a representative global value for the plains. Such a crust is thicker than previously estimated and about twice as thick as the expected thickness of crust produced at venusian spreading centers

Grimm, Robert E.↗
NASA NTRS1996

High-Performance Parallel Analysis of Coupled Problems for Aircraft Propulsion

This research program dealt with the application of high-performance computing methods to the numerical simulation of complete jet engines. The program was initiated in January 1993 by applying two-dimensional parallel aeroelastic codes to the interior gas flow problem of a bypass jet engine. The fluid mesh generation, domain decomposition and solution capabilities were successfully tested. Attention was then focused on methodology for the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by these structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field fluid elements. New partitioned analysis procedures to treat this coupled three-component problem were developed during 1994 and 1995. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. For the global steady-state axisymmetric analysis of a complete engine we have decided to use the NASA-sponsored ENG10 program, which uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor tor parallel versions of ENG10 was developed. During 1995 and 1996 we developed the capability tor the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames. Benchmark results were presented at the 1196 Computational Aeroscience meeting.

Felippa, C. A.↗
NASA NTRS1996

Far-Infrared Hydrogen Lasers in the Peculiar Star MWC 349A

Far-infrared hydrogen recombination lines H15(alpha)(169.4 micrometers), H12(alpha)(88.8 micrometers), and H10(alpha)(52.5 micrometers) were detected in the peculiar luminous star MWC 349A from the Kuiper Airborne Observatory. Here it is shown that at least H15(alpha) is strongly amplified, with the probable amplification factor being greater than or about equal to 10(exp 3) and a brightness temperature that is greater than or about equal to 10(exp 7) kelvin. The other two lines also show signs of amplification, although to a lesser degree. Beyond H10(alpha) the amplification apparently vanishes. The newly detected amplified lines fall into the laser wavelength domain. These lasers, as well as the previously detected hydrogen masers may originate in the photoionized circumstellar disk of MWC 349A and constrain the disk's physics and structure.

Strelnitski, Vladimir↗
NASA NTRS2002

Minimizing Cache Misses Using Minimum-Surface Bodies

A number of known techniques for improving cache performance in scientific computations involve the reordering of the iteration space. Some of these reorderings can be considered as coverings of the iteration space with the sets having good surface-to-volume ratio. Use of such sets reduces the number of cache misses in computations of local operators having the iteration space as a domain. First, we derive lower bounds which any algorithm must suffer while computing a local operator on a grid. Then we explore coverings of iteration spaces represented by structured and unstructured grids which allow us to approach these lower bounds. For structured grids we introduce a covering by successive minima tiles of the interference lattice of the grid. We show that the covering has low surface-to-volume ratio and present a computer experiment showing actual reduction of the cache misses achieved by using these tiles. For planar unstructured grids we show existence of a covering which reduces the number of cache misses to the level of structured grids. On the other hand, we present a triangulation of a 3-dimensional cube such that any local operator on the corresponding grid has significantly larger number of cache misses than a similar operator on a structured grid.

Frumkin, Michael↗