Engineering Papers⌕ Search

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.

At least 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Micro-wrinkling and delamination-induced buckling of stretchable electronic structures

This paper presents the results of experimental and theoretical/computational micro-wrinkles and buckling on the surfaces of stretchable poly-dimethylsiloxane (PDMS) coated with nano-scale Gold (Au) layers. The wrinkles and buckles are formed by the unloading of pre-stretched PDMS/Au structure after the evaporation of nano-scale Au layers. They are then characterized using atomic force microscopy and scanning electron microscopy. The critical stresses required for wrinkling and buckling are analyzed using analytical models. The possible interfacial cracking that can occur along with film buckling is also studied using finite element simulations of the interfacial crack growth. The implications of the results are discussed for potential applications of micro-wrinkles and micro-buckles in stretchable electronic structures and biomedical devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A Correction Scheme for Thermal Conductivity Measurement Using the Comparative Cut-bar Technique Based on a 3D Numerical Simulation

As an important factor affecting the accuracy of the thermal conductivity measurement, systematic (bias) error in the guarded comparative axial heat flow (cut-bar) method was mostly neglected by previous researches. This bias is due primarily to the thermal conductivity mismatch between sample and meter bars (reference), which is common for a sample of unknown thermal conductivity. A correction scheme, based on a finite element simulation of the measurement system, was proposed to reduce the magnitude of the overall measurement uncertainty. This scheme was experimentally validated by applying corrections on four types of sample measurements in which the specimen thermal conductivity is much smaller, slightly smaller, equal and much larger than that of the meter bar. As an alternative to the optimum guarding technique proposed before, the correction scheme can be used to minimize uncertainty contribution from the measurement system with non-optimal guarding conditions. It is especially necessary for large thermal conductivity mismatches between sample and meter bars.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2010

Effect of Microscopic Damage Events on Static and Ballistic Impact Strength of Triaxial Braid Composites

The reliability of impact simulations for aircraft components made with triaxial-braided carbon-fiber composites is currently limited by inadequate material property data and lack of validated material models for analysis. Methods to characterize the material properties used in the analytical models from a systematically obtained set of test data are also lacking. A macroscopic finite element based analytical model to analyze the impact response of these materials has been developed. The stiffness and strength properties utilized in the material model are obtained from a set of quasi-static in-plane tension, compression and shear coupon level tests. Full-field optical strain measurement techniques are applied in the testing, and the results are used to help in characterizing the model. The unit cell of the braided composite is modeled as a series of shell elements, where each element is modeled as a laminated composite. The braided architecture can thus be approximated within the analytical model. The transient dynamic finite element code LS-DYNA is utilized to conduct the finite element simulations, and an internal LS-DYNA constitutive model is utilized in the analysis. Methods to obtain the stiffness and strength properties required by the constitutive model from the available test data are developed. Simulations of quasi-static coupon tests and impact tests of a represented braided composite are conducted. Overall, the developed method shows promise, but improvements that are needed in test and analysis methods for better predictive capability are examined.

Littell, Justin D.↗
DOE OSTI2020

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Effects of Magnetic and Electric Field Uniformity on Coded Aperture Imaging Quality in a Cycloidal Mass Analyzer

Cycloidal mass analyzers are unique sector mass analyzers as they exhibit perfect double focusing, making them ideal for incorporating spatial aperture coding, which can increase the throughput of a mass analyzer without affecting the resolving power. However, the focusing properties of the cycloidal mass analyzer depend on the uniformity of the electric and magnetic fields. In this paper, finite element simulation and charged particle tracing were used to investigate the effect of field uniformity on imaging performance of a cycloidal mass analyzer. For the magnetic field, we evaluate a new permanent magnet geometry by comparing it to a traditional geometry. Results indicate that creating an aperture image in a cycloidal mass spectrometer with the same FWHM as the slit requires less than 1% variation in magnetic field strength along the ion trajectories. The new magnet design, called the opposed dipole magnet, has less than 1% field variation over an area approximately 62 × 65 mm; nearly twice the area available in a traditional design of similar size and weight. This allows ion imaging across larger detector arrays without loss of resolving power. In addition, we compare the aperture imaging quality of a traditionally used cycloidal mass spectrometer electric design with a new optimized design with improved field uniformity. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Effect of lattice orientation on compressive properties of selective laser sintered nylon lattice coupons

As AM lattices become more popular in medical devices, it is important to consider how lattice design parameters affect the mechanical integrity and performance of a device. This research investigated the effect of lattice orientation on the compressive mechanical response of five common lattice geometries: Hexagonal Honeycomb (Hex), Diamond, Voronoi Tessellation Method (VTM), sheet-based Gyroid, and Face-Centered-Cubic (FCC). Samples were tested at two relative densities and eight orientations, printed in Nylon 12 (PA2200) on an EOS P396. The mechanical response was compared to simulated finite element analysis (FEA) response for Hex and FCC lattices. Hex lattices displayed significant orientation dependance. Diamond lattices demonstrated some orientation dependance when at a lower relative density. Gyroid, FCC, and VTM all displayed minimal orientation dependance. Increasing lattice relative density appeared to reduce anisotropy among tested orientations of all lattice geometries except Hex. FEA simulations were able to produce a trend response that matched the mechanical response trends based on the investigated orientations. Furthermore, these results support the need for multi axial lattice mechanical test considerations. It also bolsters the need for lattice FEA models to have a sufficient comparator when identifying “worst-case” lattice test scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A resonant sensor composed of quartz double ended tuning fork and silicon substrate for digital acceleration measurement

Presented in this paper is a micro-resonant acceleration sensor based on the frequency shift of quartz double ended tuning fork (DETF). The structure is silicon substrate having a proof mass supported by two parallel flexure hinges as doubly sustained cantilever, with a resonating DETF located between the hinges. The acceleration normal to the chip plane induces an axial stress in the DETF beam and, in turn, a proportional shift in the beam resonant frequency. Substrate is manufactured by single-crystal silicon for stable mechanical properties and batch-fabrication processes. Electrodes on the four surfaces of DETF beam excite anti-phase vibration model, to balance inner stress and torque and imply a high quality factor. The sensor is simply packaged and operates unsealed in atmosphere for measurements. The tested natural frequency is 36.9 kHz and the sensitivity is 21 Hz/g on a nominally ±100 g device, which is in good agreement with analytical calculation and finite element simulation. The output frequency drifting is less than 0.5 Hz (0.0014% of steady output) within 1 h. The nonlinearity is 0.0019%FS and hysteresis is 0.0026%FS. The testing results confirm the feasibility of combining quartz DETF and silicon substrate to achieve a micro-resonant sensor based on simple processing for digital acceleration measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Correlation between structural and opto-electronic characteristics of crystalline Si microhole arrays for photonic light management

By employing electron paramagnetic resonance spectroscopy, transmission electron microscopy, and optical measurements, we systematically correlate the structural and optical properties with the deep-level defect characteristics of various tailored periodic Si microhole arrays, which are manufactured in an easily scalable and versatile process on nanoimprinted sol-gel coated glass. While tapered microhole arrays in a structured base layer are characterized by partly nanocrystalline features, poor electronic quality with a defect concentration of 10{sup 17} cm{sup −3} and a high optical sub-band gap absorption, planar polycrystalline Si layers perforated with periodic arrays of tapered microholes are composed of a compact crystalline structure and a defect concentration in the low 10{sup 16} cm{sup −3} regime. The low defect concentration is equivalent to the one in planar state-of-the-art solid phase crystallized Si films and correlates with a low optical sub-band gap absorption. By complementing the experimental characterization with 3-dimensional finite element simulations, we provide the basis for a computer-aided approach for the low-cost fabrication of novel high-quality structures on large areas featuring tailored opto-electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Thermal Energy Storage for Dispatchable Geothermal Power

Dispatchable generation refers to sources of electricity that can be turned on or off, or can adjust their power output accordingly to an order. Geothermal plants are usually used for base-load power rather than dispatchable power. Flexibility is a challenge for a geothermal plant because the heat source does not naturally respond quickly to changes in demand. The current state-of-the-art approach to address this challenge is to maintain the geothermal fluid flows at relatively steady rates while providing a bypass around the generation equipment as needed. Under part load, some of the geothermal fluid is pumped to the surface, bypasses the generation equipment, and is re-injected to the ground without extracting any useful enthalpy. This approach is robust but necessarily incurs high parasitic power draws at part load due to the constant, full-flow pumping power requirements. There are even more challenges and unknowns for steam and flash geothermal plants, which were the focus for this project. As a process improvement to this existing configuration of a dispatchable geothermal plant, we examined the use of a phase-change thermal-energy storage (TES) unit that will provide some “inertia” against temperature fluctuations. The TES phase change material (PCM) could be in the form of eutectic molten nitrate salts, organic sugar alcohols, or the like based on the required temperature. This configuration allowed the power plant to meet the ramp rate requirements while allowing the geothermal resource to drawn down at the duty-cycle-averaged extraction rate rather than the full-rated-output extraction rate. The project team performed an experimental and computational evaluation of a peak shaving technical approach that used thermal storage to switch between two nominal operating points: full rated power and 50% power. In the initial requirements definition phase, the team decided to examine the PCM/TES-enhanced geothermal plant in the context of a 40-MW two-stage flash plant operating in mid-Southern California. Steady-state process simulations were used to understand the equipment sizing and the energetics. Dynamic process simulations were used to understand the response rate of the system to changes in requested power output from the utility. The technical approach was also guided by a techno-economic analysis (TEA) and a high-level analysis of the reservoir. Recognizing that performance of the thermal storage was a key risk element, high-fidelity finite-element simulations of the PCM-TES heat exchanger were performed to better understand its performance and temporal response. These simulations were validated with a subscale PCM-TES experimental apparatus that mimicked the geometry of the full-scale system. The key outcomes of this Phase I project were that: 1) the PCM-TES heat exchanger met the target for gravimetric energy storage density, 2) the dispatchable geothermal plant met its target for response rate, 3) the peak shaving approach reduces the number of production wells and does not disturb the wellfield, and 4) the incentives for properly remunerating dispatchable operation were identified.

15 GEOTHERMAL ENERGY↗
NASA NTRS1997

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Qualitative measurements of adhesion or binding forces can be accomplished, for example, by using the reflection coefficient of an ultrasound or by using thermal waves (Light and Kwun, 1989, Achenbach and Parikh, 1991, and Bostrom and wickham, 1991). However, a quantitative determination of binding forces is rather difficult. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasound passes through a nonlinear material. It seems that such non-linearity can be effectively used to characterize the bond strength. Several theories have been developed to model this nonlinear effect (Adler and Nagy, 1991; Achenbach and Parikh, 1991; Parikh and Achenbach, 1992; and Hirose and Kitahara, 1992; Anastasi and Roberts, 1992). Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented by Pangraz and Arnold (1994). Recently, Tang, Cheng and Achenbach (1997) made a comparison between the experimental and simulated results based on this theoretical model. A water immersion mode-converted shear wave through-transmission setup was used by Berndt and Green (1997) to analyze the nonlinear acoustic behavior of the adhesive bond. In this project, the nonlinear responses of an adhesive joint was investigated through transmission tests of ultrasonic wave and analyzed by the finite element simulations. The higher order harmonics were obtained in the tests. It is found that the amplitude of higher harmonics increases as the aging increases, especially the 3dorder harmonics. Results from the numerical simulation show that the material nonlinearity does indeed generate higher order harmonics. In particular, the elastic-perfect plastic behavior generates significant 3rd and 5th order harmonics.

Qu, Jianmin↗
DOE OSTI2025

Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas

Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2014

Edge Fracture Prediction ofTraditional and Advanced Trimming Processes for AA6111-T4 Sheets

This work examines the traditional and advanced trimming of AA6111-T4 aluminum sheets with finite element simulations. The Rice-Tracy damage model is used for the simulation with damage parameters estimated from experimental observation of grain aspect ratio near the fracture surface of trimmed parts. Fine meshes at the shearing zone, adaptive meshing, and adaptive contact techniques are used to accurately capture the contact interactions between the sharp corner of the trimming tools and the blank to be trimmed. To the knowledge of the authors, these are the first trimming simulations that can predict the effects of shearing clearance on burr heights with quantitative accuracy for AA6111-T4 aluminum sheets. In addition, the models have also accurately reproduced the crack initiation site as well as burr and sliver formation mechanisms observed experimentally.

Adaptive Contact.↗
DOE OSTI2022

Dopant Compensation within the Intrinsic Poly-Si Isolation Region in Poly-Si/SiOx Passivated IBC Si Solar Cells

We report on the effect of dopant compensation within intrinsic poly-Si regions between p- and n-type fingers of poly-Si/SiOx passivated interdigitated back contact (IBC) solar cells using intrinsic poly-Si as the isolation region between the doped poly-Si fingers. First, we show that dopants from the doped fingers contaminate the intrinsic gap, resulting in doping of the entire intrinsic gap and overlap of the dopant tails from each finger. Next, we show that despite this doping across the gap, shunting between the doped fingers does not occur. We show that this is a result of trap-assisted compensation creating a highly resistive intrinsic region, preventing shunt. We simulate shunt resistance across the gap based on local carrier concentration and deep trap density. We show that trap defects within the poly-Si enhance compensation between the dopant tails. We experimentally confirm these predictions by scanning spreading resistance microscopy of the gap showing ~20 um domain with resistivity ~10^7 ..omega.. cm. Additionally, Kelvin probe force microscopy are compared to finite element simulations which result in the same approximate shape for potential profile, indicating diode behavior across the isolation region. These results demonstrate the powerful effect that trap defects have within the poly-Si isolation region and suggest that precision patterning is not as essential as once thought.

dopant compensation↗
DOE OSTI2019

Numerical Investigation of the Effect of Rolling on the Localized Stress and Strain Induction for Wire + Arc Additive Manufactured Structures

Cold rolling can be used in-process or post-process to improve microstructure, mechanical properties and residual stress in directed-energy-deposition techniques, such as the high deposition rate wire + arc additive manufacturing (WAAM) process. Finite element simulations of the rolling process are employed to investigate the effect of rolling parameters, in particular rolling load and roller profile radius on the residual stress field as well as plastic strain distribution for the profiled roller. The results show the response to rolling of commonly used structural metals in WAAM, i.e., AA2319, S335JR steel and Ti-6Al-4V, taking into account the presence of residual stresses. The rolling load leads to changes in the location and the maximum value of the compressive residual stresses, as well as the depth of the compressive residual stresses. However, the roller profile radius only changes the maximum value of these compressive residual stresses. Changing the rolling load influences the equivalent plastic strain close to the top surface of the wall as well as in deeper areas, whereas the influence of the roller profile radius is negligible. The plastic strain distribution is virtually unaffected by the initial residual stresses prior to rolling. Finally, design curves were generated from the simulations for different materials, suggesting ideal rolling load and roller profile combinations for microstructural improvement requiring certain plastic strains at a specific depth of the additive structure.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Application of the Quadrupole Method for Simulation of Passive Thermography

Passive thermography has been shown to be an effective method for in-situ and real time nondestructive evaluation (NDE) to measure damage growth in a composite structure during cyclic loading. The heat generation by subsurface flaw results in a measurable thermal profile at the surface. This paper models the heat generation as a planar subsurface source and calculates the resultant temperature profile at the surface using a three dimensional quadrupole. The results of the model are compared to finite element simulations of the same planar sources and experimental data acquired during cyclic loading of composite specimens.

Winfree, William P.↗
DOE OSTI2016

Modeling and simulation framework for dynamic strain localization in elasto-viscoplastic metallic materials subject to large deformations

This study describes a theoretical and computational framework for the treatment of adiabatic shear band formation in rate-sensitive polycrystalline metallic materials. From a computational perspective, accurate representation of strain localization behavior has been a long-standing challenge. In addition, the underlying physical mechanisms leading to the localization of plastic deformation are still not fully understood. The proposed framework is built around an enhanced-strain finite element formulation, designed to alleviate numerical pathologies known to arise in localization problems, by allowing a localization band of given finite width (weak discontinuity) to be embedded within individual elements. The mechanical threshold strength (MTS) model is used to represent the temperature and strain rate-dependent viscoplastic response of the material. This classical flow stress model employs an internal state variable to quantify the effect of dislocation structure evolution (work hardening and recovery). In light of growing evidence suggesting that the softening effect of dynamic recrystallization may play a significant role, alongside thermal softening, in the process of shear band formation and growth, a simple dynamic recrystallization model is proposed and cast within the context of the MTS model with the aid of the aforementioned internal state variable. An initiation criterion for shear localization in rate and temperature-sensitive materials is introduced and used in the present context of high-rate loading, where material rate-dependence is pronounced and substantial temperature increases are achieved due to the dissipative nature of viscoplastic processes. In addition, explicit time integration is adopted to facilitate treatment of the dynamic problems under consideration, where strain rates in excess of 10 4 s –1 are typically attained. Two series of experiments are conducted on AISI 316L stainless steel, employing the commonly used top-hat sample geometry and the Split-Hopkinson Pressure Bar dynamic test system. Axi-symmetric finite element simulation results are compared to cross-sectional micrographs of recovered samples and experimental load–displacement results, in order to examine the performance of the proposed framework and demonstrate its effectiveness in treating the initiation and growth of adiabatic shear banding in dynamically loaded metallic materials. These comparisons demonstrate that thermal softening alone is insufficient to induce shear localization behaviors observed in some materials, such as stainless steel, and support the hypothesis that dynamic recrystallization and/or other softening mechanisms play an essential role in this process.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗