Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite Element Method”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Novel refractive index sensor based on fiber bragg grating in nano-bore optical fiber

In this paper, we propose a novel refractive index sensor with large dynamic range based on fiber Bragg grating in the nano-bore optical fiber. The spectral responses of grating on the filling index in the bore with different bore diameters and index ranges are numerically studied by the finite element method and coupled mode theory based transmission matrix method. It is found that the Bragg wavelength of grating is sensitive to the filling index in the bore. The influences of bore diameter on the wavelength sensitivity, resolution and dynamic range of the index sensor are analyzed. There exists a trade off between the linear dynamic range and sensitivity of index when the bore diameter varies. The optimized wavelength sensitivity and resolution of filling index are 53.6923 nm/index unit and 9.3123 × 10{sup −6} index under a linear index range of 1.0–1.4 using 2 µm bore diameter. The optimized linear dynamic range of index is 1.0–1.48 with wavelength sensitivity of 13.3056 nm/index unit and resolution of 3.7578 × 10{sup −5} index using 1 µm bore diameter. Our proposed index sensor has the technical advantages of simple all-fiber structure, large index range, and potential capability to synchronously detect index, temperature and strain.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Residual Stress Asymmetry in Thin Sheets of Double-Sided Shot Peened Aluminum

The local mechanical properties and residual stresses of a double-sided shot-peened aluminum alloy AA7050-T7451 sheet were measured by nanoindentation and simulated by the finite element method (FEM). The variation of the pileup area around the indentation due to residual stresses, as well as changes in hardness as a function of cold working, was assessed. Numerical simulations were conducted to investigate the residual stress profile. When the indentation contact area was corrected to account for pileup induced from residual stresses, the hardness increased about 10% over the range of depths at which compressive residual stresses exist, indicating some strain hardening did occur. The maximum compressive residual stress for the second shot-peened surface is on the order of 10% higher than the first surface; this is ascribed to the sequencing during the peening process, in which the first side peening influenced the stress throughout the part thickness for these thin plate structures. The ability to couple FEM simulations and rapid experimental assessments of residual stresses and hardness allows for future prediction of part thickness effects in double-sided peening processes.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Measurement strategy for rectangular electrical capacitance tomography sensor

To investigate the influence of the measurement strategy for the rectangular electrical capacitance tomography (ECT) sensor, a Finite Element Method (FEM) is utilized to create the model for simulation. The simulation was carried out using COMSOL Multiphysics(trade mark, serif) and Matlab(trade mark, serif). The length-width ratio of the rectangular sensing area is 5. Twelve electrodes are evenly arranged surrounding the pipe. The covering ratio of the electrodes is 90%. The capacitances between different electrode pairs are calculated for a bar distribution. The air of the relative permittivity 1.0 and the material of the permittivity 3.0 are used for the calibration. The relative permittivity of the second phase is 3.0. The noise free and noise data are used for the image reconstruction using the Linear Back Projection (LBP). The measurement strategies with 1-, 2- and 4- electrode excitation are compared using the correlation coefficient. Preliminary results show that the measurement strategy with 2-electrode excitation outperforms other measurement strategies with 1- or 4-electrode excitation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

New developments in ALCYONE 2.0 fuel performance code

ALCYONE is a multi-dimensional (1D, 2D or 3D) fuel performance code dedicated to pressurized water reactors (PWR). This code is co-developed in the PLEIADES platform by CEA, EDF and AREVA and simulates fuel pellets and cladding evolutions during base irradiation, power ramps and accidental situations. Sophisticated chemical and physical models are used. The thermal and mechanical behavior is solved with the finite element method based code CAST3M. ALCYONE 2.0 is the next release of the code and it includes several new features: improvements in incidental and accidental situations modeling, new models available for 2D and 3D schemes and an evolution of neutron model to simulate gadolinium fuel. ALCYONE 2.0 uses a new version of PLEIADES platform which allows the parallelized computation of rod axial slices in order to decrease computation times. To better simulate accidental situations, an axial transfer model for gases and a modeling of heaters for Loss Of Coolant Accident (LOCA) experiments have been added to the code. The axial transfer model can simulate an axial internal pressure gradient in the fuel rod. For power ramps simulations in the 1D scheme, a new model for dish filling is available. This version will be validated on more than 400 rods covering a wide range of irradiation conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Metal Big Area Additive Manufacturing: Process Modeling and Validation

Metal Big Area Additive Manufacturing (mBAAM) is a new additive manufacturing (AM) technology for printing large-scale 3D objects. mBAAM is based on the gas metal arc welding process and uses a continuous feed of welding wire to manufacture an object. An electric arc forms between the wire and the substrate, which melts the wire and deposits a bead of molten metal along the predetermined path. In general, the welding process parameters and local conditions determine the shape of the deposited bead. The sequence of the bead deposition and the corresponding thermal history of the manufactured object determine the long range effects, such as thermal-induced distortions and residual stresses. Therefore, the resulting performance or final properties of the manufactured object are dependent on its geometry and the deposition path, in addition to depending on the basic welding process parameters. Physical testing is critical for gaining the necessary knowledge for quality prints, but traversing the process parameter space in order to develop an optimized build strategy for each new design is impractical by pure experimental means. Computational modeling and optimization may accelerate development of a build process strategy and saves time and resources. Because computational modeling provides these opportunities, we have developed a physics-based Finite Element Method (FEM) simulation framework and numerical models to support the mBAAM process s development and design. In this paper, we performed a sequentially coupled heat transfer and stress analysis for predicting the final deformation of a small rectangular structure printed using the mild steel welding wire. Using the new simulation technologies, material was progressively added into the FEM simulation as the arc weld traversed the build path. In the sequentially coupled heat transfer and stress analysis, the heat transfer was performed to calculate the temperature evolution, which was used in a stress analysis to evaluate the residual stresses and distortions. In this formulation, we assume that physics is directionally coupled, i.e. the effect of stress of the component on the temperatures is negligible. The experiment instrumentation (measurement types, sensor types, sensor locations, sensor placements, measurement intervals) and the measurements are presented. The temperatures and distortions from the simulations show good correlation with experimental measurements. Ongoing modeling work is also briefly discussed.

Additive Manufacturing↗
DOE OSTI2018

Numerical Simulation of Phase Transitions in Type-II Annular Superconductor Using Time-dependent Ginzburg-Landau Equations

Time-dependent Ginzburg-Landau equations were solved by finite element method in two-dimensional space for order parameter and energy components of the annular superconducting sample in steady magnetic fields. Vortices preferred to penetrate from the inner surface of the annulus due to lesser energy required at the concave surface. A transition magnetic field strength was observed in spatial averages of carrier concentration and energy components, showing small bumps and abrupt variations, indicating phase transition from a non-vortex to vortex state. These effects were observed to repeat with every subsequent entry of a set of vortices into the sample; transition magnetic field strength was found to depend inversely on the annular width of the sample. The present work gives a better understanding of energy variations during phase transition from non-vortex to vortex state and predicts that vortex state can be avoided by tuning the wire thickness in practical applications, e.g., superconducting electromagnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Multi-scale simulation of reaction, transport and deactivation in a SBA-16 supported catalyst for the conversion of ethanol to butadiene

Effective design of heterogeneous catalytic systems requires careful coordination of physiochemical phenomena that span orders or magnitude in length and time scales. Multiscale modelling and simulation tools can provide useful insight for understanding and improving the performance of such systems. Mesoporous silica catalyst support materials present a versatile platform wherein support attributes may be tuned to control transport phenomena throughout the system. In this study we develop an integrated multiscale modelling approach for the conversion of ethanol to butadiene over a SBA-16 mesoporous silica supported catalyst. We use molecular dynamics simulations to calculate domain specific diffusivities of reactants and products in the various domains of SBA-16's nanostructure. This is used to calculate a resultant effective diffusivity through its microstructure using finite element method (FEM) models of the tessellated unit cell. The resultant intraparticle effective diffusivity enables the development of a reduced-order reactor level FEM model that is able to implicitly account for the SBA-16's transport effects without the explicit consideration of its detailed geometry. Experimental bench-scale conversion data for the ethanol to butadiene process over SBA-16 is used to extract intrinsic kinetic rate constants for a simplified reaction mechanism using this reduced order bench-scale reactor FEM model. The model is able to reproduce experimental trends in conversion and activity lifetimes. The impact of varying key support attributes such as pore-size and particle size on catalytic activity lifetimes is then explored. We demonstrate the utility of this multi-scale modelling strategy in guiding catalyst design for the development of rational strategies to improve the performance of heterogeneous supported catalytic systems.

09 BIOMASS FUELS↗
DOE OSTI2014

Revealing origin of quasi-one dimensional current transport in defect rich two dimensional materials

The presence of defects in graphene have for a long time been recognized as a bottleneck for its utilization in electronic and mechanical devices. We recently showed that micro four-point probes may be used to evaluate if a graphene film is truly 2D or if defects in proximity of the probe will lead to a non-uniform current flow characteristic of lower dimensionality. In this work, simulations based on a finite element method together with a Monte Carlo approach are used to establish the transition from 2D to quasi-1D current transport, when applying a micro four-point probe to measure on 2D conductors with an increasing amount of line-shaped defects. Clear 2D and 1D signatures are observed at low and high defect densities, respectively, and current density plots reveal the presence of current channels or branches in defect configurations yielding 1D current transport. A strong correlation is found between the density filling factor and the simulation yield, the fraction of cases with 1D transport and the mean sheet conductance. The upper transition limit is shown to agree with the percolation threshold for sticks. Finally, the conductance of a square sample evaluated with macroscopic edge contacts is compared to the micro four-point probe conductance measurements and we find that the micro four-point probe tends to measure a slightly higher conductance in samples containing defects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Investigation on Microsegregation of IN718 Alloy During Additive Manufacturing via Integrated Phase-Field and Finite-Element Modeling

In this work, we apply a multi-scale model combining finite-element method (FEM) and phase-field model (PFM) to simulate the evolution of solidification microstructures at different locations within a molten pool of an additively manufactured IN718 alloy. Specifically, the FEM is used to calculate the shape of molten pool and the relative thermal gradient G at the macroscale. Then, the calculated thermal information is input into PFM for microstructure simulation. Finally, the morphology of solidification structures and formation of Laves phase at different sites are studied and compared. We found that the solidification site with a large angle between the temperature gradient and the preferred crystalline orientation could build up a high niobium (Nb) concentration in the liquid during solidification but has less possibility of forming continuous long chain morphology of Laves phase particles. This finding provides an understanding of the microstructure evolution inside the molten pool of IN718 alloy during solidification. Further, the finding indicates that the site with a large misorientation angle will have a good hot cracking resistance after solidification.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Asynchronous distributed-memory task-parallel algorithm for compressible flows on unstructured 3D Eulerian grids

Here, we discuss the implementation of a finite element method, used to numerically solve the Euler equations of compressible flows, using an asynchronous runtime system (RTS). The algorithm is implemented for distributed-memory machines, using stationary unstructured 3D meshes, combining data-, and task-parallelism on top of the Charm++ RTS. Charm++’s execution model is asynchronous by default, allowing arbitrary overlap of computation and communication. Task-parallelism allows scheduling parts of an algorithm independently of, or dependent on, each other. Built-in automatic load balancing enables continuous redistribution of computational load by migration of work units based on real-time CPU load measurement. The RTS also features automatic checkpointing, fault tolerance, resilience against hardware failure, and supports power-, and energy-aware computation. We demonstrate scalability up to 25 x 10 9 cells at $\mathscr{O}$10 4 compute cores and the benefits of automatic load balancing for irregular workloads. The full source code with documentation is available at https://quinoacomputing.org.

42 ENGINEERING↗
DOE OSTI2016

Deformation of contact surfaces in a vacuum interrupter after high-current interruptions

In a high-current interruption, the contact surface in a vacuum interrupter might be severely damaged by constricted vacuum arcs causing a molten area on it. As a result, a protrusion will be initiated by a transient recovery voltage after current zero, enhancing the local electric field and making breakdowns occur easier. The objective of this paper is to simulate the deformation process on the molten area under a high electric field by adopting the finite element method. A time-dependent Electrohydrodynamic model was established, and the liquid-gas interface was tracked by the level-set method. From the results, the liquid metal can be deformed to a Taylor cone if the applied electric field is above a critical value. This value is correlated to the initial geometry of the liquid metal, which increases as the size of the liquid metal decreases. Moreover, the buildup time of a Taylor cone obeys the power law t = k × E{sup −3}, where E is the initial electric field and k is a coefficient related to the material property, indicating a temporal self-similar characteristic. In addition, the influence of temperature has little impact on the deformation but has great impact on electron emission. Finally, the possible reason to initiate a delayed breakdown is associated with the deformation. The breakdown does not occur immediately when the voltage is just applied upon the gap but is postponed to several milliseconds later when the tip is formed on the liquid metal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Electromagnetic and Rotational Characteristics of a Superconducting Flywheel Energy Storage System Utilizing a Radial-Type High-Temperature Superconducting Bearing

A 2 kW/28.5 kJ superconducting flywheel energy storage system (SFESS) with a radial-type high-temperature superconducting (HTS) bearing was set up to study the electromagnetic and rotational characteristics. The structure of the SFESS as well as the design of its main parts was reported. A mathematical model based on the finite element method (FEM) was established to research the electromagnetic characteristics of the HTS bearing during the levitation process, which show that a part of the magnetic flux penetrates into the edge of the HTS bulks and then goes back to the opposite pole of the permanent magnet rotor (PMR). The induced current mainly distributes in the edge of the HTS bulks, indicating that larger force acts on the edge part of the HTS bulks and probably causes them to crack. The free rotations of the rotor at different steady-state speeds of 2500–5000 rpm and its radial vibration were displayed. The induced voltage of the stator winding of the motor in this process was analyzed. The rotational characteristics are related to the vibration of the rotor. Below the resonant frequency, the vibration increases significantly with the speed. Enhancing the radial stiffness to limit the vibration amplitude of the rotor is an effective approach to improve the speed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Atomistically derived cohesive zone model of intergranular fracture in polycrystalline graphene

Pristine single crystal graphene is the strongest known two-dimensional material, and its nonlinear anisotropic mechanical properties are well understood from the atomic length scale up to a continuum description. However, experiments indicate that grain boundaries in the polycrystalline form reduce the mechanical behavior of polycrystalline graphene. Herein, we perform atomistic-scale molecular dynamics simulations of the deformation and fracture of graphene grain boundaries and express the results as continuum cohesive zone models (CZMs) that embed notions of the grain boundary ultimate strength and fracture toughness. To facilitate energy balance, we employ a new methodology that simulates a quasi-static controlled crack propagation which renders the kinetic energy contribution to the total energy negligible. We verify good agreement between Griffith's critical energy release rate and the work of separation of the CZM, and we note that the energy of crack edges and fracture toughness differs by about 35%, which is attributed to the phenomenon of bond trapping. This justifies the implementation of the CZM within the context of the finite element method (FEM). To enhance computational efficiency in the FEM implementation, we discuss the use of scaled traction-separation laws (TSLs) for larger element sizes. As a final result, we have established that the failure characteristics of pristine graphene and high tilt angle bicrystals differ by less than 10%. This result suggests that one could use a unique or a few typical TSLs as a good approximation for the CZMs associated with the mechanical simulations of the polycrystalline graphene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Cracks and blisters formed close to a silicon wafer surface by He-H co-implantation at low energy

We have studied the effect of reducing the implantation energy towards low keV values on the areal density of He and H atoms stored within populations of blister cavities formed by co-implantation of the same fluence of He then H ions into Si(001) wafers and annealing. Using a variety of experimental techniques, we have measured blister heights and depth from the surface, diameter, areal density of the cracks from which they originate as functions of implantation energy and fluence. We show that there is a direct correlation between the diameters of the cracks and the heights of the associated blisters. This correlation only depends on the implantation energy, i.e., only on the depth at which the cracks are located. Using finite element method modeling, we infer the pressure inside the blister cavities from the elastic deformations they generate, i.e., from the height of the blisters. From this, we demonstrate that the gas pressure within a blister only depends on the diameter of the associated crack and not on its depth position and derive an analytical expression relating these parameters. Relating the pressure inside a blister to the respective concentrations of gas molecules it contains, we deduce the areal densities of He and H atoms contained within the populations of blisters. After low-energy implantations (8 keV He{sup +}, 3 keV H{sup +}), all the implanted He and H atoms contribute to the formation of the blisters. There is no measurable exo-diffusion of any of the implanted gases, in contrast to what was assumed at the state of the art to explain the failure of the Smart-Cut technology when using very low energy ion implantation for the fabrication of ultra-thin layers. Alternative explanations must be investigated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2013

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

Plasmons in inhomogeneously doped neutral and charged graphene nanodisks

We study plasmons in graphene nanodisks including the effect of inhomogeneity in the distribution of the doping charge. Specifically, we discuss the following two configurations: charged disks containing a fixed amount of additional carriers, which are self-consistently distributed along the surface to produce a uniform DC potential; and neutral disks exposed to a neighboring external point charge. A suitable finite-element method is elaborated to compute the charge density associated with the plasmons in the electrostatic limit. For charged disks, we find dipolar plasmons similar to those of uniformly doped graphene structures, in which the plasmon induced charge piles up near the edges. In contrast, in neutral disks placed near an external point charge, plasmons are strongly localized away from the edges. Surprisingly, a single external electron is enough to trap plasmons. The disks also display axially symmetric dark-plasmons, which can be excited through external illumination by coupling them to a neighboring metallic element. Our results have practical relevance for graphene nanophotonics under inhomogeneous doping conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Highly sensitive SPR PCF biosensors based on Ag/TiN and Ag/ZrN configurations

We numerically present and analyze a surface plasmon resonance (SPR) bimetallic photonic crystal fiber (PCF) biosensor using full vectorial finite element method. The bimetallic configuration is based on Silver and (TiN or ZrN) as alternative plasmonic materials. The reported design relies on PCF with an external channel to house the analyte. The inner surface of the channel is coated by a silver layer followed by TiN or ZrN to protect the silver layer from oxidation. The proposed biosensor can be implemented to detect the unknown analytes, organic chemicals, biological analytes, and biomolecules via the sensitive resonant peaks to the analyte refractive index variations for both polarized modes with high linearity. In this study, a comparison is made between the Ag/TiN configuration and the Ag/ZrN counterpart in terms of sensitivity, losses, and linearity. The geometrical parameters are optimized to achieve high refractive index sensitivities of 7000 nm/RIU for quasi-transverse electric (TE) mode and 6900 nm/RIU for quasi-transverse magnetic (TM) mode, respectively for the Ag/TiN configuration. However, the quasi TE and quasi TM modes achieve high sensitivity of 5300 and 5400 nm/RIU, respectively using the Ag/ZrN configuration. The standard PCF fabrication technologies can be used to fabricate the proposed SPR PCF biosensors.

36 MATERIALS SCIENCE↗