SEARCH · Engineering Papers
Results for “finite element simulation”
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.
Simulations of the effects of proppant placement on the conductivity and mechanical stability of hydraulic fractures
Here, we generate a wide range of models of proppant-packed fractures using discrete element simulations, and measure fracture conductivity using finite element flow simulations. This allows for a controlled computational study of proppant structure and its relationship to fracture conductivity and stress in the proppant pack. For homogeneous multi-layered packings, we observe the expected increase in fracture conductivity with increasing fracture aperture, while the stress on the proppant pack remains nearly constant. This is consistent with the expected behavior in conventional proppant-packed fractures, but the present work offers a novel quantitative analysis with an explicit geometric representation of the proppant particles. In single-layered packings (i.e. proppant monolayers), there is a drastic increase in fracture conductivity as the proppant volume fraction decreases and open flow channels form. However, this also corresponds to a sharp increase in the mechanical stress on the proppant pack, as measured by the maximum normal stress relative to the side crushing strength of typical proppant particles. We also generate a variety of computational geometries that resemble highly heterogeneous proppant packings hypothesized to form during channel fracturing. In some cases, these heterogeneous packings show drastic improvements in conductivity with only moderate increase in the stress on the proppant particles, suggesting that in certain applications these structures are indeed optimal. We also compare our computer-generated structures to micro computed tomography imaging of a manually fractured laboratory-scale shale specimen, and find reasonable agreement in the geometric characteristics.
Steam generator tube rupture simulation using extended finite element method
Not Available
The importance of sub-meter-scale snow roughness on conductive heat flux of Arctic sea ice
Abstract The conductive heat flux through the snow and ice is a critical component of the mass and energy budgets in the Arctic sea ice system. We use high horizontal resolution (3–15 cm) measurements of snow topography to explore the impacts of sub-meter-scale snow surface roughness on heat flux as simulated by the Finite Element method. Simulating horizontal heat flux in a variable snow cover modestly increases the total simulated heat flux. With horizontal heat flux, as opposed to simple 1D-vertical heat flux modeling, the simulated heat flux is 10% greater than that for uniform snow with the same mean snow thickness for a 31.5 × 21 m region of sea ice (the largest region we studied). Vertical-only (1D) heat flux simulates just a 6% increase for the same region. However, this is highly dependent on observation resolution. Had we measured the snow cover at 1 m horizontal spacing or greater, simulating horizontal heat flux would not have changed the net heat flux from that simulated with vertical-only heat flux. These findings suggest that measuring and modeling snow roughness at sub-meter horizontal scales may be necessary to accurately represent horizontal heat flux on level Arctic sea ice.
In Situ 3D Imaging of Catalysis Induced Strain in Gold Nanoparticles
Multielectron transfer processes are crucially important in energy and biological science but require favorable catalysts to achieve fast kinetics. Nanostructuring catalysts can dramatically improve their properties, which can be difficult to understand due to strain- and size-dependent thermodynamics, the influence of defects, and substrate-dependent activities. Here, we report three-dimensional (3D) imaging of single gold nanoparticles during catalysis of ascorbic acid decomposition using Bragg coherent diffractive imaging (BCDI). Local strains were measured in single nanoparticles and modeled using reactive molecular dynamics (RMD) simulations and finite element analysis (FEA) simulations. RMD reveals the pathway for local strain generation in the gold lattice: chemisorption of hydroxyl ions. FEA reveals that the RMD results are transferable to the nanocrystal sizes studied in the experiment. Our study probes the strain-activity connection and opens a powerful avenue for theoretical and experimental studies of nanocrystal catalysis.
Synchronous and Asynchronous Time Integration for Multiscale Simulations Using Hybridized Finite Element Methods.
Abstract not provided.
Simulation of crack induced nonlinear elasticity using the combined finite-discrete element method
Numerical simulation of nonlinear elastic wave propagation in solids with cracks is indispensable for decoding the complicated mechanisms associated with the nonlinear ultrasonic techniques in NonDestructive Testing (NDT). Here, we introduce a two-dimensional implementation of the combined finitediscrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), to explicitly simulate the crack induced nonlinear elasticity in solids with both horizontal and inclined cracks. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. The simulation results show that for cracked models, nonlinear elasticity is generated only when the excitation amplitude is large enough to trigger the contact between crack surfaces, and the nonlinear behavior is very sensitive to the crack surface contact. The simulations reveal the influence of normal and tangential contact on the nonlinear elasticity generation. Moreover, the results demonstrate the capabilities of FDEM for decoding the causality of nonlinear elasticity in cracked solid and its potential to assist in Non-Destructive Testing (NDT).
Force Limiting Vibration Tests Evaluated from both Ground Acoustic Tests and FEM Simulations of a Flight Like Vehicle System Assembly
Marshall Space Flight Center has conducted a series of ground acoustic tests with the dual goals of informing analytical judgment, and validating analytical methods when estimating vibroacoustic responses of launch vehicle subsystems. The process of repeatedly correlating finite element-simulated responses with test-measured responses has assisted in the development of best practices for modeling and post-processing. In recent work, force transducers were integrated to measure interface forces at the base of avionics box equipment. Other force data was indirectly measured using strain gauges. The combination of these direct and indirect force measurements has been used to support and illustrate the advantages of implementing the Force Limiting approach for equipment qualification tests. The comparison of force response from integrated system level tests to measurements at the same locations during component level vibration tests provides an excellent illustration. A second comparison of the measured response cases from the system level acoustic tests to finite element simulations has also produced some principles for assessing the suitability of Finite Element Models (FEMs) for making vibroacoustics estimates. The results indicate that when FEM models are employed to guide force limiting choices, they should include sufficient detail to represent the apparent mass of the system in the frequency range of interest.
Finite-element analyses and fracture simulation in thin-sheet aluminum alloy
A two-dimensional, elastic-plastic finite-element analysis was used with a critical crack-tip-opening angle (CTOA) fracture criterion to model stable crack growth in thin-sheet 2024-T3 aluminum alloy under monotonic loading after precracking at different cyclic stress levels. Tests were conducted on three types of specimens: middle-crack, three-hole-crack and blunt-notch tensile specimens. An experiment technique was developed to measure CTOA during crack growth initiation and stable tearing using a high-resolution video camera and recorder. Crack front shapes were also measured during initiation and stable tearing using a fatigue marker-load technique. Three-dimensional elastic-plastic finite-element analyses of these crack shapes for stationary cracks were conducted to study the crack-front opening displacements. Predicted load against crack extension on middle-crack tension specimens agreed well with test results even for large-scale plastic deformations. The analyses were able to predict the effects of specimen size and precracking stress history on stable tearing. Predicted load against load-line displacements agreed well with test results up to maximum load bu the analyses tended to overpredict displacements as crack grew beyond the maximum load under displacement-controlled conditions. During the initiation phase, the measured CTOA values were high but decreased and remained nearly constant after a small amount of stable tearing. The constant value of CTOA agree well with the calculated value from the finite-element analysis. The larger CTOA values measured at the sheet surface during the initiation phase may be associated with the crack tunneling observed in the tests. Three-dimensional analyses for nonstraight crack fronts predicted much higher displacements near the free surface than in the interior.
Simulating Space Capsule Water Landing with Explicit Finite Element Method
A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.
Assessment of a Hybrid Continuous/Discontinuous Galerkin Finite Element Code for Geothermal Reservoir Simulations
FALCON (“Fracturing And Liquid CONvection”) is a hybrid continuous / discontinuous Galerkin finite element geothermal reservoir simulation code based on the MOOSE (“Multiphysics Object-Oriented Simulation Environment”) framework being developed and used for multiphysics applications. In the present work, a suite of verification and validation (“V&V”) test problems for FALCON was defined to meet the design requirements, and solved to the interests of enhanced geothermal system (“EGS”) design. Furthermore, the intent for this test problem suite is to provide baseline comparison data that demonstrates the performance of the FALCON solution methods. The simulation problems vary in complexity from singly mechanical or thermo process, to coupled thermo-hydro-mechanical processes in geological porous media. Numerical results obtained by FALCON agreed well with either the available analytical solution or experimental data, indicating the verified and validated implementation of these capabilities in FALCON. Some form of solution verification has been attempted to identify sensitivities in the solution methods, where possible, and suggest best practices when using the FALCON code.
Crystal plasticity modeling of 3rd generation multi-phase AHSS with martensitic transformation
Abstract not provided
Calibration of higher eigenmodes of cantilevers
A method is presented for calibrating the higher eigenmodes (resonant modes) of atomic force microscopy cantilevers that can be performed prior to any tip-sample interaction. The method leverages recent efforts in accurately calibrating the first eigenmode by providing the higher-mode stiffness as a ratio to the first mode stiffness. A one-time calibration routine must be performed for every cantilever type to determine a power-law relationship between stiffness and frequency, which is then stored for future use on similar cantilevers. Then, future calibrations only require a measurement of the ratio of resonant frequencies and the stiffness of the first mode. This method is verified through stiffness measurements using three independent approaches: interferometric measurement, AC approach-curve calibration, and finite element analysis simulation. Power-law values for calibrating higher-mode stiffnesses are reported for several cantilever models. Once the higher-mode stiffnesses are known, the amplitude of each mode can also be calibrated from the thermal spectrum by application of the equipartition theorem.
The Application of Hydraulic and Thermal Stimulation Techniques to Create Enhanced Geothermal Systems
Geothermal energy production requires temperature at depth, a working fluid to transport heat to the surface, and a network of fluid pathways connecting injection and production wells. These requirements constrain the availability of geothermal energy resources. Enhanced Geothermal System (EGS) development expands geothermal resource availability by artificially generating networks of fluid pathways in the reservoir. This is accomplished through hydraulic and thermal stimulation of the geothermal reservoir. Hydraulic and thermal techniques were applied at the Raft River geothermal field as part of a Department of Energy (DOE) EGS development project. The test well RRG-9 ST1 was hydraulically stimulated three times between February 2012 and April 2014. Since June, 2013 continuous injection of sub-reservoir temperature plant water has been used to thermally stimulate the well. The stimulation program has resulted in a large increase in the injectivity of the well, an increase from less than 20 gpm in June 2013 to nearly 1,000 gpm in April 2016. Wellbore imaging and temperature analysis identified a northeast striking fracture zone intersecting the well between 5,640 and 5,660 ft. measured depth (MD) has nominally accepted all of the injected fluid. Microseismic activity, related to stimulation program and plant activity was used to track possible fluid pathways within the reservoir. The microseismic data indicates that the injected fluid moves primarily to the northeast through the Narrows Zone, a northeast striking fault structure that bisects the field. Tracers injected into RRG-9 ST1 and recovered at the production wells RRG-2 and RRG-4, located to the northeast of RRG-9 ST1, support this theory. It is hypothesized that water injected into RRG-9 ST1 moves through the intersecting fracture zone which connects into the Narrows Zone, and then moves along the Narrows Zone to the northeast. This conceptual model was numerical simulated using FALCON a finite element reservoir simulation code developed by Idaho National Laboratory. The modeling suggests that both thermal and hydraulic fracturing mechanisms played an important role in improving the injectivity at RRG-9 ST1. The stimulation program at Raft River has successfully demonstrated the effectiveness of hydraulic and thermal stimulation techniques. RRG-9 ST1 is now in commercial use.
Study of Electromagnetic Scattering From Material Object Doped Randomly With Thin Metallic Wires Using Finite Element Method
A new numerical simulation method using the finite element methodology (FEM) is presented to study electromagnetic scattering due to an arbitrarily shaped material body doped randomly with thin and short metallic wires. The FEM approach described in many standard text books is appropriately modified to account for the presence of thin and short metallic wires distributed randomly inside an arbitrarily shaped material body. Using this modified FEM approach, the electromagnetic scattering due to cylindrical, spherical material body doped randomly with thin metallic wires is studied.
Investigation of Cutter–Woodchip Contact Pressure in a New Biomass Comminution System
A Crumbler ® rotary shear system was recently developed to process fuel-grade woodchips into feedstocks having more uniform particle sizes for improved flowability and higher thermochemical output yield than a traditional hammer mill. It represents a significant innovation in the state of the art of biomass comminution equipment. However, the milling unit has experienced significant abrasive wear of the feeding teeth and cutters when processing hard and dusty feedstocks, such as logging residue and corn stover. Here, we present initial simulation results from a finite-element stress analysis of the rotary shear, which is used to investigate the dynamics of interaction between the cutters and woodchips. This is the first step for optimizing the cutter design to improve the system lifetime and processing efficiency, thus reducing downtime and improving overall productivity.
Statistical of a Characterization and the Simulation of a Reverberation Chamber Using Finite-Element Techniques
The statistical characterization of a simulation reverberation chamber is performed by considering a two-dimensional finite element model. This model includes a source to study the particular modal fields that couple into either a transverse electric or transverse magnetic configuration. The analysis includes a characterization of the basic field statistics, max-to-average ratio, normalized standard deviation, stirring ratio, and field uniformity. The shielding effectiveness of an aperture will be studied that will provide some insight into the nature of the fields coupled from a complex to a noncomplex environment.
Electromagnetics Simulations with Vector-valued Finite Elements in MOOSE
Vector-valued finite elements have become an essential part of modern electromagnetic (EM) numerical simulation. This talk will describe the development and initial validation of a MOOSE-based EM solver application, ELK (Electromagnetics Library for Kinetics and fluids), that utilizes the recently added vector-valued finite element system. The primary motivation for the development of ELK was to study EM wave propagation in radio-frequency (rf) driven systems, such as those encountered in the study of antenna structures, waveguides, and plasma physics. In simulating these scenarios with traditional scalar-valued finite elements, boundary conditions on tangential field components can sometimes be difficult to define, and field discontinuities encountered with certain geometric or structural features (e.g., abrupt corners and inhomogeneous materials) can often lead to non-physical results. Vector-valued elements can largely mitigate and even eliminate these difficulties. ELK verification and validation efforts thus far have focused on standard wave propagation benchmarks for antennas and waveguides, and the impact of vector-valued elements on the implementation and results of these simulations will be described. Progress using ELK for plasma physics simulation, where inhomogeneous simulation domains are commonplace, will also be discussed.