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 181 records · Page 10

Composite structural materials

A multifaceted program is described in which aeronautical, mechanical, and materials engineers interact to develop composite aircraft structures. Topics covered include: (1) the design of an advanced composite elevator and a proposed spar and rib assembly; (2) optimizing fiber orientation in the vicinity of heavily loaded joints; (3) failure mechanisms and delamination; (4) the construction of an ultralight sailplane; (5) computer-aided design; finite element analysis programs, preprocessor development, and array preprocessor for SPAR; (6) advanced analysis methods for composite structures; (7) ultrasonic nondestructive testing; (8) physical properties of epoxy resins and composites; (9) fatigue in composite materials, and (10) transverse thermal expansion of carbon/epoxy composites.

Ansell, G. S.↗

On the utility of finite element theory for computational fluid dynamics

An implicit finite element numerical solution algorithm is derived for the compressible Navier-Stokes equations expressed in generalized coordinates. The theoretical basis utilizes a Galerkin-Weighted Residuals formulation, and extremization of approximation error within the context of a multipole expansion. A von Neumann analysis for a simplified form indicates the algorithm fourth- to sixth-order phase accurate, with third-order dissipation for the elementary linear element construction. Performance is improved for the algorithm constructed using quadratic interpolation. Numerical experiments for shocked duct flows are employed to optimize the several algorithm parameters. Additional numerical solutions validate algorithm accuracy and utility for aerodynamics applications.

Baker, A. J.↗

Thermal-stress-free fasteners for joining orthotropic materials

Hot structures fabricated from orthotropic materials are an attractive design option for future high speed vehicles. Joining subassemblies of these materials with standard cylindrical fasteners can lead to loose joints or highly stressed joints due to thermal stress. A method has been developed to eliminate thermal stresses and maintain a tight joint by shaping the fastener and mating hole. This method allows both materials (fastener and structure), with different coefficients of thermal expansion (CTEs) in each of the three material directions, to expand freely with temperature yet remain in contact. For the assumptions made in the analysis, the joint will remain snug, yet free of thermal stress at any temperature. Finite element analysis was used to verify several thermal-stress-free fasteners and to show that conical fasteners, which are thermal-stress-free for isotropic materials, can reduce thermal stresses for transversely isotropic materials compared to a cylindrical fastener. Equations for thermal-stress-free shapes are presented and typical fastener shapes are shown.

Blosser, M. L.↗

A geometrically nonlinear analysis of interlaminar stresses in unsymmetrically laminated plates subjected to uniform thermal loading

An analytical study of interlaminar stresses in unsymmetrically laminated plates is presented. The study examines the linear elastic, large deflection response of square laminated composite plates subjected to uniform thermal loading. Both cross-ply and angle-ply, symmetric and unsymmetric, laminates are examined to evaluate the effects of mismatch between adjacent layers in elastic constants and coefficients of thermal expansion. A geometrically nonlinear kinematic description is used to predict the large out-of-plane (transverse) deflections. The nonlinear, three-dimensional boundary value problems are formulated from elasticity theory and approximate solutions are determined using the finite element method. A global/local analysis procedure is used to obtain improved free edge stress predictions. For the laminates and loading conditions considered, the results indicate that the out-of-plane deflections of the unsymmetric laminates reduce interlaminar shear stresses. These deflections also reduce interlaminar normal stresses in some laminates and increase these stresses for others. The results indicate that as the out-of-plane deflections become large, the differences in interlaminar stress predictions between linear and nonlinear theory can become quite large.

Norwood, D. S.↗

3-D Modeling of Directional Solidification of a Non-Dilute Alloy with Temperature and Concentration Fields Coupling via Materials Properties Dependence and via Double Diffusive Convection

Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.

Bune, Andris V.↗

The free and forced vibrations of structures using the finite dynamic element method

In addition to an extensive review of the literature on exact and corrective displacement based methods of vibration analysis, a few theorems are proven concerning the various structural matrices involved in such analyses. In particular, the consistent mass matrix and the quasi-static mass matrix are shown to be equivalent, in the sense that the terms in their respective Taylor expansions are proportional to one another, and that they both lead to the same dynamic stiffness matrix when used with the appropriate stiffness matrix.

Fergusson, Neil J.↗

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗

Linearized frequency domain Landau-Lifshitz-Gilbert equation formulation

We present a general finite element linearized Landau-Lifshitz-Gilbert equation (LLGE) solver for magnetic systems under weak time-harmonic excitation field. The linearized LLGE is obtained by assuming a small deviation around the equilibrium state of the magnetic system. Inserting such expansion into LLGE and keeping only first order terms gives the linearized LLGE, which gives a frequency domain solution for the complex magnetization amplitudes under an external time-harmonic applied field of a given frequency. We solve the linear system with an iterative solver using generalized minimal residual method. We construct a preconditioner matrix to effectively solve the linear system. The validity, effectiveness, speed, and scalability of the linear solver are demonstrated via numerical examples.

36 MATERIALS SCIENCE↗

Transverse cracking and stiffness reduction in composite laminates

A study of transverse cracking mechanism in composite laminates is presented using a singular hybrid finite element model. The model provides the global structural response as well as the precise local crack-tip stress fields. An elasticity basis for the problem is established by employing Lekhnitskii's complex variable potentials and method of eigenfunction expansion. Stress singularities associated with the transverse crack are obtained by decomposing the deformation into the symmetric and antisymmetric modes and proper boundary conditions. A singular hybrid element is thereby formulated based on the variational principle of a modified hybrid functional to incorporate local crack singularities. Axial stiffness reduction due to transverse cracking is studied. The results are shown to be in very good agreement with the existing experimental data. Comparison with simple shear lag analysis is also given. The effects of stress intensity factors and strain energy density on the increase of crack density are analyzed. The results reveal that the parameters approach definite limits when crack densities are saturated, an evidence of the existence of characteristic damage state.

Yuan, F. G.↗

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

Computational Modeling of Photovoltaic Mini-Modules Undergoing Accelerated Stress Testing

A finite element model of a four-cell photovoltaic mini-module was developed and compared to experimental results from an accelerated stress test protocol in order to validate that computational models can accurately represent their physical counterparts when subjected to mechanical loading and to assess mini-module representativeness against full scale photovoltaic modules. Deflected shapes across the simulated mini-modules were compared to measured mini-module shapes when subjected to various pressure loads. Displaced mini-module shape results constrained to the experimental protocols of 0.4 mm and 1.1 mm of displacement at the mini-module center were compared to experimental results of full-size modules subjected to module qualification test load levels of 1.0 kPa and 2.4 kPa, to assess if the bending of mini-modules was representative of full-sized modules under the load. Temperature cycling was incorporated into the model to simulate the impacts of stress due to thermal expansion of the backsheet and cells. A preliminary uncertainty analysis was performed to show how variations in material properties and geometric parameters change the simulation results.

computational modeling↗

Thermal-stress-free fasteners for joining orthotropic materials

Hot structures fabricated from orthotropic materials are an attractive design option for future high speed vehicles. Joining subassemblies of these materials with standard cylindrical fasteners can lead to loose joints or highly stressed joints due to thermal stress. A method has been developed to eliminate thermal stress and maintain a tight joint by shaping the fastener and mating hole. This method allows both materials (fastener and structure), with different coefficients of thermal expansion (CTE's) in each of the three principal material directions, to expand freely with temperature yet remain in contact. For the assumptions made in the analysis, the joint will remain snug, yet free of thermal stress at any temperature. Finite element analysis was used to verify several thermal-stress-free fasteners and to show that conical fasteners, which are thermal-stress-free for isotropic materials, can reduce thermal stresses for transversely isotropic materials compared to a cylindrical fastener. Equations for thermal-stress-free shapes are presented and typical fastener shapes are shown.

Blosser, Max L.↗

Thermal stresses in a spherical pressure vessel having temperature-dependent, transversely isotropic, elastic properties

Rayleigh-Ritz and modified Rayleigh-Ritz procedures are used to construct approximate solutions for the response of a thick-walled sphere to uniform pressure loads and an arbitrary radial temperature distribution. The thermoelastic properties of the sphere are assumed to be transversely isotropic and nonhomogeneous; variations in the elastic stiffness and thermal expansion coefficients are taken to be an arbitrary function of the radial coordinate and temperature. Numerical examples are presented which illustrate the effect of the temperature-dependence upon the thermal stress field. A comparison of the approximate solutions with a finite element analysis indicates that Ritz methods offer a simple, efficient, and relatively accurate approach to the problem.

Tauchert, T. R.↗

Acoustic transmission in non-uniform ducts with mean flow. I - The method of weighted residuals. II - The finite element method

The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.

Eversman, W.↗

Calculation of Thermally-Induced Displacements in Spherically Domed Ion Engine Grids

An analytical method for predicting the thermally-induced normal and tangential displacements of spherically domed ion optics grids under an axisymmetric thermal loading is presented. A fixed edge support that could be thermally expanded is used for this analysis. Equations for the displacements both normal and tangential to the surface of the spherical shell are derived. A simplified equation for the displacement at the center of the spherical dome is also derived. The effects of plate perforation on displacements and stresses are determined by modeling the perforated plate as an equivalent solid plate with modified, or effective, material properties. Analytical model results are compared to the results from a finite element model. For the solid shell, comparisons showed that the analytical model produces results that closely match the finite element model results. The simplified equation for the normal displacement of the spherical dome center is also found to accurately predict this displacement. For the perforated shells, the analytical solution and simplified equation produce accurate results for materials with low thermal expansion coefficients.

Soulas, George C.↗

Machine-Learning-aided Approach for Predicting the Thermal Expansion Behaviors in Advanced Test Reactor Capsules (NURETH-20 full paper)

Instrumented experiments at test reactors are essential to deploying new advanced reactor systems. Designing new experiments and generating data on specific conditions require both time and cost investment. A high-fidelity model of the experiment environment can be created using finite element analysis software to support the actual experiments, but computation time is still a concern in applying outcomes to real-time usage (e.g., a digital twin). This research proposes a machine-learning-aided approach to temperature and displacement predictions, based on the thickness of the outer gas gap on the experimental capsule used for the in-pile demonstration of a novel thermal conductivity probe in the Advanced Test Reactor. The capsule consisted of U10Zr fuel, a rodlet, sodium, and inner and outer capsules. There were gas gaps between the fuel and rodlet and between the inner and outer capsule. The learning data consisted of an experimental capsule’s radial distributions of temperature and displacement, as obtained from Abaqus and the physical features. For the first step, temperature was predicted using three positional parameters. Then the displacement was predicted using six different positional parameters. Each physical feature was normalized to be both nondimensional and standardized. The temperature and displacement predictions showed good agreement in all cases involving interpolation and extrapolation. Also, data similarity enhancement increased the similarity between training and target data increasing the predictive accuracy of machine-learning models. In some cases of extrapolation, the accuracy of the machine-learning model showed limited performance, but still data similarity enhancement improved the accuracy.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗