Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Mechanics”

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 181 records · Page 10

Trinity by the Numbers: The Computing Effort that Made Trinity Possible

Utilizing rarely used source materials and identifying points of bias among more commonly used sources, this critical review provides a more complete representation of wartime Los Alamos computing operations and personnel, including the Laboratory’s typically underrepresented human computers and how they contributed to the success of the Trinity test. This paper also identifies how the Laboratory’s unusual wartime computing demands served as a formative experience among many Los Alamos personnel and consultants who contributed significantly to the development and use of mechanized computing at and beyond Los Alamos after the war.

79 ASTRONOMY AND ASTROPHYSICS↗

Metal matrix composite micromechanics: In-situ behavior influence on composite properties

Recent efforts in computational mechanics methods for simulating the nonlinear behavior of metal matrix composites have culminated in the implementation of the Metal Matrix Composite Analyzer (METCAN) computer code. In METCAN material nonlinearity is treated at the constituent (fiber, matrix, and interphase) level where the current material model describes a time-temperature-stress dependency of the constituent properties in a material behavior space. The composite properties are synthesized from the constituent instantaneous properties by virtue of composite micromechanics and macromechanics models. The behavior of metal matrix composites depends on fabrication process variables, in situ fiber and matrix properties, bonding between the fiber and matrix, and/or the properties of an interphase between the fiber and matrix. Specifically, the influence of in situ matrix strength and the interphase degradation on the unidirectional composite stress-strain behavior is examined. These types of studies provide insight into micromechanical behavior that may be helpful in resolving discrepancies between experimentally observed composite behavior and predicted response.

Murthy, P. L. N.↗

Improving Computational Efficiency of Mechanical Finite Element Method Simulations for PV Modules: Preprint

In this work we have elucidated the trade off between structural mechanics FEM model accuracy and computation time by employing lower fidelity viscoelastic models. Results indicate that computation time can easily be cut in half while only expecting a potential maximum error of 10 % by considering lower fidelity models. A novel approach to produce the Prony Series fit for viscoelastic characterization has also been presented. Employing this approach, we were able to achieve a further 10 % reduction in computation time without further sacrificing simulation accuracy.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗

An arbitrary grid CFD algorithm for configuration aerodynamics analysis. Volume 1: Theory and validations

This report documents the user input and output data requirements for the FEMNAS finite element Navier-Stokes code for real-gas simulations of external aerodynamics flowfields. This code was developed for the configuration aerodynamics branch of NASA ARC, under SBIR Phase 2 contract NAS2-124568 by Computational Mechanics Corporation (COMCO). This report is in two volumes. Volume 1 contains the theory for the derived finite element algorithm and describes the test cases used to validate the computer program described in the Volume 2 user guide.

Baker, A. J.↗

An arbitrary grid CFD algorithm for configuration aerodynamics analysis. Volume 2: FEMNAS user guide

This report documents the user input and output data requirements for the FEMNAS finite element Navier-Stokes code for real-gas simulations of external aerodynamics flowfields. This code was developed for the configuration aerodynamics branch of NASA ARC, under SBIR Phase 2 contract NAS2-124568 by Computational Mechanics Corporation (COMCO). This report is in two volumes. Volume 1 contains the theory for the derived finite element algorithm and describes the test cases used to validate the computer program described in the Volume 2 user guide.

Manhardt, Paul D.↗

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Prediction of properties of intraply hybrid composites

Equations based on the mixtures rule are presented for predicting the physical, thermal, hygral, and mechanical properties of unidirectional intraply hybrid composites (UIHC) from the corresponding properties of their constituent composites. Bounds were derived for uniaxial longitudinal strengths, tension, compression, and flexure of UIHC. The equations predict shear and flexural properties which agree with experimental data from UIHC. Use of these equations in a composites mechanics computer code predicted flexural moduli which agree with experimental data from various intraply hybrid angleplied laminates (IHAL). It is indicated, briefly, how these equations can be used in conjunction with composite mechanics and structural analysis during the analysis/design process.

Chamis, C. C.↗

Mechanics of intraply hybrid composites - Properties, analysis and design

A mechanics theory is developed for predicting the physical thermal, hygral and mechanical properties (including various strengths) of unidirectional intraply hybrid composites (UIHC) based on unidirectional properties of the constituent composites. Procedures are described which can use this theory in conjunction with composite mechanics computer codes and general purpose structural analysis finite element programs for the analysis/design of structural components made from intraply hybrid angleplied laminates (IHAL). Comparisons with limited data show that this theory predicts mechanical properties of UIHC and flexural stiffnesses of IHAL which are in good agreement with experimental data. The theory developed herein makes it possible to design and optimize structural components from IHAL based on a large class of available constituent fibers.

Chamis, C. C.↗

Computational simulation of surface waviness in graphite/epoxy woven composites due to initial curing

Several models simulating plain weave, graphite/epoxy woven composites are presented, along with the effects that the simultaneous application of pressure and thermal loads have on their surfaces. The surface effects created by moisture absorption are also examined. The computational simulation consisted of using a two-dimensional finite element model for the composite. The properties of the finite element (FE) model are calculated by using the in-house composite mechanics computer code ICAN (Integrated Composite ANalyzer). MSC/NASTRAN is used for the FE analysis which yields the composite's top surface normalized displacements. These results demonstrate the importance of parameters such as the cure temperature (T sub o) and the resin content in the curing process of polymer-matrix composites. The modification of these parameters will help tailor the composite system to the desired requirements and applications.

Sanfeliz, Jose G.↗

Computational simulation of surface waviness in graphite/epoxy woven composites due to initial curing

Several models simulating plain weave, graphite/epoxy woven composites are presented, along with the effects that the simultaneous application of pressure and thermal loads have on their surfaces. The surface effects created by moisture absorption are also examined. The computational simulation consisted of using a two-dimensional finite element model for the composite. The properties of the finite element (FE) model are calculated by using the in-house composite mechanics computer code ICAN (Integrated Composite ANalyzer). MSC/NASTRAN is used for the FE analysis which yields the composite's top surface normalized displacements. These results demonstrate the importance of parameters such as the cure temperature (T sub o) and the resin content in the curing process of polymer-matrix composites. The modification of these parameters will help tailor the composite system to the desired requirements and applications.

Sanfeliz, Jose G.↗

Uncertainty Reduction With Multi-Model Monte Carlo for Crystal Plasticity Simulations of Additively Manufactured Metals

In this work, multi-model Monte Carlo estimators are developed to reduce uncertainty in quantities of interest (QoIs) extracted from crystal plasticity simulations of additively manufactured (AM) metals. A significant concern in AM parts is uncertainty in mechanical properties caused in part by complex microstructures that arise from the AM process. Quantifying uncertainty in microstructure-sensitive behavior using experiments alone is costly, especially when mechanical allowables must be established. Quantitative relationships among microstructure, micromechanical metrics like slip accumulation, crack initiation, and failure are also difficult to capture with limited experiments. Crystal plasticity material models instead enable computational prediction of micromechanical stress and strain fields given a discretized microstructure. However, high-fidelity finely discretized crystal plasticity simulations are computationally expensive, while lower-fidelity models are less accurate and generally biased, making uncertainty quantification and reduction computationally difficult as well. Multi-model Monte Carlo methods leverage correlations between high- and low-fidelity models to produce unbiased estimators for QoIs with reduced uncertainty relative to standard Monte Carlo. Crystal plasticity QoIs considered in this work include yield strength and the mean and extreme values of micromechanical fields that are relevant to crack initiation. Multi-model Monte Carlo estimators are developed for each individual QoI and several groups of QoIs. The results of this work establish relationships among model correlations, sample allocation, and uncertainty reduction for different combinations of QoIs and demonstrate a trend of less uncertainty reduction as QoIs become more sensitive to local microstructure. Limitations from using pilot samples to estimate model covariances and train low-fidelity models are also addressed. The uncertainty reduction achieved by multi-model Monte Carlo is an important step toward using computational mechanics models to predict microstructure-sensitive crack initiation and failure in AM parts.

uncertainty quantification↗

METCAN: The metal matrix composite analyzer

Metal matrix composites (MMC) are the subject of intensive study and are receiving serious consideration for critical structural applications in advanced aerospace systems. MMC structural analysis and design methodologies are studied. Predicting the mechanical and thermal behavior and the structural response of components fabricated from MMC requires the use of a variety of mathematical models. These models relate stresses to applied forces, stress intensities at the tips of cracks to nominal stresses, buckling resistance to applied force, or vibration response to excitation forces. The extensive research in computational mechanics methods for predicting the nonlinear behavior of MMC are described. This research has culminated in the development of the METCAN (METal Matrix Composite ANalyzer) computer code.

Hopkins, Dale A.↗

A Lanczos algorithm for vibration, suckling and termal analysis

This paper reviews an eigensolver algorithm based on the Lanczos Method for vibration, buckling and thermal analysis. The original code was written for inclusion in the Computational Mechanics Testbed (COMET), a general purpose finite element code. A portable version of the Lanczos code that is optimized for high-performance supercomputers has been developed. Special features of the algorithm include the capability to compute rigid body modes, thermal modes and Lanczos vectors that are derived from the applied load vector. The latter is necessary when using the Lanczos vectors as reduced-basis vectors in transient structural response and transient heat conduction calculations. The modularity of the code allows the user the option of including the most up-to-date utilities, such as the equation solver best suited for the application. The algorithm is discussed in detail and results of several applications are presented. Timing results for a vibration application indicate that the Lanczos algorithm is twenty times faster than the subspace iteration method which has been extensively used in the past.

Bostic, Susan W.↗

High-tip-speed fiber composite compressor blades: Vibration and strength analysis

An analytical procedure is described which couples composite mechanics computer codes with NASTRAN. This procedure was used to perform a detailed analysis of a high-tip-speed fiber composite compressor fan blade. The results indicate that the various vibration modes of this blade are highly coupled. Mechanical load ply stresses are well below the corresponding room temperature strengths. Lamination residual stresses are likely to cause transply cracks and interply delamination. Transply cracks and relaxation of root fixity decrease the vibrational frequencies whereas centrifugal stiffening increases them. Comparisons of results for various parameters are presented in tabular and graphical form.

Chamis, C. C.↗

High-tip-speed fiber composite compressor blades - Vibration and strength analyses

An analytical procedure is described which couples composite mechanics computer codes with NASTRAN. This procedure was used to perform a detailed analysis of a high-tip-speed fiber composite compressor fan blade. The results indicate that the various vibration modes of this blade are highly coupled. Mechanical load ply stresses are well below the corresponding room temperature strengths. Lamination residual stresses are likely to cause transply cracks and interply delamination. Transply cracks and relaxation of root fixity decrease the vibrational frequencies whereas centrifugal stiffening increases them. Comparisons of results for various parameters are presented in tabular and graphical form.

Chamis, C. C.↗

Nonrecursive formulations of multibody dynamics and concurrent multiprocessing

Since the late 1980's, research in recursive formulations of multibody dynamics has flourished. Historically, much of this research can be traced to applications of low dimensionality in mechanism and vehicle dynamics. Indeed, there is little doubt that recursive order N methods are the method of choice for this class of systems. This approach has the advantage that a minimal number of coordinates are utilized, parallelism can be induced for certain system topologies, and the method is of order N computational cost for systems of N rigid bodies. Despite the fact that many authors have dismissed redundant coordinate formulations as being of order N(exp 3), and hence less attractive than recursive formulations, we present recent research that demonstrates that at least three distinct classes of redundant, nonrecursive multibody formulations consistently achieve order N computational cost for systems of rigid and/or flexible bodies. These formulations are as follows: (1) the preconditioned range space formulation; (2) penalty methods; and (3) augmented Lagrangian methods for nonlinear multibody dynamics. The first method can be traced to its foundation in equality constrained quadratic optimization, while the last two methods have been studied extensively in the context of coercive variational boundary value problems in computational mechanics. Until recently, however, they have not been investigated in the context of multibody simulation, and present theoretical questions unique to nonlinear dynamics. All of these nonrecursive methods have additional advantages with respect to recursive order N methods: (1) the formalisms retain the highly desirable order N computational cost; (2) the techniques are amenable to concurrent simulation strategies; (3) the approaches do not depend upon system topology to induce concurrency; and (4) the methods can be derived to balance the computational load automatically on concurrent multiprocessors. In addition to the presentation of the fundamental formulations, this paper presents new theoretical results regarding the rate of convergence of order N constraint stabilization schemes associated with the newly introduced class of methods.

Kurdila, Andrew J.↗

Bayesian Symbolic Regression: Addressing Challenges in Estimating Fractional Bayes Factors and Application to Fatigue Crack Growth Modeling

This research pioneers advancements in computational mechanics by integrating Bayesian-based uncertainty quantification into symbolic regression, specifically focusing on the critical task of accurately estimating the fractional Bayes factor for selecting arbitrary equations. In our exploration, we rigorously study two prominent methods—sequential Monte Carlo and the Laplace approximation—employed for computing the fractional Bayes factor. Our findings underscore the limitations of the Laplace approximation, revealing its diminished accuracy in nonlinear and multimodal scenarios. Specifically, the Laplace approximation is shown to underpredict fractional Bayes factor on a wide set of equations associated with a symbolic regression benchmark. This comparative analysis sheds light on the nuanced performance of these techniques, guiding researchers toward more informed choices in uncertainty quantification within symbolic regression. Furthermore, we showcase the practical utility of these enhanced symbolic regression tools through their application to a real-world problem in fatigue crack growth modeling, emphasizing their efficacy in capturing the complexities of mechanical systems.

Geoffrey Bomarito↗