Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational homogenization”

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 19 records

MIPS: The good, the bad and the useful

Many authors are critical of the use of MIPS (Millions of Instructions per Second) as a measure of computer power. Some feel that MIPS are meaningless. While there is justification for some of the criticism of MIPS, sometimes the criticism is carried too far. MIPS can be a useful number for planning and estimating purposes when used in a homogeneous computer environmnet. Comparisons between published MIPS ratings and benchmark results reveal that there does exist a high positive correlation between MIPS and tested performance, given a homogeneous computer environment. MIPS should be understood so as not to be misused. It is not correct that the use of MIPS is always inappropriate or inaccurate

Richardson, Jerry K.↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool for the analysis of flowfields near the injector face has been developed and used to analyze the Fastrac engine. The objective was to produce a computational design tool which was detailed enough to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows. To obtain a model which could be used to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe liquid and vapor sub- and super-critical flows, the model included thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. A homogeneous model was constructed such that the local state of the flow was determined directly, i.e. the quality of the flow was calculated. Such a model does not identify drops or their distribution, but it does allow the flow along the injector face and into the acoustic cavity to be predicted. It also allows the film coolant flow to be accurately described. The initial evaluation of the injector code was made by simulating cold flow from an unlike injector element and from a like-on-like overlapping fan (LOL) injector element. The predicted mass flux distributions of these injector elements compared well to cold flow test results. These are the same cold flow tests which serve as the data base for the JANNAF performance prediction codes. The flux distributions 1 inch downstream of the injector face are very similar; the differences were somewhat larger at further distances from the faceplate. Since the cold flow testing did not achieve good mass balances when integrations across the entire fan were made, the CFD simulation appears to be reasonable alternative to future cold flow testing. To simulate the Fastrac, an RP-1/LOX combustion model must be chosen. This submodel must be relatively simple to accomplish three-dimensional, multiphase flow simulations. Single RP-1 pyrolysis and partial oxidation steps were chosen and the combustion was completed with the wet CO mechanism. Soot was also formed with a single global reaction. To validate the combustion submodel, global data from gas generator tests and from subscale motor test were used to predict qualitatively correct mean molecular weights, temperature, and soot levels. Because such tests do not provide general kinetics rates, the methodology is not necessarily appropriate for other than rocket type flows conditions. Soot predictions were made so that radiation heating to the motor walls can be made. These initial studies of the Fastrac were for a small region close to the injector face and chamber wall which included a segment of the acoustic cavity. The region analyzed includes 11 individual orifice holes to represent the LOL elements and the H2 film coolant holes. Typical results of this simulation are shown in Figure 1. At this point the only available test data to verify the predictions are temperatures measured in the acoustic cavity. These temperatures are in reasonable agreement at about 2000R (1111 K). Future work is expected to include improving the computational efficiency or the CFD model and/or using more computer capacity than the single Pentium PC with which these simulations were made.

Farmer, Richard C.↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool which emphasizes the analysis of flowfields near the injector face of liquid rocket engines has been developed and used to simulate preliminary configurations of NASA's Fastrac and vortex engines. This computational design tool is sufficiently detailed to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows and the combusting flow which results. In order to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe sub- and supercritical liquid and vapor flows, the model utilized thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. The model was constructed such that the local quality of the flow was determined directly. Since both the Fastrac and vortex engines utilize RP-1/LOX propellants, a simplified hydrocarbon combustion model was devised in order to accomplish three-dimensional, multiphase flow simulations. Such a model does not identify drops or their distribution, but it does allow the recirculating flow along the injector face and into the acoustic cavity and the film coolant flow to be accurately predicted.

Farmer, R.↗

Effects of organized turbulence structures on the phase distortion in a coherent optical beam propagating through a turbulent shear flow

Phase distortion in a coherent optical beam propagating through a turbulent shear flow is studied. The instantaneous distribution of the index refraction is represented by a passive-scalar field in a computed homogeneous shear flow. The flow contains organized vortical structures (hairpin eddies), which are characteristic of turbulent shear flows. The phase distortion induced by turbulent fluctuations is calculated from the optical path difference through the flow. A conceptual model is proposed for the distribution of scalar fluctuations produced by the hairpin vortices in the shear flow. It is shown that the phase distortion of an optical beam can be minimized by propagating the beam at an angle approximately normal to the organized vortical structures in a turbulent shear flow.

Truman, C. Randall↗

Micromechanics for particulate reinforced composites

A set of micromechanics equations for the analysis of particulate reinforced composites is developed using the mechanics of materials approach. Simplified equations are used to compute homogenized or equivalent thermal and mechanical properties of particulate reinforced composites in terms of the properties of the constituent materials. The microstress equations are also presented here to decompose the applied stresses on the overall composite to the microstresses in the constituent materials. The properties of a 'generic' particulate composite as well as those of a particle reinforced metal matrix composite are predicted and compared with other theories as well as some experimental data. The micromechanics predictions are in excellent agreement with the measured values.

Murthy, Pappu L. N.↗

Materials with periodic internal structure: Computation based on homogenization and comparison with experiment

The combination of thermal and mechanical loading expected in practice means that constitutive equations of metal matrix composites must be developed which deal with time-independent and time-dependent irreversible deformation. Also, the internal state of composites is extremely complicated which underlines the need to formulate macroscopic constitutive equations with a limited number of state variables which represent the internal state at the micro level. One available method for calculating the macro properties of composites in terms of the distribution and properties of the constituent materials is the method of homogenization whose formulation is based on the periodicity of the substructure of the composite. A homogenization procedure was developed which lends itself to the use of the finite element procedure. The efficiency of these procedures, to determine the macroscopic properties of a composite system from its constituent properties, was demonstrated utilizing an aluminum plate perforated by directionally oriented slits. The selection of this problem is based on the fact that, extensive experimental results exist, the macroscopic response is highly anisotropic, and that the slits provide very high stress gradients which severely test the effectiveness of the computational procedures. Furthermore, both elastic and plastic properties were investigated so that the application to practical systems with inelastic deformation should be able to proceed without difficulty. The effectiveness of the procedures was rigorously checked against experimental results and with the predictions of approximate calculations. Using the computational results it is illustrated how macroscopic constitutive equations can be expressed in forms of the elastic and limit load behavior.

Jansson, S.↗

Reactant conversion in homogeneous turbulence - Mathematical modeling, computational validations, and practical applications

The presently obtained closed-form analytical expressions, which predict the limiting rate of mean reactant conversion in homogeneous turbulent flows under the influence of a binary reaction, are derived via the single-point pdf method based on amplitude mapping closure. With this model, the maximum rate of the mean reactant's decay can be conveniently expressed in terms of definite integrals of the parabolic cylinder functions. The results obtained are shown to be in good agreement with data generated by direct numerical simulations.

Madnia, C. K.↗

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗

Computer program for thin-wire structures in a homogeneous conducting medium

A computer program is presented for thin-wire antennas and scatters in a homogeneous conducting medium. The anaylsis is performed in the real or complex frequency domain. The program handles insulated and bare wires with finite conductivity and lumped loads. The output data includes the current distribution, impedance, radiation efficiency, gain, absorption cross section, scattering cross section, echo area and the polarization scattering matrix. The program uses sinusoidal bases and Galerkin's method.

Richmond, J. H.↗

Improved methods for large-eddy simulations of turbulence

By using Fourier transforms for evaluating spatial derivatives, the accuracy of the large-eddy simulation of homogeneous isotropic turbulence is improved. Especially, the treatment of certain terms that arise in filtering the equations is considerably improved in both speed and accuracy. Use of vorticity as the principal variable is shown to be a viable and potentially useful alternative to the primitive variables. A method of deriving conservation properties of numerical schemes is given which is much simpler than previous methods and is widely applicable. The methods are applied to the computation of homogeneous isotropic turbulence, and it is found that the subgrid scale model is improved by using finite differences in place of 'exact' derivatives.

Mansour, N. N.↗

Graph Partitioning for Parallel Applications in Heterogeneous Grid Environments

The problem of partitioning irregular graphs and meshes for parallel computations on homogeneous systems has been extensively studied. However, these partitioning schemes fail when the target system architecture exhibits heterogeneity in resource characteristics. With the emergence of technologies such as the Grid, it is imperative to study the partitioning problem taking into consideration the differing capabilities of such distributed heterogeneous systems. In our model, the heterogeneous system consists of processors with varying processing power and an underlying non-uniform communication network. We present in this paper a novel multilevel partitioning scheme for irregular graphs and meshes, that takes into account issues pertinent to Grid computing environments. Our partitioning algorithm, called MiniMax, generates and maps partitions onto a heterogeneous system with the objective of minimizing the maximum execution time of the parallel distributed application. For experimental performance study, we have considered both a realistic mesh problem from NASA as well as synthetic workloads. Simulation results demonstrate that MiniMax generates high quality partitions for various classes of applications targeted for parallel execution in a distributed heterogeneous environment.

Bisws, Rupak↗

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗

The preparation and archiving of machine-readable astronomical data

The preparation of machine-readable astronomical catalogs for archiving and distribution from the ADC is described. Examples of specific format and data structure changes to improve compatibility with other computers, to homogenize the presentation of astronomical data among catalogs, and to increase storage efficiency are given. Various checks which are often made to insure internal consistency and for comparison with published versions are identified and the importance of permanently archiving astronomical data in the best possible condition is emphasized.

Warren, W. H., Jr.↗

Micromechanics-Based Inelastic Finite Element Analysis Accomplished Via Seamless Integration of MAC/GMC

A critical issue in the micromechanics-based analysis of composite structures becomes the availability of a computationally efficient homogenization technique: one that is 1) Capable of handling the sophisticated, physically based, viscoelastoplastic constitutive and life models for each constituent; 2) Able to generate accurate displacement and stress fields at both the macro and the micro levels; 3) Compatible with the finite element method. The Generalized Method of Cells (GMC) developed by Paley and Aboudi is one such micromechanical model that has been shown to predict accurately the overall macro behavior of various types of composites given the required constituent properties. Specifically, the method provides "closed-form" expressions for the macroscopic composite response in terms of the properties, size, shape, distribution, and response of the individual constituents or phases that make up the material. Furthermore, expressions relating the internal stress and strain fields in the individual constituents in terms of the macroscopically applied stresses and strains are available through strain or stress concentration matrices. These expressions make possible the investigation of failure processes at the microscopic level at each step of an applied load history.

Arnold, Steven M.↗

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗

Micromechanical Modeling of High-Strain Thin-Ply Composites

This paper presents a micromechanical analysis for the elastic and viscoelastic behavior of high-strain thin-ply composites. The modeling approach is based on unit cell homogenization. The geometry of the internal weave architecture and ply configuration is characterized via micrographic analysis and explicitly modeled in the unit cell. The composites are modeled as Kirchhoff plates and the homogenization analysis computes the effective relaxation ABD matrix represented by Prony series using the elastic and viscoelastic properties of the constituent fiber and matrix. The formulation of the micromechanical model and numerical implementation d. Composite laminates with 3-ply and 4-ply configurations are studied.

Hamillage, Milinda Yapa↗

Critical assessment of Reynolds stress turbulence models using homogeneous flows

In modeling the rapid part of the pressure correlation term in the Reynolds stress transport equations, extensive use has been made of its exact properties which were first suggested by Rotta. These, for example, have been employed in obtaining the widely used Launder, Reece and Rodi (LRR) model. Some recent proposals have dropped one of these properties to obtain new models. We demonstrate, by computing some simple homogeneous flows, that doing so does not lead to any significant improvements over the LRR model and it is not the right direction in improving the performance of existing models. The reason for this, in our opinion, is that violation of one of the exact properties can not bring in any new physics into the model. We compute thirteen homogeneous flows using LRR (with a recalibrated rapid term constant), IP and SSG models. The flows computed include the flow through axisymmetric contraction; axisymmetric expansion; distortion by plane strain; and homogeneous shear flows with and without rotation. Results show that for most general representation for a model linear in the anisotropic tensor, performs either better or as good as the other two models of the same level.

Shabbir, Aamir↗