SEARCH · Engineering Papers
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.
This submission contains the presentation slides and recordings from EGS Collab Modeling and Simulation Working Group (MSWG) teleconferences number 99 through 128. These teleconferences served three objectives for the project: 1) share simulation results, 2) communicate field activities and results to the simulation teams, and 3) hold open scientific discussions on EGS topics.
The core structures of future nuclear systems require tolerance to extreme irradiation, and some critical components, for example, the fuel cladding in Sodium-cooled Fast Reactors (SFRs), have to maintain mechanical integrity to very high doses of 200 -400 dpa at high temperatures up to 700 degrees C. The high Cr nanostructured ferritic alloys (NFAs) are under intense research worldwide as a candidate core material. Although the NFAs have some admirable characteristics for high-temperature applications, their crack sensitivity is very high at high temperatures. The fracture toughness of high strength NFAs is unacceptably low above 300 degrees C. The objective of this study is to develop processes and microstructures with improved high temperature fracture toughness and ductility. To optimize the afterextrusion heat treatment condition, both the computational simulation technique on phase equilibrium and the basic microstructural and mechanical characterization have been carried out. 9 Cr-NFA was produced by the mechanical alloying of pre-alloyed Fe-9Cr base metallic powder and yttria particles, and subsequent extrusion. The post-extrusion heat-treatments of various conditions were applied to the asextruded NFA. The tensile and fracture toughness tests were conducted for as-extruded and heat-treated samples at up to 700 degrees C. Fracture toughness of the NFA has increased by more than 40% at every testing temperature after heat-treatment in the inter-critical temperature range. The increment of fracture toughness of the NFA after post-extrusion heat-treatment is attributed to the increased strength at below 500 degrees C, and an increased ductility at 700 degrees C.
The resistance of polycrystalline materials to intergranular fracture depends on the toughness of the interfaces between neighboring grains. Using finite element simulations on model, two-dimensional microstructures, we demonstrate that the steady-state crack growth resistance of a polycrystal that contains exclusively high toughness interfaces is suboptimal: its crack growth resistance may be increased by introducing a minority fraction of weak interfaces at random locations within the interface network. Furthermore, the weaker are these interfaces, the greater is their toughening effect. These weak interfaces toughen the material through crack deflection and branching, nucleation of secondary cracks, and inelastic dissipation within intact interfaces. Furthermore, our findings suggest that fracture resistance-oriented material development should not aim to eliminate all weak interfaces, but rather judiciously incorporate a low fraction of them.
The peridynamic theory of solid mechanics is applied to modeling the deformation and fracture of micrometer-sized particles made of organic crystalline material. A new peridynamic material model is proposed to reproduce the elastic–plastic response, creep, and fracture that are observed in experiments. The model is implemented in a three-dimensional, meshless Lagrangian simulation code. In the small deformation, elastic regime, the model agrees well with classical Hertzian contact analysis for a sphere compressed between rigid plates. Under higher load, material and geometrical nonlinearity is predicted, leading to fracture. Finally, the material parameters for the energetic material CL-20 are evaluated from nanoindentation test data on the cyclic compression and failure of micrometer-sized grains.
Carbon nanotubes (CNTs) are promising candidates for use as reinforcements in next generation structural composite materials because of their extremely high specific stiffness and strength. They cannot, however, be viewed as simple replacements for carbon fibers because there are key differences between these materials in areas such as handling, processing, and matrix design. It is impossible to know for certain that CNT composites will represent a significant advance over carbon fiber composites before these various factors have been optimized, which is an extremely costly and time intensive process. This work attempts to place an upper bound on CNT composite mechanical properties by performing molecular dynamics simulations on idealized model systems with a reactive forcefield that permits modeling of both elastic deformations and fracture. Amorphous carbon (AC) was chosen for the matrix material in this work because of its structural simplicity and physical compatibility with the CNT fillers. It is also much stiffer and stronger than typical engineering polymer matrices. Three different arrangements of CNTs in the simulation cell have been investigated: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. The SWNT and MWNT array systems are clearly idealizations, but the SWNT bundle system is a step closer to real systems in which individual tubes aggregate into large assemblies. The effect of chemical crosslinking on composite properties is modeled by adding bonds between the CNTs and AC. The balance between weakening the CNTs and improving fiber-matrix load transfer is explored by systematically varying the extent of crosslinking. It is, of course, impossible to capture the full range of deformation and fracture processes that occur in real materials with even the largest atomistic molecular dynamics simulations. With this limitation in mind, the simulation results reported here provide a plausible upper limit on achievable CNT composite properties and yield some insight on the influence of processing conditions on the mechanical properties of CNT composites.
This report provides an analysis of engineering and geologic data collected from the Hydraulic Fracturing Test Site #1 (HTSF-1) project in the southern Midland Basin, Reagan County, Texas. The site is being studied as part of the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative at the National Energy Technology Laboratory (NETL). The data collected is intended to provide a basis of understanding of the site and to construct simulation models. This report provides a summary analysis of the engineering and geologic data collected from the project to provide a basis of understanding of the site and to construct simulation models. The report examines various possible correlations in engineering properties based on natural fracture data from the program, which consisted of 11 horizontal wells, one vertical well, and one slant well. The horizontal wells are in two horizons: 1) Upper and 2) Middle Wolfcamp formations. Program data include: fracture frequency and fracture orientation data from four core runs in a slant well; fracture spacing and orientation from a vertical pilot well; laboratory triaxial testing and mineralogical determinations; and porosity results from magnetic resonance analyses from various wells, together with observations based on the data collection. In addition, available references were reviewed on the site for additional insights. As the focus of the report is on the natural system, hydraulic fracture data from the site were not examined in detail in this report. Data variability is the chief observation in examination of the database. Fracture frequency in the slant well can range from sections with values as high as five fractures per ft to sections up to 100+ ft in length with no natural observed fractures. Fracture spacing across is typically less than 10 ft, but can range up to hundreds of feet. Apparent fracturing shows the trends in two predominate orientations, E-SW and WNW-ESE, but minor variations exist. The rock units vary across the site from siliceous mudstones to calcareous mudstones, showing a general layering with depth. The laboratory properties such as strength and modulus show no apparent trend with depth, but appear to correlate with rock mineralogy with high strength and modulus values where calcium content is high. In addition, an attempt to examine variability and mineralogy on a larger scale was made using a color-coded system based on gamma ray measurements. A staged colored approach was adopted, presuming that lower gamma ray values indicate higher value of calcium content (blue scale) and that higher gamma ray values indicate higher clay mineral content (orange scale). As provided in report appendices, the system correlated well with visual examination of the slant core and the petrofabric analyses of the vertical pilot well. The results showed large variability in mineral content along the horizontal plane across the site. The change in mineral content was also rapid, on a scale less than that of the average hydraulic fracture stage length of about 180 ft.
Fracture toughness is an important material property that plays a role in determining the in-service mechanical performance and adhesion of thin ceramic films. Unfortunately, measuring thin film fracture toughness is affected by influences from the substrate and the large residual stresses that can exist in the films. In this paper, we explore a promising new technique that potentially overcomes these problems based on nanoindentation testing of micro-pillars produced by focused ion beam milling of the films. By making the pillar diameter approximately equal to its length, the residual stress in the pillar’s upper portion is almost fully relaxed, and when indented with a sharp Berkovich indenter, the pillars fracture by splitting at reproducible loads that are readily quantified by a sudden displacement excursion in the load displacement behavior. Cohesive finite element simulations are used to analyze and develop, for a given material, a simple relation between the critical load at failure, pillar radius, and fracture toughness. The main novel aspect of this work is that neither crack geometries nor crack sizes need to be measured post test. Furthermore, the residual stress can be measured at the same time with toughness, by comparing the indentation results from the stress-free pillars and the as-deposited film. The method is tested on three different hard coatings formed by physical vapor deposition: titanium nitride, chromium nitride, and a CrAlN/Si 3 N 4 nanocomposite. Results compare well to independently measured values of fracture toughness for the three brittle films. The technique offers several benefits over existing methods.
The effects of carbon nanotube (CNT) sidewall fluorination on the interface toughness of the CNT epoxy interface have been comprehensively investigated. Nanoscale quantitative single-CNT pull-out experiments have been conducted on individual fluorinated CNTs embedded in an epoxy matrix, in situ, within a scanning electron microscope (SEM) using an InSEM{sup ®} nanoindenter assisted micro-device. Equations that were derived using a continuum fracture mechanics model have been applied to compute the interfacial fracture energy values for the system. The interfacial fracture energy values have also been independently computed by modeling the fluorinated graphene-epoxy interface using molecular dynamics simulations and adhesion mechanisms have been proposed.
Corrosion-induced concrete cracking is a significant durability problem for reinforced concrete structures. Considerable research has been carried out in the last few decades to understand and model the expansion mechanism of the corrosion products around the reinforcing bar and simulate the cracking behaviour of the concrete cover. In this paper, a new corrosion model based on non-uniform corrosion expansion is formulated and validated against experimental data. A meso-scale fracture model, consisting of aggregates, cement paste/mortar and ITZ, is established for the cases of both middle and side reinforcing bars. Under the developed corrosion and concrete fracture model, the cracking phenomena of the concrete cover are accurately simulated. It has been found that the non-uniform corrosion model can be used to express the realistic corrosion rust progression around the reinforcing bar, with the best accuracy. It has also been found that some microcracks occur before they are connected to form the dominating discrete crack which usually appears on the concrete surface. Moreover, the effects of the corrosion variables, as well as other key material and geometric parameters, on surface cracking of concrete are investigated.
n this work, a hydrocode PAGOSA with FLIP+MPM is presented and exercised to investigate the fracture in ductile material. The merit of PAGOSA with FLIP+MPM to solve the advection problem is first illustrated by a solid piston periodically moving in a sealed tube with air. Furthermore, the ability of PAGOSA with FLIP+MPM to capture the fracture in material is shown by a simple stretching fracture in ductile material. In both of two benchmark problems, the PAGOSA results and analytical solutions are also presented for comparison. Then PAGOSA with FLIP+MPM is used to model complex spall in ductile material, which is a challenging problem in engineering applications. The convergences of PAGOSA with FLIP+MPM—based on both the mesh size and the marker density—are investigated by monitoring free surface velocity. To further show the ability of PAGOSA with FLIP+MPM to predict fracture in ductile material, the numerical results are compared with the experimental results and other published numerical results. Moreover, the effect of the spall parameter in PAGOSA with FLIP+MPM on the numerical simulation is also analyzed by investigating the free surface velocity. Finally, the effects of the peak compressive stress, the tensile strain rate, and the loading rate on the spallation are further investigated. The numerical results show that PAGOSA with FLIP+MPM can improve PAGOSA's performance when applied to predicting fracture in ductile materials, and is robust enough to accurately predict spall fracture in ductile materials.
Explore the source record for details and available documents.
The purpose of this project was to develop methods and tools that will aid in safety evaluation of nuclear fuels and licensing of nuclear reactors relating to accidents.The objectives were to develop more detailed and faster computations of fission product transport and aerosol evolution as they generally relate to nuclear fuel and/or nuclear reactor accidents. The two tasks in the project related to molecular transport in nuclear fuel and aerosol transport in reactor vessel and containment. For both the tasks, explorations of coupling of Direct Simulation Monte Carlo with Navier-Stokes solvers or the Sectional method were not successful. However, Mesh free methods for the Direct Simulation Monte Carlo method were successfully explored.These explorations permit applications to porous and fractured media, and arbitrary geometries.The computations were carried out in Mathematica and are fully parallelized. The project has resulted in new computational tools (algorithms and programs) that will improve the fidelity of computations to actual physics, chemistry and transport of fission products in the nuclear fuel and aerosol in reactor primary and secondary containments.
A finite element model (FEM) was developed to simulate the mixed-mode delamination failure at the encapsulant interface within PV modules under mixed-mode loading condi- tions. Parameters to characterize the bilinear traction-separation constitutive behavior of cohesive elements were determined by testing the pure mode I critical fracture energy of the photovotaic(PV) cell interface using the width-tapered beam and fitting the simulated load-displacement curve with the tested curve. This developed FEM can be used to predict the critical temperature changes, at which delamination would initiate, of PV modules with different sizes and provide guidance for the design of PV module with better reliability.
TOUGH is a widely used general purpose numerical simulator designed for fluid and heat flows of multiphase, multicomponent fluid mixtures in porous and fractured media, which has been applied to many real-world problems such as underground geological storage, geothermal reservoir engineering or nuclear waste disposal, to name a few. `toughio` is a lightweight, object-oriented and vectorized Python library that aims to provide user-friendly routines to facilitate pre- and post-processing of a TOUGH simulation. Although its primary target is the latest version TOUGH3, `toughio` can also be used to read/write TOUGH2 input/output files. Currently, to the best of our knowledge, only PyTOUGH offers an exhaustive number of features to carry out a complete TOUGH simulation using a scripting language. `toughio` and PyTOUGH share the same objectives, yet with different approaches and philosophies. On the one hand, a PyTOUGH mesh is represented as a MULGRAPH geometry where elements can be unstructured horizontally but only layered vertically. It also checks that the model parameters are valid to ensure that a TOUGH simulation runs smoothly. On the other hand, although it provides basic meshing features, `toughio` mostly relies on common third-party softwares (e.g. Abaqus, FLAC3D, Gmsh, LaGriT) to generate the mesh by importing and converting it to a TOUGH mesh, which also conveniently facilitates the coupling of TOUGH with any other simulator that also supports the same mesh formats. Besides, `toughio` is designed to be just a TOUGH input and output file reader/writer. It assumes that the user knows what he is doing and does not verify the validity of the model parameters (`toughio` does verify the functional arguments).
Molecular dynamics simulations are used to investigate the mechanical behavior of Cu nanowires (NWs) with planar defects such as grain boundaries (GBs), twin boundaries (TBs), stacking faults (SFs), etc. To investigate how the planar defects affect the deformation and fracture mechanisms of naowires, three types of nanowires are considered in this paper: (1) polycrystalline Cu nanowire; (2) single-crystalline Cu nanowire with twin boundaries; and (3) single-crystalline Cu nanowire with stacking faults. Because of the large fraction of atoms at grain boundaries, the energy of grain boundaries is higher than that of the grains. Thus, grain boundaries are proved to be the preferred sites for dislocations to nucleate. Moreover, necking and fracture prefer to occur at the grain boundary interface owing to the weakness of grain boundaries. For Cu nanowires in the presence of twin boundaries, it is found that twin boundaries can strength nanowires due to the restriction of the movement of dislocations. The pile up of dislocations on twin boundaries makes them rough, inducing high energy in twin boundaries. Hence, twin boundaries can emit dislocations, and necking initiates at twin boundaries. In the case of Cu nanowires with stacking faults, all pre-existing stacking faults in the nanowires are observed to disappear during deformation, giving rise to a fracture process resembling the samples without stacking fault.
In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.
This report describes the results of the first of its kind monitoring of a hydrofracture operation with electromagnetic measurements. The researchers teamed with oil and gas producer Encana Corporation to design and execute a borehole to surface monitoring of three fracture stages at a well pad in central Colorado. The field project consisted of an equipment upgrade, a survey design and modeling phase, several weeks of data collection, and data processing and interpretation. Existing Depth to Surface Resistivity (DSR) instrumentation was upgraded to allow for continuous high precision recording from downhole sources. The full system can now collect data continuously for up to 72 hours, which is sufficient to measure data for 10 frac stages. Next we used numerical modeling and frac treatment data supplied by Encana to design a field survey to detect EM signal from induced fractures. Prior to modeling we developed a novel technique for using well casing as an antenna for a downhole source. Modeling shows that 1) a measurable response for an induced fracture could be achieved if the facture fluid was of high salinity 2) an optimum fracture response is created when the primary source field is parallel to the well casing but perpendicular to the fracture direction. In mid-July, 2014 we installed an array of more than 100 surface sensors, distributed above the treatment wells and extending for approximately 1 km north and 750 m eastward. We applied a 0.6 Hz square wave signal to a downhole current electrode located in a horizontal well 200 m offset from the treatment well with a return electrode on the surface. The data were transmitted to a recording trailer via Wi-Fi where we monitored receiver and transmitter channels continuously in a 72-hour period which covered 7 frac stages, three of which were high salinity. Although the background conditions were very noisy we were able to extract a clear signal from the high salinity stages. Initial data interpretation attempts consisting of trying to directly invert collected data using a simple background starting model did not produce reasonable models. We next used a simulated data set to develop a constrained inversion workflow that places the downhole fracture in the correct location. Finally, we used a combination of forward and inverse models to fit the collected data to a model that incorporated 3 frac stages. This project provided a first-of-its-kind application of EM technology to map an induced fracture. The results demonstrated a clear anomaly during the operation that can be fit to a reasonable model of an induced fracture.
Commercially pure aluminum alloy, AA1100, presents good electrical and thermal conductivity, high formability, and low cost. Those favorable characteristics have the potential to enable bipolar plates with improved economics and enhanced performance compared to current stainless steel bipolar plates for proton exchange membrane fuel cells. An accurate constitutive model is essential to develop and optimize processing parameters and effectively control the forming process. Here, the objective of this work is to develop a constitutive model of AA1100 that is able to simulate stress-strain relation, formed geometry, and predict the onset of fracture strain to avoid forming failure. Initially, a set of tensile tests at temperature between 300 and 500°C and strain rate between 0.005 and 1.0/s were conducted to examine the deformation behavior. Then, a set of damage-based unified visco-plastic constitutive equations is proposed and calibrated based on the results of stress-strain data. A genetic algorithm optimization method is applied to search for best fitting material constants in constitutive equations. The proposed model shows good predictability of both the stress-strain relation and fracture strain at low strain rate and high temperature conditions. The accuracy of proposed model is also evaluated statistically. A comparison of the proposed model with three popular models (Arrhenius-type mode, Johnson-Cook model and Zerilli-Armstrong model) was made. The proposed model shows the best experimental agreement with correlation coefficient of 0.96 in contrast to 0.25, 0.38 and 0.75 for the popular models, respectively. The proposed model can help to optimize the elevated temperature forming process and guide die design to enable optimal geometric features in the formed components.