Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

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↗
NASA NTRS1996

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.↗
DOE OSTI2015

Enabling Structured Exploration of Workflow Performance Variability in Extreme-Scale Environments

Workflows are taking an Workflows are taking an increasingly important role in orchestrating complex scientific processes in extreme scale and highly heterogeneous environments. However, to date we cannot reliably predict, understand, and optimize workflow performance. Sources of performance variability and in particular the interdependencies of workflow design, execution environment and system architecture are not well understood. While there is a rich portfolio of tools for performance analysis, modeling and prediction for single applications in homogenous computing environments, these are not applicable to workflows, due to the number and heterogeneity of the involved workflow and system components and their strong interdependencies. In this paper, we investigate workflow performance goals and identify factors that could have a relevant impact. Based on our analysis, we propose a new workflow performance provenance ontology, the Open Provenance Model-based WorkFlow Performance Provenance, or OPM-WFPP, that will enable the empirical study of workflow performance characteristics and variability including complex source attribution.

Enabling structured workflows↗
NASA NTRS1990

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.↗
DOE OSTI2020

A numerical-homogenization based phase-field fracture modeling of linear elastic heterogeneous porous media

Most porous media, such as geomaterials and biomaterials are highly heterogeneous in nature, and they contain large variations of microscopic pore structures, such as pore sizes, pore distribution, and pore shapes. The oscillation of microscopic structures is a substantial challenge in theoretical characterization and is usually ignored in continuous modeling. However, mechanical behavior of porous media such as deformation and failure, are essentially impacted by the microscopic heterogeneity which needs to be considered in modeling a porous media. Here, this research proposes a numerical modeling framework with a capability to investigate the effect of microscopic heterogeneity on the macroscopic fracture behavior in porous media by using a numerical homogenization technique, combined with the phase-field fracture modeling method. This numerical modeling strategy computes a homogenized elasticity tensor based on microscopic heterogeneous pore structures heterogenous porous domain by solving boundary value problems at microscopic domain. The strain energy and subsequent propagation of macroscopic fractures will be updated using homogenized stiffness information. Using this numerical scheme, the microscopic pore structure’s impact on the fracture behavior through the homogenized elastic tensor will be taken into account. This multiscale technique is benchmarked against classical problems. Finally, the results highlight the importance of the underlying pore structure and reveal that both fracture strength and propagation path can be influenced by the microscopic heterogeneity.

36 MATERIALS SCIENCE↗
NASA NTRS1992

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

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

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.↗
DOE OSTI2020

Computationally efficient evaluation of optimum homogeneous charge compression ignition operating range with accelerated multizone engine cycle simulation

The very intensive calculations necessary to define a performance map requiring evaluation of over a hundred individual operating points can be efficiently conducted with accelerated multizone for engine cycle simulation, leading to a definition of regions of acceptable and optimum homogeneous charge compression ignition operation. Accelerated multizone for engine cycle simulation has the virtue of enabling accurate evaluation of many operating conditions based on thermal stratification data from a single fluid mechanics run at motored conditions. This is possible because thermal stratification is more sensitive to engine geometry than to operating conditions. In this article, accuracy of accelerated multizone for engine cycle simulation is demonstrated by comparison with experimental data for iso-octane homogeneous charge compression ignition operation over a broad range of lean equivalence ratios (0.14–0.28). The validated accelerated multizone for engine cycle simulation model is then applied to generating a performance map for an engine controlled by appropriately adjusting equivalence ratio and internal exhaust gas recirculation. Regions of acceptable and optimum combustion are identified. It is finally demonstrated that while indicated mean effective pressure remains low for optimum homogeneous charge compression ignition operation (1–4 bar), this is sufficient for a large fraction of typical driving in light-duty vehicles. Much driving including idle can therefore be done in homogeneous charge compression ignition mode at high efficiency and low (essentially zero) NOx and particulate matter emissions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2018

A homogeneous model for monotone mixed horizontal linear complementarity problems

In this paper, we propose a homogeneous model for the class of mixed horizontal linear complementarity problems. The proposed homogeneous model is always solvable and provides the solution of the original problem if it exists, or a certificate of infeasibility otherwise. Our formulation preserves the sparsity of the original formulation and does not reduce to the homogeneous model of the equivalent standard linear complementarity problem. We study the properties of the model and show that interior-point methods can be used efficiently for the numerical solutions of the homogeneous problem. Finally, numerical experiments show convincingly that it is more efficient to use the proposed homogeneous model for the mixed horizontal linear complementarity problem than to use known homogeneous models for the equivalent standard linear complementarity problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Microstructure and Mechanical Behavior of High-Entropy Alloys

High entropy alloys (HEA) have generated interest in recent years due to their unique positioning within the alloy world. By incorporating a number of elements in high proportion, usually of equal atomic percent, they have high configurational entropy, and thus, they hold the promise of interesting and useful properties such as enhanced strength and alloy stability. The present study investigates the mechanical behavior, fracture characteristics and microstructure of two singlephase FCC high-entropy alloys (HEA's) CoCrFeNi and CoCrFeNiMn with some detailed attention given to melting, homogenization and thermo-mechanical processing. Ingots approaching 8 kg in mass were made by vacuum induction melting (VIM) to avoid the extrinsic factors inherent to small scale laboratory button samples. A computationally based homogenization heat treatment was given to both alloys in order to eliminate any solidification segregation. The alloys were then fabricated in the usual way (forging, followed by hot rolling) with typical thermo-mechanical processing parameters employed. Transmission electron microscopy (TEM) was subsequently used to assess the single phase nature of the alloys prior to mechanical testing. Tensile specimens (ASTM E8) were prepared with tensile mechanical properties obtained from room temperature through 800°C. Material from the gage section of selected tensile specimens was extracted to document room and elevated temperature deformation within the HEA's. Fracture surfaces were also examined to note fracture failure modes. The tensile behavior and selected tensile properties were compared with results in the literature for similar alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Porting hypre to heterogeneous computer architectures: Strategies and experiences

We report that linear systems are occurring in many applications, and solving them can take a large amount of the total simulation time. The high performance library hypre provides a variety of interfaces and linear solvers, including various multigrid methods, that have achieved good scalability on a variety of homogeneous parallel computer architectures. Heterogeneous architectures with nodes that have both CPUs and accelerators provide new challenges, since they require more fine-grained parallelism and reduced data movement between different memories on a single node as well as across nodes. We will discuss our experiences and strategies to port hypre to heterogeneous computers with accelerators, including the design of a new memory model, the use of abstractions, the BoxLoop macros in the structured and semi-structured interfaces, and the restructuring of algebraic multigrid (AMG) into modular components. We present numerical experiments comparing CPU and GPU performance for several test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Analyzing photovoltaic module mechanics using composite plate theories and finite element solutions

Deflection and stress calculated from an experimentally validated, high-fidelity finite element model (FEM) of a photovoltaic module experiencing mechanical load was compared to results from a simplified FEM treating the module laminate as a homogenized composite using a rule of mixtures approach, and further compared to analytical calculations treating the module as a Kirchoff-Love flat plate. The goal of this study was to determine the error incurred by analyzing module mechanics with varying levels of simplification, since resolving the aspect ratios of a module is computationally expensive. Homogenized FEMs were found to underpredict peak deflection under a 1.0 kPa load by between 13 and 19% for lower and upper bound application of the rule of mixtures. However, module shape was captured, implying that a useful replication of a resolved model could be achieved with a reduced, calibrated material stiffness. Homogenized stress results captured glass layer tensile stress components to within 46 to 52% at a sample location of interest, though agreement was poor through the remainder of the laminate due to the lack of material resolution. For plate theory, deflection was overpredicted by 45 to 67% for upper and lower bound homogenizations, and frame-adjacent module shapes were not adequately replicated. Stress results mirrored FEM trends but magnitudes were not well correlated to resolved model values. These results support the use of homogenized laminate models for module shape derivation, though resolved models remain necessary for stress analyses. The accuracy of plate theory was found to be inadequate for most applications.

14 SOLAR ENERGY↗
DOE OSTI2014

A spectral-Lagrangian Boltzmann solver for a multi-energy level gas

In this paper a spectral-Lagrangian method is proposed for the full, non-linear Boltzmann equation for a multi-energy level gas typical of a hypersonic re-entry flow. Internal energy levels are treated as separate species and inelastic collisions (leading to internal energy excitation and relaxation) are accounted for. The formulation developed can also be used for the case of a gas mixture made of monatomic gases without internal energy (where only elastic collisions occur). The advantage of the spectral-Lagrangian method lies in the generality of the algorithm in use for the evaluation of the elastic and inelastic collision operators, as well as the conservation of mass, momentum and energy during collisions. The latter is realized through the solution of constrained optimization problems. The computational procedure is based on the Fourier transform of the partial elastic and inelastic collision operators and exploits the fact that these can be written as weighted convolutions in Fourier space with no restriction on the cross-section model. The feasibility of the proposed approach is demonstrated through numerical examples for both space homogeneous and in-homogeneous problems. Computational results are compared with those obtained by means of the DSMC method in order to assess the accuracy of the proposed spectral-Lagrangian method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

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