Search NASA⌕ Search

SEARCH · Search NASA

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Investigation of Micro-Scale Architectural Effects on Damage of Composites

This paper presents a three-dimensional, energy based, anisotropic, stiffness reduction, progressive damage model for composite materials and composite material constituents. The model has been implemented as a user-defined constitutive model within the Abaqus finite element software package and applied to simulate the nonlinear behavior of a damaging epoxy matrix within a unidirectional composite material. Three different composite microstructures were considered as finite element repeating unit cells, with appropriate periodicity conditions applied at the boundaries. Results representing predicted transverse tensile, longitudinal shear, and transverse shear stress-strain curves are presented, along with plots of the local fields indicating the damage progression within the microstructure. It is demonstrated that the damage model functions appropriately at the matrix scale, enabling localization of the damage to simulate failure of the composite material. The influence of the repeating unit cell geometry and the effect of the directionality of the applied loading are investigated and discussed.

damage↗
DOE OSTI2020

Microwave Cavity Simulation Using Ansys HFSS

The design of microwave cavity detectors for axion dark matter research is often accomplished using advanced full-wave electromagnetic simulation software tools. These tools provide a cost-effective approach to evaluate a wide variety of cavity configurations, frequency tuning mechanisms, and conductive or dielectric materials and coatings. One simulation software package used for this application is Ansys High Frequency Structure Simulator (HFSS), which is based upon the well-established finite element method. HFSS includes numerous features useful for microwave cavity design such as parametric geometry modeling, adaptive meshing algorithm, curvilinear mesh elements, driven modal and eigenmode matrix solvers, and optimization algorithms. This paper describes the use of the HFSS software to simulate microwave cavities for axion haloscope detectors, with an example tutorial for a cylindrical cavity. Excellent agreement between the simulated and analytical results is shown for the resonant frequency, quality factor, and form factor.

Ansys HFSS, Cavity simulation, finite element mode↗
NASA NTRS1994

Orientation effects on the measurement and analysis of critical CTOA in an aluminum alloy sheet

Fracture tests were conducted on 76.2mm wide, 2.3mm thick middle crack tension (M(T)) specimens machined from 2024-T3 aluminum sheet. The specimens were tested on the T-L orientation and comparisons were made to similar tests conducted in the L-T orientation. Measurement of critical crack tip opening angle (CTOA), applied stress, and crack front shape were made as a function of crack extension. A two-dimensional, elastic-plastic finite element analysis was used to simulate the fracture behavior for both orientations. The results indicate that the T-L orientation had a 10 percent lower stress at fracture than similar tests conducted in the L-T orientation. Correspondingly, the critical CTOA in the T-L tests reached a constant value of 4.7 degrees after 2-3mm of crack extension and the L-T tests reached a value of 6 degrees. The fracture surfaces of the T-L specimens were observed to remain flat, while those of the L-T specimens transitioned to a 45 degree slant fracture after about 2-3mm of crack extension. The tunneling behavior of the two orientations also differed; the T-L specimens reached a deeply tunneled stabilized crack front shape while, the L-T specimens were observed to have only a small amount of tunneling once the crack began to grow on the 45 degree slant. The two-dimensional, elastic-plastic finite element analysis was able to simulate the fracture behavior for both the T-L and L-T orientations.

Sutton, M. A.↗
DOE OSTI2022

Measurement and Simulation of the Magnetic Fields from a 555 Timer Integrated Circuit Using a Quantum Diamond Microscope and Finite-Element Analysis

Quantum diamond microscope (QDM) magnetic field imaging is an emerging interrogation and diagnostic technique for integrated circuits (ICs). To date, the ICs measured with a QDM have been either too complex for us to predict the expected magnetic fields and benchmark the QDM performance or too simple to be relevant to the IC community. In this paper, we establish a 555 timer IC as a “model system” to optimize QDM measurement implementation, benchmark performance, and assess IC device functionality. To validate the magnetic field images taken with a QDM, we use a spice electronic circuit simulator and finite-element analysis (FEA) to model the magnetic fields from the 555 die for two functional states. Furthermore, we compare the advantages and the results of three IC-diamond measurement methods, confirm that the measured and simulated magnetic images are consistent, identify the magnetic signatures of current paths within the device, and discuss using this model system to advance QDM magnetic imaging as an IC diagnostic tool.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1997

Thermal History and Mantle Dynamics of Venus

One objective of this research proposal is to develop a 3-D thermal history model for Venus. The basis of our study is a finite-element computer model to simulate thermal convection of fluids with highly temperature- and pressure-dependent viscosities in a three-dimensional spherical shell. A three-dimensional model for thermal history studies is necessary for the following reasons. To study planetary thermal evolution, one needs to consider global heat budgets of a planet throughout its evolution history. Hence, three-dimensional models are necessary. This is in contrasts to studies of some local phenomena or local structures where models of lower dimensions may be sufficient. There are different approaches to treat three-dimensional thermal convection problems. Each approach has its own advantages and disadvantages. Therefore, the choice of the various approaches is subjective and dependent on the problem addressed. In our case, we are interested in the effects of viscosities that are highly temperature dependent and that their magnitudes within the computing domain can vary over many orders of magnitude. In order to resolve the rapid change of viscosities, small grid spacings are often necessary. To optimize the amount of computing, variable grids become desirable. Thus, the finite-element numerical approach is chosen for its ability to place grid elements of different sizes over the complete computational domain. For this research proposal, we did not start from scratch and develop the finite element codes from the beginning. Instead, we adopted a finite-element model developed by Baumgardner, a collaborator of this research proposal, for three-dimensional thermal convection with constant viscosity. Over the duration supported by this research proposal, a significant amount of advancements have been accomplished.

Hsui, Albert T.↗
NASA NTRS1995

Effect of various approximations on predicted progressive failure in plain weave composites

Three-dimensional finite element analysis was used to simulate progressive failure of a plain weave composite subjected to in-plane extension. The loading was parallel to one of the tow directions. The effects of various characteristics of the finite element model on predicted behavior were examined. The predicted behavior was found to be sensitive to quadrature order, mesh refinement, and the material degradation model. Also the sensitivity of the predictions to the tow waviness was studied. The predicted strength decreased considerably with increased waviness. More numerical studies and comparisons with experimental data are needed to establish reliable guidelines for accurate progressive failure prediction.

Whitcomb, John↗
DOE OSTI2018

Determination of Thermoelastic Properties of Carbon Nanotube/Epoxy Composites Using Finite Element Method

In the present study, three-dimensional finite element (FE) analysis has been performed for determination of thermoelastic properties of carbon nanotube (CNT)/epoxy composites considering both linear and nonlinear stress–strain relations of CNT. A square representative volume element has been considered for the FE analysis (FEA), and the effect of considering nonlinear stress–strain relation of CNT on the prediction of thermoelastic properties has been determined by comparing the results with the theoretical results. Results show that in order to obtain realistic values of thermoelastic properties, the nonlinear stress–strain relation of CNT needs to be used in the FEA. Not considering nonlinear stress–strain relation of CNTs in the FEA may lead to estimating erroneous values of thermoelastic properties of CNT composites which may have adverse effect in the design of structural components made of such composites subjected to thermomechanical loading.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Crack Growth Simulation and Residual Strength Prediction in Airplane Fuselages

This is the final report for the NASA funded project entitled "Crack Growth Prediction Methodology for Multi-Site Damage." The primary objective of the project was to create a capability to simulate curvilinear fatigue crack growth and ductile tearing in aircraft fuselages subjected to widespread fatigue damage. The second objective was to validate the capability by way of comparisons to experimental results. Both objectives have been achieved and the results are detailed herein. In the first part of the report, the crack tip opening angle (CTOA) fracture criterion, obtained and correlated from coupon tests to predict fracture behavior and residual strength of built-up aircraft fuselages, is discussed. Geometrically nonlinear, elastic-plastic, thin shell finite element analyses are used to simulate stable crack growth and to predict residual strength. Both measured and predicted results of laboratory flat panel tests and full-scale fuselage panel tests show substantial reduction of residual strength due to the occurrence of multi-site damage (MSD). Detailed comparisons of n stable crack growth history, and residual strength between the predicted and experimental results are used to assess the validity of the analysis methodology. In the second part of the report, issues related to crack trajectory prediction in thin shells; an evolving methodology uses the crack turning phenomenon to improve the structural integrity of aircraft structures are discussed, A directional criterion is developed based on the maximum tangential stress theory, but taking into account the effect of T-stress and fracture toughness orthotropy. Possible extensions of the current crack growth directional criterion to handle geometrically and materially nonlinear problems are discussed. The path independent contour integral method for T-stress evaluation is derived and its accuracy is assessed using a p- and hp-version adaptive finite element method. Curvilinear crack growth is simulated in coupon tests and in full-scale fuselage panel tests. Both T-stress and fracture toughness orthotropy are found to be essential to predict the observed crack paths. The analysis methodology and software program (FRANC3D/STAGS) developed herein allows engineers to maintain aging aircraft economically while insuring continuous airworthiness. Consequently, it will improve the technology to support the safe operation of the current aircraft fleet as well as the design of more damage-tolerant aircraft for the next generation fleet.

Chen, Chuin-Shan↗
DOE OSTI2018

Parallel Element Agglomeration Multilevel Monte Carlo Library (ParELAGMC)

ParELAGMC is a parallel distributed memory C++ library for multilevel Monte Carlo (MLMC) simulations with algebraically constructed coarse spaces for estimating statistical moments of a quantity of interest involving solutions of partial differential equations with uncertain coefficients. A primary focus of the library is the scalable generation of Gaussian random fields using a novel stochastic PDE (SPDE) sampling technique for large-scale simulations. A key feature in the SPDE sampling technique is a novel non-matching mesh embedding technique to mitigate artificial boundary artifacts that are inherent in the method. Example applications to uncertainty quantification of subsurface flow simulations in the mixed finite element setting are provided as well.

Osborn, SarahV.↗
DOE OSTI2021

CARDINAL

Cardinal is a code that automates the workflow for iteratively solving neutral particle transport via the OpenMC code and computational fluid dynamics via the NekRS code within the Multiphysics Object Oriented Simulation Environment (MOOSE), a finite element framework for solving general partial differential equations. Cardinal establishes a mapping between the Monte Carlo geometry representation used in OpenMC and the spectral element mesh used in NekRS to other physics applications built upon the MOOSE framework, allowing high-resolution particle transport and fluid dynamics physics feedback to the MOOSE application ecosystem. Heat generation rates as solved by OpenMC are used as heat sources in MOOSE, while the temperature and density fields from NekRS are used as conjugate heat transfer boundary conditions and/or source fields in MOOSE. OpenMC, Nek5000, and MOOSE are open source, community-developed simulation codes.

NOVAK, APRILJEAN RODGERS↗
DOE OSTI2018

Effect of interfacial transition zone on the Young's modulus of carbon nanofiber reinforced cement concrete

Identifying the properties of the region where the cement paste meets the aggregate surface (interfacial transition zone, ITZ) is critical to understanding the strength and fracture behavior of carbon nanofiber (CNF) reinforced cement concrete. In this study, the finite element method is employed to investigate the effect of the ITZ on the Young's modulus of cement concrete made with CNF. The numerical models for cement concrete with and without CNF are constructed based on the digital image processing technique. To consider the interface effect, the concept of ''effective aggregate'' is put forward, namely, the Young's modulus and Poisson's ratio for each aggregate particle are replaced by an effective Young's modulus and effective Poisson's ratio, in which the effect of thickness and the Young's modulus of ITZ is taken into account in an averaged manner. Then, the quantitative nanomechanical mapping (based on atomic force microscopy) technique is adopted to measure the Young's modulus and thickness of ITZ with and without CNFs, which are further used as the input parameters in the numerical model. The numerical simulation results are verified by experimental testing, which indicates that the ITZ effect should be considered when implementing the numerical simulation. In addition, this analysis shows that compared with the plain cement concrete, CNFs can greatly enhance the mechanical properties of ITZ, which will in turn improve the Young's modulus of cement concrete significantly. Finally, the effect of thickness and Young's modulus of ITZ on the Young's modulus of CNF reinforced cement concrete is also discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1998

Parallel Three-Dimensional Computation of Fluid Dynamics and Fluid-Structure Interactions of Ram-Air Parachutes

This is a final report as far as our work at University of Minnesota is concerned. The report describes our research progress and accomplishments in development of high performance computing methods and tools for 3D finite element computation of aerodynamic characteristics and fluid-structure interactions (FSI) arising in airdrop systems, namely ram-air parachutes and round parachutes. This class of simulations involves complex geometries, flexible structural components, deforming fluid domains, and unsteady flow patterns. The key components of our simulation toolkit are a stabilized finite element flow solver, a nonlinear structural dynamics solver, an automatic mesh moving scheme, and an interface between the fluid and structural solvers; all of these have been developed within a parallel message-passing paradigm.

Tezduyar, Tayfun E.↗
NASA NTRS1995

Analysis of a 170-GHz frequency doubler with an array of planar diodes

Analysis of a planar diode multiplier from 85-170 GHz is described. The doubler uses a waveguide mount with two series pairs of diodes in a balanced structure. Because of the difficulties in conventional scale model measurements, numerical electromagnetic simulation based on the finite element method was chosen for the analysis, using a commercial program. To optimize the diode design, the de-embedded diode terminal impedance was studied, as well as the power balance between the diodes. The analysis showed that the matching of the diode impedance to that of the waveguide is quite sensitive to the diode substrate thickness. Thicknesses from 25-100 micron in GaAs were studied as well as 100-micron-thick quartz. The accuracy of the theoretical analysis was verified by careful measurements using a slotted line to determine the diode terminal impedance under large signal pump, for frequencies between 80 and 90 GHz. Good agreement between the measured and simulated diode terminal impedance was observed, although full agreement requires the addition of an empirical loss term. Several options are considered for the source of this loss.

DIODE MULTIPLIER↗
DOE OSTI2020

Validation of high-fidelity ion cyclotron range of frequencies antenna coupling simulations in full 3D geometry against experiments in the ASDEX Upgrade tokamak

We address the validation of finite elements ion cyclotron range of frequencies (ICRF) antenna coupling simulation against experiments performed in the ASDEX Upgrade tokamak. Measurements of the loading resistance in ICRF-heated, magnetically-perturbed 3D plasma discharges are compared against numerical predictions of the RAPLICASOL code. To this end, the 3D induction field and the 3D density profile are modeled by concatenating the PARVMEC, BMW and EMC3-EIRENE codes. The 3D density is input to RAPLICASOL, where full-wave simulations are performed on a finite element mesh retaining full 3D geometry in the ICRF antenna model and the plasma description. The results are further compared with RAPLICASOL simulations using a 1D density profile as measured at the outboard midplane in the same experiments. We find that simulations using a 1D density profile overestimate the change in loading resistance by a factor of ~197 -248%, while simulations using the full 3D density profile are in agreement with experiments within a factor ~16%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluation of accuracy and convergence of numerical coupling approaches for poroelasticity benchmark problems

Accurate modeling of subsurface flow and transport processes is vital as the prevalence of subsurface activities such as carbon sequestration, geothermal recovery, and nuclear waste disposal increases. Computational modeling of these problems leverages poroelasticity theory, which describes coupled fluid flow and mechanical deformation. Although fully coupled monolithic schemes are accurate for coupled problems, they can demand significant computational resources for large problems. In this work, a fixed stress scheme is implemented into the Sandia Sierra Multiphysics toolkit. Two implementation methods, along with the fully coupled method, are verified with one-dimensional (1D) Terzaghi, 2D Mandel, and 3D Cryer sphere benchmark problems. The impact of a range of material parameters and convergence tolerances on numerical accuracy and efficiency was evaluated. Overall the fixed stress schemes achieved acceptable numerical accuracy and efficiency compared to the fully coupled scheme. However, the accuracy of the fixed stress scheme tends to decrease with low permeable cases, requiring the finer tolerance to achieve a desired numerical accuracy. For the fully coupled scheme, high numerical accuracy was observed in most of cases except a low permeability case where an order of magnitude finer tolerance was required for accurate results. Finally, a two-layer Terzaghi problem and an injection–production well system were used to demonstrate the applicability of findings from the benchmark problems for more realistic conditions over a range of permeability. Simulation results suggest that the fixed stress scheme provides accurate solutions for all cases considered with the proper adjustment of the tolerance. Furthermore, this work clearly demonstrates the robustness of the fixed stress scheme for coupled poroelastic problems, while a cautious selection of numerical tolerance may be required under certain conditions with low permeable materials.

58 GEOSCIENCES↗
DOE OSTI2015

STRESS AND FAILURE ANALYSIS OF RAPIDLY ROTATING ASTEROID (29075) 1950 DA

Rozitis et al. recently reported that near-Earth asteroid (29075) 1950 DA, whose bulk density ranges from 1.0 g cm{sup –3} to 2.4 g cm{sup –3}, is a rubble pile and requires a cohesive strength of at least 44-76 Pa to keep from failing due to its fast spin period. Since their technique for giving failure conditions required the averaged stress over the whole volume, it discarded information about the asteroid's failure mode and internal stress condition. This paper develops a finite element model and revisits the stress and failure analysis of 1950 DA. For the modeling, we do not consider material hardening and softening. Under the assumption of an associated flow rule and uniform material distribution, we identify the deformation process of 1950 DA when its constant cohesion reaches the lowest value that keeps its current shape. The results show that to avoid structural failure the internal core requires a cohesive strength of at least 75-85 Pa. It suggests that for the failure mode of this body, the internal core first fails structurally, followed by the surface region. This implies that if cohesion is constant over the whole volume, the equatorial ridge of 1950 DA results from a material flow going outward along the equatorial plane in the internal core, but not from a landslide as has been hypothesized. This has additional implications for the likely density of the interior of the body.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

HPC for the EEC: Industrial Trade Tools for the Aging Energy Infrastructure

The petrochemical and refining sectors are challenged with reliably delivering safer and cleaner energy to US consumers, while meeting an ever-growing global demand. The petrochemical and refining infrastructure in the USA is aging, and corrosion and damage mechanisms are constant threats to mechanical integrity, safety, and profitability. However, governments and industry stakeholders are reluctant to replace or upgrade the existing infrastructure due to the immense cost. To better understand and mitigate the risks of aging energy infrastructure, strong technical analysis capabilities, combined with optimized monitoring and decision making, is critical. Today, this is accomplished via complex simulations and data analysis. To this end, advanced High-Performance-Computing (HPC) software will be leveraged and integrated in easy-to-use industrial trade tools, to lower the barrier for new users, increase the ease of access for experienced users, and allow for smarter decisions to be made in the midstream and downstream energy sectors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗