Search NASA⌕ Search

SEARCH · Search NASA

Results for “FEM 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 19 records

Force Limiting Vibration Tests Evaluated from both Ground Acoustic Tests and FEM Simulations of a Flight Like Vehicle System Assembly

Marshall Space Flight Center has conducted a series of ground acoustic tests with the dual goals of informing analytical judgment, and validating analytical methods when estimating vibroacoustic responses of launch vehicle subsystems. The process of repeatedly correlating finite element-simulated responses with test-measured responses has assisted in the development of best practices for modeling and post-processing. In recent work, force transducers were integrated to measure interface forces at the base of avionics box equipment. Other force data was indirectly measured using strain gauges. The combination of these direct and indirect force measurements has been used to support and illustrate the advantages of implementing the Force Limiting approach for equipment qualification tests. The comparison of force response from integrated system level tests to measurements at the same locations during component level vibration tests provides an excellent illustration. A second comparison of the measured response cases from the system level acoustic tests to finite element simulations has also produced some principles for assessing the suitability of Finite Element Models (FEMs) for making vibroacoustics estimates. The results indicate that when FEM models are employed to guide force limiting choices, they should include sufficient detail to represent the apparent mass of the system in the frequency range of interest.

Smith, Andrew↗

An efficient and accurate model of the coax cable feeding structure for FEM simulations

An efficient and accurate coax cable feed model is proposed for microstrip or cavity-backed patch antennas in the context of a hybrid finite element method (FEM). A TEM mode at the cavity-cable junction is assumed for the FEM truncation and system excitation. Of importance in this implementation is that the cavity unknowns are related to the model fields by enforcing an equipotential condition rather than field continuity. This scheme proved quite accurate and may be applied to other decomposed systems as a connectivity constraint. Comparisons of our predictions with input impedance measurements are presented and demonstrate the substantially improved accuracy of the proposed model.

Gong, Jian↗

Arbitrarily shaped dual-stacked patch antennas: A hybrid FEM simulation

A dual-stacked patch antenna is analyzed using a hybrid finite element - boundary integral (FE-BI) method. The metallic patches of the antenna are modeled as perfectly electric conducting (PEC) plates stacked on top of two different dielectric layers. The antenna patches may be of any shape and the lower patch is fed by a coaxial cable from underneath the ground plane or by an aperture coupled microstrip line. The ability of the hybrid FEM technique for the stacked patch antenna characterization will be stressed, and the EM coupling mechanism is also discussed with the aid of the computed near field patterns around the patches.

Gong, Jian↗

Occupant Response Analysis of a Full-Scale Crash Test of a Fokker F28 Fellowship Aircraft

In this report, the capability of component fuselage testing and finite element models (FEMs) to predict occupant injury risk when applied to full-scale aircraft crash testing was evaluated. Component level and full-scale crash tests of a Fokker-F28 aircraft were performed in conjunction with representative FEM simulations. Lumbar spine injury risk, calculated through Anthropomorphic Test Device (ATD) outputs measured in the full-scale test, was compared with predictions made by component level testing and simulations. A quantitative assessment of FEM prediction accuracy was made using the International Organization for Standardization (ISO) Technical Report (TR) 16250 curve comparison methodology. Lumbar load injury risk measured in the full-scale test configuration was found to be significantly higher than in the component tests performed. The FEM provided a closer prediction of occupant injury risk than component testing but exhibited limitations in the multi-axis load environment. This work will help inform the application of component testing and FEM simulation when used to evaluate aircraft crashworthiness.

Jacob B Putnam↗

Transformer Leakage Inductance Design Methodology

The leakage inductance exhibited by a transformer depends on its winding geometry, which generally involves the selection of several key design parameters in addition to the winding structure and the interleaving configuration. With few resources explaining the effects of these design choices on the observed leakage inductance, numerous trial-and-error iterations become necessary to realize the desired leakage inductance. This paper explores more than a hundred winding geometries feasible in a 2-winding transformer comprising the same magnetic core, number of turns, and wire gauge, and finds the leakage inductance for each unique design using 2-D finite element method (FEM) simulations in association with the semi-analytical double- 2-D model. These leakage inductances are plotted and further analyzed to understand the effects of different design parameters on the effective leakage inductance. The results presented herein and the conclusions drawn from this research can serve as a valuable resource for future design practitioners from both industry and academia.

24 POWER TRANSMISSION AND DISTRIBUTION↗

New Hybrid Model for Evaluating the Frequency-Dependent Leakage Inductance of a Variable Inductance Transformer (VIT)

Skin and proximity effects can cause a significant drop in the effective leakage inductance of a transformer when the operating frequency is increased. Although the magnetic image method-based double-2-D model can calculate the lowfrequency leakage inductance with sufficient accuracy, it is inherently a frequency-independent model. While Dowell’s 1-D model uses frequency-dependent relations to account for both skin and proximity effects, its accuracy is severely affected by the assumed winding geometry. In this paper, a hybrid model is proposed that uses superposition to combine a modified Dowell’s model with the double-2-D model. The proposed model is investigated on a variable inductance transformer (VIT)—a partially-filled transformer whose leakage inductance can be varied by moving one of the windings mechanically. The frequency-dependent leakage inductances of the VIT evaluated using the hybrid model are in excellent agreement with the corresponding finite element method (FEM) simulated and experimentally measured values, thereby validating the proposed hybrid model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗

Measurement of amplification and absorption of a THz quantum-cascade metasurface free-space amplifier

An active amplifying metasurface based on a quantum-cascade gain material at 2.7 THz is studied. The metasurface is first evaluated as the active component of an external cavity laser with excellent beam quality and frequency tunability from 2.55–2.8 THz. Amplification and absorption of the metasurface alone are then separately measured at a single frequency using a probe signal from a CO 2 -pumped gas laser operating at 2.743 THz. The metasurface reflectance vs bias is measured and compared with expectations from non-equilibrium Green’s function simulations of the quantum-cascade gain material and FEM simulations of the metasurface reflectance. A peak amplification on the order of 0.3 dB is measured. Design strategies are discussed for increasing single-reflection amplification (upward of 10 dB) and reducing power dissipation. Further increased amplification could be achieved by cascading multiple metasurfaces.

Curwen, Christopher A.↗

Stability Analysis of Intertank Formed Skin/Stringer Compression Panel with Simulated Damage

The External Tank (ET) is a component of the Space Shuttle launch vehicle that contains fuel and oxidizer. During launch, the ET supplies the space shuttle main engines with liquid hydrogen and liquid oxygen. In addition to supplying fuel and oxidizer, it is the backbone structural component of the Space Shuttle. It is comprised of a liquid hydrogen (LH2) tank and a liquid oxygen (LOX) tank, which are separated by an Intertank. The Intertank is a stringer-stiffened cylindrical structure with hat-section stringers that are roll formed from aluminum-lithium alloy Al-2090. Cracks in the Intertank stringers of the STS-133 ET were noticed after a November 5, 2010 launch attempt. The cracks were approximately nine inches long and occurred on the forward end of the Intertank (near the LOX tank), along the fastener line, and were believed to have occurred while loading the ET with the cryogenic propellants. These cracks generated questions about the structural integrity of the Intertank. In order to determine the structural capability of the Intertank with varying degrees of damage, a finite element model (FEM) simulating a 1995 compression panel test was analyzed and correlated to test data. Varying degrees of damage were simulated in the FEM, and non-linear stability analyses were performed. The high degree of similarity between the compression panel and the Intertank provided confidence that the ET Intertank would have similar capabilities.

Harper, David W.↗

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High Throughput Computational Framework of Materials Properties for Extreme Environments

This project aims to establish a framework capable of efficiently predicting the properties of structural materials for service in harsh environments over a wide range of temperatures and over long periods of time. The approach is to develop and integrate high throughput first-principles calculations in combination with machine learning (ML) methods, perform high throughput CALPHAD (calculations of phase diagrams) modeling, and carry out finite element method (FEM) simulations. Relevant to high temperature service in fossil power system, nickel-based superalloys such as Inconel 740 and Haynes 282 as well as the associated (Ni-Cr-Co)-Al-C-Fe-Mn-Mo-Nb-Si-Ti system, were investigated. The present framework was built on the concept of phase-based property data, in which properties of individual phases are modeled as a function of internal and external independent variables. This project established an open-source infrastructure with the following capabilities: (1) High throughput implementation of first-principles calculations at finite temperatures and variable compositions using both accurate phonon calculations and the efficient Debye model for thermodynamic properties, elastic constants, diffusion coefficients, vacancy formation, stacking and twin faults, and dislocation mobility; i.e., using the developed code DFTTK; (2) Machine learning capabilities to predict the above properties so that the number of first-principles calculations can be significantly reduced; e.g., using the developed code SIPFENN; (3) High throughput CALPHAD modeling of the above properties as a function of temperature and composition using our unique capability based on ESPEI and PyCalphad; (4) New capabilities to predict the stress-strain behavior of individual phases; and (5) New models for tensile strength prediction in common FEM software with the crystal plasticity finite element simulations (CPFEM).

, Ni-based superalloys↗

Adaptive Data-Driven Deep-Learning Surrogate Model for Frontal Polymerization in Dicyclopentadiene

Frontal polymerization (FP) is a self-sustaining curing process that enables rapid and energy-efficient manufacturing of thermoset polymers and composites. Computational methods conventionally used to simulate the FP process are time-consuming, and repeating simulations are required for sensitivity analysis, uncertainty quantification, or optimization of the manufacturing process. Here, in this work, we develop an adaptive surrogate deep-learning model for FP of dicyclopentadiene (DCPD), which predicts the evolution of temperature and degree of cure orders of magnitude faster than the finite-element method (FEM). The adaptive algorithm provides a strategy to select training samples efficiently and save computational costs by reducing the redundancy of FEM-based training samples. The adaptive algorithm calculates the residual error of the FP governing equations using automatic differentiation of the deep neural network. A probability density function expressed in terms of the residual error is used to select training samples from the Sobol sequence space. The temperature and degree of cure evolution of each training sample are obtained by a 2D FEM simulation. The adaptive method is more efficient and has a better prediction accuracy than the random sampling method. With the well-trained surrogate neural network, the FP characteristics (front speed, shape, and temperature) can be extracted quickly from the predicted temperature and degree-of-cure fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel processing route for design and manufacturing of metal toughened nanoceramics: Al–Al 2 O 3 nanocermets

Monolithic alumina is a highly brittle ceramic with a very low fracture toughness, making it an unreliable material for mechanical and structural applications. Ceramic-metal composites (cermets) are the subject of intensive interest to obtain materials which can address the intrinsic brittleness of ceramics without compromising their high strength and light weight. In the present investigation, aluminium-alumina (Al–Al 2 O 3 ) nanocermets were developed. The nanocermets possessed the toughness of Al phase and the strength of reinforced nano Al 2 O 3 particles. Here, the nanocermets green bodies were prepared such that upon sintering each adjacent Al 2 O 3 grain was conformally coated with a continuous Al film. Mechanical milling via a combination of cryo-milling and ball-milling was used to prepare Al–Al 2 O 3 nano-powder mixture. Spark plasma sintering was used to consolidate the green nanocermets. Results show that specimens prepared from 35 h milled powders exhibited superior mechanical properties as compared to other nanocermets, with an improvement in microhardness, elastic modulus, compressive strength, and flexural strength by 88%, 78%, 69% and 40% with respect to sintered pure Al. The indentation toughness of 35 h milled nanocermet showed an increment of 12.5% with respect to sintered pure Al 2 O 3 . Characterization techniques such as XRD, SEM, HR-TEM, X-ray tomography were employed to investigate the phase evolution and structural morphology, while FEM simulations were performed for understanding the crack propagation mechanisms in nanocermets.

36 MATERIALS SCIENCE↗

Evaluation of adsorption and mechanical strength of 13X zeolite mixtures with phyllosilicate binders using molecular dynamics simulation and positron annihilation spectroscopy

There is growing interest in developing zeolites with novel internal structures that have optimal adsorptive capacity and high mechanical strength, while offering advantages, such as being light weight. We integrate computational and experimental methods to explore the effect of binder/zeolite types, and weight percentages on the mechanical strength of 13X zeolite and adsorption capacities of N 2 , H 2 O, and CO 2 for additive manufacturing (AM) applications with the goal of maximizing both adsorption and strength. Zeolite 13X mixtures and phyllosilicate binders (either bentonite or kaolin) are combined using molecular dynamics (MD) simulations to create structures with various binder/zeolite weight percentages. Further, adsorption capabilities and mechanical strength are assessed using the grand canonical Monte Carlo (GCMC) and ReaxFF modules, respectively. Our modeling shows that an optimized zeolite/binder ratio for N 2 adsorption is around 15 wt% for kaolin and roughly 10 wt% for bentonite. The resulting parameters can be applied to facilitate macro-scale computational fluid dynamics (CFD) and finite element method (FEM) simulations of an AM zeolite structure. We also performed Positron Annihilation Lifetime Spectroscopy (PALS) measurements on zeolite samples to explore the effect of changes in the internal volume. The results show an inverse relationship between the free volume and the solid loading and adding a binder changes the morphology of the zeolite-binder compound and decreases open-volume area significantly.

36 MATERIALS SCIENCE↗

Generalizing synthetic data-trained acoustic predictive models to real-world measurements

Acoustic Resonance Spectroscopy (ARS) is highly sensitive to structural properties such as material, geometry, and environmental conditions; as a consequence, it can noninvasively measure internal properties that are unobservable by most other methods. Because of its sensing capabilities and low implementation cost and complexity, ARS has potential as a paradigm shift in noninvasive sensing, characterization, and monitoring applications. However, extracting specific properties from ARS measurements, comprising the vibration spectrum of a test object, is challenging due to the sensitivity of the spectra to other structural changes not being measured, e.g. manufacturing tolerances, component coupling, environmental variation, etc. Neural Networks are promising tools for identifying trends in ARS measurements, but their training typically requires large datasets, which are often impractical to obtain for real-world systems. Synthetic data can be simulated efficiently, but discrepancies between synthetic and real-world data frequently lead to poor generalization when testing on the real-world data. We propose a novel ARS model training framework that enables networks trained exclusively on synthetic ARS data to generalize effectively to real-world measurements. Our approach leverages the Correlation Alignment (CORAL) technique to enforce the extraction of features common to both synthetic and real-world domains. As a case study, we demonstrate noninvasive ARS-based pressure measurements in sealed systems. Finite element method (FEM) simulations were used to generate synthetic training data across diverse vessel configurations and pressure conditions, and model performance was then tested on real-world measurements. We demonstrate that robust machine learning models for ARS can be developed without large real-world datasets, significantly broadening the applicability of ARS for noninvasive sensing. Moreover, the approach is extensible to other sensing modalities where synthetic data are abundant but real-world data are limited.

36 MATERIALS SCIENCE↗

A numerical and experimental approach to compare the effect of sample thickness in small in-situ SEM and large ex-situ tensile testing in Alloy 709

ASTM standards for tensile tests define specific sample size requirements regardless of grain size of the material. However, sample size requirements for testing should be considered in conjunction with the number of grains within its reduced cross-sectional area. This is particularly important for in-situ Scanning Electron Microscope (SEM) tensile tests, as they have to be conducted on smaller samples. In this study, a comprehensive experimental and numerical evaluation of the effect of specimen thickness (and the resulting number of grains within the cross-section) were conducted using in-situ SEM tensile test (on submillimeter thick samples) and ex-situ tensile tests (on samples of 0.68–5.9 mm thickness) and the results are compared with FEM simulations outcome. All tests were conducted at room temperature and the results are correlated to the number of grains within the thickness of specimens. The tensile test results indicated that even though the 0.2% proof stress and the tensile strength do not vary for the current range of thicknesses (with number of grains in the cross-section varying from 13 to 118), a difference in necking mechanism exists. Post tensile strength, thinner samples undergo more shear failure and diagonal localized necking whereas thicker samples experience more diffuse necking indicated by decreased area of shear failure at the edges and increased area of dimpled tensile failure at the center of the sample. Further, FEM results complement the experimental findings by showing the formation of conjugated localized shear bands on the upper and lower surface in the 0.68 mm thick sample and shear bands combining to form singular shear bands in thicker samples. These results also confirm the validity of the in-situ SEM tensile tests conducted on thinner samples as long as the required minimum number of grains exist within the cross-section.

36 MATERIALS SCIENCE↗

Integrated Model of the Eye/Optic Nerve Head Biomechanical Environment

Visual Impairment and Intracranial Pressure (VIIP) syndrome is a concern for long-duration space flight. Previously, it has been suggested that ocular changes observed in VIIP syndrome are related to the cephalad fluid shift that results in altered fluid pressures [1]. We are investigating the impact of changes in intracranial pressure (ICP) using a combination of numerical models, which simulate the effects of various environment conditions, including finite element (FE) models of the posterior eye. The specific interest is to understand how altered pressures due to gravitational changes affect the biomechanical environment of tissues of the posterior eye and optic nerve sheath. METHODS: Additional description of the numerical modeling is provided in the IWS abstract by Nelson et al. In brief, to simulate the effects of a cephalad fluid shift on the cardiovascular and ocular systems, we utilized a lumped-parameter compartment model of these systems. The outputs of this lumped-parameter model then inform boundary conditions (pressures) for a finite element model of the optic nerve head (Figure 1). As an example, we show here a simulation of postural change from supine to 15 degree head-down tilt (HDT), with primary outcomes being the predicted change in strains at the optic nerve head (ONH) region, specifically in the lamina cribrosa (LC), retrolaminar optic nerve, and prelaminar neural tissue (PLNT). The strain field can be decomposed into three orthogonal components, denoted as the first, second and third principal strains. We compare the peak tensile (first principal) and compressive (third principal) strains, since elevated strain alters cell phenotype and induces tissue remodeling. RESULTS AND CONCLUSIONS: Our lumped-parameter model predicted an IOP increase of c. 7 mmHg after 21 minutes of 15 degree HDT, which agreed with previous reports of IOP in HDT [1]. The corresponding FEM simulations predicted a relative increase in the magnitudes of the peak tensile and compressive strains in the lamina cribrosa of 42 and 43, respectively (Fig. 2). The corresponding changes in the optic nerve strains were 17 and 39, while in the PLNT they were 47 and 43. These magnitudes of relative elevations in peak strains may induce a phenotypic response in resident mechano-responsive resident cells [2]. This approach may be expanded to investigate other environmental changes (e.g. parabolic flight). Through our VIIP SCHOLAR project, we will validate and improve these integrated models by measuring patient-specific changes in optic nerve sheath geometry in patients with idiopathic intracranial hypertension before and after lumbar puncture and CSF removal.

biomechanics↗