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.
A geometrical interpretation of the grouping and variation rate terms of Reynolds stresses in compressible turbulent flow is presented. Examples of applications of the method are also given.
The history-dependent behaviors of classical plasticity models are often driven by internal variables evolved according to phenomenological laws. The difficulty to interpret how these internal variables represent a history of deformation, the lack of direct measurement of these internal variables for calibration and validation, and the weak physical underpinning of those phenomenological laws have long been criticized as barriers to creating realistic models. In this work, geometric machine learning on graph data (e.g. finite element solutions) is used as a means to establish a connection between nonlinear dimensional reduction techniques and plasticity models. Geometric learning-based encoding on graphs allows the embedding of rich time-history data onto a low-dimensional Euclidean space such that the evolution of plastic deformation can be predicted in the embedded feature space. Finally, a corresponding decoder can then convert these low-dimensional internal variables back into a weighted graph such that the dominating topological features of plastic deformation can be observed and analyzed.
The thrust of advanced nuclear reactor demonstrations demands the accelerated qualification of in-core materials to enable licensing processes and developing the performance data. Due to the high-operating temperatures of such reactors, long-term mechanical behavior under constant load is essential to determine the geometrical and mechanical integrity of in-core components during operation and off-normal conditions. Therefore, the thermal creep behavior of neutron-irradiated advanced reactor materials must be determined. The feasibility of using subsize specimens for the irradiation campaigns and the limited available infrastructure challenge the assessment of thermal creep behavior of advanced reactor materials. Therefore, the U.S. Department of Energy Office of Nuclear Energy (DOE-NE) National Reactor Innovation Center (NRIC) prioritizes the development of a thermal creep testing infrastructure for multiple subsize specimens to accelerate the demonstration and deployment of advanced reactor concepts. This report describes the activities for the construction of a thermal creep testing capability at Idaho National Laboratory (INL). The overall project consists of conceptual design, out-of-cell demonstration, and in-cell demonstration phases. During fiscal year (FY)-2021, the team finished the conceptual design of a thermal creep test facility that can test multiple subsize specimens. This conceptual design consisted of the determination of technical and functional requirements, the determination of the design space, and the preparation of the technical drawings. Technical and functional requirements were categorized as required and desired capabilities and the conceptual design was performed to meet all the required capabilities with the flexibility to achieve the desired capabilities. The design space identified the operational capacity of the creep frame for different advanced reactor relevant materials with the consideration of the feasible operation in the hotcell at the Fuels and Applied Science Building (FASB) at INL. Based on the requirements and design space, the multiple specimen creep frame was designed. The official INL engineering drawing process was started and the procurement of materials for construction was initiated. For FY-2022, the out-of-cell demonstration and final installation of the multiple creep frame is planned.
Laser Powder Bed Fusion (LPBF) enables the additive manufacturing of complex lattice structures. However, the fabrication of lattice structures via LPBF poses challenges in achieving the intended geometrical accuracy due to their inherent complexity. This study investigates the effects of LPBF processing parameters, specifically laser power and scanning speed, on the geometrical characteristics, surface quality, and mechanical behavior of Inconel 718 lattices structures. The results reveal that processing parameters required for the fabrication of near-full dense structures do not translate effectively to lattice configurations, as variations in energy input influence lattice geometry and surface quality. In this work, strut thickness, open-pore size, open-cell porosity, and surface roughness were measured, and the mechanical properties of the lattices were evaluated under shear loading. The findings indicate that lower energy inputs, achieved by reducing laser power and increasing scanning speed, yield porous structures but lead to mechanical degradation. In contrast, high energy inputs lead to lattices with enhanced strength but result in undesirable open-pore blockage and dimensional inaccuracies. These findings provide insights into tailoring LPBF parameters for dimensional accuracy in lattices and correlating the processing parameters to mechanical performance and surface roughness.
We report the negative Poisson’s (NPR) ratio in a two-dimensional (2D) material is a counterintuitive mechanical property that facilitates the development of nanoscale devices with sophisticated functionality. Inspired by the peculiar buckled lower-symmetric, trilayered geometry of pentagonal monolayers, we theoretically predict penta-SiCN, a ternary auxetic metallic monolayer with highly tunable NPR. The penta-SiCN is structurally, thermally, dynamically, and mechanically stable, and sustainable at and beyond room temperature with experimental feasibility. It possesses nontrivial geometrical and mechanical isotropy and relatively moderate thickness. Remarkably, the shorter and quasi sp3-hybridized C–N bond and the rigidity against the strain allow the monolayer to possess a high value of NPR (-0.136), even higher than that of black phosphorene, extendable up to -0.639 by 4% of biaxial stretching. On the other hand, the 2D Young’s modulus of 129.88 N/m decreases to 41.34 N/m at equivalent stretching, indicating relative softening and flexibility. Interestingly, a buckled-to-planar phase transition is identified at 10% biaxial strain before it suffers the fracture at 16%. Additionally, the strong optical anisotropy, absorbance (up to 6.51 x 105 cm -1 ), and presence of plasmon frequency demonstrate its potential application in optomechanical and plasmonics.
Accurately predicting redeposition is vital for high-Z plasma-facing component (PFC) survivability in magnetic confinement fusion. In this study, we categorize high-Z redeposition into three mechanisms: geometric-driven (prompt), sheath-driven (local), and scrape-off-layer-driven (far) redeposition. To investigate these mechanisms, we employ Monte Carlo transport codes to simulate azimuthally symmetric tungsten source erosion and redeposition in a tokamak. By iteratively analyzing critical parameters, we evaluate redeposition scaling for each mechanism. Specifically, we investigate the impact of magnetic-field-to-PFC pitch angle assumptions on PFC losses into the scrape-off layer. Our findings reveal significant pitch angle sensitivity due to an asymmetric prompt vs local redeposition trade-off. These results enhance our understanding of redeposition phenomena in fusion plasma environments.
Planning for manipulation in complex environments often requires reasoning about the geometric and mechanical constraints which are posed by the task. In planning assembly operations, the automatic generation of operations sequences depends on the geometric feasibility of paths which permit parts to be joined into subassemblies. Feasible locations and collision-free paths must be present for part motions, robot and grasping motions, and fixtures. This paper describes an approach to reasoning about the feasibility of straight-line paths among three-dimensional polyhedral parts using an algebra of polyhedral cones. A second method recasts the feasibility conditions as constraints in a nonlinear optimization framework. Both algorithms have been implemented and results are presented.
In this work, we propose a special relativistic framework for quantum mechanics. It is based on introducing a Hilbert space for events. Events are taken as primitive notions (as customary in relativity), whereas quantum systems (e.g. fields and particles) are emergent in the form of joint probability amplitudes for position and time of events. Textbook relativistic quantum mechanics and quantum field theory can be recovered by dividing the event Hilbert spaces into space and time (a foliation) and then conditioning the event states onto the time part. Our theory satisfies the full Lorentz symmetry as a ‘geometric’ unitary transformation, and possesses relativistic observables for space (location of an event) and time (position in time of an event).
We present a family of preconditioning strategies for the contact problem in fractured and faulted porous media. We combine low-order continuous finite elements to simulate the bulk deformation with piecewise constant Lagrange multipliers to impose the frictional contact constraints. This formulation is not uniformly inf-sup stable and requires stabilization. We improve previous work by Franceschini et al. (2020) by introducing a novel jump stabilization technique that requires only local geometrical and mechanical properties. We then design scalable preconditioning strategies that take advantage of the block structure of the Jacobian matrix using a physics-based partitioning of the unknowns by field type, namely displacement and Lagrange multipliers. The key to the success of the proposed preconditioners is a pseudo-Schur complement obtained by eliminating the Lagrange multiplier degrees of freedom, which can then be efficiently solved using an optimal multigrid method. Numerical results, including complex real-world problems, are presented to illustrate theoretical properties, scalability and robustness of the preconditioner. A comparison with other approaches available in the literature is also provided.
A two-phase approach and a computational procedure are used for predicting the variability of the response of stiffened composite panels associated with variations in the geometric and material parameters of the structures. In the first phase, hierarchical sensitivity analysis is used to identify the major parameters that have the most effect on the response quantities of interest. In the second phase, the major parameters are taken to be fuzzy parameters, and a fuzzy set analysis is used to determine the range of variation of the response, associated with preselected variations in the major parameters. Numerical results are presented showing the variability of the response of panels with both continuous and terminated stiffeners associated with variations in the micro mechanical and geometric parameters. Both flat and curved panels are considered.
This study focused on the analysis of geometrically necessary dislocation (GND) densities for five selected fine-grained magnesium samples. Among the samples, three were tested under different fatigue-loading conditions at 0 °C, one experienced quasi-static tensile loading at 0 °C, and one represented the as-rolled state. The fatigue-tested samples were chosen according to the relationship between the maximum loading stress of a test and the material’s yield strength. This study provides new insights on the deformation mechanism of fine-grained magnesium at 0 °C. It is observed that the average GND densities were increased by 95~111% for the tested samples when compared with the as-rolled sample. It is especially interesting that there is a significant increase in the average GND density for the sample that experienced the fatigue loading with a low-maximum applied stress, and the maximum applied stress was lower than the material’s yield strength. This observation implies that the grain boundary mediated the dislocation-emission mechanism.
Theoretically, the strength of glass processed in vacuum should be higher due to outgassing of contaminants normally present in the glass, such as bulk water in the form of OH bonds that tends to weaken the glass structure. In this research, small discs of a few types of glass have been subjected to various temperatures for extended periods of time in vacuum. Their strength was then tested using a standard flexure technique, facilitated by a custom-designed test fixture, and the results were compared to glass tested in air using the same fixture. The purpose of the glass blowing investigation was to prove the basic feasibility of a high-level concept for in-space manufacture of optical elements. The central requirement was that the glass bubble had to be blown into a support structure such that the bubble could be handled by manipulation of the structure. The blown bubble attached itself to a mullite ring geometrically and mechanically, as a demonstration in the initial experiments described here, by expanding through and around it. The vacuum system used was custom made, as were most of the components of the system, such as the heating element, the glass and ring support structure, and the gas inlet system that provided the pressure needed to blow the glass.
This article presents a survey of the core-congruential formulation (CCF) for geometrically nonlinear mechanical finite elements based on the total Lagrangian (TL) kinematic description. Although the key ideas behind the CCF can be traced back to Rajasekaran and Murray in 1973, it has not subsequently received serious attention. The CCF is distinguished by a two-phase development of the finite element stiffness equations. The initial phase developed equations for individual particles. These equations are expressed in terms of displacement gradients as degrees of freedom. The second phase involves congruential-type transformations that eventually binds the element particles of an individual element in terms of its node-displacement degrees of freedom. Two versions of the CCF, labeled direct and generalized, are distinguished. The direct CCF (DCCF) is first described in general form and then applied to the derivation of geometrically nonlinear bar, and plane stress elements using the Green-Lagrange strain measure. The more complex generalized CCF (GCCF) is described and applied to the derivation of 2D and 3D Timoshenko beam elements. Several advantages of the CCF, notably the physically clean separation of material and geometric stiffnesses, and its independence with respect to the ultimate choice of shape functions and element degrees of freedom, are noted. Application examples involving very large motions solved with the 3D beam element display the range of applicability of this formulation, which transcends the kinematic limitations commonly attributed to the TL description.
To investigate unsaturated soil-water interaction at micro-scale, this work extends the numerical manifold method (NMM) by incorporating a soil-water coupling model considering specific capillary water distribution and capillary force calculation. The soil skeleton is constructed by a soil skeleton generation algorithm with random polygons. To more realistically capture the interaction between soil grains and capillary water, a capillary mechanics-based geometric algorithm is proposed to iteratively calculate the capillary water distribution. The capillary forces corresponding to the capillary water distribution are calculated based on the Young-Laplace equation. The proposed capillary water solving framework is first verified by reproducing the soil-water characteristic curve and the capillary water distribution of an ideal contact-disk model against analytical solutions. To further validate the ability of the capillary water solving framework to predict hydraulic behavior of the real soil, a laboratory test on the Toyoura sand is reproduced numerically. Then an ideal direct shear test is performed to further validate the two-way soil-water coupling procedure, in which a comparison between the numerical and analytical results regarding the shear strength and matric suction is presented. Finally, microscopic hydraulic and compression tests are conducted on two soil specimens with the same porosity and mean grain diameter but different uniformity coefficients. The results elucidate that the extended method is a potential tool to explore unsaturated soil behaviors at micro-scale.
The United States High Performance Research Reactor project is tasked with fuel development and qualification leading to conversion of higher power research and test reactors in the US from high-enriched uranium (HEU) to low-enriched uranium (LEU) fuels. Here, this manuscript identifies the functional and operational design requirements of the first miniature test plate (mini-plate [MP]) irradiation campaign (MP-1) of commercially fabricated LEU U-10Mo monolithic plate-type fuel and is the precursor to a large parametric mini-plate test (MP-2) aimed at producing the data to support regulatory qualification of the LEU U-10Mo monolithic fuel. The manuscript (1) provides a general description of the selected U-10Mo LEU fuel and (2) defines the overall experiment design and functional requirements to accomplish the specific test objective of MP-1, which is to confirm that the commercially manufactured LEU U-10Mo monolithic fuel meets the established requirements of geometric stability, mechanical integrity and stable and predictable behavior. The fuel testing parameters are established by the need to bound performance behavior within the operational envelope of the reactors being converted.
High velocity collisions can cause melting, a situation involving elasticity, plasticity, and phase transformations in tandem, which is very challenging for mechanical modeling. Here, this paper advances the use of smoothed particle hydrodynamics (SPH) simulations for this complex situation, investigating stainless steel microspheres impacting a tin substrate. After extensive validation against a range of mechanical and geometrical measurements of experimental impacts, the model provides a detailed, quantitative view of melt volume production over a broad range of impact velocities. It also enables direct decomposition of the impact energy partitioning, including analysis of the energy absorbed by plasticity in the impactor, the system kinetic energy, and the adiabatic heating of the substrate in addition to the melting itself. Since these quantities are not available experimentally, SPH modeling substantially enhances the quantitative understanding of such impacts. Finally, the model speaks to the phase characteristics of the ejecta formed during impact, confirming the experimental inference that the ejected material likely contains a significant fraction of solid material.
Fractures are conduits for fluid flow in low-permeability geological formations. Multiphase flow properties of fractures are important in natural processes and in engineering applications such as the evaluation of the sealing capacity of caprocks and productivity of hydrocarbon-bearing tight rocks. Investigations of flow and transport through fractures typically focus on the effects of fracture geometric and mechanical factors such as aperture, roughness, and compressibility. The wettability of the fracture surfaces and its influence on microscale interfacial phenomena and macroscale effective transport properties are seldom studied. In this study, we investigated the effect of heterogeneous wetting properties on the displacement of water by supercritical CO2 through a series of lattice Boltzmann method simulations. The results show the evolution of the CO2 plume within a fracture is controlled by both the roughness of the aperture field and the wetting distribution. We combined these factors into a capillary pressure map that can be related to the macroscopic flow behavior of the fracture. We observed that heterogeneous wetting distributions promote the residual trapping of water where lower capillary pressures allowed for isolated water pockets in higher capillary pressure zones. Analysis of fracture unsteady relative permeability shows the effect of wetting on permeability evolution and provides support for the viscous-coupling relative permeability model. Finally, analysis of the steady-state relative permeability and saturation demonstrates a strong correlation between permeability and the standard deviation of the capillary pressure field. Thus, characterizing the distribution of wetting properties of fractures is crucial to understanding multiphase fracture flow and transport properties.
Flows and rotation, particularly E×B rotation, are critical to improving plasma performance, and waves are a primary tool of plasma control. Thus, it is paramount to understand under what conditions waves can drive E×B flows in plasmas. In this didactic review, an invited paper accompanying the 2023 Marshall N. Rosenbluth Doctoral Thesis Award, this question is answered in the context of momentum-conserving quasilinear theory. There are two primary frameworks for momentum-conserving quasilinear theories that can handle both resonant and nonresonant particles: Eulerian averaging theories and oscillation-center Hamiltonian theories. There are also two different paradigmatic wave problems: plane-wave initial value problems, and steady-state boundary value problems. Here, it is shown that each of these frameworks “naturally” works better with a different problem type. By using these theories, one finds a great difference in the behavior of time- vs space-dependent waves. A time-evolving plane wave can only drive flow if the electromagnetic momentum of the wave, given by the Poynting flux, changes. This result precludes flow drive by any planar electrostatic wave. In contrast, a steady-state spatially evolving wave can drive flow whenever there is divergence in the flux of Minkowski momentum, a completely different physical quantity. This review aims to provide a high-level, intuitive understanding of the very different behaviors observed for these two types of problem.