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 127 records · Page 7

A Hydrogen Peroxide Hot-Jet Simulator for Wind-Tunnel Tests of Turbojet-Exit Models

A turbojet-engine-exhaust simulator which utilizes a hydrogen peroxide gas generator has been developed for powered-model testing in wind tunnels with air exchange. Catalytic decomposition of concentrated hydrogen peroxide provides a convenient and easily controlled method of providing a hot jet with characteristics that correspond closely to the jet of a gas turbine engine. The problems associated with simulation of jet exhausts in a transonic wind tunnel which led to the selection of a liquid monopropellant are discussed. The operation of the jet simulator consisting of a thrust balance, gas generator, exit nozzle, and auxiliary control system is described. Static-test data obtained with convergent nozzles are presented and shown to be in good agreement with ideal calculated values.

Runckel, Jack F.↗

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗

Numerical simulation of a compressible homogeneous, turbulent shear flow

A direct, low Reynolds number, numerical simulation was performed on a homogeneous turbulent shear flow. The full compressible Navier-Stokes equations were used in a simulation on the ILLIAC IV computer with a 64,000 mesh. The flow fields generated by the code are used as an experimental data base, to examine the behavior of the Reynols stresses in this simple, compressible flow. The variation of the structure of the stresses and their dynamic equations as the character of the flow changed is emphasized. The structure of the tress tensor is more heavily dependent on the shear number and less on the fluctuating Mach number. The pressure-strain correlation tensor in the dynamic uations is directly calculated in this simulation. These correlations are decomposed into several parts, as contrasted with the traditional incompressible decomposition into two parts. The performance of existing models for the conventional terms is examined, and a model is proposed for the 'mean fluctuating' part.

Feiereisen, W. J.↗

Parallel DSMC Solution of Three-Dimensional Flow Over a Finite Flat Plate

This paper describes a parallel implementation of the direct simulation Monte Carlo (DSMC) method. Runtime library support is used for scheduling and execution of communication between nodes, and domain decomposition is performed dynamically to maintain a good load balance. Performance tests are conducted using the code to evaluate various remapping and remapping-interval policies, and it is shown that a one-dimensional chain-partitioning method works best for the problems considered. The parallel code is then used to simulate the Mach 20 nitrogen flow over a finite-thickness flat plate. It is shown that the parallel algorithm produces results which compare well with experimental data. Moreover, it yields significantly faster execution times than the scalar code, as well as very good load-balance characteristics.

Nance, Robert P.↗

Smectite clays in Mars soil - Evidence for their presence and role in Viking biology experimental results

Evidence for the presence of smectite clays in Martian soils is reviewed and results of experiments with certain active clays simulating the Viking biology experiments are reported. Analyses of Martian soil composition by means of X-ray fluorescence spectrometry and dust storm spectroscopy and Martian geological history strongly suggest the presence of a mixture of weathered ferro-silicate minerals, mainly nontronite and montmorillonite, accompanied by soluble sulphate salts, as major constituents. Samples of montmorillonite and nontronite incubated with (C-14)-formate or the radioactive nutrient medium solution used in the Viking Labeled Release experiment, were found to produce patterns of release of radioactive gas very similar to those observed in the Viking experiments, indicating the iron-catalyzed decomposition of formate as the reaction responsible for the Viking results. The experimental results of Hubbard (1979) simulating the results of the Viking Pyrolytic Release experiment using iron montmorillonites are pointed out, and it is concluded that many of the results of the Viking biology experiments can be explained in terms of the surface activity of smectite clays in catalysis and adsorption.

Banin, A.↗

Alternative Modal Basis Selection Procedures for Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of the three reduced-order analyses are compared with the results of the computationally taxing simulation in the physical degrees of freedom. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

An optimized index of human cardiovascular adaptation to simulated weightlessness

Prolonged exposure to weightlessness is known to produce a variety of cardiovascular changes, some of which may influence the astronaut's performance during a mission. In order to find a reliable indicator of cardiovascular adaptation to weightlessness, we analyzed data from nine male subjects after a 24-hour period of normal activity and after a period of simulated weightlessness produced by two hours in a launch position followed by 20 hours of 6 degrees head-down tilt plus pharmacologically induced diuresis (furosemide). Heart rate, arterial pressure, thoracic fluid index, and radial flow were analyzed. Autoregressive spectral estimation and decomposition were used to obtain the spectral components of each variable from the subjects in the supine position during pre- and post-simulated weightlessness. We found a significant decrease in heart rate power and an increase in thoracic fluid index power in the high frequency region (0.2-0.45 Hz) and significant increases in radial flow and arterial pressure powers in the low frequency region (<0.2 Hz) in response to simulated weightlessness. However, due to the variability among subjects, any single variable appeared limited as a dependable index of cardiovascular adaptation to weightlessness. The backward elimination algorithm was then used to select the best discriminatory features from these spectral components. Fisher's linear discriminant and Bayes' quadratic discriminant were used to combine the selected features to obtain an optimal index of adaptation to simulated weightlessness. Results showed that both techniques provided improved discriminant performance over any single variable and thus have the potential for use as an index to track adaptation and prescribe countermeasures to the effects of weightlessness.

NASA Discipline Number 14-10↗

Lumley's PODT definition of large eddies and a trio of numerical procedures

Lumley's 1967 Moscow paper provided, for the first time, a completely rational definition of the physically-useful term 'large eddy', popular for a half-century. The numerical procedures based upon his results are: (1) PODT (Proper Orthogonal Decomposition Theorem), which extracts the Large Eddy structure of stochastic processes from physical or computer simulation two-point covariances, and 2) LEIM (Large-Eddy Interaction Model), a predictive scheme for the dynamical large eddies based upon higher order turbulence modeling. Earlier Lumley's work (1964) forms the basis for the final member of the triad of numerical procedures: this predicts the global neutral modes of turbulence which have surprising agreement with both structural eigenmodes and those obtained from the dynamical equations. The ultimate goal of improved engineering design tools for turbulence may be near at hand, partly due to the power and storage of 'supermicrocomputer' workstations finally becoming adequate for the demanding numerics of these procedures.

Payne, Fred R.↗

Upward Flame Spread Over Thin Solids in Partial Gravity

The effects of partial-gravity, reduced pressure, and sample width on upward flame spread over a thin cellulose fuel were studied experimentally and the results were compared to a numerical flame spread simulation. Fuel samples 1-cm, 2-cm, and 4-cm wide were burned in air at reduced pressures of 0.2 to 0.4 atmospheres in simulated gravity environments of 0.1-G, 0.16-G (Lunar), and 0.38-G (Martian) onboard the NASA KC-135 aircraft and in normal-gravity tests. Observed steady flame propagation speeds and pyrolysis lengths were approximately proportional to the gravity level. Flames spread more quickly and were longer with the wider samples and the variations with gravity and pressure increased with sample width. A numerical simulation of upward flame spread was developed including three-dimensional Navier-Stokes equations, one-step Arrhenius kinetics for the gas phase flame and for the solid surface decomposition, and a fuel-surface radiative loss. The model provides detailed structure of flame temperatures, the flow field interactions with the flame, and the solid fuel mass disappearance. The simulation agrees with experimental flame spread rates and their dependence on gravity level but predicts a wider flammable region than found by experiment. Some unique three-dimensional flame features are demonstrated in the model results.

Feier, I. I.↗

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

Fire Safety↗

Determining the Cause of Reduced Concurrent Flame Spread over Thin Solid Fuels in Low Pressure and Low Gravity

The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.

microgravity↗

Computer Simulation of the far Infrared Collision Induced Absorption Spectra of Gaseous CO2

Far infrared collision induced absorption spectra of gaseous CO2 were computed using molecular dynamics simulations. The quadrupole and hexadecapole multipolar induction, through the trace, and the anisotropy of the molecular polarizability were found to be insufficient to represent properly the dipole induction mechanism. For a detailed analysis of the induction process the spectra obtained were decomposed into components resulting from different terms of the induced dipole. Based on this decomposition, all additional overlap contribution for each term was proposed. When spectra were recomputed including such overlap, good agreement between experiment and simulation was achieved over the temperature range at which measurements exist (233-400 K). The use of an anisotropic intermolecular potential was found to be of critical importance for obtaining the right shape of the far wings of the spectra.

Gruszka, Marcin↗

A model for the simulation of atmospheric turbulence

A method that produces realistic simulations of atmospheric turbulence is developed and analyzed. The procedure makes use of a generalized spectral analysis, often called a proper orthogonal decomposition or the Karhunen-Loeve expansion. A set of criteria, emphasizing a realistic appearance, a correct spectral shape, and non-Gaussian statistics, is selected in order to evaluate the model turbulence. An actual turbulence record is analyzed in detail, providing both a standard for comparison and input statistics for the generalized spectral analysis. The simulation utilizes the eigenfunction-expansion procedure to produce preliminary time histories of the three velocity components simultaneously. It is found that important cross-statistical features are reasonably well-behaved.

Petersen, E. L.↗

Polar decomposition for attitude determination from vector observations

This work treats the problem of weighted least squares fitting of a 3D Euclidean-coordinate transformation matrix to a set of unit vectors measured in the reference and transformed coordinates. A closed-form analytic solution to the problem is re-derived. The fact that the solution is the closest orthogonal matrix to some matrix defined on the measured vectors and their weights is clearly demonstrated. Several known algorithms for computing the analytic closed form solution are considered. An algorithm is discussed which is based on the polar decomposition of matrices into the closest unitary matrix to the decomposed matrix and a Hermitian matrix. A somewhat longer improved algorithm is suggested too. A comparison of several algorithms is carried out using simulated data as well as real data from the Upper Atmosphere Research Satellite. The comparison is based on accuracy and time consumption. It is concluded that the algorithms based on polar decomposition yield a simple although somewhat less accurate solution. The precision of the latter algorithms increase with the number of the measured vectors and with the accuracy of their measurement.

Bar-Itzhack, Itzhack Y.↗

The Thermal Stability of Unsymmetrical Dimethylhydrazine

The thermal stability of unsymmetrical dimethylhydrazine was investigated in a static system simulating conditions in an almost-empty fuel tank. The self-ignition temperature and spontaneous decomposition temperature of the pure fuel were determined at atmospheric pressure to be 454 and 740 F. respectively, with the larger (740 F) value, obtained in an inert atmosphere of nitrogen, representing the minimum temperature that would cause a rapid exothermic reaction. The addition of 40 weight percent dimethylenetriamine to unsymmetrical dimethylhydrazine did not significantly affect these properties.

Spakowski, Adolph E.↗

Alternative Modal Basis Selection Procedures For Reduced-Order Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of a computationally taxing full-order analysis in physical degrees of freedom are taken as the benchmark for comparison with the results from the three reduced-order analyses. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗