SEARCH · Engineering Papers
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.
An approach is proposed for calculating the transient torque imbalance of rotating shafts in the diesel generators of autonomous power units. The maximum shear strain of the rotor shaft in the generator is assessed by the finite-element method.
This report details the development and demonstration of an entirely new capability in Pronghorn, namely the ability to model liquid-metal fast reactor (LMFR) flow conditions on the engineering scale. We developed two modeling approaches for LMFR that can be used separately or be combined into a hybrid simulation: (1) a modern subchannel capability called Pronghorn-Subchannel for square and hexagonal lattices, and (2) a porous flow capability for LMFR geometries. The report emphasizes the novel aspects of the developed subchannel capability and the interoperability of the subchannel capability, porous flow capability, and multiphysics tools within the multiphysics object oriented simulation environment (MOOSE). Here we demonstrate the ability to: (1) Accurately model subchannel flow in hexagonal lattices; (2) Couple the subchannel flow model to multidimensional finite-element method (FEM) or finite-volume method (FVM) heat conduction models; (3) Model LMFRs using Pronghorn’s porous media FVM approach; (4) Couple porous flow FVM and subchannel models in a single simulation; (5) Explicitly model inter-wrapper flows along with conjugate heat transfer from the intra-element flow; and (6) Demonstrate the numerical robustness of the subchannel algorithm by simulating intra-element flow recirculation in a high-buoyancy, low-flow fuel element.
The implementation of the finite element method on unstructured triangular grids is described and the development of centered finite element schemes for the solution of the compressible Euler equation on general triangular and tetrahedral grids is discussed. Explicit and implicit Lax-Wendroff type methods and a method based upon the use of explicit multistep timestepping are considered. In the latter case, the convergence behavior of the method is accelerated by the incorporation of a fully unstructured multigrid procedure. The advancing front method for generating unstructured grids of triangles and tetrahedra is described and the application of adaptive mesh techniques to both steady and transient flow analysis is illustrated.
We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.
Low-enriched uranium alloyed with 10 wt% molybdenum (U-10Mo) has been identified as a promising alternative to highly enriched uranium as fuel for the nation’s high-performance research reactors. Manufacturing the U-10Mo alloy consists of multiple complex thermomechanical processes, which is highly challenging for computational modeling. The integrated computational materials engineering (ICME) concept supports building a microstructure-based framework to investigate the effect of manufacturing processes on the material microstructure. In this report, an ICME model that combines homogenization, hot rolling, annealing, and cold rolling is presented as an example of implementing the ICME concept for modeling the manufacturing processes of the U-10Mo alloy. The integrated model enables information to be passed between the different model components, leading to an improved prediction and a better understanding of microstructure across multiple processes. The integrated model was also implemented to investigate the variation of chemical composition, thickness of the zirconium interlayer, grain growth, and carbide redistribution and fracture of U-10Mo during the thermomechanical processes. With the homogenization model, a molybdenum concentration map can be reconstructed from any given microstructure, and the homogenization time needed to reach a desired homogenization level can also be approximated. The microstructure-based finite-element rolling model, shows that the homogenization process reduces the nonuniformity in the thickness of the zirconium interlayers on the U-10Mo. With the output from the rolling model, a statistical analysis of second-phase hard particles was conducted to quantify the particle redistribution and interparticle spacing between carbides during the deformation. Also, new stringer identification criteria were proposed to predict the correlation between stringer volume fractions and hot-rolling reduction. A macro-finite-element method hot-rolling model showed that higher reduction passes can be achieved within the mill separation force range, and a new “aggressive” rolling schedule was proposed. The model also showed the importance of can material on the quality of rolled foils, with stronger can materials resulting in hot-rolled U-10Mo free of defects (waviness, dog-boning). The deformation-induced recrystallization model demonstrated capability to use inputs from deformation simulations of strain accumulation in U-10Mo samples to predict grain growth behavior of the fuel during thermal treatment stages. Finally, a plane-strain compression mode simulated particle fracture and calculates the void volume fraction generated during cold rolling. The macro-finite-element method cold-rolling model showed that smaller diameter rolls are necessary to achieve the desired foil thickness.
Continuum resonance damping is an important factor in determining the stability of certain global modes in fusion plasmas. A number of analytic and numerical approaches have been developed to compute this damping, particularly, in the case of the toroidicity-induced shear Alfvén eigenmode. This paper compares results obtained using an analytical perturbative approach with those found using resistive and complex contour numerical approaches. It is found that the perturbative method does not provide accurate agreement with reliable numerical methods for the range of parameters examined. This discrepancy exists even in the limit where damping approaches zero. When the perturbative technique is implemented using a standard finite element method, the damping estimate fails to converge with radial grid resolution. The finite elements used cannot accurately represent the eigenmode in the region of the continuum resonance, regardless of the number of radial grid points used.
It is well-known that stresses (and strains) calculated by a displacement-based finite element analysis are generally not as accurate as the displacements. In addition, the calculated stress field is typically discontinuous at element interfaces. Because the stresses are typically of more interest than the displacements, several procedures have been proposed to obtain a smooth stress field, given the finite element stresses, and to improve the accuracy. Hinton and Irons introduced global least squares smoothing of discrete data defined on a plane using a finite element formulation. Tessler and co-workers recently developed a conceptually similar formulation for smoothing of two-dimensional data based on a discrete least square approximation with a penalty constraint. The penalty constraint results in a stress field which is C(exp 1)-continuous, a result not previously obtained. The approach requires additional, 'smoothing' finite element analysis and for their two-dimensional application, they used a conforming C(exp 0)-continuous triangular finite element based on a conforming plate element. This paper presents the results of a detailed investigation into the application of Tessler's smoothing procedure to the smoothing of finite element stresses from one-dimensional problems. Although the one-dimensional formulation has some practical applicability, such as in truss, beam, axisymmetric mechanics, and one-dimensional heat conduction, the primary motivation for developing the one-dimensional smoothing case is to explore the characteristics of the general smoothing strategy. In particular, it is used to describe the behavior of the method and to explore the suitability of criteria proposed for the smoothing analysis. Prior to presenting numerical results, the variational formulation of the smoothing strategy is presented and a criterion for the smoothing analysis is described.
Flaws resulting from improper welding and forging are usually modeled as cracks in flat plates, hollow cylinders or spheres. The stress intensity factor solutions for these crack cases are of great practical interest. This report describes some recent efforts at improving the stress intensity factor solutions for cracks in such geometries with emphasis on hollow cylinders. Specifically, two crack configurations for cylinders are documented. One is that of a surface crack in an axial plane and the other is a part-through thumb-nail crack in a circumferential plane. The case of a part-through surface crack in flat plates is used as a limiting case for very thin cylinders. A combination of the two cases for cylinders is used to derive a relation for the case of a surface crack in a sphere. Solutions were sought which cover the entire range of the geometrical parameters such as cylinder thickness, crack aspect ratio and crack depth. Both the internal and external position of the cracks are considered for cylinders and spheres. The finite element method was employed to obtain the basic solutions. Power-law form of loading was applied in the case of flat plates and axial cracks in cylinders and uniform tension and bending loads were applied in the case of circumferential (thumb-nail) cracks in cylinders. In the case of axial cracks, the results for tensile and bending loads were used as reference solutions in a weight function scheme so that the stress intensity factors could be computed for arbitrary stress gradients in the thickness direction. For circumferential cracks, since the crack front is not straight, the above technique could not be used. Hence for this case, only the tension and bending solutions are available at this time. The stress intensity factors from the finite element method were tabulated so that results for various geometric parameters such as crack depth-to-thickness ratio (a/t), crack aspect ratio (a/c) and internal radius-to-thickness ratio (R/t) or the crack length-to-width ratio (2c/W) could be obtained by interpolation and extrapolation. Such complete tables were then incorporated into the NASA/FLAGRO computer program which is widely used by the aerospace community for fracture mechanics analysis.
Skin and proximity effects can cause a significant drop in the effective leakage inductance of a transformer when the operating frequency is increased. Although the magnetic image method-based double-2-D model can calculate the lowfrequency leakage inductance with sufficient accuracy, it is inherently a frequency-independent model. While Dowell’s 1-D model uses frequency-dependent relations to account for both skin and proximity effects, its accuracy is severely affected by the assumed winding geometry. In this paper, a hybrid model is proposed that uses superposition to combine a modified Dowell’s model with the double-2-D model. The proposed model is investigated on a variable inductance transformer (VIT)—a partially-filled transformer whose leakage inductance can be varied by moving one of the windings mechanically. The frequency-dependent leakage inductances of the VIT evaluated using the hybrid model are in excellent agreement with the corresponding finite element method (FEM) simulated and experimentally measured values, thereby validating the proposed hybrid model.
A finite element method to model dynamic structural systems undergoing large rotations is presented. The dynamic systems are composed of rigid joint bodies and flexible beam elements. The configurations of these systems are subject to change due to the relative motion in the joints among interconnected elastic beams. A body fixed reference is defined for each joint body to describe the joint body's displacements. Using the finite element method and the kinematic relations between each flexible element and its corotational reference, the total displacement field of an element, which contains gross rigid as well as elastic effects, can be derived in terms of the translational and rotational displacements of the two end nodes. If one end of an element is hinged to a joint body, the joint body's displacements and the hinge degree of freedom at the end are used to represent the nodal displacements. This results in a highly coupled system of differential equations written in terms of hinge degrees of freedom as well as the rotational and translational displacements of joint bodies and element nodes.
The depth-dependent strain partitioning across the interfaces in the growth direction of the NiAl/Cr(Mo) nanocomposite between the Cr and NiAl lamellae was directly measured experimentally and simulated using a finite element method (FEM). Depth-resolved X-ray microdiffraction demonstrated that in the as-grown state both Cr and NiAl lamellae grow along the 111 direction with the formation of as-grown distinct residual ~0.16% compressive strains for Cr lamellae and ~0.05% tensile strains for NiAl lamellae. Three-dimensional simulations were carried out using an implicit FEM. First simulation was designed to study residual strains in the composite due to cooling resulting in formation of crystals. Strains in the growth direction were computed and compared to those obtained from the microdiffraction experiments. Second simulation was conducted to understand the combined strains resulting from cooling and mechanical indentation of the composite. Numerical results in the growth direction of crystal were compared to experimental results confirming the experimentally observed trends.
This paper introduces the concept of a novel variable inductance transformer (VIT) whose leakage and magnetizing inductances can be varied dynamically and independently to aid the advancement of research in isolated power electronic converters. In this transformer, the leakage inductance can be varied by varying the length of overlap between the primary and secondary windings, while the magnetizing inductance can be varied by varying the air gap between the core legs. Individual controls for the two inductances ensure that one can be varied independent of the other, as needed to satisfy different power conversion objectives. In this paper, the analytical models of the variable leakage and magnetizing inductances are presented. Analytical results are compared with those obtained from Finite Element Method (FEM) and an experimental prototype of a VIT. Analytical solutions of the variable leakage inductance present an error of less than 2.35 % when compared to FEM solutions, while the analytical solutions of the variable magnetizing inductance present a maximum error of 1.82 % when compared to experimental measurements.
Finite element models were used to study the effects of notch angle variations on the stress state within an Iosipescu shear test speciment. These analytical results were also studied to determine the feasibility of using strain gage rosettes and a modified extensometer to measure shear strains in this test specimen. Analytical results indicate that notch angle variations produced only small differences in simulated shear properties. Both strain gage rosettes and the modified extensometer were shown to be feasible shear strain transducers for the test method. The Iosipoescu shear test fixture was redesigned to incorporate several improvements. These improvements include accommodation of a 50 percent larger specimen for easier measurement of shear train, a clamping mechanism to relax strict tolerances on specimen width, and a self contained alignment tool for use during specimen installation. A set of in-plane and interlaminar shear properties were measured for three graphite fabric/epoxy composites of T300/934 composite material. The three weave patterns were Oxford, 5-harness satin, and 8-harness satin.
Finite Element Methods (FEM) and Spectral Element Methods (SEM) are crucial for solving partial differential equations (PDEs) on complex geometries. SEM offers superior accuracy due to potential superconvergence for simple domains. Challenges persist for domains with curved boundaries, restricting SEM’s advantages in real-world applications. A proposed solution is the introduction of a novel strategy to enhance accuracy and maintain superconvergence of SEM in curved domains. The strategy includes a mesh-generation procedure with geometrically refined elements near curved boundaries and a post-processing phase using the Adaptive Extended Stencil Finite Element Method (AES-FEM). The method, named AES-FEM post-processed Spectral Element Method (ApSEM), aligns the accuracy of non-tensor-product elements with superconvergent spectral elements.
Recently, the piezotronic effect has been observed in two-dimensional single-layer MoS{sub 2} materials, which have potential applications in force and pressure triggered or controlled electronic devices, sensors, and human-machine interfaces. However, classical theory faces the difficulty in explaining the mechanism of the piezotronic effect for the top- and enclosed-contacted MoS{sub 2} transistors, since the piezoelectric charges are assumed to exist only at the edge of the MoS{sub 2} flake that is far from the electronic transport pathway. In the present study, we identify the piezoelectric charges at the MoS{sub 2}/metal-MoS{sub 2} interface by employing both the density functional theory and finite element method simulations. This interface is on the transport pathway of both top- and enclosed-contacted MoS{sub 2} transistors, thus it is capable of controlling their transport properties. This study deepens the understanding of piezotronic effect and provides guidance for the design of two-dimensional piezotronic devices.
Gas turbine engines experience “rub” when the rotating blades come in contact with a static abradable coating. This results in extreme strain rates and dynamics inside a high-temperature/high-pressure environment. Current rub models are phenomenological and do not reflect the underlying microstructures, thus limiting their prediction accuracy. In this work, a microstructure-informed, reduced order modeling framework is introduced for simulating abradable coating “rub" behavior. This framework comprises a microscale model constructed based on digitized abradable microstructure and explicitly simulates the mechanical behavior of each constituent phases and their interactions. After calibration and validation with experiment data, the calibrated microscale model is used to generate data across a vast range of applied strain rates and temperature with various load paths. Then, the virtually generated data are used to fit the macroscopic-reduced order model, which enables fast component scale rub simulation without compromising the integrity of the complex material behavior. In conclusion, the proposed effort will address the technical challenge of predicting abradable material behavior during rub through the application of multiscale modeling from microstructure to engines behavior, effectively reducing the development costs and time of new abradable material for better “rub” properties.
The ability to accurately predict the thermomechanical deformation response of advanced composite materials continues to play an important role in the development of these strategic materials. Analytical models that predict the effective behavior of composites are used not only by engineers in performing structural analysis of large-scale composite components but also by material scientists in developing new material systems. For an analytical model to fulfill these two distinct functions, it must be based on a micromechanics approach that uses physically based deformation and life constitutive models, and it must allow one to generate the average (macro) response of a composite material given the properties of the individual constituents and their geometric arrangement. Only then can such a model be used by a material scientist to investigate the effect of different deformation mechanisms on the overall response of the composite and, thereby, identify the appropriate constituents for a given application. However, if a micromechanical model is to be used in a large-scale structural analysis it must be (1) computationally efficient, (2) able to generate accurate displacement and stress fields at both the macro and micro level, and (3) compatible with the finite element method. In addition, new advancements in processing and fabrication techniques now make it possible to engineer the architectures of these advanced composite systems. Full utilization of these emerging manufacturing capabilities require the development of a computationally efficient micromechanics analysis tool that can accurately predict the effect of microstructural details on the internal and macroscopic behavior of composites. Computational efficiency is required because (1) a large number of parameters must be varied in the course of engineering (or designing) composite materials and (2) the optimization of a material's microstructure requires that the micromechanics model be integrated with optimization algorithms. From this perspective, analytical approaches that produce closed-form expressions which describe the effect of a material's internal architecture on the overall material behavior are preferable to numerical methods such as the finite element or finite difference schemes.