Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid 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 253 records · Page 14

The Considere Condition and Rapid Stretching of Linear and Branched Polymer Melts

We analyze the onset of "necking" and subsequent filament failure during the transient uniaxial elongation of viscoelastic fluid samples in extensional rheometers. In the limit of rapid elongation (such that no molecular relaxation occurs), the external work applied is all stored elastically and the Considere criterion originally developed in solid mechanics can be used to quantitatively predict the critical Hencky strain to failure. By comparing the predictions of the Doi-Edwards model for linear homopolymer melts with those of the "Pom-Pom" model for prototypical branched melts we show that the critical strain to failure in rapid elongation of a rubbery material is intimately linked to the molecular topology of the chain, especially the degree of chain branching. The onset of necking instability is monotonically shifted to larger Hencky strains as the number of branches is increased. Numerical computations at finite Deborah numbers also show that there is an optimal range of deformation rates over which homogeneous extensions can be maintained to large strain. We also consider other rapid homogeneous stretching deformations, such as biaxial and planar stretching, and show that the degree of stabilization afforded by inclusion of material with long-chain branching is a sensitive function of the imposed mode of deformation.

McKinley, Gareth H.↗

Equivalent-Continuum Modeling of Nano-Structured Materials

A method has been developed for modeling structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with an equivalent-continuum model. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As an important example with direct application to the development and characterization of single-walled carbon nanotubes, the model has been applied to determine the effective continuum geometry of a graphene sheet. A representative volume element of the equivalent-continuum model has been developed with an effective thickness. This effective thickness has been shown to be similar to, but slightly smaller than, the interatomic spacing of graphite.

Odegard, Gregory M.↗

Enhancing the Remote Variable Operations in NPSS/CCDK

Many scientific applications in aerodynamics and solid mechanics are written in Fortran. Refitting these legacy Fortran codes with distributed objects can increase the code reusability. The remote variable scheme provided in NPSS/CCDK helps programmers easily migrate the Fortran codes towards a client-server platform. This scheme gives the client the capability of accessing the variables at the server site. In this paper, we review and enhance the remote variable scheme by using the operator overloading features in C++. The enhancement enables NPSS programmers to use remote variables in much the same way as traditional variables. The remote variable scheme adopts the lazy update approach and the prefetch method. The design strategies and implementation techniques are described in details. Preliminary performance evaluation shows that communication overhead can be greatly reduced.

Sang, Janche↗

High-End Computing for Incompressible Flows

The objective of the First MIT Conference on Computational Fluid and Solid Mechanics (June 12-14, 2001) is to bring together industry and academia (and government) to nurture the next generation in computational mechanics. The objective of the current talk, 'High-End Computing for Incompressible Flows', is to discuss some of the current issues in large scale computing for mission-oriented tasks.

Kwak, Dochan↗

End Effects and Load Diffusion in Composite Structures

The research carried out here builds on our previous NASA supported research on the general topic of edge effects and load diffusion in composite structures. Further fundamental solid mechanics studies were carried out to provide a basis for assessing the complicated modeling necessary for large scale structures used by NASA. An understanding of the fundamental mechanisms of load diffusion in composite subcomponents is essential in developing primary composite structures. Specific problems recently considered were focussed on end effects in sandwich structures and for functionally graded materials. Both linear and nonlinear (geometric and material) problems have been addressed. Our goal is the development of readily applicable design formulas for the decay lengths in terms of non-dimensional material and geometric parameters. Analytical models of load diffusion behavior are extremely valuable in building an intuitive base for developing refined modeling strategies and assessing results from finite element analyses. The decay behavior of stresses and other field quantities provides a significant aid towards this process. The analysis is also amenable to parameter study with a large parameter space and should be useful in structural tailoring studies.

Horgan, Cornelius O.↗

Fidelity of the Integrated Force Method Solution

The theory of strain compatibility of the solid mechanics discipline was incomplete since St. Venant's 'strain formulation' in 1876. We have addressed the compatibility condition both in the continuum and the discrete system. This has lead to the formulation of the Integrated Force Method. A dual Integrated Force Method with displacement as the primal variable has also been formulated. A modest finite element code (IFM/Analyzers) based on the IFM theory has been developed. For a set of standard test problems the IFM results were compared with the stiffness method solutions and the MSC/Nastran code. For the problems IFM outperformed the existing methods. Superior IFM performance is attributed to simultaneous compliance of equilibrium equation and compatibility condition. MSC/Nastran organization expressed reluctance to accept the high fidelity IFM solutions. This report discusses the solutions to the examples. No inaccuracy was detected in the IFM solutions. A stiffness method code with a small programming effort can be improved to reap the many IFM benefits when implemented with the IFMD elements. Dr. Halford conducted a peer-review on the Integrated Force Method. Reviewers' response is included.

Hopkins, Dale↗

Equivalent-Continuum Modeling With Application to Carbon Nanotubes

A method has been proposed for developing structure-property relationships of nano-structured materials. This method serves as a link between computational chemistry and solid mechanics by substituting discrete molecular structures with equivalent-continuum models. It has been shown that this substitution may be accomplished by equating the vibrational potential energy of a nano-structured material with the strain energy of representative truss and continuum models. As important examples with direct application to the development and characterization of single-walled carbon nanotubes and the design of nanotube-based devices, the modeling technique has been applied to determine the effective-continuum geometry and bending rigidity of a graphene sheet. A representative volume element of the chemical structure of graphene has been substituted with equivalent-truss and equivalent continuum models. As a result, an effective thickness of the continuum model has been determined. This effective thickness has been shown to be significantly larger than the interatomic spacing of graphite. The effective thickness has been shown to be significantly larger than the inter-planar spacing of graphite. The effective bending rigidity of the equivalent-continuum model of a graphene sheet was determined by equating the vibrational potential energy of the molecular model of a graphene sheet subjected to cylindrical bending with the strain energy of an equivalent continuum plate subjected to cylindrical bending.

Odegard, Gregory M.↗

Two- and Three-Dimensional Numerical Experiments Representing Two Limiting Cases of an In-Line Pair of Finger Seal Components

The work presented here concerns the numerical development and simulation of the flow, pressure patterns and motion of a pair of fingers arranged behind each other and axially aligned in-line. The fingers represent the basic elemental component of a Finger Seal (FS) and form a tight seal around the rotor. Yet their flexibility allows compliance with rotor motion and in a passive-adaptive mode complies also with the hydrodynamic forces induced by the flowing fluid. While the paper does not treat the actual staggered configuration of a finger seal, the inline arrangement represents a first step towards that final goal. The numerical 2-D (axial-radial) and 3-D results presented herein were obtained using a commercial package (CFD-ACE+). Both models use an integrated numerical approach, which couples the hydrodynamic fluid model (Navier-Stokes based) to the solid mechanics code that models the compliance of the fingers.

Braun, M. J.↗

On the Formulation of Weakly Singular Displacement/Traction Integral Equations; and Their Solution by the MLPG Method

In this paper, a very simple method is used to derive the weakly singular traction boundary integral equation based on the integral relationships for displacement gradients. The concept of the MLPG method is employed to solve the integral equations, especially those arising in solid mechanics. A moving Least Squares (MLS) interpolation is selected to approximate the trial functions in this paper. Five boundary integral Solution methods are introduced: direct solution method; displacement boundary-value problem; traction boundary-value problem; mixed boundary-value problem; and boundary variational principle. Based on the local weak form of the BIE, four different nodal-based local test functions are selected, leading to four different MLPG methods for each BIE solution method. These methods combine the advantages of the MLPG method and the boundary element method.

Atluri, Satya N.↗

Aviation Safety Modeling and Simulation (ASMM) Propulsion Fleet Modeling: A Tool for Semi-Automatic Construction of CORBA-based Applications from Legacy Fortran Programs

Within NASA's Aviation Safety Program, NASA GRC participates in the Modeling and Simulation Project called ASMM. NASA GRC s focus is to characterize the propulsion systems performance from a fleet management and maintenance perspective by modeling and through simulation predict the characteristics of two classes of commercial engines (CFM56 and GE90). In prior years, the High Performance Computing and Communication (HPCC) program funded, NASA Glenn in developing a large scale, detailed simulations for the analysis and design of aircraft engines called the Numerical Propulsion System Simulation (NPSS). Three major aspects of this modeling included the integration of different engine components, coupling of multiple disciplines, and engine component zooming at appropriate level fidelity, require relatively tight coupling of different analysis codes. Most of these codes in aerodynamics and solid mechanics are written in Fortran. Refitting these legacy Fortran codes with distributed objects can increase these codes reusability. Aviation Safety s modeling and simulation use in characterizing fleet management has similar needs. The modeling and simulation of these propulsion systems use existing Fortran and C codes that are instrumental in determining the performance of the fleet. The research centers on building a CORBA-based development environment for programmers to easily wrap and couple legacy Fortran codes. This environment consists of a C++ wrapper library to hide the details of CORBA and an efficient remote variable scheme to facilitate data exchange between the client and the server model. Additionally, a Web Service model should also be constructed for evaluation of this technology s use over the next two- three years.

Sang, Janche↗

Applications of Automation Methods for Nonlinear Fracture Test Analysis

Using automated and standardized computer tools to calculate the pertinent test result values has several advantages such as: 1. allowing high-fidelity solutions to complex nonlinear phenomena that would be impractical to express in written equation form, 2. eliminating errors associated with the interpretation and programing of analysis procedures from the text of test standards, 3. lessening the need for expertise in the areas of solid mechanics, fracture mechanics, numerical methods, and/or finite element modeling, to achieve sound results, 4. and providing one computer tool and/or one set of solutions for all users for a more "standardized" answer. In summary, this approach allows a non-expert with rudimentary training to get the best practical solution based on the latest understanding with minimum difficulty.Other existing ASTM standards that cover complicated phenomena use standard computer programs: 1. ASTM C1340/C1340M-10- Standard Practice for Estimation of Heat Gain or Loss Through Ceilings Under Attics Containing Radiant Barriers by Use of a Computer Program 2. ASTM F 2815 - Standard Practice for Chemical Permeation through Protective Clothing Materials: Testing Data Analysis by Use of a Computer Program 3. ASTM E2807 - Standard Specification for 3D Imaging Data Exchange, Version 1.0 The verification, validation, and round-robin processes required of a computer tool closely parallel the methods that are used to ensure the solution validity for equations included in test standard. The use of automated analysis tools allows the creation and practical implementation of advanced fracture mechanics test standards that capture the physics of a nonlinear fracture mechanics problem without adding undue burden or expense to the user. The presented approach forms a bridge between the equation-based fracture testing standards of today and the next generation of standards solving complex problems through analysis automation.

Allen, Phillip A.↗

Advanced Computational Modeling Approaches for Shock Response Prediction

Motivation: (1) The activation of pyroshock devices such as explosives, separation nuts, pin-pullers, etc. produces high frequency transient structural response, typically from few tens of Hz to several hundreds of kHz. (2) Lack of reliable analytical tools makes the prediction of appropriate design and qualification test levels a challenge. (3) In the past few decades, several attempts have been made to develop methodologies that predict the structural responses to shock environments. (4) Currently, there is no validated approach that is viable to predict shock environments overt the full frequency range (i.e., ~100 Hz to 10 kHz). Scope: (1) Model, analyze, and interpret space structural systems with complex interfaces and discontinuities, subjected to shock loads. (2) Assess the viability of a suite of numerical tools to simulate transient, non-linear solid mechanics and structural dynamics problems, such as shock wave propagation.

Derkevorkian, Armen↗

Order-Reduced Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique for accurately simulating nonlinear composite material behavior. The HFGMC uses a higher-order approximation for the subcell displacement field that allows for a more accurate determination of the subcell stressstrain fields at the cost of some computational efficiency. In order to reduce computational costs associated with the solution of the ensuing system of simultaneous equations, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents was reduced in size through the use of a Petrov-Galerkin-based Proper Orthogonal Decomposition order-reduction scheme. A number of cases were presented that address the computational feasibility of using order-reduction techniques to solve solid mechanics problems involving complex microstructures.

order-reduction↗

Material Response Modeling of Ablative Thermal Protection Systems using PATO

Developing thermal protection systems (TPS) for future space vehicles involves an extensive design and test cycle. A key component of the cycle is determining margins for a safe design that rely on confidence in the performance of the TPS. Predicting the complicated multiphysics phenomena that occur during atmospheric entry requires high-fidelity modeling tools. To this end, the Porous-material Analysis Toolbox based on OpenFOAM (PATO) has been developed. PATO is an open-source software for computing material response of reactive porous materials submitted to high-temperature environments. Recent applications of PATO include computation of the full heatshield 3D material response from the Mars Science Laboratory atmospheric entry and ablation during arc jet testing. Current efforts are underway to loosely couple PATO with other discipline specialized codes including hypersonic computational fluid dynamics (CFD) to assess the effects of pyrolysis-gas blowing into the boundary layer and computational solid mechanics to address modeling of mechanical erosion. Surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation are also being added.

Thermal Protection Systems↗

Material Response Modeling of Ablative Thermal Protection Systems using PATO

Developing thermal protection systems (TPS) for future space vehicles involves an extensive design and test cycle. A key component of the cycle is determining margins for a safe design that rely on confidence in the performance of the TPS. Predicting the complicated multiphysics phenomena that occur during atmospheric entry requires high-fidelity modeling tools. To this end, the Porous-material Analysis Toolbox based on OpenFOAM (PATO) has been developed. PATO is an open-source software for computing material response of reactive porous materials submitted to high-temperature environments. Recent applications of PATO include computation of the full heatshield 3D material response from the Mars Science Laboratory atmospheric entry and ablation during arc jet testing. Current efforts are underway to loosely couple PATO with other discipline specialized codes including hypersonic computational fluid dynamics (CFD) to assess the effects of pyrolysis-gas blowing into the boundary layer and computational solid mechanics to address modeling of mechanical erosion. Surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation are also being added.

Thermal Protection Systems↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Benchmarking and Exploring Parameter Space of the 2-Phase Bubble Tracking Model for Liquid Mercury Target Simulation

High intensity proton pulses strike the Spallation Neutron Source (SNS)’s mercury target to provide bright neutron beams. These strikes deposit extensive energy into the mercury and its steel vessel. Prediction of the resultant loading on the target is difficult when helium gas is intentionally injected into the mercury to reduce the loading and to mitigate the pitting damage on the vessel. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. We present a study comparing the measured target strains in the SNS target station with the simulation results of the solid mechanics simulation framework. We investigate a wide range of various physical model parameters, including the number of bubble families, bubble size distribution, viscosity, surface tension, etc. to understand their impact on simulation accuracy. Our initial findings reveal that using 8-10 bubble families in the model renders a simulation strain envelope that covers the experimental ones. Further optimization studies are planned to predict the strain response more accurately.

Lin, Lianshan↗