Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Chemical Kinetics of the TPS and Base Bleeding During Flight Test

The present research deals with thermal degradation of polyurethane foam (PUF) during flight test. Model of thermal decomposition was developed that accounts for polyurethane kinetics parameters extracted from thermogravimetric analyses and radial heat losses to the surrounding environment. The model predicts mass loss of foam, the temperature and kinetic of release of the exhaust gases and char as function of heat and radiation loads. When PUF is heated, urethane bond break into polyol and isocyanate. In the first stage, isocyanate pyrolyses and oxidizes. As a result, the thermo-char and oil droplets (yellow smoke) are released. In the second decomposition stage, pyrolysis and oxidization of liquid polyol occur. Next, the kinetics of chemical compound release and the information about the reactions occurring in the base area are coupled to the CFD simulations of the base flow in a single first stage motor vertically stacked vehicle configuration. The CFD simulations are performed to estimate the contribution of the hot out-gassing, chemical reactions, and char oxidation to the temperature rise of the base flow. The results of simulations are compared with the flight test data.

Osipov, Viatcheslav↗
NASA NTRS1994

Shortwave feedbacks and El Nino-Southern Oscillation: Forced ocean coupled ocean-atmosphere experiments

Changes in tropical sea surface temperature (SST) can produce changes in cloudiness that modify incoming solar shortwave (SW) radiation, which in turn affects SST. The effects of this negative feedback on Pacific interannual variability are examined in forced ocean model and hybrid coupled ocean-atmosphere model simulations. Two empirical schemes are used to model the large-scale, low-frequency response of surface SW to SST anomalies. The first scheme attempts to account for the nonlocal nature of the atmospheric response to SST based patterns of covariability analyzed through singular value decomposition. In the observations the primary coupled mode of variability is composed of a SW anomaly in the central Pacific that covaries with anomalous SST in the eastern Pacific. This is applied in the model as a nonlocal feedback. The second scheme examines the effects of a purely local feedback with a spatially varying coefficient of magnitude chosen similar to the first scheme. In almost all cases the second scheme behaved similarly to the first, presumably because the correlation scale of SST is large enough for El Nino-Southern Oscillation (ENSO) dynamics that there is little sensitivity to the local approximation in the SW feedback. In simulations forced by time series of observed wind stress the SW feedback induced very minor SST damping. Results for a simpified heat budget analysis showed that while the SW feedback increased the local heat flux damping on SST, it also induced a mean shallowing of the mixed layer. The resulting changes in both the local mean vertical temperature gradient and the zonal velocity response to the wind stress acted to oppose the local heat flux damping effects. When the observed SW anomalies were applied to forced simulations, the simulated SST anomalies were modified as expected, and agreement with observed SST improved. In coupled simulations the SW feedbacks had greater impact than in the case of specified stress. The main effects were to decrease the magnitude of the warm and cold SST anomalies in the central Pacific, while leaving the pattern and evolution of ENSO anomalies essentially unchanged elsewhere. The SW feedbacks thus produce a modest improvement of the model ENSO SST pattern compared with observations, although they tended to shorten the period of the model ENSO cycle. Overall the results suggest that large-scale SW feedbacks are of quantitative importance to simulating some aspects of the ENSO cycle but are not critical to the overall occurrence of the phenomenon.

Waliser, Duane E.↗
NASA NTRS1991

Proper orthogonal decomposition based turbulence modeling

The proper orthogonal decomposition is utilized to provide one mode estimates of various turbulence quantities in the near-wall region of the turbulent boundary layer. The experimental data of Herzog (1986) and the direct numerical simulation (DNS) results of Zang (1991) proved the two-point correlation tensor. The correlation tensor is then used as the kernel in the integral eigenvalue problem from which the proper orthogonal modes are obtained. The one mode estimates of the turbulent kinetic energy, dissipation, and triple moments all exhibit the proper asymtotic behavior near the wall. As a validation test, a model for the triple moments is developed and its near-wall behavior is studied. This information can then be used to obtain an estimate of the damping function near the wall, which is then compared to known existing behavior. Since this is an initial investigation, the streamwise and spanwise dependence is surpressed and only the wall normal direction (y+) behavior of various turbulence statistics near the wall is studied. Before examining one mode estimates of the triple moments and the dissipation (on which the model is based), various second order moments are examined. If these do not exhibit the proper behavior near the wall (as discussed by Speziale et al 1990) it does not seem reasonable to expect the one mode estimate of the triple moment and dissipation to have the right behavior. Only the results obtained using the experimental data of Herzog are discussed since the results from the simulations of Zang exhibit similar behavior as discussed above.

Glauser, Mark N.↗
NASA NTRS1996

Concurrent Probabilistic Simulation of High Temperature Composite Structural Response

A computational structural/material analysis and design tool which would meet industry's future demand for expedience and reduced cost is presented. This unique software 'GENOA' is dedicated to parallel and high speed analysis to perform probabilistic evaluation of high temperature composite response of aerospace systems. The development is based on detailed integration and modification of diverse fields of specialized analysis techniques and mathematical models to combine their latest innovative capabilities into a commercially viable software package. The technique is specifically designed to exploit the availability of processors to perform computationally intense probabilistic analysis assessing uncertainties in structural reliability analysis and composite micromechanics. The primary objectives which were achieved in performing the development were: (1) Utilization of the power of parallel processing and static/dynamic load balancing optimization to make the complex simulation of structure, material and processing of high temperature composite affordable; (2) Computational integration and synchronization of probabilistic mathematics, structural/material mechanics and parallel computing; (3) Implementation of an innovative multi-level domain decomposition technique to identify the inherent parallelism, and increasing convergence rates through high- and low-level processor assignment; (4) Creating the framework for Portable Paralleled architecture for the machine independent Multi Instruction Multi Data, (MIMD), Single Instruction Multi Data (SIMD), hybrid and distributed workstation type of computers; and (5) Market evaluation. The results of Phase-2 effort provides a good basis for continuation and warrants Phase-3 government, and industry partnership.

Abdi, Frank↗
NASA NTRS1998

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.↗
NASA NTRS2016

Effects of Different Heave Motion Components on Pilot Pitch Control Behavior

The study described in this paper had two objectives. The first objective was to investigate if a different weighting of heave motion components decomposed at the center of gravity, allowing for a higher fidelity of individual components, would result in pilot manual pitch control behavior and performance closer to that observed with full aircraft motion. The second objective was to investigate if decomposing the heave components at the aircraft's instantaneous center of rotation rather than at the center of gravity could result in additional improvements in heave motion fidelity. Twenty-one general aviation pilots performed a pitch attitude control task in an experiment conducted on the Vertical Motion Simulator at NASA Ames under different hexapod motion conditions. The large motion capability of the Vertical Motion Simulator also allowed for a full aircraft motion condition, which served as a baseline. The controlled dynamics were of a transport category aircraft trimmed close to the stall point. When the ratio of center of gravity pitch heave to center of gravity heave increased in the hexapod motion conditions, pilot manual control behavior and performance became increasingly more similar to what is observed with full aircraft motion. Pilot visual and motion gains significantly increased, while the visual lead time constant decreased. The pilot visual and motion time delays remained approximately constant and decreased, respectively. The neuromuscular damping and frequency both decreased, with their values more similar to what is observed with real aircraft motion when there was an equal weighting of the heave of the center of gravity and heave due to rotations about the center of gravity. In terms of open- loop performance, the disturbance and target crossover frequency increased and decreased, respectively, and their corresponding phase margins remained constant and increased, respectively. The decomposition point of the heave components only had limited effects on pilot manual control behavior and performance.

manual control↗
NASA NTRS1990

Numerical simulation of spatially-evolving instability

A computational study of the spatial stability of plane Poiseuille flow is presented. The numerical scheme employs a time-splitting method to integrate the full Navier-Stokes equations using spectral collocation/finite-difference discretization on a non-staggered mesh. The eigenvalue decomposition procedure is applied for the solution of the Poisson equations using the capacitance matrix technique. The buffer domain method is incorporated for the outflow boundary conditions. The input perturbation velocities are obtained by solving the Orr-Sommerfeld equation for the nonlinear eigenvalue problem employing the companion matrix method. Computational results are compared with the linear theory for two-dimensional disturbances.

Danabasoglu, G.↗
NASA NTRSDate not supplied

Dynamic Analysis of Reynolds Number Effects on Trailing Edge Transonic Vortex Shedding and Its Impact on Turbine Blade Aerodynamic Performance

Time-resolved, high-speed self-aligned focusing schlieren images were acquired in the NASA Glenn Research Center Transonic Turbine Blade Cascade facility to help understand the aerodynamic behavior of high-pressure, thick trailing edge turbine blades. The trailing edge thickness of 9% of axial chord tested represents simulated ceramic matrix composite fabrication constraints, which was verified previously to possess a high-loss flow regime at high inlet turbulence conditions over a narrow range of Reynolds numbers and at a fixed design exit Mach number of 0.74. Our high-speed images, which were acquired at 10 distinct Reynolds numbers, show a significant increase in energy from flow oscillations due to transonic vortex shedding at Reynolds numbers corresponding to the high loss conditions. For those conditions, strong acoustic waves turn into shock waves. Spectral Proper Orthogonal Decomposition of the high-speed images shows acoustic waves from trailing edge vortex shedding at all conditions, with increased spectral energy at the high-loss conditions and slightly increasing frequency (about 6%) as a function of Reynolds number. Analysis of potential feedback timing is performed using velocity fields from a previous LES simulation, considering different feedback mechanisms. Most noteworthy is the acoustic/shock-boundary layer interaction mechanism on the suction surface at the blade geometric throat, which likely plays an important role in realistic curved blade passages.

Trailing Edge↗
NASA NTRS1994

Optical Implementation of Matching Pursuit for Image Representation

We have developed a technique for image analysis, representation, and decomposition. This technique was motivated by Stephane Mallat's matching-pursuit algorithm. We've altered and simplified the mechanics of his algorithm to enable an extremely fast implementation via optical processing. Initial computer simulations show that our algorithm is capable of decomposing and representing a 2-D image as a linear combination of basis images with both high speed and high fidelity.

matching pursuit adaptive wavelet transform image ↗
NASA NTRS1992

A real-time robot arm collision avoidance system

A data structure and update algorithm are presented for a prototype real-time collision avoidance safety system simulating a multirobot workspace. The data structure is a variant of the octree, which serves as a spatial index. An octree recursively decomposes 3D space into eight equal cubic octants until each octant meets some decomposition criteria. The N-objects octree, which indexes a collection of 3D primitive solids is used. These primitives make up the two (seven-degrees-of-freedom) robot arms and workspace modeled by the system. As robot arms move, the octree is updated to reflect their changed positions. During most update cycles, any given primitive does not change which octree nodes it is in. Thus, modification to the octree is rarely required. Cycle time for interpreting current arm joint angles, updating the octree to reflect new positions, and detecting/reporting imminent collisions averages 30 ms on an Intel 80386 processor running at 20 MHz.

Shaffer, Clifford A.↗
NASA NTRS2024

Assessing the Impact of Climate Change on Atmospheric Rivers: A Modeling Perspective

Atmospheric rivers (ARs) play a crucial role in the poleward transport of water vapor, and the AR-associated precipitation is a critical component of global water supplies, making it critical that we understand how ARs may change in the future. To approach this issue, integrations of the NASA Goddard Institute for Space Studies global climate model ModelE version 2.1 (GISSE2.1) are employed. Multiple configurations of the model simulating different climates are analyzed: (1) the last-glacial maximum; (2) present day; (3) the end of the 21st century. The thermodynamic and dynamic components of changes to AR frequency are analyzed using a decomposition method. This method utilizes differences in distinct AR seasonal climatology frequencies derived from various vertically integrated water vapor transport (IVT) thresholds to resolve AR frequency into its components. Global mean state changes in poleward AR frequency for different climates are dominated by precipitable water vapor (PWV) changes. A set of idealized cold and warm climates in which present day sea surface temperatures are uniformly changed are considered for a targeted analysis of the south Pacific Ocean basin. For this analysis, frequency and distribution of AR events in the model runs are analyzed by comparing them to changes in the jet stream as well as the Eulerian storm tracks and low-level baroclinicity. Latitudinal shifts in the ARs in the south Pacific Ocean basin using our integrations are not as tightly coupled to these two storm-related climatological metrics in the midlatitudes but fare better on the poleward side of the storm tracks.

Atmospheric rivers↗
NASA NTRS1980

Handling Qualities of Large Flexible Aircraft

The effects on handling qualities of elastic modes interaction with the rigid body dynamics of a large flexible aircraft are studied by a mathematical computer simulation. An analytical method to predict the pilot ratings when there is a severe modes interactions is developed. This is done by extending the optimal control model of the human pilot response to include the mode decomposition mechanism into the model. The handling qualities are determined for a longitudinal tracking task using a large flexible aircraft with parametric variations in the undamped natural frequencies of the two lowest frequency, symmetric elastic modes made to induce varying amounts of mode interaction.

Poopaka, S.↗
NASA NTRS1993

A study of CFD algorithms applied to complete aircraft configurations

Research aimed at increasing the efficiency and productivity of an overall CFD simulation process, from geometry definition to post processing of results as applied to complete aircraft configurations is reviewed. The progress was achieved through enhancements in geometry processing, surface grid generation, grid refinement, grid adaptation, and domain decomposition strategies. Algorithms for surface grid redistribution, automatic remapping, elliptic refinement and grid adaptation utilize the direct and inverse nonuniform rational B spline based surface representation.

Soni, Bharat K.↗
NASA NTRS2006

Processes Influencing Ozone Levels in Alaskan Forest Fires Plumes during Long-Range Transport over the North Atlantic

A case of long-range transport of a biomass burning plume from Alaska to Europe is analyzed using a Lagrangian approach. This plume was sampled several times in the free troposphere over North America, the North Atlantic and Europe by 3 different aircraft during the IGAC Lagrangian 2K4 experiment which was part of the ICARTT/ITOP measurement intensive in summer 2004. Measurements in the plume showed enhanced values of CO, VOCs and NOy, mainly in form of PAN. Observed O3 levels increased by 17 ppbv over 5 days. A photochemical trajectory model, CiTTyCAT, is used to examine processes responsible for the chemical evolution of the plume. The model was initialized with upwind data, and compared with downwind measurements. The influence of high aerosol loading on photolysis rates in the plume is investigated using in-situ aerosol measurements in the plume and lidar retrievals of optical depth as input into a photolysis code (Fast-J), run in the model. Significant impacts on photochemistry are found with a decrease of 18 percent in O3 production and 24 percent in O3 destruction over 5 days when including aerosols. The plume is found to be chemically active with large O3 increases attributed primarily to PAN decomposition during descent of the plume towards Europe. The predicted O3 changes are very dependent on temperature changes during transport, and also, on water vapor levels in the lower troposphere which can lead to O3 destruction. Simulation of mixing/dilution was necessary to reproduce observed pollutants level in the plume. Mixing was simulated using background concentrations from measurements in air masses in close proximity to the plume, and mixing timescales (averaging 6.25 days) were derived from CO changes. Observed and simulated O3/CO correlations in the plume are also compared in order to evaluate the photochemistry in the model. Observed slopes changed from negative to positive over 5 days. This change, which can be attributed largely to photochemistry, is well reproduced by multiple model runs even if slope values are slightly underestimated suggesting small underestimation of photochemical processes. The possible impact of this biomass burning plume on O3 levels in the Europe boundary layer is also examined by running the model for a further 5 days, and comparing with data collected at surface sites, such as Jungfraujoch, which showed small O3 increases and elevated CO levels. The model predicts significant changes in O3 over this 10 days due to photochemistry but the signal is largely lost due to the effects of dilution. However, measurement in several others BB plumes over Europe show that O3 impact of Alaskan fires can be punctually significant over Europe.

Real, E.↗
NASA NTRS1989

Vortex breakdown simulation - A circumspect study of the steady, laminar, axisymmetric model

The incompressible axisymmetric steady Navier-Stokes equations are written using the streamfunction-vorticity formulation. The resulting equations are discretized using a second-order central-difference scheme. The discretized equations are linearized and then solved using an exact LU decomposition, Gaussian elimination, and Newton iteration. Solutions are presented for Reynolds numbers (based on vortex core radius) 100-1800 and swirl parameter 0.9-1.1. The effects of inflow boundary conditions, the location of farfield and outflow boundaries, and mesh refinement are examined. Finally, the stability of the steady solutions is investigated by solving the time-dependent equations.

Salas, M. D.↗
NASA NTRS1986

Control of elastic robotic systems by nonlinear inversion and modal damping

Energy efficient, lightweight robot arms for space applications have considerable structural flexibility. For large and fast motions, both the nonlinear coupled dynamics and the elastic behavior of the robots must be considered in control system designs. This paper presents an approach to the control of a class of flexible robotic systems. A control law is derived which decouples the joint-angle motion from the flexible motion and asymptotically decomposes the elastic dynamics into two subsystems, representing the transverse vibrations of the elastic link in two orthogonal planes. This decomposition allows the design of an elastic mode stabilizer independently based on lower order models representing structural flexibility. The closed-loop system is shown to be globally asymptotically stable and robust to uncertainty in system parameters. Simulation results are presented to show that large, fast control of joint angles can be performed in spite of space vehicle motion and uncertainty in the payload.

Singh, S. N.↗
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 NTRS1991

General Monte Carlo reliability simulation code including common mode failures and HARP fault/error-handling

A Monte Carlo Fortran computer program was developed that uses two variance reduction techniques for computing system reliability applicable to solving very large highly reliable fault-tolerant systems. The program is consistent with the hybrid automated reliability predictor (HARP) code which employs behavioral decomposition and complex fault-error handling models. This new capability is called MC-HARP which efficiently solves reliability models with non-constant failures rates (Weibull). Common mode failure modeling is also a specialty.

Platt, M. E.↗