Engineering Papers⌕ Search

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.

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

At least 91 records · Page 5

A high-frequency, low-power resonant radio-frequency neutron spin flipper for high-resolution spectroscopy

Here, we present a resonant-mode, transverse-field, radio-frequency (rf) neutron spin flipper design that uses high-temperature superconducting films to ensure sharp transitions between uniform magnetic field regions. Resonant mode allows for low-power, high-frequency operation but requires strict homogeneity of the magnetic fields inside the device. This design was found to efficiently flip neutrons at 96.6 ± 0.6% at an effective frequency of 4 MHz in bootstrap configuration with a beam size of 2.4 × 2.5 cm 2 and a wavelength of 0.4 nm. The high frequency and efficiency enable this device to perform high-resolution neutron spectroscopy with comparable performance with currently implemented rf flipper designs. The limitation of the maximum frequency was found due to the field homogeneity of the device. We numerically analyze the maximum possible efficiency of this design using a Bloch solver simulation with magnetic fields generated from finite-element simulations. We also discuss future improvements of the efficiency and frequency to the design based on the experimental and simulation results.

47 OTHER INSTRUMENTATION↗

Machine learning–assisted prediction of heat fluxes through thermally anisotropic building envelopes

Thermally anisotropic building envelope (TABE) is a novel active building envelope that can save energy use to maintain thermal comfort in buildings by redirecting heat and coolness from building envelopes to thermal loops. Finite element models (FEMs) can be used to compute the heat fluxes through TABEs, but the high computational cost of finite element simulations has prevented parametric studies and design optimizations. This paper proposes a domain knowledge–informed, finite element–based machine learning framework to reduce the computation cost for the energy management of buildings installed with TABE that uses a ground thermal loop. First, the training heat flux data set was generated by FEM simulations with different thermal loop schedules. Then, both shallow learning models (i.e., multivariate linear regression and eXtreme Gradient Boost, or XGBoost) and a deep learning model (i.e., deep neural network, or DNN) were trained to predict the heat fluxes. Domain knowledge was used for data preprocessing and feature selection. Finally, the suitability of the selected machine learning model was tested under different thermal loop schedules. Herein, the case study results showed that: (1) XGBoost can be as accurate as DNN (coefficient of determination equal to 0.81) with much less training time; (2) the annual energy cost savings for different thermal loop schedules obtained by the XGBoost-predicted and FEM-calculated heat fluxes are consistent, having a difference of only 4%; and (3) XGBoost can reduce the computation time for the annual energy analysis of the case study building with a given thermal loop schedule from around 12 h by using FEM to less than 1 min.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Toughening of interface networks through the introduction of weak links

The resistance of polycrystalline materials to intergranular fracture depends on the toughness of the interfaces between neighboring grains. Using finite element simulations on model, two-dimensional microstructures, we demonstrate that the steady-state crack growth resistance of a polycrystal that contains exclusively high toughness interfaces is suboptimal: its crack growth resistance may be increased by introducing a minority fraction of weak interfaces at random locations within the interface network. Furthermore, the weaker are these interfaces, the greater is their toughening effect. These weak interfaces toughen the material through crack deflection and branching, nucleation of secondary cracks, and inelastic dissipation within intact interfaces. Furthermore, our findings suggest that fracture resistance-oriented material development should not aim to eliminate all weak interfaces, but rather judiciously incorporate a low fraction of them.

36 MATERIALS SCIENCE↗

Nanoscopic Imaging of Self-Propelled Ultrasmall Catalytic Nanomotors

Ultrasmall nanomotors (<100 nm) are highly desirable nanomachines for their size-specific advantages over their larger counterparts in applications spanning nanomedicine, directed assembly, active sensing, and environmental remediation. While there are extensive studies on motors larger than 100 nm, the design and understanding of ultrasmall nanomotors have been scant due to the lack of high-resolution imaging of their propelled motions with orientation and shape details resolved. Here, we report the imaging of the propelled motions of catalytically powered ultrasmall nanomotors─hundreds of them─at the nanometer resolution using liquid-phase transmission electron microscopy. These nanomotors are Pt nanoparticles of asymmetric shapes (“tadpoles” and “boomerangs”), which are colloidally synthesized and observed to be fueled by the catalyzed decomposition of NaBH4 in solution. Statistical analysis of the orientation and position trajectories of fueled and unfueled motors, coupled with finite element simulation, reveals that the shape asymmetry alone is sufficient to induce local chemical concentration gradient and self-diffusiophoresis to act against random Brownian motion. Our work elucidates the colloidal design and fundamental forces involved in the motions of ultrasmall nanomotors, which hold promise as active nanomachines to perform tasks in confined environments such as drug delivery and chemical sensing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Automated design of an additive manufactured compact broadband antenna for plasma reflectometry

Broadband antennas operating in the gigahertz frequency range are regularly used for plasma reflectometry diagnostics. Due to a lack of space and unique diagnostic constraints, these antennas are often custom in design and frequency range. Recent advances in additive manufacturing of high temperature copper alloys allow for expanded freedom in design of these diagnostic antennas. In this work, a heuristic simulated annealing algorithm is used alongside 3-D finite element simulation to automate the design of a double ridged rectangular horn antenna for a reflectometry diagnostic on the DIII-D tokamak. Optimization of antenna performance given the design constraints results in a compact broadband (6-20+ GHz) antenna design. Measured transmission from the additively manufactured antenna matches simulation within reasonable error, and experimental plasma electron density profiles from the DIII-D high-field side scrape-off layer are shown.

Additive manufacturing↗

Phase Shift Control of a Three-Phase Inverter for Balanced Secondary Currents in Misaligned Three-Phase Inductive Power Transfer Systems

This paper presents a strategy for generating fundamental reference voltage commands for a double-sided LCC tuned three-phase wireless power transfer (WPT) system to achieve balanced output currents when the receiver and transmitter are not perfectly aligned. Formulas for translating the normalized fundamental current commands into phase-shift commands for two of the phase legs in a two-level three-phase inverter are derived. The methodology is validated through circuit simulations performed using coupler models derived from finite-element simulations. Simulation results for a 1.5 k W system design are included, which indicate near perfect correction of imbalance in the secondary currents for a typical misalignment case.

Pries, Jason↗

Crushing behaviors of unidirectional carbon fiber reinforced plastic composites under dynamic bending and axial crushing loading

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

42 ENGINEERING↗

TRUST Sensors in Environments: Thermocouples (SE-TC), Release FY25

The Delivery Environments Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is a broad project intended to analyze simplified problems experimentally and with modeling and simulation. The purpose of analyzing these simplified problems is to extend solution methods to more complex problems, as well as understand deficiencies and gaps in knowledge of methods currently used in more complex analyses. The TRUST project encompasses several smaller testbeds intended to isolate individual phenomena. The testbed under consideration in this report is the Sensors in Environments: Thermocouples testbed. In previous years, the purpose of this testbed was to quantify uncertainty of thermocouple sensors. To accomplish this, an aluminum plate was placed in a thermal chamber and subject to various types of thermal loading. Thermocouples were placed in various locations on the aluminum plate in various configurations (e.g., embedded in the plate, placed under Kapton tape), and an effort was made to quantify uncertainty in these measurements. Finite element simulations were performed to investigate how sensitive these measurements were to parameters such as the boundary conditions on the plate and material properties. However, a fundamental source of uncertainty in this analysis was the convective heat transfer from the plate. Convective heat transfer is a complex physical phenomenon comprised of a number of interacting sub-processes and is difficult to predict accurately a priori. As such, the main purpose of this testbed in FY25 was to better understand, both experimentally and numerically, the convective heat transfer from the plate. This is a highly applicable problem to several more complex problems, as convective heat transfer occurs in nearly all problems where a body is moving through air. Numerically, this required a two-step approach. First, the air flow in the thermal chamber was in vestigated using computational fluid dynamics. The commercial solver Fluent was used to perform these simulations. From these simulations, a heat transfer coefficient over the surface of the plate was calculated. This heat transfer was then used as boundary conditions for finite element heat transfer simulations within the plate, which were performed using Abaqus. Significant effort was devoted to automating the handoff between these two solvers. Experimentally, previous thermocouple results in the plate were used to validate the time-dependent thermal profiles produced from Abaqus. Further experimental efforts were performed both to help validate the Fluent simulations and to inform its boundary conditions. For example, hot-wire anemometers were used to measure the velocity in the chamber, which would be particularly useful in understanding the chamber inlet velocity. Thermocouple measurements were also taken in the chamber, instead of only on the plate, to serve as validation evidence for the Fluent simulations. Numerical results showed that the Fluent to Abaqus workflow matched previous plate thermocouple measurements well. This type of handoff is useful for more complex experiments, or those that are not able to be examined in as great of detail as this testbed, as it was performed without any experimental input. Experimental results, however, were more mixed. The anemometers proved unreliable, with inconsistent measurements across all anemometers, even at locations that were nearly identical. On the other hand, the thermocouples provided a relatively rich view of the temperature field in the chamber.

42 ENGINEERING↗

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗

Parametric Finite Element Analysis of Naturally Corroded Steel Specimens Using 3D Surface Laser Scans

Corrosion is considered a uniform thickness reduction design guideline of the maritime industry. However, additionally, the corroded and irregular morphology of the surface affects the steel's load-bearing capacity and its impact on the strength and elongation behaviour of the steel is not yet fully understood. These effects on the local behaviour of steel structures under tensile loading were investigated with tensile tests on naturally corroded steel specimens and nonlinear finite element simulations including the corroded surface morphology with a uniform surface idealation. The models also include the deformed specimen shape. The developed approach led to highly accurate parametric finite element models predicting the ultimate tensile strength and longitudinal position of fracture. The results show that all included aspects are essential for accurate simulations, while solely the maximum available surface resolution was not as decisive.

corrosion↗

Strain localization criteria for viscoplastic geomaterials

This report presents a viscoplastic localization criterion to detect quasi-instantaneous (i.e., load-induced) and delayed (creep-induced) strain localization in rate-dependent solids. The study is based on the theory of controllability and a viscoplastic description of the mechanical response. Analytical precursors of unstable states are defined through systems of ordinary differential equations (OEDs). The use of the proposed criteria is illustrated at the material point level through a set of strain localization analyses simulating active strain localization of a porous rock. In addition, full-field finite element simulations of compression tests conducted under various pressures are reported to demonstrate the role of local unstable viscoplasticity in the spontaneous propagation of deformation bands under stationary boundary conditions. The study shows that the viscoplastic localization criterion maintains a negative sign as long as the behavior is unstable, that is, the rate of deformation is accelerating. The sign switch coincides with the transition to decelerating deformation. The analyses revealed that pulses of overstress always emerge in correspondence with the growth of unstable behavior, and the peak matches the transition to stable behavior. The local responses recovered from full-field analyses were consistent with those observed in analyses at material point level and the predictions of the presented theory.

36 MATERIALS SCIENCE↗

A Comparative Study of the Efficacy of Local/Global and Parametric/Nonparametric Machine Learning Methods for Establishing Structure–Property Linkages in High-Contrast 3D Elastic Composites

Reduced-order structure–property (S-P) linkages play a pivotal role in the tailored design of materials for advanced engineering components. There is a critical need to distill these from the simulation datasets aggregated using sophisticated, computationally expensive, physics-based numerical tools (e.g., finite element methods). The recent emergence of materials data science approaches has opened new avenues for addressing this challenge. In this paper, we critically compare the relative merits of the application of four distinct machine learning approaches for their efficacy in extracting microstructure-property linkages from the finite element simulation data aggregated on high-contrast elastic composites with different microstructures. The machine learning approaches selected for the study have included different combinations of local/global and parametric/nonparametric approaches. Furthermore, the nonparametric approaches selected for this study are based on Gaussian Process (GP) models that allow for a formal treatment of uncertainty quantification in the predicted values. The predictive performances of these different approaches have been compared against each other using rigorous cross-validation error metrics. Furthermore, their sensitivity to both the dataset size and dimensionality has been investigated.

36 MATERIALS SCIENCE↗

FEA-aided investigation of the effective thermal conductivity in a medium with embedded spheres

For multiple applications in nuclear energy, the ability to accurately represent material behavior with a simplified model is important to facilitate practical engineering-scale simulations. Here, we focus on the homogenized thermal response of a medium containing spherical inclusions, similar to a fuel form (compact or pebble) containing (TRISO) particles. A review on effective thermal conductivity (ETC) modeling is performed considering a random distribution of mono-sized spherical inclusions in a continuous matrix, with a primary focus on the analytical models. Finite element simulations are performed to evaluate each analytical model in varying material conditions. The model predictions are compared with the expected results obtained from the finite element predictions in addition to the Wiener and Hashin-Shtrikman bounds. Lastly, we included a practical discussion on the homogenization applied to a TRISO fuel pebble.

42 ENGINEERING↗

Thermal-Mechanical Elastic-Plastic and Ductile Failure Model Calibrations for 304L Stainless Steel Alloy

Numerical simulations of metallic structures undergoing rapid loading into the plastic range require material models that accurately represent the response. In general, the material response can be seen as having four interrelated parts: the baseline response under slow loading, the effect of strain rate, the conversion of plastic work into heat and the effect of temperature. In essence, the material behaves in a thermal-mechanical manner if the loading is fast enough so when heat is generated by plastic deformation it raises the temperature and therefore influences the mechanical response. In these cases, appropriate models that can capture the aspects listed above are necessary. The material of interest here is 304L stainless steel, and the objective of this work is to calibrate thermal-mechanical models: one for the constitutive behavior and another for failure. The work was accomplished by first designing and conducting a material test program to provide data for the calibration of the models. The test program included uniaxial tension tests conducted at room temperature, 150 and 300 C and at strain rates between 10 –4 and 10 3 1/s. It also included notched tension and shear-dominated compression hat tests specifically designed to calibrate the failure model. All test specimens were extracted from a single piece of plate to maintain consistency. The constitutive model adopted was a modular $J_2$ plasticity model with isotropic hardening that included rate and temperature dependence. A criterion for failure initiation based on a critical value of equivalent plastic strain fitted the failure data appropriately and was adopted. Possible ranges of the values of the parameters of the models were determined partially on historical data from calibrations of the same alloy from other lots and are given here. The calibration of the parameters of the models were based on finite element simulations of the various material tests using relatively ne meshes and hexahedral elements. When using the model in structural finite element calculations, however, element formulations and sizes different from those in the calibration are likely to be used. A brief investigation demonstrated that the failure initiation predictions can be particularly sensitive to the element selection and provided an initial guide to compensate for the effect of element size in a specific example.

36 MATERIALS SCIENCE↗

Bayesian calibration of strength model parameters from Taylor impact data

Materials strength plays a key role in determining the mechanical response of engineered structures. As such, accurate strength models are crucial in simulations involving complex loading conditions, particularly when deformation in the plastic regime is deemed important. In this work, a Gaussian process based surrogate for a finite element simulation of a Taylor impact test is developed and used for Bayesian calibration of the Preston–Tonks–Wallace strength model. Here, the surrogate model is shown to closely approximate the salient features of the Taylor cylinder deformation and is validated against simulation output before being used in the strength model calibration routine. The results show that Taylor impact test data can be used in the calibration of constitutive equations through the use of a combination of data science techniques, namely Gaussian processes and Bayesian inference.

36 MATERIALS SCIENCE↗

Far-field thermal imaging below diffraction limit

Non-uniform self-heating and temperature hotspots are major concerns compromising the performance and reliability of submicron electronic and optoelectronic devices. At deep submicron scales where effects such as contact-related artifacts and diffraction limits accurate measurements of temperature hotspots, non-contact thermal characterization can be extremely valuable. In this work, we use a Bayesian optimization framework with generalized Gaussian Markov random field (GGMRF) prior model to obtain accurate full-field temperature distribution of self-heated metal interconnects from their thermoreflectance thermal images (TRI) with spatial resolution 2.5 times below Rayleigh limit for 530nm illumination. Finite element simulations along with TRI experimental data were used to characterize the point spread function of the optical imaging system. In addition, unlike iterative reconstruction algorithms that use ad hoc regularization parameters in their prior models to obtain the best quality image, we used numerical experiments and finite element modeling to estimate the regularization parameter for solving a real experimental inverse problem.

42 ENGINEERING↗

LAROMance Grade 91 Model Integration in NEML2

New reactor designs are targeting higher operating temperatures for increased thermal efficiency when compared to the current fleet of light water reactors. Designing structural components for these high temperature environments with reliable long-term operations requires material models that can accurately capture the deformation mechanisms active in these environments. The LAROMance surrogate material models are based on a database of mechanistic crystal plasticity simulations for high-temperature conditions. Inputs to the LAROMance models reflect the microstructural pedigree of the material, like dislocation densities and precipitate contents. Based on the evolution of these microstructural features, the LAROMance model provides the engineering scale constitutive model response. The LAROMance model was recently parameterized for Grade 91, a high temperature alloy. In the present work, the Grade 91 LAROMance model is implemented in the New Material Model Library, version 2 (NEML2). NEML2 provides a modular way to build material models from smaller blocks and was developed to vectorize the material update to efficiently run on modern computational architectures with graphics processing unit accelerators. NEML2 constitutive models can be used in simulations based on the multiphysics object-oriented simulation environment (MOOSE). This report provides details on the implementation of the Grade 91 LAROMance model in NEML2 and its verification of engineering scale finite element simulations in MOOSE.

42 - ENGINEERING↗

Evolution of Microstructure and Creep Behavior in an Fe-Ni-Cr-Nb-C Alloy during Service in Hydrocarbon Cracker Tubes

Heat resistant, cast austenitic Fe-Ni-Cr-Nb-C alloy tubes are used by the petrochemical industry to transport hydrocarbons at temperatures between 850 and 1050 °C. Exposure to carbon from the process stream and nitrogen from the air causes severe carburization and nitridation in the tube walls, which influences creep resistance and degrades their ductility. A better understanding of these processes would enable their service life to be extended. To this end, we report the results of an experimental and computational study of the changes in microstructure and creep behavior in a representative alloy during service. The microstructure of the fresh material is compared to that of 8- and 10-year-old tubes. The results confirm the presence of M{sub 23}C{sub 6}, M{sub 7}C{sub 3}, and M{sub 2}(CN) phases reported in prior studies, and also reveal the nucleation and growth of the carbo-nitride phase M{sub 6}(CN). We show that the spatial variations of carbide and nitride particle volume fractions in aged tubes can be predicted by a simple reaction–diffusion equation. In addition, the uniaxial compressive creep response of the fresh and aged tubes was obtained as functions of temperature and strain rate. The materials deform by power-law creep, with strain rate versus stress relation as ε-dot =Aσ{sup n}. Aging reduces the stress exponent from n = 8 in fresh tubes to n = 5 in aged tubes, and also reduces the pre-exponential coefficient A. Tests on model alloys with composition identical to that of the matrix combined with finite element simulations of creep in model microstructures reveal that Cr depletion is responsible for the change in stress exponent, while the change in A is caused by the increase in particle volume fraction. By combining the microstructure-scale simulations of creep with the predictions of the reaction–diffusion model of microstructure evolution, the creep behavior of tubes left in service beyond 10 years is estimated.

36 MATERIALS SCIENCE↗