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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Tuning the local temperature during feedback controlled electromigration in gold nanowires

Feedback controlled electromigration (FCE) in metallic nanowires has been widely used for various purposes. However, the control of the local temperature during FCE remains a challenge. Here, we report that the environment temperature can be used as a knob to tune the local temperature during FCE in gold nanowires. FCE was performed in gold nanowires at various environment temperatures ranging from 4.2 K to 300 K. We find that the dissipated power normalized by the cross section area of the nano constriction is linearly proportional to the environment temperature. Interestingly, the estimated local maximum temperature parabolically depends on the environment temperature. A minimum in the local temperature can be reached if an appropriate environment temperature is chosen. Our findings are well supported by the finite element simulation. Moreover, the data indicates the coupling between FCE triggering current density and local temperature.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2016

Composition and Manufacturing Effects on Electrical Conductivity of Li/FeS 2 Thermal Battery Cathodes

The electrical conductivity is key to the performance of thermal battery cathodes. In this work we present the effects of manufacturing and processing conditions on the electrical conductivity of Li/FeS2 thermal battery cathodes. Finite element simulations were used to compute the conductivity of three-dimensional microcomputed tomography cathode microstructures and compare results to experimental impedance spectroscopy measurements. A regression analysis reveals a predictive relationship between composition, processing conditions, and electrical conductivity; a trend which is largely erased after thermally-induced deformation. Moreover, the trend applies to both experimental and simulation results, although is not as apparent in simulations. This research is a step toward a more fundamental understanding of the effects of processing and composition on thermal battery component microstructure, properties, and performance.

25 ENERGY STORAGE↗
DOE OSTI2021

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↗
DOE OSTI2018

Generalized radial-return mapping algorithm for anisotropic von Mises plasticity framed in material eigenspace

A computationally efficient integration algorithm for anisotropic plasticity is proposed, which is identified as a generalization of the radial-return mapping algorithm to anisotropy. The algorithm is based upon formulation within the eigenspace of a material anisotropy tensor associated with anisotropic quadratic von Mises (J 2 ) plasticity (also called Hill plasticity), for which it is shown to ensure that the flow rule remains associative, ie, the normality condition is satisfied. Extension of the algorithm to include anisotropic elasticity (anisotropic elastoplasticity) is further provided, made possible by the identification of a certain fourth-order material tensor dependent on both the elastic and plastic anisotropy. The derivation of the fully elastoplastically anisotropic algorithm involves further complexity, but the resulting algorithm is shown to closely resemble the purely plastically anisotropic one once the appropriate eigenspace is identified. The proposed generalized radial-return algorithm is compared to a classical closest-point projection algorithm, for which it is shown to provide considerable advantage in computational cost. In conclusion, the efficiency, accuracy, and robustness of the algorithm are demonstrated through various illustrative test cases and in the finite element simulation of Taylor impact tests on tantalum.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Millimeter wave detection via Autler-Townes splitting in rubidium Rydberg atoms

In this paper, we demonstrate the detection of millimeter waves via Autler-Townes splitting in {sup 85}Rb Rydberg atoms. This method may provide an independent, atom-based, SI-traceable method for measuring mm-wave electric fields, which addresses a gap in current calibration techniques in the mm-wave regime. The electric-field amplitude within a rubidium vapor cell in the WR-10 wave guide band is measured for frequencies of 93.71 GHz and 104.77 GHz. Relevant aspects of Autler-Townes splitting originating from a four-level electromagnetically induced transparency scheme are discussed. We measured the E-field generated by an open-ended waveguide using this technique. Experimental results are compared to a full-wave finite element simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1979

Two-dimensional finite-element analyses of simulated rotor-fragment impacts against rings and beams compared with experiments

Finite element modeling alternatives as well as the utility and limitations of the two dimensional structural response computer code CIVM-JET 4B for predicting the transient, large deflection, elastic plastic, structural responses of two dimensional beam and/or ring structures which are subjected to rigid fragment impact were investigated. The applicability of the CIVM-JET 4B analysis and code for the prediction of steel containment ring response to impact by complex deformable fragments from a trihub burst of a T58 turbine rotor was studied. Dimensional analysis considerations were used in a parametric examination of data from engine rotor burst containment experiments and data from sphere beam impact experiments. The use of the CIVM-JET 4B computer code for making parametric structural response studies on both fragment-containment structure and fragment-deflector structure was illustrated. Modifications to the analysis/computation procedure were developed to alleviate restrictions.

Stagliano, T. R.↗
DOE OSTI2015

Crystal plasticity analysis of stress partitioning mechanisms and their microstructural dependence in advanced steels

Two-phase advanced steels contain an optimized combination of high yield strength and large elongation strain at failure, as a result of stress partitioning between a hard phase (martensite) and a ductile phase (ferrite or austenite). Provided with strong interfaces between the constituent phases, the failure in the brittle martensite phase will be delayed by the surrounding geometric constraints, while the rule of mixture will dictate a large strength of the composite. To this end, the microstructural design of these composites is imperative especially in terms of the stress partitioning mechanisms among the constituent phases. Based on the characteristic microstructures of dual phase and multilayered steels, two polycrystalline aggregate models are constructed to simulate the microscopic lattice strain evolution of these materials during uniaxial tensile tests. By comparing the lattice strain evolution from crystal plasticity finite element simulations with advanced in situ diffraction measurements in literature, this study investigates the correlations between the material microstructure and the micromechanical interactions on the intergranular and interphase levels. Finally, it is found that although the applied stress will be ultimately accommodated by the hard phase and hard grain families, the sequence of the stress partitioning on grain and phase levels can be altered by microstructural designs. Implications of these findings on delaying localized failure are also discussed.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Interaction of Delaminations and Matrix Cracks in a CFRP Plate, Part II: Simulation Using an Enriched Shell Finite Element Model

Numerical simulations are presented of a recently developed test which creates multiple delaminations in a CFRP laminate specimen that grow and interact via transverse matrix cracks. A novel shell element enriched with the Floating Node Method, and a damage algorithm based on the Virtual Crack Closure Technique, were used to successfully simulate the tests. Additionally, a 3D high mesh fidelity model based on cohesive zones and continuum damage mechanics was used to simulate the tests and act as a representative of other similar state-of-the-art high mesh fidelity modeling techniques to compare to the enriched shell element. The enriched shell and high mesh fidelity models had similar levels of accuracy and generally matched the experimental data. With runtimes of 36 minutes for the shell model and 55 hours for the high mesh fidelity model, the shell model is 92 times faster than the high- fidelity simulation.

Mark W McElroy↗
DOE OSTI2016

Thermal conductivity measurement of amorphous dielectric multilayers for phase-change memory power reduction

In this work, we investigate the temperature-dependent thermal conductivities of few nanometer thick alternating stacks of amorphous dielectrics, specifically SiO{sub 2}/Al{sub 2}O{sub 3} and SiO{sub 2}/Si{sub 3}N{sub 4}. Experiments using steady-state Joule-heating and electrical thermometry, while using a micro-miniature refrigerator over a wide temperature range (100–500 K), show that amorphous thin-film multilayer SiO{sub 2}/Si{sub 3}N{sub 4} and SiO{sub 2}/Al{sub 2}O{sub 3} exhibit through-plane room temperature effective thermal conductivities of about 1.14 and 0.48 W/(m × K), respectively. In the case of SiO{sub 2}/Al{sub 2}O{sub 3}, the reduced conductivity is attributed to lowered film density (7.03 → 5.44 × 10{sup 28 }m{sup –3} for SiO{sub 2} and 10.2 → 8.27 × 10{sup 28 }m{sup –3} for Al{sub 2}O{sub 3}) caused by atomic layer deposition of thin-films as well as a small, finite, and repeating thermal boundary resistance (TBR) of 1.5 m{sup 2} K/GW between dielectric layers. Molecular dynamics simulations reveal that vibrational mismatch between amorphous oxide layers is small, and that the TBR between layers is largely due to imperfect interfaces. Finally, the impact of using this multilayer dielectric in a dash-type phase-change memory device is studied using finite-element simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Multi-mode interference revealed by two photon absorption in silicon rich SiO{sub 2} waveguides

Photoluminescence (PL) from Si nanocrystals (NCs) excited by two-photon absorption (TPA) has been observed in Si nanocrystal-based waveguides fabricated by plasma enhanced chemical vapor deposition. The TPA excited photoluminescence emission resembles the one-photon excited photoluminescence arising from inter-band transitions in the quantum confined Si nanocrystals. By measuring the non-linear transmission of waveguides, a large TPA coefficient of β up to 10{sup −8 }cm/W has been measured at 1550 nm. These values of β depend on the Si NCs size and are two orders of magnitude larger than the bulk silicon value. Here, we propose to use the TPA excited visible PL emission as a tool to map the spatial intensity profile of the 1550 nm propagating optical modes in multimode waveguides. In this way, multimode interference has been revealed experimentally and confirmed through a finite element simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

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↗
DOE OSTI2015

Measurement and simulation of thermoelectric efficiency for single leg

Thermoelectric efficiency measurements were carried out on n-type bismuth telluride legs with the hot-side temperature at 100 and 150°C. The electric power and heat flow were measured individually. Water coolant was utilized to maintain the cold-side temperature and to measure heat flow out of the cold side. Leg length and vacuum pressure were studied in terms of temperature difference across the leg, open-circuit voltage, internal resistance, and heat flow. Finite-element simulation on thermoelectric generation was performed in COMSOL Multiphysics, by inputting two-side temperatures and thermoelectric material properties. The open-circuit voltage and resistance were in good agreement between the measurement and simulation. Much larger heat flows were found in measurements, since they were comprised of conductive, convective, and radiative contributions. Parasitic heat flow was measured in the absence of bismuth telluride leg, and the conductive heat flow was then available. Finally, the maximum thermoelectric efficiency was derived in accordance with the electric power and the conductive heat flow.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Direct current dielectric barrier assistant discharge to get homogeneous plasma in capacitive coupled discharge

In this paper, we propose a method to get more homogeneous plasma in the geometrically asymmetric capacitive coupled plasma (CCP) discharge. The dielectric barrier discharge (DBD) is used for the auxiliary discharge system to improve the homogeneity of the geometrically asymmetric CCP discharge. The single Langmuir probe measurement shows that the DBD can increase the electron density in the low density volume, where the DBD electrodes are mounted, when the pressure is higher than 5 Pa. By this manner, we are able to improve the homogeneity of the plasma production and increase the overall density in the target volume. At last, the finite element simulation results show that the DC bias, applied to the DBD electrodes, can increase the homogeneity of the electron density in the CCP discharge. The simulation results show a good agreement with the experiment results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Stability Considerations Using a Microscopic Stability Model Applied to a 2G Thin Film Coated Superconductor

A numerical, finite element simulation of local temperature fields, critical current and transport current distributions is applied to a 2G thin film coated high-temperature superconductor. The focus is on simulation of quenches originating from transport current locally exceeding critical current density. As in previous reports, a statistical treatment of superconductor parameters (critical temperature, current and magnetic field) is applied for this analysis; it shall take into account uncertainties possibly arising from shortcomings in conductor manufacture or during measurements of their properties. The results of the calculations are used as quasi-boundary (driving) conditions for a subsequent transient microscopic numerical stability analysis. Traditionally, stability analysis usually applies analytic stability functions that incorporate conventional, phonon-related timescales, t and disturbances. The question is whether decay of electron pairs and subsequent relaxation of the excited state to a new dynamic equilibrium carrier concentration, the “electron aspect” of the stability problem, under the same disturbances, proceeds on another timescale, t{sup ′}. Is this timescale identical to the traditional (phonon) timescale, t? If not, how large are the differences? The recently reported microscopic stability model, now applied to a thin film conductor, is consulted to find an answer to these questions. A time limit, t{sub Quench}, can be extracted from the simulations, as the time of immediate onset of a quench. This is a new approach in stability considerations since, conventionally, a temperature limit, T < T{sub Crit}, is set as stability criterion. If electrical operation or cooling conditions cannot immediately respond to this challenge, the thin film superconductor, if this time limit is exceeded, will hardly be able to return to zero loss current transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Time Scale Effects on the Temperature Coefficient of Reactivity

For a nuclear reactor, the temperature coefficient of reactivity (TCR) is of significance for both power control and safety analysis, but it is not always a constant. For example, the TCR may be dependent upon temperature due to Doppler broadening, and may be also dependent upon time scale due to heat transfer and thermal expansion. A fast burst reactor is a favorable nuclear facility to simulate the time scale effects on the TCR. Fast burst reactors are generally unmoderated research reactors of enriched uranium metal or uranium alloys. The neutrons are fast and the Doppler broadening effects can be neglected, and the TCR is mainly contributed by fuel thermal expansion. Fast burst reactors can operate on both steady state power mode and burst power mode. In steady state power mode, the systems are assumed to be in thermal equilibrium. In burst mode, the full width at half maximum (FWHM) of the power is from tens of microseconds to a few milliseconds, and the heat transfer is generally neglected. If FWHM is short enough, thermal expansion will have a time delay due to inertial effect. This paper is aimed to study the exact time scale effects on the TCR due to thermal expansion. The TCR contributed by the Doppler broadening effects was neglected, because it is relatively small (The k{sub eff} as a function of temperature was also calculated using the Monte Carlo particle transport code JMCT by our group, and no definite difference due to the temperature change was obtained near room temperature), and there is no method to evaluate the time scale dependence of the Doppler broadening effects yet. In the time scale from 1 x 10{sup -6} s to 1 x 10{sup 5} s, the TCR of a Godiva-like spherical fast burst reactor was calculated with perturbation theory. A Godiva-like spherical fast burst reactor was studied, the calculation of the time scale effects on TCR involved thermal conduction, thermoelasticity and reactivity variation analysis. Both the thermal and the thermoelastic problems were solved using ANSYS software. ANSYS is software implementation of finite element simulations, which is widely used in engineering field. The reactivity variation due to thermal expansion was calculated with perturbation theory.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

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↗
NASA NTRS2016

Microgravity-Induced Physiological Fluid Redistribution: Computational Analysis to Assess Influence of Physiological Parameters

Space flight impacts human physiology in many ways, the most immediate being the marked cephalad (headward) shift of fluid upon introduction into the microgravity environment. This physiological response to microgravity points to the redistribution of blood and interstitial fluid as a major factor in the loss of venous tone and reduction in heart muscle efficiency which impact astronaut performance. In addition, researchers have hypothesized that a reduction in astronaut visual acuity, part of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, is associated with this redistribution of fluid. VIIP arises within several months of beginning space flight and includes a variety of ophthalmic changes including posterior globe flattening, distension of the optic nerve sheath, and kinking of the optic nerve. We utilize a suite of lumped parameter models to simulate microgravity-induced fluid redistribution in the cardiovascular, central nervous and ocular systems to provide initial and boundary data to a 3D finite element simulation of ocular biomechanics in VIIP. Specifically, the lumped parameter cardiovascular model acts as the primary means of establishing how microgravity, and the associated lack of hydrostatic gradient, impacts fluid redistribution. The cardiovascular model consists of 16 compartments, including three cerebrospinal fluid (CSF) compartments, three cranial blood compartments, and 10 thoracic and lower limb blood compartments. To assess the models capability to address variations in physiological parameters, we completed a formal uncertainty and sensitivity analysis that evaluated the relative importance of 42 input parameters required in the model on relative compartment flows and compartment pressures. Utilizing the model in a pulsatile flow configuration, the sensitivity analysis identified the ten parameters that most influenced each compartment pressure. Generally, each compartment responded appropriately to parameter variations associated with itself and adjacent compartments. However, several unexpected interactions between components, such as between the choroid plexus and the lower capillaries, were found, and are due to simplifications in the formulation of the model. The analysis illustrates that highly influential parameters and those that have unique influences within the model formulation must be tightly controlled for successful model application.

gravitational physiology↗