Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Domain decomposition”

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 271 records · Page 15

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗

Global energetics analysis using 3-dimensional normal mode decomposition

The Goddard Laboratory for Atmospheric Sciences (GLAS) analysis of the FGGE observations for a 25-day period in January 1979 is examined using the normal mode energetics scheme. The results from the energetics analysis are compared with the data of Tanaka (1985). Kinetic energy spectra of the barotropic mode in the meridional-mode and frequency domains, and the total diabatic process are studied. It is observed that there are significant differences between GLAS and GFDL analyses results in the barotropic energy of gravity modes and the normal energetics scheme is applicable as a diagnostic tool.

Tanaka, Hiroshi↗

Constraint-based integration of planning and scheduling for space-based observatory management

Progress toward the development of effective, practical solutions to space-based observatory scheduling problems within the HSTS scheduling framework is reported. HSTS was developed and originally applied in the context of the Hubble Space Telescope (HST) short-term observation scheduling problem. The work was motivated by the limitations of the current solution and, more generally, by the insufficiency of classical planning and scheduling approaches in this problem context. HSTS has subsequently been used to develop improved heuristic solution techniques in related scheduling domains and is currently being applied to develop a scheduling tool for the upcoming Submillimeter Wave Astronomy Satellite (SWAS) mission. The salient architectural characteristics of HSTS and their relationship to previous scheduling and AI planning research are summarized. Then, some key problem decomposition techniques underlying the integrated planning and scheduling approach to the HST problem are described; research results indicate that these techniques provide leverage in solving space-based observatory scheduling problems. Finally, more recently developed constraint-posting scheduling procedures and the current SWAS application focus are summarized.

Muscettola, Nicola↗

Stability of constant gain systems with vector feedback

The state space, the controllability, and the observability concepts are discussed in connection with the proposed stability analysis which permits drastic dimensional reductions for a vector feedback problem. Any constant gain system's stability can thus be analyzed in the frequency domain with a single Nyquist plot. The analysis considers the total system with all loops closed, a disturbance vector as input, and the feedback vector as output. All constant gain systems are shown to be decomposable into stable subsystems where the degree of the decomposition determines the dimensions. The maximum decomposition results in the state-space approach which is the limit case. The method is demonstrated with the stability analysis of the pogo phenomenon, an oscillatory interaction between the propulsion and the structure of a space vehicle. This problem, with eigenvalues over a hundred, was drastically but rigorously reduced to a stability analysis of a 4x4 matrix.

Vonpragenau, G. L.↗

Refinement and application of acoustic impulse technique to study nozzle transmission characteristics

An improved acoustic impulse technique was developed and was used to study the transmission characteristics of duct/nozzle systems. To accomplish the above objective, various problems associated with the existing spark-discharge impulse technique were first studied. These included (1) the nonlinear behavior of high intensity pulses, (2) the contamination of the signal with flow noise, (3) low signal-to-noise ratio at high exhaust velocities, and (4) the inability to control or shape the signal generated by the source, specially when multiple spark points were used as the source. The first step to resolve these problems was the replacement of the spark-discharge source with electroacoustic driver(s). These included (1) synthesizing on acoustic impulse with acoustic driver(s) to control and shape the output signal, (2) time domain signal averaging to remove flow noise from the contaminated signal, (3) signal editing to remove unwanted portions of the time history, (4) spectral averaging, and (5) numerical smoothing. The acoustic power measurement technique was improved by taking multiple induct measurements and by a modal decomposition process to account for the contribution of higher order modes in the power computation. The improved acoustic impulse technique was then validated by comparing the results derived by an impedance tube method. The mechanism of acoustic power loss, that occurs when sound is transmitted through nozzle terminations, was investigated. Finally, the refined impulse technique was applied to obtain more accurate results for the acoustic transmission characteristics of a conical nozzle and a multi-lobe multi-tube supressor nozzle.

Salikuddin, M.↗

Planning paths through a spatial hierarchy - Eliminating stair-stepping effects

Stair-stepping effects are a result of the loss of spatial continuity resulting from the decomposition of space into a grid. This paper presents a path planning algorithm which eliminates stair-stepping effects induced by the grid-based spatial representation. The algorithm exploits a hierarchical spatial model to efficiently plan paths for a mobile robot operating in dynamic domains. The spatial model and path planning algorithm map to a parallel machine, allowing the system to operate incrementally, thereby accounting for unexpected events in the operating space.

Slack, Marc G.↗

Thin film processing of photorefractive BaTiO3

During the period covered by this report, October 11, 1991 through October 10, 1992, the research has progressed in a number of different areas. The sol-gel technique was initially studied and experimentally evaluated for depositing films of BaTiO3. The difficulties with the precursors and the poor quality of the films deposited lead to the investigation of pulsed laser deposition as an alternative approach. The development of the pulsed laser deposition technique has resulted in continuous improvements to the quality of deposited films of BaTiO3. The initial depositions of BaTiO3 resulted in amorphous films, however, as the pulsed laser deposition technique continued to evolve, films were deposited in the polycrystalline state, then the textured polycrystalline state, and most recently heteroepitaxial films have also been successfully deposited on cubic (100) oriented SrTiO3 substrates. A technique for poling samples at room temperature and in air is also undergoing development with some very preliminary but positive results. The analytical techniques, which include x-ray diffraction, ferroelectric analysis, UV-Vis spectrophotometry, scanning electron microscopy with x-ray compositional analysis, optical and polarized light microscopy, and surface profilometry have been enhanced to allow for more detailed evaluation of the samples. In the area of optical characterization, a pulsed Nd:YAG laser has been incorporated into the experimental configuration. Now data can also be acquired within various temporal domains resulting in more detailed information on the optical response of the samples and on their photorefractive sensitivity. The recent establishment of collaborative efforts with two departments at Johns Hopkins University and the Army Research Lab at Fort Belvoir has also produced preliminary results using the metallo-organic decomposition technique as an alternative method for thin film processing of BaTiO3. RF and DC sputtering is another film deposition approach that should be initiated in the near future. Other techniques for optical characterization, which may even allow for intragrannular (within single grains) investigations, are also being considered.

Schuster, Paul R.↗

Improving NASA's Multiscale Modeling Framework for Tropical Cyclone Climate Study

One of the current challenges in tropical cyclone (TC) research is how to improve our understanding of TC interannual variability and the impact of climate change on TCs. Recent advances in global modeling, visualization, and supercomputing technologies at NASA show potential for such studies. In this article, the authors discuss recent scalability improvement to the multiscale modeling framework (MMF) that makes it feasible to perform long-term TC-resolving simulations. The MMF consists of the finite-volume general circulation model (fvGCM), supplemented by a copy of the Goddard cumulus ensemble model (GCE) at each of the fvGCM grid points, giving 13,104 GCE copies. The original fvGCM implementation has a 1D data decomposition; the revised MMF implementation retains the 1D decomposition for most of the code, but uses a 2D decomposition for the massive copies of GCEs. Because the vast majority of computation time in the MMF is spent computing the GCEs, this approach can achieve excellent speedup without incurring the cost of modifying the entire code. Intelligent process mapping allows differing numbers of processes to be assigned to each domain for load balancing. The revised parallel implementation shows highly promising scalability, obtaining a nearly 80-fold speedup by increasing the number of cores from 30 to 3,335.

tropical cyclone interannual variability↗

Sensitivity to Snapshot Frequency in the POD-Based Reduced-Order Modeling of Flow over a Gaussian Bump

The flow over a Gaussian bump at a Reynolds number of 10^6, based on the bump length and freestream conditions, and a Mach number of 0.2, was simulated using wall-modeled large eddy simulation (WMLES). The incoming turbulent flow exhibits relaminarization upstream of the bump apex, due to the favourable pressure gradient acceleration induced by the bump curvature, as well as separation, reattachment and a detached shear layer downstream of the apex. A Galerkin reduced-order model (ROM) was then constructed using a proper orthogonal decomposition (POD) modal basis. The complex, three-dimensional flow features are captured by spatial POD modes differently based on the frequency of the high-fidelity snapshots originally provided. POD modes constructed from differently sampled snapshot matrices were then truncated to retain 90% of the turbulence kinetic energy in the domain of interest. A parametric study is conducted to evaluate the effect of the frequency of snapshots used to construct the POD-ROM.

turbulence↗

Multiresolution Representation Using Biorthogonal Multiwavelets

We generalize Harten's multiresolution representation to biorthogonal multiwavelets. Several variants are considered. For example, a given array of discrete point values is transformed to point values and derivatives or point 'values and cell averages'. Compact Hermite interpolation is used in the decomposition and reconstruction algorithm. The resulting basis functions that are symmetric or skewsymmetric, compact, and smooth with optimal order accuracy. Harten's approach has several advantages: the multiresolution scheme is inherently discrete, non-periodic boundary conditions are easy to implement, and the representation can be extended to unstructured grids in bounded domains. We demonstrate the compression features of the new mutliwavelets by application to variable scale piecewise smooth functions with jump discontinuities typical of numerical solutions of nonlinear hyperbolic conservation laws.

Warming, Robert F.↗

Experimental Shock Decomposition of Siderite and the Origin of Magnetite in Martian Meteorite ALH84001

Shock recovery experiments to determine whether magnetite could be produced by the decomposition of iron-carbonate were initiated. Naturally occurring siderite was first characterized by electron microprobe (EMP), transmission electron microscopy (TEM), Mossbauer spectroscopy, and magnetic susceptibility measurements to be sure that the starting material did not contain detectable magnetite. Samples were shocked in tungsten-alloy holders (W=90%, Ni=6%, Cu=4%) to further insure that any iron phases in the shock products were contributed by the siderite rather than the sample holder. Each sample was shocked to a specific pressure between 30 to 49 GPa. Previously reported results of TEM analyses on 49 GPa experiments indicated the presence of nano-phase spinel-structured iron oxide. Transformation of siderite to magnetite as characterized by TEM was found in the 49 GPa shock experiment. Compositions of most magnetites are greater than 50% Fe sup(+2) in the octahedral site of the inverse spinel structure. Magnetites produced in shock experiments display the same range of single-domain, superparamagnetic sizes (approx. 50 100 nm), compositions (100% magnetite to 80% magnetite-20% magnesioferrite), and morphologies (equant, elongated, euhedral to subhedral) as magnetites synthesized by Golden et al. (2001) or magnetites grown naturally by MV1 magnetotactic bacteria, and as the magnetites in Martian meteorite ALH84001. Fritz et al. (2005) previously concluded that ALH84001 experienced approx. 32 GPa pressure and a resultant thermal pulse of approx. 100 - 110 C. However, ALH84001 contains evidence of local temperature excursions high enough to 1 melt feldspar, pyroxene, and a silica-rich phase. This 49 GPa experiment demonstrates that magnetite can be produced by the shock decomposition of siderite as a result of local heating to greater than 470 C. Therefore, magnetite in the rims of carbonates in Martian meteorite ALH84001 could be a product of shock devolatilization of siderite as well.

Bell, Mary Sue↗

Multiphysics Computational Analysis of a Solid-Core Nuclear Thermal Engine Thrust Chamber

The objective of this effort is to develop an efficient and accurate computational heat transfer methodology to predict thermal, fluid, and hydrogen environments for a hypothetical solid-core, nuclear thermal engine - the Small Engine. In addition, the effects of power profile and hydrogen conversion on heat transfer efficiency and thrust performance were also investigated. The computational methodology is based on an unstructured-grid, pressure-based, all speeds, chemically reacting, computational fluid dynamics platform, while formulations of conjugate heat transfer were implemented to describe the heat transfer from solid to hydrogen inside the solid-core reactor. The computational domain covers the entire thrust chamber so that the afore-mentioned heat transfer effects impact the thrust performance directly. The result shows that the computed core-exit gas temperature, specific impulse, and core pressure drop agree well with those of design data for the Small Engine. Finite-rate chemistry is very important in predicting the proper energy balance as naturally occurring hydrogen decomposition is endothermic. Locally strong hydrogen conversion associated with centralized power profile gives poor heat transfer efficiency and lower thrust performance. On the other hand, uniform hydrogen conversion associated with a more uniform radial power profile achieves higher heat transfer efficiency, and higher thrust performance.

Wang, Ten-See↗

Heuristic Area Cost Estimation for Observational Coverage Schedulers

This paper presents a comparison of heuris- tics used to estimate the amount of time it would take for a spacecraft to image an area using Boustrophedon decomposition (Choset and Pignon 1998). Machine learning tech- niques are used to characterize algorithmic performance of coverage algorithms. It is shown that an ordinary least-squares linear model is among the most accurate in a set of constant and linear order regression models both in terms of memory consumption and schedule duration. These are demonstrated using the ASPEN planning system (Fukunaga et al. 1997) on the Eagle Eye domain.

Knight, Russell↗

Monitoring Composite Material Pressure Vessels with a Fiber-Optic/Microelectronic Sensor System

We discuss the concept of an integrated, fiber-optic/microelectronic distributed sensor system that can monitor composite material pressure vessels for Air Force space systems to provide assessments of the overall health and integrity of the vessel throughout its entire operating history from birth to end of life. The fiber optic component would include either a semiconductor light emitting diode or diode laser and a multiplexed fiber optic sensing network incorporating Bragg grating sensors capable of detecting internal temperature and strain. The microelectronic components include a power source, a pulsed laser driver, time domain data acquisition hardware, a microprocessor, a data storage device, and a communication interface. The sensing system would be incorporated within the composite during its manufacture. The microelectronic data acquisition and logging system would record the environmental conditions to which the vessel has been subjected to during its storage and transit, e.g., the history of thermal excursions, pressure loading data, the occurrence of mechanical impacts, the presence of changing internal strain due to aging, delamination, material decomposition, etc. Data would be maintained din non-volatile memory for subsequent readout through a microcomputer interface.

Klimcak, C.↗

Global models: Robot sensing, control, and sensory-motor skills

Robotics research has begun to address the modeling and implementation of a wide variety of unstructured tasks. Examples include automated navigation, platform servicing, custom fabrication and repair, deployment and recovery, and science exploration. Such tasks are poorly described at onset; the workspace layout is partially unfamiliar, and the task control sequence is only qualitatively characterized. The robot must model the workspace, plan detailed physical actions from qualitative goals, and adapt its instantaneous control regimes to unpredicted events. Developing robust representations and computational approaches for these sensing, planning, and control functions is a major challenge. The underlying domain constraints are very general, and seem to offer little guidance for well-bounded approximation of object shape and motion, manipulation postures and trajectories, and the like. This generalized modeling problem is discussed, with an emphasis on the role of sensing. It is also discussed that unstructured tasks often have, in fact, a high degree of underlying physical symmetry, and such implicit knowledge should be drawn on to model task performance strategies in a methodological fashion. A group-theoretic decomposition of the workspace organization, task goals, and their admissible interactions are proposed. This group-mechanical approach to task representation helps to clarify the functional interplay of perception and control, in essence, describing what perception is specifically for, versus how it is generically modeled. One also gains insight how perception might logically evolve in response to needs of more complex motor skills. It is discussed why, of the many solutions that are often mathematically admissible to a given sensory motor-coordination problem, one may be preferred over others.

Schenker, Paul S.↗

An Algorithm for Integrated Subsystem Embodiment and System Synthesis

Consider the statement,'A system has two coupled subsystems, one of which dominates the design process. Each subsystem consists of discrete and continuous variables, and is solved using sequential analysis and solution.' To address this type of statement in the design of complex systems, three steps are required, namely, the embodiment of the statement in terms of entities on a computer, the mathematical formulation of subsystem models, and the resulting solution and system synthesis. In complex system decomposition, the subsystems are not isolated, self-supporting entities. Information such as constraints, goals, and design variables may be shared between entities. But many times in engineering problems, full communication and cooperation does not exist, information is incomplete, or one subsystem may dominate the design. Additionally, these engineering problems give rise to mathematical models involving nonlinear functions of both discrete and continuous design variables. In this dissertation an algorithm is developed to handle these types of scenarios for the domain-independent integration of subsystem embodiment, coordination, and system synthesis using constructs from Decision-Based Design, Game Theory, and Multidisciplinary Design Optimization. Implementation of the concept in this dissertation involves testing of the hypotheses using example problems and a motivating case study involving the design of a subsonic passenger aircraft.

Lewis, Kemper↗

A Hybrid Procedural/Deductive Executive for Autonomous Spacecraft

The New Millennium Remote Agent (NMRA) will be the first AI system to control an actual spacecraft. The spacecraft domain places a strong premium on autonomy and requires dynamic recoveries and robust concurrent execution, all in the presence of tight real-time deadlines, changing goals, scarce resource constraints, and a wide variety of possible failures. To achieve this level of execution robustness, we have integrated a procedural executive based on generic procedures with a deductive model-based executive. A procedural executive provides sophisticated control constructs such as loops, parallel activity, locks, and synchronization which are used for robust schedule execution, hierarchical task decomposition, and routine configuration management. A deductive executive provides algorithms for sophisticated state inference and optimal failure recover), planning. The integrated executive enables designers to code knowledge via a combination of procedures and declarative models, yielding a rich modeling capability suitable to the challenges of real spacecraft control. The interface between the two executives ensures both that recovery sequences are smoothly merged into high-level schedule execution and that a high degree of reactivity is retained to effectively handle additional failures during recovery.

Pell, Barney↗

Development of Fast Algorithms Using Recursion, Nesting and Iterations for Computational Electromagnetics

In the first phase of our work, we have concentrated on laying the foundation to develop fast algorithms, including the use of recursive structure like the recursive aggregate interaction matrix algorithm (RAIMA), the nested equivalence principle algorithm (NEPAL), the ray-propagation fast multipole algorithm (RPFMA), and the multi-level fast multipole algorithm (MLFMA). We have also investigated the use of curvilinear patches to build a basic method of moments code where these acceleration techniques can be used later. In the second phase, which is mainly reported on here, we have concentrated on implementing three-dimensional NEPAL on a massively parallel machine, the Connection Machine CM-5, and have been able to obtain some 3D scattering results. In order to understand the parallelization of codes on the Connection Machine, we have also studied the parallelization of 3D finite-difference time-domain (FDTD) code with PML material absorbing boundary condition (ABC). We found that simple algorithms like the FDTD with material ABC can be parallelized very well allowing us to solve within a minute a problem of over a million nodes. In addition, we have studied the use of the fast multipole method and the ray-propagation fast multipole algorithm to expedite matrix-vector multiplication in a conjugate-gradient solution to integral equations of scattering. We find that these methods are faster than LU decomposition for one incident angle, but are slower than LU decomposition when many incident angles are needed as in the monostatic RCS calculations.

Chew, W. C.↗