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

Statistical analysis of near-wall structures in turbulent channel flow

The two-point correlation tensor is used to analyze near-wall structures in turbulent channel flow. Simulation results have been utilized to obtain the correlation tensor with sufficient spatial resolution to permit accurate differentiation for computing the vorticity correlation, and the resolution of sharp features such as the shear layer and its corresponding region of spanwise vorticity. Proper orthogonal decomposition and stochastic estimation were observed to yield similar results when appropriate conditions were employed for the estimation. Resulting structures were such that they would be detected by standard conditional sampling methods. The results of decomposing the vorticity field and the velocity field were significantly different, indicating that if a coherent structure is found that dominates the velocity fields, the curl of that velocity structure will not dominate the vorticity field.

Moser, Robert D.↗

Development of Multiobjective Optimization Techniques for Sonic Boom Minimization

A discrete, semi-analytical sensitivity analysis procedure has been developed for calculating aerodynamic design sensitivities. The sensitivities of the flow variables and the grid coordinates are numerically calculated using direct differentiation of the respective discretized governing equations. The sensitivity analysis techniques are adapted within a parabolized Navier Stokes equations solver. Aerodynamic design sensitivities for high speed wing-body configurations are calculated using the semi-analytical sensitivity analysis procedures. Representative results obtained compare well with those obtained using the finite difference approach and establish the computational efficiency and accuracy of the semi-analytical procedures. Multidisciplinary design optimization procedures have been developed for aerospace applications namely, gas turbine blades and high speed wing-body configurations. In complex applications, the coupled optimization problems are decomposed into sublevels using multilevel decomposition techniques. In cases with multiple objective functions, formal multiobjective formulation such as the Kreisselmeier-Steinhauser function approach and the modified global criteria approach have been used. Nonlinear programming techniques for continuous design variables and a hybrid optimization technique, based on a simulated annealing algorithm, for discrete design variables have been used for solving the optimization problems. The optimization procedure for gas turbine blades improves the aerodynamic and heat transfer characteristics of the blades. The two-dimensional, blade-to-blade aerodynamic analysis is performed using a panel code. The blade heat transfer analysis is performed using an in-house developed finite element procedure. The optimization procedure yields blade shapes with significantly improved velocity and temperature distributions. The multidisciplinary design optimization procedures for high speed wing-body configurations simultaneously improve the aerodynamic, the sonic boom and the structural characteristics of the aircraft. The flow solution is obtained using a comprehensive parabolized Navier Stokes solver. Sonic boom analysis is performed using an extrapolation procedure. The aircraft wing load carrying member is modeled as either an isotropic or a composite box beam. The isotropic box beam is analyzed using thin wall theory. The composite box beam is analyzed using a finite element procedure. The developed optimization procedures yield significant improvements in all the performance criteria and provide interesting design trade-offs. The semi-analytical sensitivity analysis techniques offer significant computational savings and allow the use of comprehensive analysis procedures within design optimization studies.

Chattopadhyay, Aditi↗

Design and Performance Analysis of a Massively Parallel Atmospheric General Circulation Model

In the 1990's computer manufacturers are increasingly turning to the development of parallel processor machines to meet the high performance needs of their customers. Simultaneously, atmospheric scientists study weather and climate phenomena ranging from hurricanes to El Nino to global warming that require increasingly fine resolution models. Here, implementation of a parallel atmospheric general circulation model (GCM) which exploits the power of massively parallel machines is described. Using the horizontal data domain decomposition methodology, this FORTRAN 90 model is able to integrate a 0.6 deg. longitude by 0.5 deg. latitude problem at a rate of 19 Gigaflops on 512 processors of a Cray T3E 600; corresponding to 280 seconds of wall-clock time per simulated model day. At this resolution, the model has 64 times as many degrees of freedom and performs 400 times as many floating point operations per simulated day as the model it replaces.

Schaffer, Daniel S.↗

Decomposition of Electron Ionization Mass Spectra for Space Application Using a Monte-Carlo approach

Rationale: Quadrupole mass spectrometers equipped with an electron ionization (EI) sources have been widely used in space exploration to investigate the composition of planetary surfaces and atmospheres. However, the complexity of the samples and the minimal calibration for the fragmentation of molecules in the ionization chambers have prevented the deconvolution of the majority of the mass spectra obtained at different targets, thus limiting the determination of the exact composition of the samples analyzed. We propose a Monte‐Carlo approach to solve this issue mathematically. Methods: We decomposed simulated mass spectra of mixtures acquired with unit resolving power mass spectrometers and EI sources into the sum of the single components fragmentation patterns weighted by their relative concentration using interior‐point least‐square fitting. To fit compounds with poorly known fragmentation patterns, we used a Monte‐Carlo method to vary the intensity of individual fragment ions. We then decomposed the spectrum thousands of times to obtain a statistical distribution. Results: By performing the deconvolution on a mixture of seven different molecules with interfering fragmentation patterns (H2O, O2, CH4, Ar, N2, C2H4, and C2H6) we show that this approach retrieves the mixing ratio of the individual components more accurately than regular mass spectra decomposition methods that rely on fragmentation patterns from general databases. It also provides the probability density function for each species's mixing ratio. Conclusions: By removing the solution degeneracy in the decomposition of mass spectra, the method described herein could significantly increase the scientific retrieval from archived space flight mass spectrometry data, where calibration of the ionization source is no longer an option.

Thomas Gautier↗

Design of the hybrid automated reliability predictor

The design of the Hybrid Automated Reliability Predictor (HARP), now under development at Duke University, is presented. The HARP approach to reliability prediction is characterized by a decomposition of the overall model into fault-occurrence and fault-handling sub-models. The fault-occurrence model is a non-homogeneous Markov chain which is solved analytically, while the fault-handling model is a Petri Net which is simulated. HARP provides automated analysis of sensitivity to uncertainties in the input parameters and in the initial state specifications. It then produces a predicted reliability band as a function of mission time, as well as estimates of the improvement (narrowing of the band) to be gained by a specified amount of reduction in uncertainty.

Geist, R.↗

Flexible spacecraft maneuver - Inverse attitude control and modal stabilization

A control law is presented for three-axis rotational maneuvers of a spacecraft (orbiter)-beam-tip body (antenna or a reflector) configuration based on nonlinear inversion and modal velocity feedback. Using invertibility and functional reproducibility results, a decoupling attitude control law is presented such that, in the closed-loop system, the attitude angles of the spacecraft are independently controlled using the control moments acting on the space vehicle. This controller asymptotically decouples the flexible dynamics from the rigid one and also allows the decomposition of the elastic dynamics into two subsystems representing the transverse deflections of the beam in two orthogonal planes. These low-order subsystems are used for derivation of a modal velocity feedback stabilizer using the force and moment actuators at the end body. Simulation results are presented to show that, in the closed-loop system, attitude control and elastic mode stabilization are accomplished in spite of the parameter uncertainty and disturbance torque input in the system.

Singh, Sahjendra N.↗

Optical matrix-vector processing for computational fluid dynamics

An optical processor to solve partial differential equations for computational fluid dynamics applications is considered. This application is new and original for optical processors. The algorithms that are used are optical realizations of the Newton-Raphson method for nonlinear equations and a new optical LU direct decomposition and Gauss-Seidel iterative solution to the resultant linear algebraic equations. These algorithms are used to solve Burger's equation (a specific form of the momentum equation in fluid dynamics). The nonlinear equations provide 1-D velocity data at each time step. Simulation results of optical processing with these algorithms on computational fluid dynamics data is included.

Perlee, Caroline J.↗

Emerging CFD technologies and aerospace vehicle design

With the recent focus on the needs of design and applications CFD, research groups have begun to address the traditional bottlenecks of grid generation and surface modeling. Now, a host of emerging technologies promise to shortcut or dramatically simplify the simulation process. This paper discusses the current status of these emerging technologies. It will argue that some tools are already available which can have positive impact on portions of the design cycle. However, in most cases, these tools need to be integrated into specific engineering systems and process cycles to be used effectively. The rapidly maturing status of unstructured and Cartesian approaches for inviscid simulations makes suggests the possibility of highly automated Euler-boundary layer simulations with application to loads estimation and even preliminary design. Similarly, technology is available to link block structured mesh generation algorithms with topology libraries to avoid tedious re-meshing of topologically similar configurations. Work in algorithmic based auto-blocking suggests that domain decomposition and point placement operations in multi-block mesh generation may be properly posed as problems in Computational Geometry, and following this approach may lead to robust algorithmic processes for automatic mesh generation.

Aftosmis, Michael J.↗

Developing an Accurate CFD Based Gust Model for the Truss Braced Wing Aircraft

The increased flexibility of long endurance aircraft having high aspect ratio wings necessitates attention to gust response and perhaps the incorporation of gust load alleviation. The design of civil transport aircraft with a strut or truss-braced high aspect ratio wing furthermore requires gust response analysis in the transonic cruise range. This requirement motivates the use of high fidelity nonlinear computational fluid dynamics (CFD) for gust response analysis. This paper presents the development of a CFD based gust model for the truss braced wing aircraft. A sharp-edged gust provides the gust system identification. The result of the system identification is several thousand time steps of instantaneous pressure coefficients over the entire vehicle. This data is filtered and downsampled to provide the snapshot data set from which a reduced order model is developed. A stochastic singular value decomposition algorithm is used to obtain a proper orthogonal decomposition (POD). The POD model is combined with a convolution integral to predict the time varying pressure coefficient distribution due to a novel gust profile. Finally the unsteady surface pressure response of the truss braced wing vehicle to a one-minus-cosine gust, simulated using the reduced order model, is compared with the full CFD.

Bartels, Robert E.↗

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↗

Three axis rotational maneuver and vibration stabilization of elastic spacecraft

A control law for three-axis rotational maneuvers of a spacecraft beam-tip body configuration based on non-linear inversion and modal velocity feedback is presented. A decoupling attitude control law is presented such that in the closed-loop system the attitude angles of the spacecraft are independently controlled, using the control moments acting on the space vehicle. This controller asymptotically decouples the flexible dynamics from the rigid one and also allows the decomposition of the elastic dynamics into two subsystems representing the transverse deflections of the beam in two orthogonal planes. These low-order subsystems are used for the derivation of a modal velocity feedback stabilizer using the force and moment actuators at the end body. Simulation results are presented to show the capability of the controller.

Singh, Sahjendra N.↗

Flexible Ablators Char Depths LHMEL Test Results

Char and pyrolysis zone depths give physical evidence of peak temperature reached in depth: The pyrolyzing material acts as a temperature indicator within its characteristic thermal decomposition range. A matrix of novel flexible ablators were laser tested in one component of material screening for NASA Entry, Descent and Landing research for future Mars missions. LHMEL tests were run both on virgin materials, and on previously charred materials for a dual pulse simulation of the heating due to aerocapture followed by atmospheric entry. The test models were machined to expose the cross-sections. Char measurements were made at three locations near the center of the exposed area. Data are presented showing the char depths developed in these flexible materials, grouped by reinforcing fiber and pyrolyzing material type.

White, Susan↗

Experimental Study of the Stability of Aircraft Fuels at Elevated Temperatures

An experimental study of fuel stability was conducted in an apparatus which simulated an aircraft gas turbine fuel system. Two fuels were tested: Jet A and Number 2 Home Heating oil. Jet A is an aircraft gas turbine fuel currently in wide use. No. 2HH was selected to represent the properties of future turbine fuels, particularly experimental Reference Broad Specification, which, under NASA sponsorship, was considered as a possible next-generation fuel. Tests were conducted with varying fuel flow rates, delivery pressures and fuel pretreatments (including preheating and deoxygenation). Simulator wall temperatures were varied between 422K and 672K at fuel flows of 0.022 to 0.22 Kg/sec. Coking rate was determined at four equally-spaced locations along the length of the simulator. Fuel samples were collected for infrared analysis. The dependence of coking rate in Jet A may be correlated with surface temperature via an activation energy of 9 to 10 kcal/mole, although the results indicate that both bulk fluid and surface temperature affect the rate of decomposition. As a consequence, flow rate, which controls bulk temperature, must also be considered. Taken together, these results suggest that the decomposition reactions are initiated on the surface and continue in the bulk fluid. The coking rate data for No. 2 HH oil are very highly temperature dependent above approximately 533K. This suggests that bulk phase reactions can become controlling in the formation of coke.

Vranos, A.↗

Use of laboratory simulated pyrolysis in tracing the history of sedimentary organic matter

Results from laboratory simulated pyrolyses experiments show that in addition to depth of burial, preservation of kerogen, and hence any morphologic structure in it, is also dependent on the mineral matrix with which it is associated. In the presence of clay minerals, and especially under dry conditions, extractable lipids released during kerogen decomposition are more rapidly destroyed than in the presence of calcite or chert matrices. The result is production of gas, polar bitumen and pyrobitumen and destruction of biomarkers. During such an early reorganization of the kerogen, the biomarker constituents can be destroyed, or unrecognizably altered. The above process of organic residues maturation appears to be inhibited in the presence of water and is significantly reduced where kerogen is hosted in limestones, dolomites or cherts. These minerals have been characteristically found to be the most reliable in yielding morphological fossils and small quantities of extractable bitumen in Archean and Proterozoic rocks. To understand the validity of chemical and morphological fossils, in the early geologic record, it will be necessary to understand the process of kerogen in sedimentary rocks. To test the role of various minerals on the preservation process, kerogen extracted from a variety of rocks has been heated together with montmorillonite, illite and calcite. The kinetics of the process has been monitored and the products quantitatively identified.

Kaplan, I. R.↗

Chemical Characterization and Reactivity of Fuel-Oxidizer Reaction Product

Fuel-oxidizer reaction product (FORP), the product of incomplete reaction of monomethylhydrazine and nitrogen tetroxide propellants prepared under laboratory conditions and from firings of Shuttle Reaction Control System thrusters, has been characterized by chemical and thermal analysis. The composition of FORP is variable but falls within a limited range of compositions that depend on three factors: the fuel-oxidizer ratio at the time of formation; whether the composition of the post-formation atmosphere is reducing or oxidizing; and the reaction or post-reaction temperature. A typical composition contains methylhydrazinium nitrate, ammonium nitrate, methylammonium nitrate, and trace amounts of hydrazinium nitrate and 1,1-dimethylhydrazinium nitrate. Thermal decomposition reactions of the FORP compositions used in this study were unremarkable. Neither the various compositions of FORP, the pure major components of FORP, nor mixtures of FORP with propellant system corrosion products showed any unusual thermal activity when decomposed under laboratory conditions. Off-limit thruster operations were simulated by rapid mixing of liquid monomethylhydrazine and liquid nitrogen tetroxide in a confined space. These tests demonstrated that monomethylhydrazine, methylhydrazinium nitrate, ammonium nitrate, or Inconel corrosion products can induce a mixture of monomethylhydrazine and nitrogen tetroxide to produce component-damaging energies. Damaging events required FORP or metal salts to be present at the initial mixing of monomethylhydrazine and nitrogen tetroxide.

David, Dennis D.↗

Control system software, simulation, and robotic applications

All essential existing capabilities needed to create a man-machine interaction dynamics and performance (MMIDAP) capability are reviewed. The multibody system dynamics software program Order N DISCOS will be used for machine and musculo-skeletal dynamics modeling. The program JACK will be used for estimating and animating whole body human response to given loading situations and motion constraints. The basic elements of performance (BEP) task decomposition methodologies associated with the Human Performance Institute database will be used for performance assessment. Techniques for resolving the statically indeterminant muscular load sharing problem will be used for a detailed understanding of potential musculotendon or ligamentous fatigue, pain, discomfort, and trauma. The envisioned capacity is to be used for mechanical system design, human performance assessment, extrapolation of man/machine interaction test data, biomedical engineering, and soft prototyping within a concurrent engineering (CE) system.

Frisch, Harold P.↗

Evaluation of the Dornier Gmbh interactive grid generation system

An interactive grid generation program, INGRID, is investigated and evaluated. A description of the task and work performed, a description and evaluation of INGRID, and a discussion of the possibilities for bringing INGRID into the NASA and Numerical Aerodynamic Simulator (NAS) computing environments is included. The interactive grid generation program was found to be a viable approach for grid generation and determined that it could be converted to work in the NAS environment but that INGRID does not solve the fundamentally hard problems associated with grid generation, specifically, domain decomposition.

Brown, Robert L.↗

Investigation of Transitional Shock-Wave/Boundary Layer Interactions Using Direct Numerical Simulations

Interaction of a transitional boundary layer with a normal-shock is investigated using unstructured tetrahedral meshes under the numerical framework of the space-time conservation element, solution element (CESE) method. The computations mimic recent experimental efforts at the University of Tennessee Space Institute, where a Mach 2.0 flow interacts with a tall cylinder attached to a flat plate. The location of the cylinder with respect to the flat plate leading edge determines if the incoming boundary layer is laminar, transitional or fully turbulent. Four representative flow conditions exemplifying laminar and transitional boundary layers are analyzed by direct numerical simulations. Similar to what was observed in the experiments for the case of transitional interaction, the computations reveal an intermittent upstream influence (UI) shock that repeatedly travels upstream from the lambda-foot toward the leading edge before vanishing. Through detailed unsteady flow analysis obtained using Fourier analysis and dynamic mode decomposition techniques, the presence of disturbances with similar frequencies as those measured in experiments were identified in the flow along with locations that appear to influence the dynamics of the flow.

Venkatachari, Balaji Shankar↗