Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element expansion”

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 199 records · Page 11

Toward an In-Depth Material Model for Cermet Nuclear Thermal Rocket Fuel Elements

The development and qualification of nuclear thermal propulsion (NTP) fuel element technologies would be aided by an in-depth model of material response and failure modes at operating conditions. Integrated computational materials engineering techniques have the potential to provide such a model, as demonstrated here through three case studies focused on a tungsten-uranium mononitride cermet fuel. The first case focuses on the erosion of tungsten (W, also named wolfram), a nominal coating/cladding material, in hot hydrogen. Ab initio techniques are used to calculate erosion rates and thermal expansion at NTP operating conditions. The second focuses on the stability of uranium mononitride (UN) fuels at high temperature and in the presence of hydrogen. Phase diagram techniques reveal potential instabilities and decomposition pathways at high hydrogen concentrations. The third focuses on using microstructure information to predict high temperature mechanical response and failure of tungsten, used in refractory cermet materials. Combined finite element and discrete dislocation dynamics techniques provide mechanical properties in agreement with experimental methods. The integration of these techniques for an all-encompassing material model is discussed.

cermet↗

Thermal expansion accommodation in a jet engine frame

Design advancements to enhance stress accommodation in gas turbine engine frames are described. Consideration is given to mechanical stiffness to maintain an adequate spring rate, and to thermal stresses in both ambient and transient modes. Noting that thermal stresses occur due to differing temperatures at different parts of the frame, the matching of the thermal expansion rates of different structural materials is emphasized. Adequate stiffness is necessary to avoid dynamic reactions leading to case cracking as a result of normal imbalances of moving parts. Stress and deflection analysis of a total frame concept using a three-dimensional finite element stress analysis computer program, with modelling of all frame components, is presented, and illustrations are provided.

Schneider, M. H.↗

Single Axis Piezoceramic Gimbal

This paper describes the fabrication, testing, and analysis of a single axis piezoceramic gimbal. The fabrication process consist of pre-stressing a piezoceramic wafer using a high-temperature thermoplastic polyimide and a metal foil. The differential thermal expansion between the ceramic and metal induces a curvature. The pre-stressed, curved piezoceramic is mounted on a support mechanism and a mirror is attached to the piezoceramic. A plot of gimbal angle versus applied voltage to the piezoceramic is presented. A finite element analysis of the piezoceramic gimbal is described. The predicted gimbal angle versus applied voltage is compared to experimental results.

Horner, Garnett C.↗

Single Axis Piezoceramic Gimbal

This paper describes the fabrication, testing, and analysis of a single axis piezoceramic gimbal. The fabrication process consists of pre-stressing a piezoceramic wafer using a high-temperature thermoplastic polyimide and a metal foil. The differential thermal expansion between the ceramic and metal induces a curvature. The pre-stressed, curved piezoceramic is mounted on a support mechanism and a mirror is attached to the piezoceramic. A plot of gimbal angle versus applied voltage to the piezoceramic is presented. A finite element analysis of the piezoceramic gimbal is described. The predicted gimbal angle versus applied voltage is compared to experimental results.

Horner, Garnett↗

Dimensional stability of curved panels with cocured stiffeners and cobonded frames

Closed form and finite element analyses are presented for axial direction and transverse direction dimensional stability of skin/stringer panels. Several sensitivity studies are presented to illustrate the influence of various design parameters on the dimensional stability of these panels. Panel geometry, material properties (stiffness and coefficient of thermal expansion), restraint conditions and local details, such as resin fillets, all combine to influence dimensional stability, residual and assembly forces.

Mabson, G. E.↗

The Effect of Temperature Dependent Material Nonlinearities on the Response of Piezoelectric Composite Plates

Previously developed analytical formulations for piezoelectric composite plates are extended to account for the nonlinear effects of temperature on material properties. The temperature dependence of the composite and piezoelectric properties are represented at the material level through the thermopiezoelectric constitutive equations. In addition to capturing thermal effects from temperature dependent material properties, this formulation also accounts for thermal effects arising from: (1) coefficient of thermal expansion mismatch between the various composite and piezoelectric plies and (2) pyroelectric effects on the piezoelectric material. The constitutive equations are incorporated into a layerwise laminate theory to provide a unified representation of the coupled mechanical, electrical, and thermal behavior of smart structures. Corresponding finite element equations are derived and implemented for a bilinear plate element with the inherent capability to model both the active and sensory response of piezoelectric composite laminates. Numerical studies are conducted on a simply supported composite plate with attached piezoceramic patches under thermal gradients to investigate the nonlinear effects of material property temperature dependence on the displacements, sensory voltages, active voltages required to minimize thermal deflections, and the resultant stress states.

Lee, Ho-Jun↗

Direct Calculation of the Scattering Amplitude Without Partial Wave Analysis

Two new developments in scattering theory are reported. We show, in a practical way, how one can calculate the full scattering amplitude without invoking a partial wave expansion. First, the integral expression for the scattering amplitude f(theta) is simplified by an analytic integration over the azimuthal angle. Second, the full scattering wavefunction which appears in the integral expression for f(theta) is obtained by solving the Schrodinger equation with the finite element method (FEM). As an example, we calculate electron scattering from the Hartree potential. With minimal computational effort, we obtain accurate and stable results for the scattering amplitude.

Shertzer, J.↗

Comprehensive Material Characterization and Simultaneous Model Calibration for Improved Computational Simulation Credibility

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10 ⁻5 and 10 ⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

36 MATERIALS SCIENCE↗

Vibrations of incompressible liquid-plastic tank combinations using a series representation of the liquid

A series representation is presented of the dynamic behavior of incompressible non-viscous liquids contained in partially filled elastic shells of revolution. Using a series in which each term satisfies exactly Laplace's equation, the liquid portion of the problem becomes an integral equation on the wetted wall of the shell. Physical insight required in selecting the terms of the series is such that selection may be readily programed on a computer. Liquid-shell vibration interaction is investigated by modifying an existing shell of revolution vibration program with the series expansion of the liquid pressure loading on the wetted shell wall. Application of the modified program is demonstrated on a hemispherical tank and a realistic scaled model of a launch vehicle tank having a more complex geometry. Comparison of the present analysis with a NASTRAN finite element analysis for a hemispherical tank, and comparison with test results for a launch vehicle tank shows excellent agreement and rapid convergence of the series representation of the liquid to the true solution.

Housner, J. M.↗

Thermal stress in high temperature cylindrical fasteners

Uninsulated structures fabricated from carbon or silicon-based materials, which are allowed to become hot during flight, are attractive for the design of some components of hypersonic vehicles. They have the potential to reduce weight and increase vehicle efficiency. Because of manufacturing contraints, these structures will consist of parts which must be fastened together. The thermal expansion mismatch between conventional metal fasteners and carbon or silicon-based structural materials may make it difficult to design a structural joint which is tight over the operational temperature range without exceeding allowable stress limits. In this study, algebraic, closed-form solutions for calculating the thermal stresses resulting from radial thermal expansion mismatch around a cylindrical fastener are developed. These solutions permit a designer to quickly evaluate many combinations of materials for the fastener and the structure. Using the algebraic equations developed, material properties and joint geometry were varied to determine their effect on thermal stresses. Finite element analyses were used to verify that the closed-form solutions derived give the correct thermal stress distribution around a cylindrical fastener and to investigate the effect of some of the simplifying assumptions made in developing the closed-form solutions for thermal stresses.

Blosser, Max L.↗

Thermoelastic Formulation of Stiffened, Unsymmetric Composite Panels for Finite Element Analysis of High Speed Aircraft

An emerging technology need for capturing 3-D panel thermoelastic response with 2-D planar finite element models (FEMs) is aided with an equivalent plate stiffness and thermal coefficient formulation. The formulation is general and applies to all panel concepts. Included with the formulation is the ability to provide membrane-bending coupling of unsymmetric sections and calculation of all thermal expansion and bending responses from in-plane and through-the-thickness temperature gradients. Thermal residual strains for both the laminates and plies are included. The general formulation is defined and then applied to a hat-shaped, corrugated stiffened panel. Additional formulations are presented where required to include all of the hat's unique characteristics. Each formulation is validated independently with 3-D FEA.

Collier, Craig S.↗

A finite element technique for non-deterministic thermal deformation analyses including temperature dependent material properties

A technique utilizing the finite element displacement method is developed for the static analysis of structures subjected to non-deterministic thermal loading in which the material properties, assumed isotropic, are temperature dependent. Matrix equations are developed for the first two statistical moments of the displacements using a third order series expansion for the displacements in terms of the random temperatures. Sample problems are included to demonstrate the range of applicability of the third order series solutions. These solutions are compared with results from Monte Carlo analyses and also, for some problems, with solutions obtained by numerically integrating equations for the statistical properties of the displacements. In general, it is shown that the effect of temperature dependent material properties can have a significant effect on the covariances of the displacements.

Case, W. R., Jr.↗

Applications of perturbation techniques

Two perturbation techniques were applied to two singular perturbation problems in heat transfer to obtain uniformly valid solutions which can serve as benchmarks for finite difference and finite element techniques. In the first problem, the method of strained parameters coupled with the application of a solvability condition is used to obtain a uniform solution for the problem of unsteady heat conduction in a long nearly circular cylinder. In the second problem, the method of matched asymptotic expansion coupled with Van Dyke's matching principle is used to obtain a uniform solution for the problem of one dimensional conduction-convection heat transfer of a uniform fluid flow.

Kandil, O. A.↗

Microstructural Simulation of Ion-Irradiated Natural Rocks and Minerals

This report details the contribution of Oak Ridge National Laboratory to the Nuclear Energy University Program (NEUP) project, "Rapid Characterization of Concrete Mineralogy using Multi-Scale Tools" led by the University of Illinois at Urbana-Champaign. To support the subsequent license renewal (SLR) of the US nuclear power plants (NPPs) fleet, the performance of concrete-forming aggregates against neutron irradiation needs to be assessed. Ion irradiation is proposed as a rapid, cost-effective surrogate method to full neutron irradiation testing. To complement the characterization of ion-irradiated rock specimens, ORNL has run finite-elements and fast-Fourier transform (FFT-based) simulations using the codes MARS and Microstructure-Oriented Scientific Analysis of Irradiated Concrete (MOSAIC). The main conclusion of this analysis is that the apparent post–ion-irradiation step-height underestimates the accumulated free radiation-induced volumetric expansion (RIVE) in the ion-implanted depth by about 15% at full amorphization and about 25% toward the beginning of the ion irradiation experiment. The discrepancy is explained by the fact that the step height is proportional to the sum of the RIVE (i.e., one third of the RIVE) and the irradiation-assisted plastic strains in the vertical direction (lower than two thirds of the RIVE). Because of the large lateral compressive stresses caused by the restraining effect of the unirradiated substrate, the stress field in the mineral grains and at the grain boundary (GB) in the ion-implanted layer is much different than the comparable stress field occurring during neutron irradiation. Hence, the mismatch RIVE causing cracks in the rock-forming minerals leads to different cracking patterns. Ion irradiation continues to be used as a rapid technique to assess the RIVE potential of rocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rapid iterative reanalysis for automated design

A method for iterative reanalysis in automated structural design is presented for a finite-element analysis using the direct stiffness approach. A basic feature of the method is that the generalized stiffness and inertia matrices are expressed as functions of structural design parameters, and these generalized matrices are expanded in Taylor series about the initial design. Only the linear terms are retained in the expansions. The method is approximate because it uses static condensation, modal reduction, and the linear Taylor series expansions. The exact linear representation of the expansions of the generalized matrices is also described and a basis for the present method is established. Results of applications of the present method to the recalculation of the natural frequencies of two simple platelike structural models are presented and compared with results obtained by using a commonly applied analysis procedure used as a reference. In general, the results are in good agreement. A comparison of the computer times required for the use of the present method and the reference method indicated that the present method required substantially less time for reanalysis. Although the results presented are for relatively small-order problems, the present method will become more efficient relative to the reference method as the problem size increases. An extension of the present method to static reanalysis is described, ana a basis for unifying the static and dynamic reanalysis procedures is presented.

Bhatia, K. G.↗

Recent advances in error estimation and adaptive improvement of finite element calculations

Recent advances in adaptive finite elements are summarized. General concepts behind a posteriori error estimation, h-method data management, and algorithms for fluid-mechanics applications are then examined, and some results of numerical experiments with new adaptive codes are given. Numerical examples include supersonic flow over a 20-deg ramp, supersonic flow in expansion corners, the reflecting-shock problem, and the rotating-cone problem. Finally, future directions of research in the field are outlined.

Oden, J. T.↗

STOP Analysis and Optimization of a Very-Low-Distortion Space Instrument: HST WFC3 Case Study

New generation optical instruments with very demanding stability requirements are being proposed and developed for space applications. STOP (Structural-Thermal-Optical Performance) analysis and optimization is crucial in meeting the very tight distortion budgets of these instruments. This presentation outlines STOP analysis and optimization approach in the context of WFC3 (Wide-Field Camera 3), which is a radial instrument designed to replace the Wide-Field Planetary Camera 2 (WFPC2) of the Hubble Space Telescope (HST). WFC3 houses two separate channels, UVIS and IR, and will have greater throughput and sensitivity than WFPC2. WFC3 line-of-sight alignment budget for the UVIS and IR channels are as small as 10 and 20 milli-arcsec, respectively. Its optical bench is the most critical subsystem effecting the optical stability of WFC3 hence our effort concentrates on the design and analysis of the bench and its interfaces. Structural analysis has accompanied the mechanical design of the bench since the initial concept study. A high fidelity structural Finite Element Model (FEM) of the bench has been developed and used for minimizing its thermally induced distortions as well as sizing it to meet the stiffness and strength requirements of a Shuttle launch. The bench is a composite honeycomb panel box structure with a very low planar Coefficient of Thermal Expansion (CTE) of approximately 0.1 ppm/C. Optic components are mounted to super-INVAR inserts bonded into the panels. The bench is kinematically supported on three HST latches via interface struts, which are tailored to exhibit negative CTE to cancel out the thermal motions of the latches. The interface struts also incorporate flexure elements to minimize the mechanical distortions coming into the bench from its enclosure. Bench FEM is coupled with the enclosure FEM to quantify these effects. Short term or on-orbit STOP analysis includes distortion due to the temperature variations of the bench, the struts, and the enclosure. Long term or ground-to-orbit STOP analysis includes distortional effects of gravity release, desorption, and assembly in addition to the ground-to-orbit temperature variations. A rigorous testing program has been implemented for verifying the material properties and the analysis predictions. STOP analysis results demonstrate that both the short-term and the long-term alignment budgets will be met. Presentation will cover design and analysis details that are critical to a successful implementation of the STOP analysis and optimization process.

Kunt, Cengiz↗