Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Homogenization theory”

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 253 records · Page 14

Correlation between strength properties in standard test specimens and molded phenolic parts

This report describes an investigation of the tensile, flexural, and impact properties of 10 selected types of phenolic molding materials. The materials were studied to see in what ways and to what extent their properties satisfy some assumptions on which the theory of strength of materials is based: namely, (a) isotropy, (b) linear stress-strain relationship for small strains, and (c) homogeneity. The effect of changing the dimensions of tensile and flexural specimens and the span-depth ratio in flexural tests were studied. The strengths of molded boxes and flexural specimens cut from the boxes were compared with results of tests on standard test specimens molded from the respective materials. The nonuniformity of a material, which is indicated by the coefficient of variation, affects the results of tests made with specimens of different sizes and tests with different methods of loading. The strength values were found to depend on the relationship between size and shape of the molded specimen and size and shape of the fillers. The most significant variations observed within a diversified group of materials were found to depend on the orientation of fibrous fillers. Of secondary importance was the dependence of the variability of test results on the pieces of filler incorporated into the molding powder as well as on the size of the piece. Static breaking strength tests on boxes molded from six representative phenolic materials correlated well with falling-ball impact tests on specimens cut from molded flat sheets. Good correlation was obtained with Izod impact tests on standard test specimens prepared from the molding materials. The static breaking strengths of the boxes do not correlate with the results of tensile or flexural tests on standard specimens.

MATERIALS, MONSANTO RESINOX↗

Dipole Alignment in Rotating MHD Turbulence

We present numerical results from long-term CPU and GPU simulations of rotating, homogeneous, magnetohydrodynamic (MHD) turbulence, and discuss their connection to the spherically bounded case. We compare our numerical results with a statistical theory of geodynamo action that has evolved from the absolute equilibrium ensemble theory of ideal MHD turbulence, which is based on the ideal MHD invariants are energy, cross helicity and magnetic helicity. However, for rotating MHD turbulence, the cross helicity is no longer an exact invariant, although rms cross helicity becomes quasistationary during an ideal MHD simulation. This and the anisotropy imposed by rotation suggests an ansatz in which an effective, nonzero value of cross helicity is assigned to axisymmetric modes and zero cross helicity to non-axisymmetric modes. This hybrid statistics predicts a large-scale quasistationary magnetic field due to broken ergodicity , as well as dipole vector alignment with the rotation axis, both of which are observed numerically. We find that only a relatively small value of effective cross helicity leads to the prediction of a dipole moment vector that is closely aligned (less than 10 degrees) with the rotation axis. We also discuss the effect of initial conditions, dissipation and grid size on the numerical simulations and statistical theory.

Shebalin, John V.↗

The Dirichlet problem for the two-dimensional Helmholtz equation for an open boundary

Development of a complete theory of the two-dimensional Dirichlet problem for an open boundary. It is shown that the solution of the Dirichlet problem for an open boundary requires the solution of a Fredholm integral equation of the first kind. Although a Fredholm integral equation of the first kind usually has no solution if the kernel is continuous, owing to the logarithmic singularity of the kernel, the equation in this case is converted to a singular integral equation with a Cauchy kernel. It is proven that the homogeneous adjoint equation of the singular integral equation has no nonzero solution. By virtue of this result, and with the aid of an existence theorem known in the theory of singular integral equations, the existence of solutions of the singular integral equation, and then of the unique solution of the Fredholm integral equation of the first kind is proved.

Hayashi, Y.↗

On the instabilities of supersonic mixing layers - A high-Mach-number asymptotic theory

The stability of a family of tanh mixing layers is studied at large Mach numbers using perturbation methods. It is found that the eigenfunction develops a multilayered structure, and the eigenvalue is obtained by solving a simplified version of the Rayleigh equation (with homogeneous boundary conditions) in one of these layers which lies in either of the external streams. This analysis leads to a simple hypersonic similarity law which explains how spatial and temporal phase speeds and growth rates scale with Mach number and temperature ratio. Comparisons are made with numerical results, and it is found that this similarity law provides a good qualitative guide for the behavior of the instability at high Mach numbers. In addition to this asymptotic theory, some fully numerical results are also presented (with no limitation on the Mach number) in order to explain the origin of the hypersonic modes (through mode splitting) and to discuss the role of oblique modes over a very wide range of Mach number and temperature ratio.

Balsa, Thomas F.↗

Nearly incompressible magnetohydrodynamics, pseudosound, and solar wind fluctuations

The results of previous investigations of low Mach number MHD density fluctuations and associated dynamical structure of nearly compressible polytropic MHD are reviewed. The recipe for computing 'pseudosound' density fluctuations is generalized to give rise to formulas expressing the density spectrum in terms of the incompressible field variables in homogeneous MHD turbulence with arbitrary rotational symmetry and with arbitrary magnetic and cross helicity spectra. The known observational consequences of the polytropic theory are summarized and several tests of the theory using Voyager data are presented. Finally, compressible two-dimensional MHD simulation results are presented which appear to have a direct bearing on interpretation of the observations.

Matthaeus, William H.↗

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗

On the generation of magnetohydrodynamic waves in a stratified and magnetized fluid. II - Magnetohydrodynamic energy fluxes for late-type stars

Magnetohydrodynamic (MHD) wave energy fluxes for late-type stars are calculated, using previously obtained formulae for the source functions for the generation of MHD waves in a stratified, but otherwise uniform, turbulent atmosphere; the magnetic fields in the wave generation region are assumed to be homogeneous. In contradiction to previous results, it is shown that in this uniform magnetic field case there is no significant increase in the efficiency of MHD wave generation, at least within the theory's limits of applicability. The major results are that the MHD energy fluxes calculated for late-type stars are less than those obtained for compressible modes in the magnetic field-free case, and that these MHD energy fluxes do not vary enough for a given spectral type to explain the observed range of UV and X-ray fluxes from such stars. It is therefore concluded that MHD waves in stellar atmospheres with homogeneous magnetic fields in the wave generation region cannot explain the observed stellar coronal emissions; if such MHD waves are responsible for a significant component of stellar coronal heating, then nonuniform fields within the generation region must be appealed to.

Musielak, Z. E.↗

Homogeneous crystal nucleation in binary metallic melts

A method for calculating the homogeneous crystal nucleation frequency in binary metallic melts is developed. The free energy of crystallization is derived from regular solution models for the liquid and solid and is used, together with model-based estimates of the interfacial tension, to calculate the nucleation frequency from the classical theory. The method can account for the composition dependence of the maximum undercooling observed in a number of experiments on small droplet dispersions. It can also be used to calculate the driving force for crystal growth and to obtain more precise estimates of the homogeneous crystal nucleation frequency in glass-forming alloys. This method, although approximate, is simple to apply, and requires only knowledge of the phase diagram and a few readily available thermodynamic quantities as input data.

Thompson, C. V.↗

Kinetics of Nucleation and Crystal Growth in Glass Forming Melts in Microgravity

This flight definition project has the specific objective of investigating the kinetics of nucleation and crystal growth in high temperature inorganic oxide, glass forming melts in microgravity. It is related to one1 of our previous NASA projects that was concerned with glass formation for high temperature containerless melts in microgravity. The previous work culminated in two experiments which were conducted aboard the space shuttle in 1983 and 1985 and which consisted of melting (at 1500 C) and cooling levitated 6 to 8 mm diameter spherical samples in a Single Axis Acoustic Levitator (SAAL) furnace. Compared to other types of materials, there have been relatively few experiments, 6 to 8, conducted on inorganic glasses in space. These experiments have been concerned with mass transport (alkali diffusion), containerless melting, critical cooling rate for glass formation, chemical homogeneity, fiber pulling, and crystallization of glass forming melts. One of the most important and consistent findings in all of these experiments has been that the glasses prepared in microgravity are more resistant to crystallization (better glass former) and more chemically homogeneous than equivalent glasses made on earth (1g). The chemical composition of the melt appears relatively unimportant since the same general results have been reported for oxide, fluoride and chalcogenide melts. These results for space-processed glasses have important implications, since glasses with a higher resistance to crystallization or higher chemical homogeneity than those attainable on earth can significantly advance applications in areas such as fiber optics communications, high power laser glasses, and other photonic devices where glasses are the key functional materials. The classical theories for nucleation and crystal growth for a glass or melt do not contain any parameter that is directly dependent upon the g-value, so it is not readily apparent why glasses prepared in microgravity should be more resistant to crystallization than equivalent glasses prepared on earth. Similarly, the gravity-driven convection in a fluid melt is believed to be the primary force field that is responsible for melt homogenization on earth. Thus, it is not obvious why a glass prepared in space, where gravity-driven convection is ideally absent, would be more chemically homogeneous than a glass identically prepared on earth. The primary objective of the present research is to obtain experimental data for the nucleation rate and crystal growth rate for a well characterized silicate melt (lithium disilicate) processed entirely in space (low gravity) and compare these rates with the nucleation and crystal growth rates for a similar glass prepared identically on earth (1g).

Day, Delbert E.↗

On the oscillation of the laterally heterogeneous earth, 1

The perturbative effects, as cause by lateral inhomogeneities in the earth structure and by Coriolis force, contaminate the originally toroidal and spheroidal earth's oscillations, making them of mixed type. For this reason, in order to make the computation of the perturbations more uniform and homogeneous, it was suggested that the earth's free oscillations be expanded into a series in terms of generalized harmonics familiar from the theory of angular momentum in quantum mechanics. Making use of Gibbsian symbolism and of some operators from the theory of angular momentum, explicit expressions were deduced for the perturbative terms in the differential equation of the earth's free oscillations. Decomposition of the strain tensor in terms of canonical vectors was also obtained.

Musen, P.↗

New versions of old flow visualization systems

A small, modified Mach-Zehnder interferometer placed in series with a much larger schlieren optical system spanning the test section is examined. In one arm of the interferometer, light from the schlieren is focused through a pinhole and recollimating lens to produce a reference beam which interferes with the remaining object beam from the other arm. The object and reference beams are separated only over a small interval following the test section, and differential vibrations are greatly reduced. Color schlieren has technical, as well as aesthetic, advantages over black-and-white schlieren. Since each color is associated with a specific amount of refraction, quantitative evaluation of certain refractive-index fields becomes possible using very simple equations derived from ray trace theory. Rainbow schlieren of an acetylene flame, and the evaluated refractive index distribution are shown. Root-mean-square refractive index fluctuations in homogeneous, isotropic turbulence were determined using the rainbow schlieren, since these fluctuations determine the root-mean-square refraction, which is indicated by the overall color of the image.

Howes, W. L.↗

Estimation of homogeneous nucleation flux via a kinetic model

The proposed kinetic model for condensation under homogeneous conditions, and the onset of unidirectional cluster growth in supersaturated gases, does not suffer from the conceptual flaws that characterize classical nucleation theory. When a full set of simultaneous rate equation is solved, a characteristic time emerges, for each cluster size, at which the production rate, and its rate of conversion to the next size (n + 1) are equal. Procedures for estimating the essential parameters are proposed; condensation fluxes J(kin) exp ss are evaluated. Since there are practical limits to the cluster size that can be incorporated in the set of simultaneous first-order differential equations, a code was developed for computing an approximate J(th) exp ss based on estimates of a 'constrained equilibrium' distribution, and identification of its minimum.

Wilcox, C. F.↗

Single and double point modeling of homogeneous turbulence

Investigations carried out for revisiting homogeneous turbulent flows in the presence of mean shear, rotation, or external compression are summarized. The simplest and most concise RDT (Rapid Distortion Theory) formulation, which includes a comprehensive linear stability analysis, is used for this purpose. Such a linear approach could be extended by a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) to two point closure in order to model non-linear interactions, especially when pure Coriolis effects are present. The results are discussed in connection with databases obtained by DNS (Direct Numerical Simulations), including previous Center for Turbulence Research (CTR) results and new calculations in progress. The main goal is to contribute to and significantly improve on the rational one-point closure models in progress at the CTR and at Ecole Centrale de Lyon (ECL).

Cambon, C.↗

Adaptation of Combustion Principles to Aircraft Propulsion: Basic Considerations in the Combustion of Hydrocarbon Fuels with Air - Volume I

The report summarizes source material on combustion for flight-propulsion engineers. First, several chapters review fundamental processes such as fuel-air mixture preparation, gas flow and mixing, flammability and ignition, flame propagation in both homogenous and heterogenous media, flame stabilization, combustion oscillations, and smoke and carbon formation. The practical significance and the relation of these processes to theory are presented. A second series of chapters describes the observed performance and design problems of engine combustors of the principal types. An attempt is made to interpret performance in terms of the fundamental processes and theories previously reviewed. Third, the design of high-speed combustion systems is discussed. Combustor design principles that can be established from basic considerations and from experience with actual combustors are described. Finally, future requirements for aircraft engine combustion systems are examined.

THEORIES - COMBUSTION↗

Time Evolution of Modeled Reynolds Stresses in Planar Homogeneous Flows

The analytic expression of the time evolution of the Reynolds stress anisotropy tensor in all planar homogeneous flows is obtained by exact integration of the modeled differential Reynolds stress equations. The procedure is based on results of tensor representation theory, is applicable for general pressure-strain correlation tensors, and can account for any additional turbulence anisotropy effects included in the closure. An explicit solution of the resulting system of scalar ordinary differential equations is obtained for the case of a linear pressure-strain correlation tensor. The properties of this solution are discussed, and the dynamic behavior of the Reynolds stresses is studied, including limit cycles and sensitivity to initial anisotropies.

Jongen, T.↗

High-Fidelity Micromechanics Model Developed for the Response of Multiphase Materials

A new high-fidelity micromechanics model has been developed under funding from the NASA Glenn Research Center for predicting the response of multiphase materials with arbitrary periodic microstructures. The model's analytical framework is based on the homogenization technique, but the method of solution for the local displacement and stress fields borrows concepts previously employed in constructing the higher order theory for functionally graded materials. The resulting closed-form macroscopic and microscopic constitutive equations, valid for both uniaxial and multiaxial loading of periodic materials with elastic and inelastic constitutive phases, can be incorporated into a structural analysis computer code. Consequently, this model now provides an alternative, accurate method.

Aboudi, Jacob↗

Multiscale Prediction of Yarn Pullout Failure Mode in Unreinforced Textile Fabrics

Unreinforced woven fabrics have been implemented in a variety of high performance applications, including body armor, deployable structures, and as the reinforcement material in composites. Multiscale modeling techniques have significantly improved the capabilities of simulation-based tools to capture fabric mechanics efficiently and accurately, but often lack in their prediction of failure and require pairing with finite element analysis (FEA) software, limiting their application to the design of ‘fit-for-purpose’ materials. NASA’s Multiscale Analysis Tool (NASMAT) is a standalone multiscale program that has been traditionally used in the analysis of reinforced composites materials. More recently, it has been amended to simulate unreinforced fabric behavior by allowing the geometric state of the tows to change with applied loading due to the lack of a reinforcement material, such as the matrix seen in composites. Previous work has shown the ability of NASMAT to capture nonlinear macroscale behavior by predicting geometric changes in the state of each subcell as a function of the applied loading and allowing each subcell in the analysis to rotate according to these predicted changes, as well as predict nonlinear behavior due to the fiber breakage failure mode. In this work, the capability of predicting the onset and propagation of failure in plain woven fabrics in NASMAT is presented for the yarn pullout failure mode, which occurs when a fabric is loaded at an off-axis angle relative to the warp of weft tow direction. Yarn pullout behavior is initiated by determining the applied load in which the shear resistance of the contact area between yarn families is overcome. When failure is initiated, yarn pullout is determined to have occurred when the applied displacement, calculated from global strain, exceeds the deformed position of a given contact points between yarn families, determined from pin-joint kinematics. Contact points where pullout has occurred contribute to a global damage parameter used to modify the homogenized stiffness of the fabric, resulting in nonlinear behavior observed at the macroscale. The off-axis loading behavior and yarn pullout failure theory have been developed and implemented into NASMAT such that users can simulate off-axis tensile behavior of fabrics in a single, standalone multiscale tool. Simulations are compared to uniaxial tensile tests at various off-axis angles to demonstrate the capability of the tool in its prediction of both the onset of failure at each off-axis angles and the stress-strain behavior as failure progresses.

Materials↗