Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Elements”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Illumination-induced modulation of conductivity and Gunn oscillation properties in epitaxial GaAs

In this work, we illuminate a gallium arsenide (GaAs) Gunn device and study the light-induced changes of Gunn oscillation properties. We observe that illumination leads to the modulation of the Gunn threshold voltage, the Gunn oscillation magnitude, and the coherency of Gunn oscillation, with the nature of the modulation being closely related to the position of illumination on the device. These effects are attributed to the generation of optically excited carriers, which results in the modulation of conductivity and the electric field profile along the device. The finite element method is used to simulate the change of the field profile of the Gunn device caused by illumination. We also report an unexpected phenomenon of Gunn oscillation property manipulation with an optical chopper. In addition, wavelength-dependent, power-dependent, and pulsed illumination measurements are performed to help with further understanding the observations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Probing the properties and mechanisms of failure waves in soda-lime glass

Soda-lime glass (SLG) and other silica glasses exhibit the failure wave phenomenon under shock compression. The mechanism responsible for this peculiar behavior of glasses is still unresolved. In this study, a series of plate impact experiments was performed at three different impact stresses of 6.4, 8.3, and 10.8 GPa to better understand the mechanisms underlying the failure wave phenomenon. Specifically, spall experiments were conducted to probe the speed and existence of failure waves at different stresses in SLG. A layered glass target was used to probe the possibility of a “renucleation” of the wave at the SLG–SLG interface. When it existed, the failure wave was inferred to propagate at a speed of 1.3 km/s. However, it was observed that the failure wave phenomenon ceases to exist for impact stresses higher than 10 GPa. In experiments with a 6.4 GPa impact stress, the peak free surface velocity was significantly less than what is predicted by stress-Hugoniot calculations. This velocity deficit and other important features of the measured free surface velocity profiles were simulated using finite element analysis by incorporating an abrupt densification of SLG at a critical stress in the equation of state. Furthermore, this densification feature is similar to what would be expected of a phase transition. Although unable to unambiguously reveal the mechanism causing the failure wave phenomenon, the results of the present work clearly indicate that the failure wave causes a secondary compression and densification in SLG.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Error Estimates for the Optimal Control of a Parabolic Fractional PDE

In this work, we consider the integral definition of the fractional Laplacian and analyze a linear-quadratic optimal control problem for the so-called fractional heat equation; control constraints are also considered. We derive existence and uniqueness results, first order optimality conditions, and regularity estimates for the optimal variables. To discretize the state equation we propose a fully discrete scheme that relies on an implicit finite difference discretization in time combined with a piecewise linear finite element discretization in space. We derive stability results and a novel $L^2(0,T;L^2(\Omega))$ a priori error estimate. On the basis of the aforementioned solution technique, we propose a fully discrete scheme for our optimal control problem that discretizes the control variable with piecewise constant functions, and we derive a priori error estimates for it. We illustrate the theory with one- and two-dimensional numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Characterization of latent image of electron beam resist via critical-dimension resonant soft X-ray scattering

Metrology plays a crucial role in semiconductor manufacturing by providing accurate and precise measurement and characterization of critical parameters. With the development of high-resolution extreme ultraviolet lithography (EUVL) processes, critical dimensions are shrinking to sub-10 nm. Resist materials encounter the challenge of providing heightened sensitivity and a handle on exacerbating stochastic variations. A comprehensive understanding of the chemical profile of the latent image is pivotal for mitigating stochastic effects and optimizing pattern quality. However, the subtle differences in chemistry between the exposed and unexposed regions of the resists make it extremely challenging to characterize the latent images with sub-nanometer precision. Here, we develop the metrology with critical-dimension resonant soft X-ray scattering (CD-RSoXS) to probe the chemical profiles of latent images stored in resist after exposure. The combination of absorption spectroscopy and enhanced scattering contrast makes it possible to characterize the subtle structural and chemical variations in the latent image. Moreover, the results of the measurements are compared with the simulations with a finite element method–based Maxwell solver to extract a detailed profile of the latent and developed images. We demonstrate that the CD-RSoXS technique can provide valuable insights into the high spatial resolution and local chemical sensitivity simultaneously, which is crucial to understanding the resolution limits and stochastic effects in EUVL processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Design, performance, and cost savings of using GF-PC additively manufactured tooling for stamping of HSS 590 sheet metal

Conventional sheet metal forming tooling in the automotive industry is made up of hardened steel and used for mass-production. Prototype tooling made of metal is durable, but it is only used for a small number of parts despite its high cost, contributing heavily to the vehicle development cost. Additive Manufacturing (AM) offers a low-cost alternative for fabricating tooling suitable for low-volume production. Here, this study investigates the feasibility of using AM polymer composite tooling for the stamping of 1.5mm thickness HSS 590 steel sheets through a two-pronged approach – experimental and numerical analysis. Sheet metal stamping experiments were conducted with glass fiber reinforced polycarbonate (GF-PC) AM polymer tooling and their performance was evaluated based on various metrics such as tool deformation and part accuracy. Finite element simulations of the stamping process were performed, which accurately captured the tool deformations seen in experiments with the use of anisotropic material models. The effective use of simulations in optimizing process parameters to achieve the desired final part geometry is also demonstrated. The GF-PC AM polymer composite tooling was found to be well suited for low volume production applications with potential cost savings.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Machine tool cross beam design, fabrication, and testing using metal big area additive manufacturing

This paper describes the application of metal Big Area Additive Manufacturing (mBAAM) to the fabrication of a machine tool cross beam. The replacement of a traditional box design weldment with a new design printed by wire arc additive manufacturing using the MedUSA system at Oak Ridge National Laboratory (ORNL) is detailed. This requires a new design strategy based on the unique mBAAM capabilities. The intent of the new design is to reduce mass, while maintaining the dynamic stiffness. To compare the two designs, the natural frequencies and mode shapes are measured using impact testing and predicted using finite element analysis. It is confirmed that the printed structure dynamics agreed with the numerical model predictions, which demonstrates that it is feasible to model a large-scale mBAAM part and understand its behavior prior to printing. Another notable outcome of this study is that the significant residual stress and distortion in the print indicate that knowledge gaps remain for widespread implementation of mBAAM.

42 ENGINEERING↗
DOE OSTI2021

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Measurement of a radial flow profile with eddy current flow meters and deep neural networks

Eddy current flow meters (ECFMs) measure flows of conductive fluids. Recent interest in ECFMs has increased due to applications in advanced nuclear reactors. ECFMs are well suited for such applications, as they can provide non-invasive measurements of flow in fluids that are often difficult to measure. Traditionally, ECFMs are operated using an alternating current at a single frequency, limiting ECFMs to measure average fluid velocities, blockages, or voids. Here, we expand the capabilities of ECFMs by measuring the fluid radial velocity profile of liquid mercury. To accomplish this, we made several ECFM sensitivity measurements at a range of frequencies. Different frequencies vary the electromagnetic skin depth of the device. By adjusting frequencies, we probed the fluid velocity at various radial locations and constructed a flow-velocity profile. The relationship between the ECFM measurements and velocity profile is nonlinear and requires solving an inverse problem. Using electromagnetic finite-element simulations to train a deep neural network (DNN), we created a model that provides a stable general relationship between the sensitivity measurements of an ECFM and the fluid velocity profile. Using ECFM measurements of liquid mercury, our DNN model calculates a flow profile that agrees well with computational fluid dynamics (CFD) simulations. This technique has potential to improve flow monitoring for optimization, safe operation of conductive fluid loops, and/or validating complex CFD models.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Design, Analysis, Comparison, and Experimental Validation of Insulated Metal Substrates for High-Power Wide-Bandgap Power Modules

Direct bonded copper (DBC) substrates used in power modules have limited heat spreading and manufacturing capability due to ceramic properties and manufacturing technology. The ceramic and copper bonding is also subject to high mechanical stress due to coefficient of thermal expansion mismatch between the copper and the ceramic. For wide-bandgap (WBG) devices, it is of interest exploring new substrate technologies that can overcome some of the challenges of direct bonded copper substrates. In this technical paper, the design, analysis, and comparison of insulated metal substrates (IMSs) for high-power wide-bandgap semiconductor-based power modules are discussed. This paper starts with technical description and discussion of state-of-the-art DBC substrates with different ceramic insulators such as aluminum nitride (AlN), Al 2 O 3 , and Si 3 N 4 . Next, an introduction of IMSs and their material properties, and a design approach for SiC (silicon carbide) metal-oxide-semiconductor field-effect transistor (MOSFET)-based power modules for high-power applications is provided. The influence of dielectric thickness on the power handling capability of the substrate are also discussed. The designed IMS and DBC substrates were characterized in terms of steady-state and transient thermal performance using finite element simulation. Finally, the performance of the IMS and DBC are validated in an experimental setup under different loading and cooling temperature conditions. We report the simulation and experimental results showed that the IMS can provide high steady-state thermal performance for high-power modules based on SiC MOSFETs. Furthermore, the IMS provided enhanced transient thermal performance, which provided a reduced junction temperature when the module is operated at low fundamental output frequencies in traction drive systems.

42 ENGINEERING↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Tailored Graphene Micropatterns by Wafer-Scale Direct Transfer for Flexible Chemical Sensor Platform

2D materials, such as graphene, exhibit great potential as functional materials for numerous novel applications due to their excellent properties. The grafting of conventional micropatterning techniques on new types of electronic devices is required to fully utilize the unique nature of graphene. However, the conventional lithography and polymer-supported transfer methods often induce the contamination and damage of the graphene surface due to polymer residues and harsh wet-transfer conditions. In this work, a novel strategy to obtain micropatterned graphene on polymer substrates using a direct curing process is demonstrated. Employing this method, entirely flexible, transparent, well-defined self-activated graphene sensor arrays, capable of gas discrimination without external heating, are fabricated on 4 in. wafer-scale substrates. Finite element method simulations show the potential of this patterning technique to maximize the performance of the sensor devices when the active channels of the 2D material are suspended and nanoscaled. This study contributes considerably to the development of flexible functional electronic devices based on 2D materials.

2D materials, Graphene, Microscale patterning, Che↗
DOE OSTI2024

Design optimization of lightweight automotive seatback through additive manufacturing compression overmolding of metal polymer composites

With the growing demand for enhanced automotive fuel efficiency and environmental sustainability, there is a need for lightweighting automotive components through innovative design and manufacturing processes. Here, this study leverages a combination of numerical iterative design optimization and hybrid additive manufacturing–compression molding (AM-CM) technique for metal polymer composites to lightweight an automotive seatback. The AM-CM process enables robust mechanical interlocking between metals and composites, boasting high stiffness and strength with low overall density. Replacing metallic components with such metal polymer composites allows for comparable mechanical performance while significantly reducing the overall weight. First, the automotive seatback design space is reduced to critical load carrying regions using topology optimization and high stress concentration areas are identified using finite element analysis. Next, a lightweight metal polymer subcomponent is designed for a high stress concentration region. The full seatback frame with spatially heterogeneous material-specific design is then iteratively optimized to enable enhanced stiffness with minimal weight. Overall, the automotive seatback frame designed with location-specific metal, polymer, and metal polymer composite materials weighs 20% less than the metal-only design while exhibiting similar stiffness.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Development of Lightweight Structural Materials with Improved Properties for Fission Batteries

The notion of a “fission battery” conveys a vision focused on realizing very simple “plug-and-play” nuclear systems that can be integrated into a variety of applications requiring affordable, reliable energy in the form of electricity and/or heat and function without operations and maintenance staff. Fission batteries require lightweight structural materials to increase their mobility, and the lightweight materials must demonstrate structural resilience under various conditions. The objective of this work is to develop lattice structured lightweight structural material featuring a good combination of mechanical properties using advanced modeling and simulation together with an advanced additive manufacturing technique such as laser powder bed fusion. The preliminary results show that different lattice structures and types can be successfully meshed using nTopology software, and the lattice structure data can be successfully transformed to Multiphysics Object-Oriented Simulation (MOOSE) Environment input. Finite Element Analysis (FEA) displays that, at macro/engineering scale simulation, the weight saving design has an obvious effect on tensile behavior such as effective elastic modulus and yield stress. The novel approaches of this work are (1) development of lattice structures for improved mechanical properties using advanced simulation and modeling techniques; and (2) model predictions of the mechanical properties (e.g., strength and stress distribution) of macroscopic materials in order to preliminarily select a lattice structure for additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sensing depths in frequency domain thermoreflectance

In this work, a method is developed to calculate the length into a sample to which a Frequency Domain Thermoreflectance (FDTR) measurement is sensitive. Sensing depth and sensing radius are defined as limiting cases for the spherically spreading FDTR measurement. A finite element model for FDTR measurements is developed in COMSOL multiphysics and used to calculate sensing depth and sensing radius for silicon and silicon dioxide samples for a variety of frequencies and laser spot sizes. The model is compared to experimental FDTR measurements. Design recommendations for sample thickness are made for experiments where semi-infinite sample depth is desirable. For measurements using a metal transducer layer, the recommended sample thickness is three thermal penetration depths, as calculated from the lowest measurement frequency.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multi-scale characterization and simulation of impact welding between immiscible Mg/steel alloys

Vaporizing foil actuator spot welding method is used in this paper to join magnesium alloy AZ31 and uncoated high-strength steel DP590, which are typically considered as un-weldable due to their high physical property disparities, low mutual solubility, and the lack of any intermetallic phases. Characterization results from scanning electron microscopy (SEM) and high-resolution transmission electron microscopy (HRTEM) of the weld interface indicate that the impact creates an Mg nanocrystalline interlayer with abundant Fe particles. The interlayer exhibits intact bonding with both DP590 and AZ31 substrates. To investigate the fundamental bond formation mechanisms at the interface, a finite element (FE)-based process simulation is first performed to calculate the local temperature and deformation at the interface under the given macroscopic experimental condition. Finally, taking the FE results at the interface as inputs, molecular dynamics (MD) simulations are conducted to study the interlayer formation at the Mg/Fe interface during the impact and cooling. The results found a high velocity shearing-induced mechanical mixing mechanism that mixes Mg/Fe atoms at the interface and creates the interlayer, leading to the metallurgical bond between Mg/steel alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Approach to nonlinear magnetohydrodynamic simulations in stellarator geometry

The capability to model the nonlinear magnetohydrodynamic (MHD) evolution of stellarator plasmas is developed by extending the M3D-C 1 code to allow non-axisymmetric domain geometry. We introduce a set of logical coordinates, in which the computational domain is axisymmetric, to utilize the existing finite-element framework of M3D-C 1 . A C 1 coordinate mapping connects the logical domain to the non-axisymmetric physical domain, where we use the M3D-C 1 extended MHD models essentially without modifications. We present several numerical verifications on the implementation of this approach, including simulations of the heating, destabilization, and equilibration of a stellarator plasma with strongly anisotropic thermal conductivity, and of the relaxation of stellarator equilibria to integrable and non-integrable magnetic field configurations in realistic geometries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗