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 1,369 records · Page 76

P1 Nonconforming Finite Element Method for the Solution of Radiation Transport Problems

The simulation of radiation transport in the optically thick flux-limited diffusion regime has been identified as one of the most time-consuming tasks within large simulation codes. Due to multimaterial complex geometry, the radiation transport system must often be solved on unstructured grids. In this paper, we investigate the behavior and the benefits of the unstructured P(sub 1) nonconforming finite element method, which has proven to be flexible and effective on related transport problems, in solving unsteady implicit nonlinear radiation diffusion problems using Newton and Picard linearization methods. Key words. nonconforrning finite elements, radiation transport, inexact Newton linearization, multigrid preconditioning

Kang, Kab S.↗

An internal variable constitutive model for the large deformation of metals at high temperatures

The advent of large deformation finite element methodologies is beginning to permit the numerical simulation of hot working processes whose design until recently has been based on prior industrial experience. Proper application of such finite element techniques requires realistic constitutive equations which more accurately model material behavior during hot working. A simple constitutive model for hot working is the single scalar internal variable model for isotropic thermal elastoplasticity proposed by Anand. The model is recalled and the specific scalar functions, for the equivalent plastic strain rate and the evolution equation for the internal variable, presented are slight modifications of those proposed by Anand. The modified functions are better able to represent high temperature material behavior. The monotonic constant true strain rate and strain rate jump compression experiments on a 2 percent silicon iron is briefly described. The model is implemented in the general purpose finite element program ABAQUS.

Brown, Stuart↗

A stochastic approach to uncertainty in the equations of MHD kinematics

The magnetohydrodynamic (MHD) kinematics model describes the electromagnetic behavior of an electrically conducting fluid when its hydrodynamic properties are assumed to be known. In particular, the MHD kinematics equations can be used to simulate the magnetic field induced by a given velocity field. While prescribing the velocity field leads to a simpler model than the fully coupled MHD system, this may introduce some epistemic uncertainty into the model. If the velocity of a physical system is not known with certainty, the magnetic field obtained from the model may not be reflective of the magnetic field seen in experiments. Additionally, uncertainty in physical parameters such as the magnetic resistivity may affect the reliability of predictions obtained from this model. By modeling the velocity and the resistivity as random variables in the MHD kinematics model, we seek to quantify the effects of uncertainty in these fields on the induced magnetic field. We develop stochastic expressions for these quantities and investigate their impact within a finite element discretization of the kinematics equations. We obtain mean and variance data through Monte Carlo simulation for several test problems. Toward this end, we develop and test an efficient block preconditioner for the linear systems arising from the discretized equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Finite Element Analysis and Machine Learning Guided Design of Carbon Fiber Organosheet-Based Battery Enclosures for Crashworthiness

Carbon fiber composite can be a potential candidate for replacing metal-based battery enclosures of current electric vehicles (E.V.s) owing to its better strength-to-weight ratio and corrosion resistance. However, the strength of carbon fiber-based structures depends on several parameters that should be carefully chosen. Here, in this work, we implemented high throughput finite element analysis (FEA) based thermoforming simulation to virtually manufacture the battery enclosure using different design and processing parameters. Subsequently, we performed virtual crash simulations to mimic a side pole crash to evaluate the crashworthiness of the battery enclosures. This high throughput crash simulation dataset was utilized to build predictive models to understand the crashworthiness of an unknown set. Our machine learning (ML) models showed excellent performance (R 2 > 0.97) in predicting the crashworthiness metrics, i.e., crush load efficiency, absorbed energy, intrusion, and maximum deceleration during a crash. We believe that this FEA-ML work framework will be helpful in down select process parameters for carbon fiber-based component design and can be transferrable to other manufacturing technologies.

36 MATERIALS SCIENCE↗

Multi-Terrain Impact Testing and Simulation of a Composite Energy Absorbing Fuselage Section

Comparisons of the impact performance of a 5-ft diameter crashworthy composite fuselage section were investigated for hard surface, soft soil, and water impacts. The fuselage concept, which was originally designed for impacts onto a hard surface only, consisted of a stiff upper cabin, load bearing floor, and an energy absorbing subfloor. Vertical drop tests were performed at 25-ft/s onto concrete, soft-soil, and water at NASA Langley Research Center. Comparisons of the peak acceleration values, pulse durations, and onset rates were evaluated for each test at specific locations on the fuselage. In addition to comparisons of the experimental results, dynamic finite element models were developed to simulate each impact condition. Once validated, these models can be used to evaluate the dynamic behavior of subfloor components for improved crash protection for hard surface, soft soil, and water impacts.

Fasanella, Edwin L.↗

Multi-Terrain Impact Testing and Simulation of a Composite Energy Absorbing Fuselage Section

Comparisons of the impact performance of a 5-ft diameter crashworthy composite fuselage section were investigated for hard surface, soft soil, and water impacts. The fuselage concept, which was originally designed for impacts onto a hard surface only, consisted of a stiff upper cabin, load bearing floor, and an energy absorbing subfloor. Vertical drop tests were performed at 25-ft/s onto concrete, soft-soil, and water at NASA Langley Research Center. Comparisons of the peak acceleration values, pulse durations, and onset rates were evaluated for each test at specific locations on the fuselage. In addition to comparisons of the experimental results, dynamic finite element models were developed to simulate each impact condition. Once validated, these models can be used to evaluate the dynamic behavior of subfloor components for improved crash protection for hard surface, soft soil, and water impacts.

Fasanella, Edwin L.↗

Applicability of Loads Estimation Techniques Using Sparse Acceleration Sensor Data to Spacecraft Structural Health Monitoring

The use of structural health monitoring systems on spacecraft structures can play a crucial role in ensuring the safety, reliability, and longevity of the structure by gathering and analyzing onboard sensor data. Of specific importance is monitoring for excessive loading at critical interfaces as any off-nominal structural excitations experienced by spacecraft structures can cause early unpredicted high structural life consumption or damage. The availability and cost of flight-certified sensors along with the size of spacecraft structures and allowable payload mass drives the need for a method to estimate loads using sparsely-located sensors. Numerous approaches such as physics-based, statistical learning, and physics-enhanced statistical learning algorithms have gained popularity among structural prognostics applications. However, developing noise-robust prediction models to assess loads and structural life predictions from a sparse multi-sensor data acquisition system can be a challenging task. This paper discusses the evaluation of physics-based versus machine-learning algorithms for predicting loads and structural life at mission critical locations on the spacecraft structure using a finite element loads analysis with the application of simulated noise and noise reduction techniques. To estimate the loads from accelerations, the physics-based algorithm leverages a loads transformation matrix from a Craig-Bampton reduced finite element model. A System Equivalent Reduction Expansion Process (SEREP) and a pseudo-inverse approach are considered to expand from the onboard sensor degrees of freedom to the Craig-Bampton model degrees of freedom. The machine learning algorithm provides a data driven solution/mapping of the sensor accelerations to the loads at the mission critical locations using a high dimensionality analysis. Although these strategies produce comparable loads prediction without noise, the limitations of these strategies with incorporating simulated noise and noise reduction techniques with low signal to noise ratio signals are evaluated. The study demonstrates the immense potential of statistical learning algorithms for sparse structural prognostic models and enhancing signal denoising techniques. These findings also highlight the need for noise-resilient prognostic models and low-noise data acquisition systems onboard spacecraft structures.

Spacecraft Structural Health Monitoring↗

Influences of Quench Cooling Rate on Microstructure and Corrosion Resistance of Al-Cu-Mg Alloy Based on the End-Quenching Test

To investigate the effects of the quench cooling rate on corrosion resistance of Al-Cu-Mg alloy, an end-quenching test was conducted and the microstructures at different cooling rates were observed by SEM and TEM. Additionally, the corrosion resistance was characterized by an intergranular corrosion test and electrochemical test. Moreover, the finite element method was applied to simulate the end quenching process. The results indicate that the actual end quenching process can be approximated as one-dimensional heat transfer, and the cooling rate varies at different cooling distances. By affecting the microstructures, decreasing the cooling rate leads to a decline in the corrosion properties. Low cooling rates coarsen the constituent particles and grain boundary particles, resulting in a wide precipitation-free zone and an increase in the intensity of corrosion reactions. A high cooling rate concentrates on the intragranular precipitant, which can reduce the pitting depth and represents a conversion from localized corrosion to general corrosion.

36 MATERIALS SCIENCE↗

Compact spherical neutron polarimeter using high-T{sub c} YBCO films

We describe a simple, compact device for spherical neutron polarimetry measurements at small neutron scattering angles. The device consists of a sample chamber with very low (<0.01 G) magnetic field flanked by regions within which the neutron polarization can be manipulated in a controlled manner. This allows any selected initial and final polarization direction of the neutrons to be obtained. We have constructed a prototype device using high-T{sub c} superconducting films and mu-metal to isolate regions with different magnetic fields and tested device performance in transmission geometry. Finite-element methods were used to simulate the device’s field profile and these have been verified by experiment using a small solenoid as a test sample. Measurements are reported using both monochromatic and polychromatic neutron sources. The results show that the device is capable of extracting sample information and distinguishing small angular variations of the sample magnetic field. As a more realistic test, we present results on the characterization of a 10 μm thick Permalloy film in zero magnetic field, as well as its response to an external magnetic field.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Analysis of new composite architectures

Global/local method in conjunction with a special macro finite element is used to gain computational efficiency in the simulation of textile composites behavior. Results are included to demonstrate the effectiveness of the method. Also, 2D finite elements are used to study boundary effects in plain weave composite specimens subjected to extension, shear, and flexure loads. Effective extension, shear, and flexural moduli were found to be quite sensitive to specimen size. For extension and flexure loads stress distributions were affected by a free surface, but the free surface boundary effect did not appear to propagate very far into the interior. For shear load, the boundary effect appeared to propagate much further into the interior. The report is in the form of two technical articles, the first describes the global/local method while the second describes the use of the 2D finite elements.

Whitcomb, John D.↗

Incorporating physically-based microstructures in materials modeling: Bridging phase field and crystal plasticity frameworks

Here, the mechanical properties of materials systems are highly influenced by various features at the microstructural level. The ability to capture these heterogeneities and incorporate them into continuum-scale frameworks of the deformation behavior is considered a key step in the development of complex non-local models of failure. In this study, we present a modeling framework that incorporates physically-based realizations of polycrystalline aggregates from a phase field (PF) model into a crystal plasticity finite element (CP-FE) framework. Simulated annealing via the PF model yields ensembles of materials microstructures with various grain sizes and shapes. With the aid of a novel FE meshing technique, FE discretizations of these microstructures are generated, where several key features, such as conformity to interfaces, and triple junction angles, are preserved. The discretizations are then used in the CP-FE framework to simulate the mechanical response of polycrystalline α-iron. It is shown that the conformal discretization across interfaces reduces artificial stress localization commonly observed in non-conformal FE discretizations. The work presented herein is a first step towards incorporating physically-based microstructures in lieu of the overly simplified representations that are commonly used. In broader terms, the proposed framework provides future avenues to explore bridging models of materials processes, e.g. additive manufacturing and microstructure evolution of multi-phase multi-component systems, into continuum-scale frameworks of the mechanical properties.

36 MATERIALS SCIENCE↗

A Statistical Approach for the Concurrent Coupling of Molecular Dynamics and Finite Element Methods

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, increasing the size of the MD domain quickly presents intractable computational demands. A robust approach to surmount this computational limitation has been to unite continuum modeling procedures such as the finite element method (FEM) with MD analyses thereby reducing the region of atomic scale refinement. The challenging problem is to seamlessly connect the two inherently different simulation techniques at their interface. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the typical boundary value problem used to define a coupled domain. The method uses statistical averaging of the atomistic MD domain to provide displacement interface boundary conditions to the surrounding continuum FEM region, which, in return, generates interface reaction forces applied as piecewise constant traction boundary conditions to the MD domain. The two systems are computationally disconnected and communicate only through a continuous update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM) as opposed to a direct coupling method where interface atoms and FEM nodes are individually related. The methodology is inherently applicable to three-dimensional domains, avoids discretization of the continuum model down to atomic scales, and permits arbitrary temperatures to be applied.

Saether, E.↗

Modeling of contact theories for the manipulation of biological micro/nanoparticles in the form of circular crowned rollers based on the atomic force microscope

This article has dealt with the development and modeling of various contact theories for biological nanoparticles shaped as cylinders and circular crowned rollers for application in the manipulation of different biological micro/nanoparticles based on Atomic Force Microscope. First, the effective contact forces were simulated, and their impact on contact mechanics simulation was investigated. In the next step, the Hertz contact model was simulated and compared for gold and DNA nanoparticles with the three types of spherical, cylindrical, and circular crowned roller type contact geometries. Then by reducing the length of the cylindrical section in the circular crowned roller geometry, the geometry of the body was made to approach that of a sphere, and the results were compared for DNA nanoparticles. To anticipatory validate the developed theories, the results of the cylindrical and the circular crowned roller contacts were compared with the results of the existing spherical contact simulations. Following the development of these contact models for the manipulation of various biological micro/nanoparticles, the cylindrical and the circular crowned roller type contact theories were modeled based on the theories of Lundberg, Dowson, Nikpur, Heoprich, and Hertz for the manipulation of biological micro/nanoparticles. Then, for a more accurate validation, the results obtained from the simulations were compared with those obtained by the finite element method and with the experimental results available in previous articles. The previous research works on the simulation of nanomanipulation have mainly investigated the contact theories used in the manipulation of spherical micro/nanoparticles. However since in real biomanipulation situations, biological micro/nanoparticles of more complex shapes need to be displaced in biological environments, this article therefore has modeled and compared, for the first time, different contact theories for use in the biomanipulation of cylindrical and circular crowned roller shaped micro/nanoparticles. The results of models indicate that the contact model of Hertz achieves the largest amount of deformation for the DNA nanoparticle in cylindrical form and the contact model of Heoprich achieves the largest deformation for the circular crowned roller shaped DNA. Of course, this finding is not always true for the other nanoparticles; and considering the mechanical and environmental characteristics, different results can be obtained. Also, by comparing the deformations of different types of nanoparticles, it was determined that the platelet type nanoparticles display the highest degree of deformation in all the considered models, due to their particular mechanical characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development of a Thin Film Primary Surface Heat Exchanger for Advanced Power Cycles

This project objective is to develop a high-temperature design upgrade for an existing primary surface heat exchanger so that the redesigned hardware is capable of operation in CO 2 at temperatures up to 1,510°F (821°C) and pressure differentials up to 130 psi (9 bar). The heat exchanger is proposed for use as a recuperator in an advanced low-pressure oxy-fuel Brayton cycle that is predicted to achieve over 50% thermodynamic efficiency, although the heat exchanger could also be used in other high-temperature, low-differential pressure cycles. This report describes the progress to date, which includes continuing work performed to select and test new candidate materials for the recuperator redesign, final mechanical and thermal performance analysis results of various redesign concepts, and the preliminary design of a test loop for the redesigned recuperator including a budgetary estimate for detailed test loop design, procurement, and test operation. A materials search was performed in order to investigate high-temperature properties of many candidate materials, including high-temperature strength and nickel content. These properties were used to rank the candidate materials, resulting in a reduced list of nine materials for corrosion testing. Multiple test rigs were considered and analyzed for short-term corrosion testing and Thermal Gravimetric Analysis (TGA) was selected as the most cost-effective option for evaluating corrosion resistance of the candidate materials. In addition, tantalum, niobium, and chromium coatings were identified as potential options for increased corrosion resistance. The test results show that many materials exhibit relatively low weight gain rates, and that niobium and tantalum coatings may improve corrosion resistance for many materials, while chromium coatings appear to oxidize and debond quickly. Metallurgical analysis of alloys was also performed, showing evidence of intergranular attack in 282 that may cause long-term reliability problems in CO 2 service at these temperatures. However, long-term testing in a flowing environment is recommended in order to understand accurately the severity of the attack. Detailed economic modeling of the existing air cycle recuperator and CO 2 cycle recuperator options was also completed, including costs for material, fabrication, fuel, maintenance, and operation. The analysis results show that the increased capital cost for high-temperature materials may be offset by higher cycle efficiencies, decreasing the overall lifetime cost of the system. The economic analysis also examines costs associated with increased pressure drop and material changes for two redesign options. These results show that, even with slightly reduced performance and/or higher material costs, the lifetime cost per energy production may still be reduced by over 12%. The existing recuperator design information was provided by Solar Turbines, Inc. via several models, drawings, and design handoff meetings. Multiple fluid/thermal and structural models were created in order to analyze critical recuperator performance and mechanical strength in critical areas throughout the redesign process. These models were analyzed for a baseline condition (consistent with current Mercury 50 operation) for validation purposes. Results are presented for heat transfer coefficients and pressure drops, matching well with the existing operational data. Simulation of higher-temperature CO 2 conditions was also performed, showing a slight expected increase in both heat transfer and pressure drop. Mechanical analysis results for critical areas on the cross-flow and counter-flow sheets have also been obtained for air and CO 2 cases. These results show similar stresses in both cases but significantly reduced safety factors for the CO 2 case due to reduced yield and creep rupture strengths of alloy 625 at the higher temperatures. A concept brainstorm session and initial down-selection were completed in order to identify promising redesign options for further analysis. Detailed analysis of all promising redesign options was performed via finite element and computational fluid dynamic simulations in order to characterize mechanical and thermal-fluid performance of each option. These options included material change, various sheet thickness configurations, pitch and phasing of cross-flow and counter-flow sheets, and separator sheets. The analysis results have identified two viable redesign options that maintain existing safety margins optimally through a material change to Haynes 282 and (A) sheet thickness increases of 40% on the counter-flow sheet and 75% on the hot side cross-flow corrugation sheet or (B) addition of a separator sheet in the counter-flow section while maintaining the original counter-flow sheet thickness and increasing the cross-flow corrugation sheet thickness by 90% to account for the increase in cell height. While both options satisfy mechanical stress constraints, the separator sheet design has a higher part count, slightly reduced heat transfer, and slightly higher pressure drop than the first option and is not preferred. Finally, several test loop concepts have been developed for different full-scale and reduced-scale recuperator testing options. For each option, various loop components, such as heat exchangers, valves, heaters, and compressors, were evaluated in an effort to maximize utilization of existing resources. All concepts utilize an existing 3-MW CO 2 compressor, heater, and loop coolers, but the concepts vary by incorporating different amounts of new equipment for achieving various flow rates (all concepts operate at design pressure and temperature). The third concept achieves a 1 kg/s test without purchasing any costly equipment (coolers, heaters, blowers, etc.). Since the stacked cell design of the recuperator results in the same flow conditions at each core cell (even for a reduced-scale test). Thus, test loop Concept #3 was selected for the preliminary design. This loop design is detailed within the report, culminating in a budgetary estimate of $1,013,000.00 for the detailed design, construction, commissioning, and operation of a high-temperature recuperator test loop.

03 NATURAL GAS↗

Simulations of SITI Cookoff Experiments Carried Out with Different Lots of PBX 9502

A pressure dependent cookoff model for PBX 9502 was developed by Hobbs’ et. al. PBX 9502 is composed of 95% by mass triaminotrinitrobenzene (TATB) and a 5% by mass chlorotrifluoroethylene/vinylidine fluoride binder. The objective in this study is to implement this cookoff model in Aria to simulate Sandia Instrumented Thermal Ignition (SITI) experiments that were carried out with different manufacturing lots of PBX 9502. The SITI design consists of solid cylinders (1" diameter × 1" height) of insensitive high explosive (IHE) confined by a cylindrical aluminum case. An electric heater is wrapped around the lateral surface of the case. This heater produces a temperature heating ramp on the outer surface of the case. Internal thermocouples measure the IHE temperature rise from the center to locations close to the IHE-aluminum interface. The energetic material is heated until thermal ignition occurs. Pressure is measured with a static pressure transducer installed on top of the confinement case. Two–dimensional axisymmetric heat conduction finite element models were implemented to simulate these experiments using four options of the PBX 9502 cookoff model. In addition, the predictive ability of this thermal decomposition model is evaluated using Latin Hypercube Sampling (LHS) techniques.

42 ENGINEERING↗

Solution of a few nonlinear problems in aerodynamics by the finite elements and functional least squares methods

The numerical simulation of the transonic flows of idealized fluids and of incompressible viscous fluids, by the nonlinear least squares methods is presented. The nonlinear equations, the boundary conditions, and the various constraints controlling the two types of flow are described. The standard iterative methods for solving a quasi elliptical nonlinear equation with partial derivatives are reviewed with emphasis placed on two examples: the fixed point method applied to the Gelder functional in the case of compressible subsonic flows and the Newton method used in the technique of decomposition of the lifting potential. The new abstract least squares method is discussed. It consists of substituting the nonlinear equation by a problem of minimization in a H to the minus 1 type Sobolev functional space.

Periaux, J.↗

Simulation of thin slot spirals and dual circular patch antennas using the finite element method with mixed elements

This semi-annual report describes progress up to mid-January 1995. The report contains five sections all dealing with the modeling of spiral and patch antennas recessed in metallic platforms. Of significance is the development of decomposition schemes which separate the different regions of the antenna volume. Substantial effort was devoted to improving the feed model in the context of the finite element method (FEM). Finally, an innovative scheme for truncating finite element meshes is presented.

Gong, Jian↗