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.

SourceResearch connectionFrom yearThrough year

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

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Simulation of the temperature distribution in the selective beam melting process for polymer material

In the present contribution the temperature distribution in the selective beam melting process for polymer materials is simulated to better understand the influence of process parameters on the properties of the produced part. The basis for the developed simulation tool is the nonlinear heat equation including temperature dependent functions for the heat capacity and the heat conduction which were obtained by experimental measurements. The effect of latent heat occurring in the process is also taken into account. The heat equation is discretized in time and space where a Runge-Kutta method of Radau IIA type is used for time integration. An adaptive finite element method is applied for the discretization in space and the model is implemented into the finite element library deal.II. The heat and cooling rate as important process parameters are simulated for different beam velocities. The ability for computing these process parameters makes the simulation tool suited for optimizing the process management of selective beam melting plants.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2002

Occupant Responses in a Full-Scale Crash Test of the Sikorsky ACAP Helicopter

A full-scale crash test of the Sikorsky Advanced Composite Airframe Program (ACAP) helicopter was performed in 1999 to generate experimental data for correlation with a crash simulation developed using an explicit nonlinear, transient dynamic finite element code. The airframe was the residual flight test hardware from the ACAP program. For the test, the aircraft was outfitted with two crew and two troop seats, and four anthropomorphic test dummies. While the results of the impact test and crash simulation have been documented fairly extensively in the literature, the focus of this paper is to present the detailed occupant response data obtained from the crash test and to correlate the results with injury prediction models. These injury models include the Dynamic Response Index (DRI), the Head Injury Criteria (HIC), the spinal load requirement defined in FAR Part 27.562(c), and a comparison of the duration and magnitude of the occupant vertical acceleration responses with the Eiband whole-body acceleration tolerance curve.

Jackson, Karen E.↗

NASA NTRS2002

Occupant Responses in a Full-Scale Crash Test of the Sikorsky ACAP Helicopter

Jackson, Karen E.↗

DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

NASA NTRS1992

Computer Simulation For Design Of TWT's

A three-dimensional finite-element analytical technique facilitates design and fabrication of traveling-wave-tube (TWT) slow-wave structures. Used to perform thermal and mechanical analyses of TWT designed with variety of configurations, geometries, and materials. Using three-dimensional computer analysis, designer able to simulate building and testing of TWT, with consequent substantial saving of time and money. Technique enables detailed look into operation of traveling-wave tubes to help improve performance for future communications systems.

Bartos, Karen F.↗

DOE OSTI2022

Pronghorn Porous Media Model Validation with Pressure Drop Measurements

The verification and validation (V&V) of Pronghorn is imperative to assert its accuracy when predicting the fluid velocity, temperature, and pressure in high temperature gas-cooled reactors. Pronghorn is a coarse-mesh, intermediate-fidelity, and multidimensional thermal-hydraulics (TH) code developed by the Idaho National Laboratory (INL). New pebble bed experiments are used to observe the details of the fluid motion and pressure drop in the porous bed under the reactor normal operation. This paper focuses on the validation of the Pronghorn compressible and incompressible Navier-Stokes equations using the pressure drop measurements performed at the engineering-scale pebble bed facility at the Texas A&M university (TAMU). Various pressure drop correlations and porosity functions are implemented in both Pronghorn and STAR-CCM+ to compare the pressure drop due to the combined viscous and inertial resistances in the porous bed. The correlations accounting for the near-wall effect are also utilized to observe if the pressure drop estimates can be improved. Pronghorn porous media models predict the pressure drop well relative to the STAR-CCM+ simulation results and 1D correlations, and both the finite element method (FEM) and finite volume method (FVM) perform accurately. Pronghorn models are also validated with the experimental measurements given the different Reynolds number ranges and specific aspect ratios. The likelihood of the statistical significance between the pressure drop measurements and specific correlations or simulations is low provided that the overlap of their confidence intervals is more than the half of a single arm. Several validation metrics are reasonable in regard to the similar studies from other literature. The precise average pebble bed porosity estimation has much impact on the pressure drop, and the Foumeny and Montillet (dense packing) models carry out the accurate pressure drop prediction by considering the near-wall effect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2015

Near field fluid coupling between internal motion of the organ of Corti and the basilar membrane

The pressure distribution in each of the fluid chambers of the cochlea can be decomposed into a 1D, or plane wave, component and a near field component, which decays rapidly away from the excitation point. The transverse motion of the basilar membrane, BM, for example, generates both a 1D pressure field, which couples into the slow wave, and a local near field pressure, proportional to the BM acceleration, that generates an added mass on the BM due to the fluid motion. When the organ of Corti, OC, undergoes internal motion, due for example to outer hair cell activity, this motion will not itself generate any 1D pressure if the OC is incompressible and the BM is constrained not to move volumetrically, and so will not directly couple into the slow wave. This motion will, however, generate a near field pressure, proportional to the OC acceleration, which will act on the OC and thus increases its effective mass. The near field pressure due to this OC motion will also act on the BM, generating a force on the BM proportional to the acceleration of the OC, and thus create a “coupling mass” effect. By reciprocity, this coupling mass is the same as that acting on the OC due to the motion of the BM. This near field fluid coupling is initially observed in a finite element model of a slice of the cochlea. These simulations suggest a simple analytical formulation for the fluid coupling, using higher order beam modes across the width of the cochlear partition. It is well known that the added mass due to the near field pressure dominates the overall mass of the BM, and thus significantly affects the micromechanical dynamics. This work not only quantifies the added mass of the OC due its own motion in the fluid, and shows that this is important, but also demonstrates that the coupling mass effect between the BM and OC significantly affects the dynamics of simple micromechanical models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1998

Free-Suspension Residual Flexibility Testing of Space Station Pathfinder: Comparison to Fixed-Base Results

Application of the free-suspension residual flexibility modal test method to the International Space Station Pathfinder structure is described. The Pathfinder, a large structure of the general size and weight of Space Station module elements, was also tested in a large fixed-base fixture to simulate Shuttle Orbiter payload constraints. After correlation of the Pathfinder finite element model to residual flexibility test data, the model was coupled to a fixture model, and constrained modes and frequencies were compared to fixed-base test. modes. The residual flexibility model compared very favorably to results of the fixed-base test. This is the first known direct comparison of free-suspension residual flexibility and fixed-base test results for a large structure. The model correlation approach used by the author for residual flexibility data is presented. Frequency response functions (FRF) for the regions of the structure that interface with the environment (a test fixture or another structure) are shown to be the primary tools for model correlation that distinguish or characterize the residual flexibility approach. A number of critical issues related to use of the structure interface FRF for correlating the model are then identified and discussed, including (1) the requirement of prominent stiffness lines, (2) overcoming problems with measurement noise which makes the antiresonances or minima in the functions difficult to identify, and (3) the use of interface stiffness and lumped mass perturbations to bring the analytical responses into agreement with test data. It is shown that good comparison of analytical-to-experimental FRF is the key to obtaining good agreement of the residual flexibility values.

Tinker, Michael L.↗

DOE OSTI2019

Yearly report for the period Jan. 2018 – Feb. 2019 IC Project: w17_faultprediction “Critical Stress in Earth”

The consequence of earthquake is usually catastrophic. Currently, we cannot determine when a fault is approaching failure due to either anthropogenic activities or tectonic movement. In an anthropogenically perturbed setting, the rate of physical changes on faults far exceeds most natural processes with currently unpredictable and unmeasurable results. In a tectonic setting we generally do not know where a fault is in the earthquake cycle beyond noting how long has passed since the last rupture, which is itself sometimes unknown. The Laboratory has therefore funded a multi million multi-year LDRD-DR program that comprises experiments, numerical modeling, and machine learning in order to develop and test a highly novel approach to probe the earthquake cycle and in particular, to detect and locate stressed faults that are approaching failure. For this program, modeling of a series of sheared granular layers undergoing stick-slip behavior using the EES-17’s Combined Finite-Discrete Element Method (FDEM) based code HOSS (Hybrid Optimization Software Suite) plays a pivotal role since numerical simulations of granular fault gouge allow for analysis of the mechanical behavior of the system at a level of spatial and temporal resolution not accessible experimentally, and also offer the possibility of identifying the most important parameters governing earthquake dynamics. The simulated data also provide very detailed stress, strain and motion information for machine learning purpose to predict earthquake. The code – HOSS – has been developed in Geophysics Group in LANL during the past ten years and is used to simulate the stick-slip behavior of granular fault gouge. One of the key feature of HOSS is that it is superior to pure finite element method and discrete element method and can provide very detailed evolvement information of granular fault gouge. Additionally, its parallel computing capability makes large scale earthquake gouge simulation possible. Taking these facts into account, and also because few other fully parallel FDEM code exists, this IC project allows LANL to stay at the forefront of earthquake gouge and granular material simulation in the world.

42 ENGINEERING↗

DOE OSTI2023

Neural network ensembles and uncertainty estimation for predictions of inelastic mechanical deformation using a finite element method-neural network approach

The finite element method (FEM) is widely used to simulate a variety of physics phenomena. Approaches that integrate FEM with neural networks (NNs) are typically leveraged as an alternative to conducting expensive FEM simulations in order to reduce the computational cost without significantly sacrificing accuracy. However, these methods can produce biased predictions that deviate from those obtained with FEM, since these hybrid FEM-NN approaches rely on approximations trained using physically relevant quantities. In this work, an uncertainty estimation framework is introduced that leverages ensembles of Bayesian neural networks to produce diverse sets of predictions using a hybrid FEM-NN approach that approximates internal forces on a deforming solid body. The uncertainty estimator developed herein reliably infers upper bounds of bias/variance in the predictions for a wide range of interpolation and extrapolation cases using a three-element FEM-NN model of a bar undergoing plastic deformation. This proposed framework offers a powerful tool for assessing the reliability of physics-based surrogate models by establishing uncertainty estimates for predictions spanning a wide range of possible load cases.

42 ENGINEERING↗

NASA NTRS2018

High-Order Stabilized Finite Elements on Dynamic Meshes

The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.

Anderson, W. Kyle↗

DOE OSTI2017

Simulation of Laser Ultrasonics for Detection of Surface-Connected Rail Defects

Laser ultrasonic produces frequencies in the MHz range, enabling high accuracy and a strong ability to detect rail surface defects. This paper mainly studied on the simulation of detecting surface-connected rail defects on 60 kg rails with laser ultrasonic, established the finite element model of laser-excited ultrasonic Rayleigh wave, carried out the simulation, and verified the effectiveness of the technology through experiments. To solve the problem that laser ultrasonic is insensitive to the width of defects in actual detection, and unable to make quantitative detection of defects, this paper established a new model on the basis of improving the original model that has been verified, exciting ultrasonic at the two sides at the same time of a rail with two staggered beams of laser separately to detect irregular scratch defects on rail surface, and two groups of signal data were received through two probes. Each group of data can present the half-profile information of defects, and further form two detection images of the defect. At last, the two detection images were combined into a complete image through image processing. The results of the experiment indicate that the technology studied offers a new method for the effective quantitative detection of surface-connected defects on rail.

42 ENGINEERING↗

NASA NTRS1999

Adaptive-Grid Methods for Phase Field Models of Microstructure Development

In this work the authors show how the phase field model can be solved in a computationally efficient manner that opens a new large-scale simulational window on solidification physics. Our method uses a finite element, adaptive-grid formulation, and exploits the fact that the phase and temperature fields vary significantly only near the interface. We illustrate how our method allows efficient simulation of phase-field models in very large systems, and verify the predictions of solvability theory at intermediate undercooling. We then present new results at low undercoolings that suggest that solvability theory may not give the correct tip speed in that regime. We model solidification using the phase-field model used by Karma and Rappel.

Provatas, Nikolas↗

NASA NTRS2010

Predicting Failure Progression and Failure Loads in Composite Open-Hole Tension Coupons

Failure types and failure loads in carbon-epoxy [45n/90n/-45n/0n]ms laminate coupons with central circular holes subjected to tensile load are simulated using progressive failure analysis (PFA) methodology. The progressive failure methodology is implemented using VUMAT subroutine within the ABAQUS(TradeMark)/Explicit nonlinear finite element code. The degradation model adopted in the present PFA methodology uses an instantaneous complete stress reduction (COSTR) approach to simulate damage at a material point when failure occurs. In-plane modeling parameters such as element size and shape are held constant in the finite element models, irrespective of laminate thickness and hole size, to predict failure loads and failure progression. Comparison to published test data indicates that this methodology accurately simulates brittle, pull-out and delamination failure types. The sensitivity of the failure progression and the failure load to analytical loading rates and solvers precision is demonstrated.

Arunkumar, Satyanarayana↗

NASA NTRS2023

Acoustics of a Supersonic Mach 1.4 Axisymmetric Spike Inlet

An inlet for a supersonic commercial aircraft is analyzed for noise radiation in low-speed configuration, with auxiliary inlet open. Duct mode amplitudes generated by the fan are determined by post-processing a URANS simulation and become a boundary condition to the inlet. Far-field sound is modeled using commercial acoustics finite element software. The effect of inlet geometry and acoustic liners is simulated and discussed. Absolute levels for inlet tone noise are predicted.

Supersonics↗

DOE OSTI2014

Redwing: A MOOSE application for coupling MPACT and BISON

Fuel performance and whole core neutron transport programs are often used to analyze fuel behavior as it is depleted in a reactor. For fuel performance programs, internal models provide the local intra-pin power density, fast neutron flux, burnup, and fission rate density, which are needed for a fuel performance analysis. The fuel performance internal models have a number of limitations. These include effects on the intra-pin power distribution by nearby assembly elements, such as water channels and control rods, and the further limitation of applicability to a specified fuel type such as low enriched UO2. In addition, whole core neutron transport codes need an accurate intra-pin temperature distribution in order to calculate neutron cross sections. Fuel performance simulations are able to model the intra-pin fuel displacement as the fuel expands and densifies. These displacements must be accurately modeled in order to capture the eventual mechanical contact of the fuel and the clad; the correct radial gap width is needed for an accurate calculation of the temperature distribution of the fuel rod. Redwing is a MOOSE-based application that enables coupling between MPACT and BISON for transport and fuel performance coupling. MPACT is a 3D neutron transport and reactor core simulator based on the method of characteristics (MOC). The development of MPACT began at the University of Michigan (UM) and now is under the joint development of ORNL and UM as part of the DOE CASL Simulation Hub. MPACT is able to model the effects of local assembly elements and is able calculate intra-pin quantities such as the local power density on a volumetric mesh for any fuel type. BISON is a fuel performance application of Multi-physics Object Oriented Simulation Environment (MOOSE), which is under development at Idaho National Laboratory. BISON is able to solve the nonlinearly coupled mechanical deformation and heat transfer finite element equations that model a fuel element as it is depleted in a nuclear reactor. Redwing couples BISON and MPACT in a single application. Redwing maps and transfers the individual intra-pin quantities such as fission rate density, power density, and fast neutron flux from the MPACT volumetric mesh to the individual BISON finite element meshes. For a two-way coupling Redwing maps and transfers the individual pin temperature field and axially dependent coolant densities from the BISON mesh to the MPACT volumetric mesh. Details of the mapping are given. Redwing advances the simulation with the MPACT solution for each depletion time step and then advances the multiple BISON simulations for fuel performance calculations. Sub-cycle advancement can be applied to the individual BISON simulations and allows multiple time steps to be applied to the fuel performance simulations. Currently, only loose coupling where data from a previous time step is applied to the current time step is performed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

A matrix-free hyperviscosity formulation for high-order ALE hydrodynamics

The numerical approximation of compressible hydrodynamics is at the core of high-energy density (HED) multiphysics simulations as shocks are the driving force in experiments like inertial confinement fusion (ICF). In this work, we describe our extension of the hyperviscosity technique, originally developed for shock treatment in finite difference simulations, for use in arbitrarily high-order finite element methods for Lagrangian hydrodynamics. Hyperviscosity enables shock capturing while preserving the high-order properties of the underlying discretization away from the shock region. Specifically, we compute a high-order term based on a product of the mesh length scale to a high power scaled by a hyper-Laplacian operator applied to a scalar field. We then form the total artificial viscosity by taking a non-linear blend of this term and a traditional artificial viscosity term. We also present a matrix-free formulation for computing the finite element based hyper-Laplacian operator. Such matrix-free methods have superior performance characteristics compared to traditional full matrix assembly approaches and offer advantages for GPU based HPC hardware. We demonstrate the numerical convergence of our method and its application to complex, multi-material ALE simulations on high-order (curved) meshes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2013

Numerical-Model Investigation of the Hydrothermal Regime of a Straight-Through Shallow Cooling Pond

A mathematic model based on solution of hydrodynamics and heat-transfer equations by the finite-element method is constructed to predict the hydrothermal regime of a straight-through shallow cooling pond, which provides cooling circulating water to a repository of spent nuclear fuel. Numerical experiments made it possible to evaluate the influence exerted by wind conditions and flow rate of water in the river on the temperature of the circulating water.

42 ENGINEERING↗