SEARCH · Search NASA

Results for “Discrete element models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Weibull-Based Stochastic Simulation of Mud-Crack Damage Formation in an Environmental Barrier Coating

The integrated Finite Element Analysis Micromechanics Analysis Code/Ceramics Analysis and Reliability Evaluation of Structures (FEAMAC/CARES) program was used to simulate the formation of mudflat-cracks from thermomechanical loading on a multi-layered Environmental Barrier Coating (EBC) system deposited on a ceramic substrate. FEAMAC/CARES combines MAC/GMC (Micromechanics Analysis Code/Generalized Method of Cells) multiscale composite micromechanics code with CARES/Life probabilistic multiaxial failure criteria code and Abaqusfinite element analysis. In this work, step function elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio finite elements enabled the depiction of crack boundaries and formation of mudflat patterned damage. Demonstrated examples include finite element models of button-sized disk-shaped 3-D specimen, and a 2-D model of through-the-thickness cross-section. All models were subjected to a progressive cool down from 1300 degrees Centigrade to room temperature. Mudflat crack damage in the coating system resulted from the buildup of residual tensile stresses between the individual material constituents from thermal expansion mismatch. A 2-parameter Weibull distribution characterized the coating layer stochastic strength response and the effect of the Weibull modulus on the formation of damage was studied here.

Residual Tensile Stress↗

NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗

DOE OSTI2023

Modeling fission dynamics at the barrier in a discrete-basis formalism

A configuration-interaction model is presented for the barrier region of induced fission. The configuration space is composed of seniority-zero configurations constructed from self-consistent mean-field wave functions. The Hamiltonian matrix elements between configurations include diabatic and pairing interactions between particles. Other aspects of the Hamiltonian are treated statistically, guided by phenomenological input of compound-nucleus transmission coefficients. In this exploratory study the configuration space is restricted to neutron excitations only. A key observable calculated in the model is the fission-to-capture branching ratio. We find that both pairing and diabatic interactions are important for achieving large branching to the fission channels. In accordance with the transition-state theory of fission, the calculated branching ratio is found to be quite insensitive to the fission decay widths of the pre-scission configurations. Furthermore, the barrier-top dynamics appear to be quite different from transition-state theory in that the transport is distributed over many excited configurations at the barrier top.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting, during flight, the residual strength of aircraft structures that sustain discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. Two ductile fracture simulations are presented to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data does, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗

NASA NTRS1998

Buckling and Damage Resistance of Transversely-Loaded Composite Shells

Experimental and numerical work was conducted to better understand composite shell response to transverse loadings which simulate damage-causing impact events. The quasi-static, centered, transverse loading response of laminated graphite/epoxy shells in a [+/-45(sub n)/O(sub n)](sub s) layup having geometric characteristics of a commercial fuselage are studied. The singly-curved composite shell structures are hinged along the straight circumferential edges and are either free or simply supported along the curved axial edges. Key components of the shell response are response instabilities due to limit-point and/or bifurcation buckling. Experimentally, deflection-controlled shell response is characterized via load-deflection data, deformation-shape evolutions, and the resulting damage state. Finite element models are used to study the kinematically nonlinear shell response, including bifurcation, limit-points, and postbuckling. A novel technique is developed for evaluating bifurcation from nonlinear prebuckling states utilizing asymmetric spatial discretization to introduce numerical perturbations. Advantages of the asymmetric meshing technique (AMT) over traditional techniques include efficiency, robustness, ease of application, and solution of the actual (not modified) problems. The AMT is validated by comparison to traditional numerical analysis of a benchmark problem and verified by comparison to experimental data. Applying the technique, bifurcation in a benchmark shell-buckling problem is correctly identified. Excellent agreement between the numerical and experimental results are obtained for a number of composite shells although predictive capability decreases for stiffer (thicker) specimens which is attributed to compliance of the test fixture. Restraining the axial edge (simple support) has the effect of creating a more complex response which involves unstable bifurcation, limit-point buckling, and dynamic collapse. Such shells were noted to bifurcate into asymmetric deformation modes but were undamaged during testing. Shells in this study which were damaged were not observed to bifurcate. Thus, a direct link between bifurcation and atypical damage could not be established although the mechanism (bifurcation) was identified. Recommendations for further work in these related areas are provided and include extensions of the AMT to other shell geometries and structural problems.

Wardle, Brian L.↗

NASA NTRS2001

X-38 Experimental Aeroheating at Mach 10

This report provides an update of the hypersonic aerothermodynamic wind tunnel test program conducted at the NASA Langley Research Center in support of the X-38 program. Global surface heat transfer distributions were measured on 0.0177 and 0.0236 scale models of the proposed X-38 configuration at Mach 10 in air. The parametrics that were investigated primarily include freestream unit Reynolds numbers of 0.6 to 2.2 million per foot and body flap deflections of 15, 20, and 25 deg for an angle-of-attack of 40 deg. The model-scale variance was tested to obtain laminar, transitional, and turbulent heating levels on the defected bodyflaps. In addition, a limited investigation of forced boundary layer transition through the use of discrete roughness elements was performed. Comparisons of the present experimental results to computational predictions and previous experimental data were conducted Laminar, transitional, and turbulent heating levels were observed on the deflected body flap, which compared favorably to the computational results and to the predicted heating based on the flight aerothermodynamic database.

Berry, Scott A.↗

DOE OSTI2015

A conductive heat transfer model for particle flows over immersed surfaces

Not Available

Morris, A. B.↗

DOE OSTI2017

CFD-DEM modeling the effect of column size and bed height on minimum fluidization velocity in micro fluidized beds with Geldart B particles

The fluidization behavior of Geldart B particles in micro fluidized beds is investigated numerically using Computational Fluid Dynamics coupled with Discrete Element Method (CFD-DEM) available in the open-source Multiphase Flow with Interphase eXchanges (MFIX) code. The effects of different bed inner diameters (D) of 8 mm, 12 mm, 16 mm and various initial static bed heights (H) were examined. It is found that both decreasing the column diameter and increasing the bed height in a micro fluidized bed increases the minimum fluidization velocity (Umf). The observed overshoot in pressure drop that occurs before the onset of fluidization decreases in magnitude with increasing column diameter, however there is less sensitivity to bed height. Overall, the numerical results agree qualitatively with existing theoretical correlations and experimental studies. The simulations show that both column diameter and particle-wall friction contribute to the variation in minimum fluidization velocity. Finally, these two factors are coupled and hard to separate. The detailed influences of wall friction on minimum fluidization velocity are then investigated for a prescribed column diameter of 8 mm by varying the wall friction from 0 to 0.4.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Kinetic Plasma Simulation in the MOOSE Framework: Verification of Electrostatic Particle In Cell Capabilities

In magnetic confinement nuclear fusion reactors, the interaction between the plasma edge and plasma facing components is extremely important. At the plasma edge, a kinetic representation such as particle-in-cell (rather than a fluid representation) is required to accurately capture the plasma behavior. General purpose particle-in-cell plasma simulation capabilities have been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. This new capability is a part of the development of a new MOOSE-based framework for modeling plasma facing components, the Fusion ENergy Integrated multiphys-X (FENIX) framework. In this work, the verification of foundational particle-in-cell capabilities in FENIX is presented. This new plasma simulation capability has three main components: moving particles in discrete steps on the finite element mesh, mapping charge density from the particle's location to the finite element mesh, and solving for the electrostatic potential based on the charge density mapped from particles to the mesh. In this paper, simple verification problems demonstrating each of these new capabilities are presented, and future work includes electromagnetic capabilities and Monte Carlo collisions with neutral gas particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Computational simulation of damping in composite structures

A computational methodology is developed for the prediction of passive damping in composite structures. The method involves multiple levels of damping modeling by integrating micromechanics, laminate, and structural damping theories. The effects of temperature and moisture on structural damping are included. The simulation of damping in the structural level is accomplished with finite-element discretization. Applications are performed on graphite/epoxy composite beams, plates, and shells to illustrate the methodology. Additional parametric studies demonstrate the variation of structural modal damping with respect to ply angles, fiber volume ratio, and temperature.

Saravanos, D. A.↗

NASA NTRS1991

Computational simulation of damping in composite structures

Saravanos, D. A.↗

NASA NTRS1992

A macro-micromechanics analysis of a notched metal matrix composite

Macro- and micromechanics analysis were conducted to determine the matrix and fiber behaviors near the notch in a center-notched metal-matrix composite. In this approach, the macrolevel analysis models the entire notched specimen using a 3D finite element program that uses the vanishing-fiber-diameter model to simulate the elastic-plastic behavior of the matrix and the elastic behavior of the fiber. The microlevel behavior is analyzed using a discrete fiber-matrix model containing one fiber and the surrounding matrix. The viability of this analysis is demonstrated using results for a boron/aluminum monolayer.

Bigelow, Catherine A.↗

NASA NTRS2019

Photogrammetry-Based Analysis of the On-Orbit Structural Dynamics of the Roll-Out Solar Array

The Roll-Out Solar Array (ROSA) flight experiment was launched to the International Space Station (ISS) on June 3rd, 2017. ROSA is an innovative, lightweight solar array with a flexible substrate that makes use of the stored strain energy in its composite structural members to provide deployment without the use of motors. This paper will discuss the results of various structural dynamics experiments conducted on the ISS during the weeks following launch. Data gathered from instrumentation on the solar array wing during the experiments was previously compared with pre-flight predictions from two different Finite Element Modeling (FEM) efforts. In this paper, data generated from photogrammetry is compared with accelerometer data and used to extend previous conclusions. Whereas previous analyses were only able to track the accelerations of six discrete points on the structure and photovoltaic (PV) blanket of ROSA, the photogrammetry analysis makes available displacements for dozens of points distributed throughout the array. This larger data set makes it possible to compare higher-order PV blanket modes with FEM predictions, in addition to verifying conclusions reached using accelerometer data. The goal in this effort was to better understand the performance of ROSA and to improve modeling efforts for future designs of similar solar arrays.

Chamberlain, Matthew K.↗

DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2011

Surrogate Modeling of High-Fidelity Fracture Simulations for Real-Time Residual Strength Predictions

A surrogate model methodology is described for predicting in real time the residual strength of flight structures with discrete-source damage. Starting with design of experiment, an artificial neural network is developed that takes as input discrete-source damage parameters and outputs a prediction of the structural residual strength. Target residual strength values used to train the artificial neural network are derived from 3D finite element-based fracture simulations. A residual strength test of a metallic, integrally-stiffened panel is simulated to show that crack growth and residual strength are determined more accurately in discrete-source damage cases by using an elastic-plastic fracture framework rather than a linear-elastic fracture mechanics-based method. Improving accuracy of the residual strength training data would, in turn, improve accuracy of the surrogate model. When combined, the surrogate model methodology and high-fidelity fracture simulation framework provide useful tools for adaptive flight technology.

Spear, Ashley D.↗

DOE OSTI2018

Fourier amplitude sensitivity test applied to dynamic combined finite-discrete element methods-based simulations

Fracture propagation plays a key role for a number of applications of interest to the scientific community, from dynamic fracture processes like spallation and fragmentation in metals to failure of ceramics, airplane wings, etc. Simulations of material deformation and fracture propagation rely on accurate knowledge of material characteristics such as material strength and the amount of energy being dissipated during the fracture process. Within the combined finite-discrete element method (FDEM) framework material fracture behavior is typically described through a parametrized softening curve, which defines a stress-strain relationship unique to each material. We apply the Fourier amplitude sensitivity test to explore how each of these parameters influences the simulated damage processes and to determine the key input parameters that have the most impact on the model response. We present several sensitivity numerical experiments for the simulation of a split Hopkinson pressure bar (SHPB) test for weathered granite samples using different combinations of model parameters. We validate the obtained results against SHPB experimental data. The experiments show that the model is mostly sensitive to parameters related to tensile and shear strengths, even in the presence of other parameter perturbations. Lastly, the results suggest that the specification of tensile and shear strengths at the interfaces dominate the stress-time history of the FDEM simulation of SHPB test.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Solvable model of quantum phase transitions and the symbolic-manipulation-based study of its multiply degenerate exceptional points and of their unfolding

It is known that the practical use of non-Hermitian (i.e., typically, PT-symmetric) phenomenological quantum Hamiltonians H≠H{sup †} requires an efficient reconstruction of an ad hoc Hilbert-space metric Θ=Θ(H) which would render the time-evolution unitary. Once one considers just the N-dimensional matrix toy models H=H{sup (N)}, the matrix elements of Θ(H) may be defined via a coupled set of N{sup 2} polynomial equations. Their solution is a typical task for computer-assisted symbolic manipulations. The feasibility of such a model-completion construction is illustrated here via a discrete square well model H=p{sup 2}+V endowed with a k-parametric close-to-the-boundary interaction V. The model is shown to possess (possibly, multiply degenerate) exceptional points marking the phase transitions which are attributable, due to the exact solvability of the model at any N<∞, to the loss of the regularity of the metric. In the parameter-dependence of the energy spectrum near these singularities one encounters a broad variety of alternative, topologically non-equivalent scenarios. -- Highlights: •New elementary non-Hermitian quantum Hamiltonians with real spectra proposed. •Exceptional points found and studied. •Non-equivalent stability-loss patterns of phase transition identified. •Hermitization matrices of metrics Θ constructed via symbolic manipulations and extrapolations at all N and k.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

A mesh-independent method for planar three-dimensional crack growth in finite domains: A mesh-independent method for crack growth in 3D finite domains

The discrete crack mechanics (DCM) method is a dislocation-based crack modeling technique where cracks are constructed using Volterra dislocation loops. The method allows for the natural introduction of displacement discontinuities, avoiding numerically expensive techniques. Mesh dependence in existing computational modeling of crack growth is eliminated by utilizing a superposition procedure. The elastic field of cracks in finite bodies is separated into two parts: the infinite-medium solution of discrete dislocations and an finite element method solution of a correction problem that satisfies external boundary conditions. In the DCM, a crack is represented by a dislocation array with a fixed outer loop determining the crack tip position encompassing additional concentric loops free to expand or contract. Solving for the equilibrium positions of the inner loops gives the crack shape and stress field. The equation of motion governing the crack tip is developed for quasi-static growth problems. Additionally, convergence and accuracy of the DCM method are verified with two- and three-dimensional problems with well-known solutions. Crack growth is simulated under load and displacement (rotation) control. In the latter case, a semicircular surface crack in a bent prismatic beam is shown to change shape as it propagates inward, stopping as the imposed rotation is accommodated.

42 ENGINEERING↗