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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Mechanics of cantilever beam: Implementation and comparison of FEM and MLPG approach

Two weak form solution approaches for partial differential equations, the well known meshbased finite element method and the newer meshless local Petrov Galerkin method are described and compared on a standard test case - mechanics of cantilever beam. The implementation, solution accuracy and calculation complexity are addressed for both approaches. We found out that FEM is superior in most standard criteria, but MLPG has some advantages because of its flexibility that results from its general formulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2014

Numerical Modeling of Ocular Dysfunction in Space

Upon introduction to microgravity, the near-loss of hydrostatic pressure causes a marked cephalic (headward) shift of fluid in an astronaut's body. The fluid shift, along with other factors of spaceflight, induces a cascade of interdependent physiological responses which occur at varying time scales. Long-duration missions carry an increased risk for the development of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, a spectrum of ophthalmic changes including posterior globe flattening, choroidal folds, distension of the optic nerve sheath, kinking of the optic nerve and potentially permanent degradation of visual function. In the cases of VIIP found to date, the initial onset of symptoms occurred after several weeks to several months of spaceflight, by which time the gross bodily fluid distribution is well established. We are developing a suite of numerical models to simulate the effects of fluid shift on the cardiovascular, central nervous and ocular systems. These models calculate the modified mean volumes, flow rates and pressures that are characteristic of the altered quasi-homeostatic state in microgravity, including intracranial and intraocular pressures. The results of the lumped models provide initial and boundary data to a 3D finite element biomechanics simulation of the globe, optic nerve head and retrobulbar subarachnoid space. The integrated set of models will be used to investigate the evolution of the biomechanical stress state in the ocular tissues due to long-term exposure to microgravity.

intraocular pressure↗
DOE OSTI2014

WE-G-17A-08: Electron Gun Operation for in Line MRI-Linac Configurations: An Assessment of Beam Fidelity and Recovery Techniques for Different SIDs and Magnetic Field Strengths

Purpose: To test the functionality of medical electron guns within the fringe field of a purpose built superconducting MRI magnet, and to test different recovery techniques for a variety of imaging field strengths and SIDs. Methods: Three different electron guns were simulated using Finite Element Modelling; a standard diode gun, a standard triode gun, and a novel diode gun designed to operate within parallel magnetic fields. The approximate working regime of each gun was established by assessing exit current in constant magnetic fields of varying strength and defining ‘working’ as less than 10% change in injection current. Next, the 1.0T MRI magnet was simulated within Comsol Multiphysics. The coil currents in this model were also scaled to produce field strengths of .5, 1, 1.5 and 3T. Various magnetic shield configurations were simulated, varying the SID from 800 to 1300mm. The average magnetic field within the gun region was assessed together with the distortion in the imaging volume - greater than 150uT distortion was considered unacceptable. Results: The conventional guns functioned in fields of less than 7.5mT. Conversely, the redesigned diode required fields greater than .1T to function correctly. Magnetic shielding was feasible for SIDS of greater than 1000mm for field strengths of .5T and 1T, and 1100mm for 1.5 and 3.0T. Beyond these limits shielding resulted in unacceptable MRI distortion. In contrast, the redesigned diode could perform acceptably for SIDs of less than 812, 896, 931, and 974mm for imaging strengths of 0.5, 1.0, 1.5, 3.0T. Conclusions: For in-line MRIlinac configurations where the electron gun is operating in low field regions, shielding is a straight forward option. However, as magnetic field strength increases and the SID is reduced, shielding results in too great a distortion in the MRI and redesigning the electron optics is the preferable solution. The authors would like to acknowledge funding from the National Health and Research Council (AUS), National Institute of Health (NIH), and Cancer Institute NSW.

43 PARTICLE ACCELERATORS↗
NASA NTRS1987

Initial dynamic load estimates during configuration design

This analysis includes the structural response to shock and vibration and evaluates the maximum deflections and material stresses and the potential for the occurrence of elastic instability, fatigue and fracture. The required computations are often performed by means of finite element analysis (FEA) computer programs in which the structure is simulated by a finite element model which may contain thousands of elements. The formulation of a finite element model can be time consuming, and substantial additional modeling effort may be necessary if the structure requires significant changes after initial analysis. Rapid methods for obtaining rough estimates of the structural response to shock and vibration are presented for the purpose of providing guidance during the initial mechanical design configuration stage.

Schiff, Daniel↗
DOE OSTI2023

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗
DOE OSTI2025

Quantum State Transfer between Heterogeneous 3D QPUs

We present detailed analyses and microwave simulations of a quantum interconnect architecture that couples two distinct three-dimensional quantum processing units (QPUs): (i) a two-cell superconducting radio-frequency (SRF) cavity coupled to multiple transmon qubits and (ii) a waveguide-based cavity hosting a planar circuit with transmons and resonators. The modules are linked via a low-loss superconducting coaxial cable, which serves as a quantum bus. Finite-element electromagnetic simulations are performed to optimize coupling strengths and mode frequencies for high-fidelity quantum state transfer between cavities. These simulation results provide a comprehensive blueprint for the experimental realization of a flexible quantum interconnect architecture, coupling heterogeneous QPUs.

Yaker, Joey [Fermilab]↗
DOE OSTI2016

Laser-excited optical emission response of CdTe quantum dot/polymer nanocomposite under shock compression

Laser-driven shock compression experiments and corresponding finite element method simulations are carried out to investigate the blueshift in the optical emission spectra under continuous laser excitation of a dilute composite consisting of 0.15% CdTe quantum dots by weight embedded in polyvinyl alcohol polymer. This material is a potential candidate for use as internal stress sensors. The analyses focus on the time histories of the wavelength blue-shift for shock loading with pressures up to 7.3 GPa. The combined measurements and calculations allow a relation between the wavelength blueshift and pressure for the loading conditions to be extracted. It is found that the blueshift first increases with pressure to a maximum and subsequently decreases with pressure. This trend is different from the monotonic increase of blueshift with pressure observed under conditions of quasistatic hydrostatic compression. Additionally, the blueshift in the shock experiments is much smaller than that in hydrostatic experiments at the same pressure levels. The differences in responses are attributed to the different stress states achieved in the shock and hydrostatic experiments and the time dependence of the mechanical response of the polymer in the composite. The findings offer a potential guide for the design and development of materials for internal stress sensors for shock conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Thermal behavior characterization of a kilowatt-power-level cryogenically cooled Yb:YAG active mirror laser amplifier

We present an analysis of the thermal behavior of a high-energy kilowatt-average-power diode-pumped cryogenically cooled Yb:YAG active mirror laser amplifier based on measurements and simulations. Maps of the temperature distribution of the laser material at pump powers up to 1 kW were obtained for the first time by spatially and spectrally resolving the fluorescence induced by a scanning probe beam. The wavefront distortion resulting from the front surface deformation and the overall deformation of the gain medium assembly were measured using a Mach–Zehnder interferometer. The measured deformations agree well with the results of thermomechanical modeling using finite element method simulations, and with the results of focal length shift measurements. The relative contributions to the optical path difference (OPD) of the mechanical deformations, refractive index changes, and electronic contribution are discussed. We propose that the Cr 4+ :YAG cladding plays a significant role in both the temperature distribution and the overall OPD changes. Moreover, the pump-induced mechanical deformations of the assembly dominate the OPD changes in this kilowatt-average-pump-power cryogenically cooled Yb:YAG active mirror laser.

36 MATERIALS SCIENCE↗
NASA NTRS2000

Simulation of Aircraft Landing Gears with a Nonlinear Dynamic Finite Element Code

Recent advances in computational speed have made aircraft and spacecraft crash simulations using an explicit, nonlinear, transient-dynamic, finite element analysis code more feasible. This paper describes the development of a simple landing gear model, which accurately simulates the energy absorbed by the gear without adding substantial complexity to the model. For a crash model, the landing gear response is approximated with a spring where the force applied to the fuselage is computed in a user-written subroutine. Helicopter crash simulations using this approach are compared with previously acquired experimental data from a full-scale crash test of a composite helicopter.

Lyle, Karen H.↗
DOE OSTI2019

Strain localization and dynamic recrystallization in polycrystalline metals: Thermodynamic theory and simulation framework

We describe a theoretical and computational framework for adiabatic shear banding (ASB) and dynamic recrystallization (DRX) in polycrystalline materials. The Langer-Bouchbinder-Lookman (LBL) thermodynamic theory of polycrystalline plasticity, which we recently reformulated to describe DRX via the inclusion of the grain boundary density or the grain size as an internal state variable, provides a convenient and self-consistent way to represent the viscoplastic and thermal behavior of the material, with minimal ad-hoc assumptions regarding the initiation of yielding or onset of shear banding. Here, we implement the LBL-DRX theory in conjunction with a finite-element computational framework. Favorable comparison to experimental measurements on a top-hat AISI 316L stainless steel sample compressed with a split-Hopkinson pressure bar suggests the accuracy and usefulness of the LBL-DRX framework, and demonstrates the crucial role of DRX in strain localization.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Analysis of the static magnetic field-dependent optical transmission of Ni nanorod colloidal suspensions

The magnetic field-dependent optical transmission of dilute Ni nanorod aqueous suspensions was investigated. A series of four samples of nanorods were synthesized using the AAO template method and processed to stable colloids. The distributions of their length and diameter were characterized by analysis of TEM images and revealed average diameters of ∼25 nm and different lengths in the range of 60 nm–1100 nm. The collinear magnetic and optical anisotropy was studied by static field-dependent transmission measurements of linearly polarized light parallel and perpendicular to the magnetic field direction. The experimental results were modelled assuming the field-dependent orientation distribution function of a superparamagnetic ensemble for the uniaxial ferromagnetic nanorods in liquid dispersion and extinction cross sections for longitudinal and transversal optical polarization derived from different approaches, including the electrostatic approximation and the separation of variables method, both applied to spheroidal particles, as well as finite element method simulations of spheroids and capped cylindrical particles. The extinction cross sections were compared to reveal the differences associated with the approximations of homogeneous polarization and/or particle shape. The consequences of these approximations for the quantitative analysis of magnetic field-dependent optical transmission measurements were investigated and a reliable protocol derived. Furthermore, the changes in optical cross sections induced by electromagnetic interaction between two nanorods in parallel end-to-end and side-by-side configuration as a function of their separation were studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2008

Efficient Modal Basis Selection Criteria for Reduced-Order Nonlinear Simulation

A modal basis selection technique for a reduced-order nonlinear numerical simulation with application to two-dimensional structures is presented as a two-step procedure. A system identification analysis is first performed using proper orthogonal decomposition. Using these results, a set of load-invariant bases consisting of the normal modes is next selected. Two criteria for making the basis selection are offered; one using the modal assurance criterion and the other using the modal expansion theorem. The quality of the subsequent reduced-order analyses are examined through comparison with computationally intensive finite element nonlinear simulations in physical degrees-of-freedom. A clamped flat isotropic plate under a random acoustic loading is considered to demonstrate the procedure. It is found that the subject procedure enables formation of an accurate and computationally efficient reduced-order system applicable to a broad range of loading conditions.

Przekop, Adam↗
DOE OSTI2019

Influence of Single-Point Incremental Force Process Parameters on Forming Characteristics and Microstructure Evolution of AA-6061 Alloy Sheet

Single-point incremental forming (SPIF) process comprises a set of process variables, such as forming tool diameter, vertical pitch, spindle rotation speed, and tool velocity, which may affect the forming behavior of the sheet to be deformed. The objective of this work is to study the effect of SPIF process parameters on forming characteristics and microstructure development for the AA-6061 (T6) aluminum alloy sheet. The SPIF experiments and finite element (FE) simulations were performed at different process parameters to achieve conical shapes from the AA-6061 blanks. The effect of process parameters on forming forces, thickness uniformity in wall region of the cone and surface roughness of the blank was analyzed. A detailed microstructure study was performed to analyze the effect of process variables on microstructure and texture evolution during the SPIF process. This study reveals that the process parameters are likely to influence the texture development especially at high tool diameter and vertical pitch values. Therefore, suitability and consequences of using different combinations of tool diameter and vertical pitch values are discussed for AA-6061 alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Modelling grain-scattered ultrasound in austenitic stainless-steel welds: A hybrid model

The ultrasonic inspection of austenitic stainless steel welds can be challenging due to their coarse grain structure, charaterised by preferentially oriented, elongated grains. The anisotropy of the weld is manifested as both a ‘steering’ of the beam and the back-scatter of energy due to the macroscopic granular structure of the weld. However, the influence of weld properties, such as mean grain size and orientation distribution, on the magnitude of scattered ultrasound is not well understood. A hybrid model has been developed to allow the study of grain-scatter effects in austenitic welds. An efficient 2D Finite Element (FE) method is used to calculate the complete scattering response from a single elliptical austenitic grain of arbitrary length and width as a function of the specific inspection frequency. A grain allocation model of the weld is presented to approximate the characteristic structures observed in austenitic welds and the complete scattering behaviour of each grain calculated. This model is incorporated into a semi-analytical framework for a single-element inspection of a typical weld in immersion. Experimental validation evidence is demonstrated indicating excellent qualitative agreement of SNR as a function of frequency and a minimum SNR difference of 2 dB at a centre frequency of 2.25 MHz. Additionally, an example Monte-Carlo study is presented detailing the variation of SNR as a function of the anisotropy distribution of the weld, and the application of confidence analysis to inform inspection development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Solution of the C5G7 Benchmark using the MOCUM Transport Code with ANSYS Unstructured Mesh

The objective of this research is to use finite element analysis (FEA) software to perform the unstructured meshing for reactor geometries and port the mesh information to MOCUM code. This workflow takes the advantage of the mature meshing module of the mainstream FEA; reduces the development effort for nuclear codes; and increases the nuclear code accuracy by using high quality meshes. In this work, ANSYS is used to build the reactor model and perform the unstructured meshing. The mesh information is then imported to the MOCUM code, a method of characteristic (MOC) transport program supporting unstructured meshes, using a developed MATLAB script. This procedure is verified by modeling the C5G7 benchmark problem. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Effect of laser intensity and exposure time on photothermal therapy with nanoparticles heated by a 793-nm diode laser and tissue optical clearing

Laser-induced thermotherapy is a promising method for cancer treatment, the outcome of which is affected by the exposure time. An inappropriate exposure time and laser intensity cause incomplete tumour destruction, tumour regrowth, and metastasis. Also possible is irreversible damage, i.e. death of healthy cells, and so numerical models are necessary to provide an optimised laser intensity and exposure time for different cancerous tumours. In this study, a model based on finite element method (FEM) is used for solving the bio-heat transfer equation and the Arrhenius equation describing tissue damage. The cancerous tumour is considered as a perfect cylinder with a diameter of 20 mm and a thickness of 2, 3, 4, and 5 mm, filled up by highly absorbing nanoparticles and surrounded by healthy cylindrical tissue with a diameter of 40 mm and a length of 10 mm, which ahs a low scattering coefficient due to optical clearing. The results show that 243 s is a safe and appropriate exposure time when a diode laser with a wavelength of 793 nm and intensity of 0.75 W cm{sup −2} together with gold nanorods of concentration 0.0001 % is used for the treatment of a 3-mm-thick tumour. Then, the results are developed and extensive numerical simulations are used to reveal mathematical relationships between two critical parameters, input power and optimised exposure time, for a series of different tumour thicknesses. Treatment protocols are presented.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Ultrasonic fingerprint sensor using a piezoelectric micromachined ultrasonic transducer array integrated with complementary metal oxide semiconductor electronics

This paper presents an ultrasonic fingerprint sensor based on a 24 × 8 array of 22 MHz piezoelectric micromachined ultrasonic transducers (PMUTs) with 100 μm pitch, fully integrated with 180 nm complementary metal oxide semiconductor (CMOS) circuitry through eutectic wafer bonding. Each PMUT is directly bonded to a dedicated CMOS receive amplifier, minimizing electrical parasitics and eliminating the need for through-silicon vias. The array frequency response and vibration mode-shape were characterized using laser Doppler vibrometry and verified via finite element method simulation. The array's acoustic output was measured using a hydrophone to be ∼14 kPa with a 28 V input, in reasonable agreement with predication from analytical calculation. Pulse-echo imaging of a 1D steel grating is demonstrated using electronic scanning of a 20 × 8 sub-array, resulting in 300 mV maximum received amplitude and 5:1 contrast ratio. Because the small size of this array limits the maximum image size, mechanical scanning was used to image a 2D polydimethylsiloxane fingerprint phantom (10 mm × 8 mm) at a 1.2 mm distance from the array.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Simulation Analysis of Geometry and Material Effects for Dropkinson Bar

The reported research is motivated by the need to address a key issue affecting the Dropkinson bar apparatus. This unresolved issue is the interference of the stress wave reflected from the bar-beam boundary with the measurement of the stress-strain response of a material tested in the apparatus. The purpose of the wave beam that is currently connected to the bar is to dissipate the stress wave, but the portion of the wave reflected from the bar-beam boundary is still significant. First, we focused on understanding which parameters affect the reflected wave's arrival time at a strain gauge. Specifically, we used finite-element numerical simulations with the Sierra/SM module to study the effects of various bar-beam connection fixities, alternative wave beam materials, and alternative geometries of the Dropkinson bar system based on a monolithic design. The conclusion of this study is that a partial reflection always occurs at the bar-beam boundary (or, for a monolithic design, at a point where the bar geometry changes). Therefore, given a fixed total length of the bar, it is impossible to increase the reflected wave's arrival time by any significant amount. After reaching this conclusion, we focused instead on trying to minimize the energy of the reflected stress wave circulating up and down through the bar over a relatively long period of time (10 ms). Once again, we used numerical simulations with the Sierra/SM module to investigate the effects of various bar-beam connection fixities, alternative wave beam materials, and parameters of an asymmetric monolithic design of the bar-and-beam system. This study demonstrated that various parameters can significantly affect the energy of the wave reflections, with the difference between best and worst configurations being about one order of magnitude in terms of energy. Based on the obtained results, we conclude with concrete takeaways for Dropkinson bar users and propose potential directions for future research and optimization.

36 MATERIALS SCIENCE↗