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 793 records · Page 44

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↗

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↗

Optimization of ultrasonic array inspections using an efficient hybrid model and real crack shapes

Models which simulate the interaction of ultrasound with cracks can be used to optimize ultrasonic array inspections, but this approach can be time-consuming. To overcome this issue an efficient hybrid model is implemented which includes a finite element method that requires only a single layer of elements around the crack shape. Scattering Matrices are used to capture the scattering behavior of the individual cracks and a discussion on the angular degrees of freedom of elastodynamic scatterers is included. Real crack shapes are obtained from X-ray Computed Tomography images of cracked parts and these shapes are inputted into the hybrid model. The effect of using real crack shapes instead of straight notch shapes is demonstrated. An array optimization methodology which incorporates the hybrid model, an approximate single-scattering relative noise model and the real crack shapes is then described.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Extended isogeometric analysis of multi-material and multi-physics problems using hierarchical B-splines

Here, this paper presents an immersed, isogeometric finite element framework to predict the response of multi-material, multi-physics problems with complex geometries using locally refined discretizations. To circumvent the need to generate conformal meshes, this work uses an extended finite element method (XFEM) to discretize the governing equations on non-conforming, embedding meshes. A flexible approach to create truncated hierarchical B-splines discretizations is presented. This approach enables the refinement of each state variable field individually to meet field-specific accuracy requirements. To obtain an immersed geometry representation that is consistent across all hierarchically refined B-spline discretizations, the geometry is immersed into a single mesh, the XFEM background mesh, which is constructed from the union of all hierarchical B-spline meshes. An extraction operator is introduced to represent the truncated hierarchical B-spline bases in terms of Lagrange shape functions on the XFEM background mesh without loss of accuracy. The truncated hierarchical B-spline bases are enriched using a generalized Heaviside enrichment strategy to accommodate small geometric features and multi-material problems. The governing equations are augmented by a formulation of the face-oriented ghost stabilization enhanced for locally refined B-spline bases. We present examples for two- and three-dimensional linear elastic and thermo-elastic problems. The numerical results validate the accuracy of our framework. The results also demonstrate the applicability of the proposed framework to large, geometrically complex problems.

42 ENGINEERING↗

Characterization and performance of a field aligned ion cyclotron range of frequency antenna in Alcator C-Mod

Ion cyclotron range of frequency (ICRF) heating is expected to provide auxiliary heating for ITER and future fusion reactors where high Z metallic plasma facing components (PFCs) are being considered. Impurity contamination linked to ICRF antenna operation remains a major challenge particularly for devices with high Z metallic PFCs. Here, we report on an experimental investigation to test whether a field aligned (FA) antenna can reduce impurity contamination and impurity sources. We compare the modification of the scrape of layer (SOL) plasma potential of the FA antenna to a conventional, toroidally aligned (TA) antenna, in order to explore the underlying physics governing impurity contamination linked to ICRF heating. The FA antenna is a 4-strap ICRF antenna where the current straps and antenna enclosure sides are perpendicular to the total magnetic field while the Faraday screen rods are parallel to the total magnetic field. In principle, alignment with respect to the total magnetic field minimizes integrated E|| (electric field along a magnetic field line) via symmetry. A finite element method RF antenna model coupled to a cold plasma model verifies that the integrated E|| should be reduced for all antenna phases. Monopole phasing in particular is expected to have the lowest integrated E||. Consistent with expectations, we observed that the impurity contamination and impurity source at the FA antenna are reduced compared to the TA antenna. In both L and H-mode discharges, the radiated power is 20%–30% lower for a FA-antenna heated discharge than a discharge heated with the TA-antennas. However, inconsistent with expectations, we observe RF induced plasma potentials (via gas-puff imaging and emissive probes to be nearly identical for FA and TA antennas when operated in dipole phasing). Moreover, the highest levels of RF-induced plasma potentials are observed using monopole phasing with the FA antenna. Thus, while impurity contamination and sources are indeed reduced with the FA antenna configuration, the mechanism determining the SOL plasma potential in the presence of ICRF and its impact on impurity contamination and sources remains to be understood.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Analysis of the Effect of Infrared Radiation on Cryogenic Inertial Confinement Fusion Targets

The present work applies the finite element method to calculate the maximum allowable time that cryogenic inertial confinement fusion (ICF) targets can be exposed to infrared radiation (IR). Hence, a 3-D numerical model integrated with discrete coordinate radiation model was developed to investigate the influence of transmittance of the laser entrance holes (LEHs) and boundary conditions on the temperature field distribution and the maximum DT layer deterioration time for CH, Be, and diamond capsules. Our study shows that introducing such a radiation model can accurately obtain more detailed spatial and temporal distribution information in the ICF targets. The simulation results demonstrate that the Be and diamond capsules provided much better temperature field homogenization than the CH capsule under equivalent boundary conditions, but the CH capsule was heated more by IR radiation than the Be and diamond. In addition, the maximum DT layer deterioration time was significantly increased to 3 s when decreasing the transmittance of the LEH from 0.2 to 0.01. However, either reducing the capsule IR absorption or increasing the inner hohlraum IR absorption demonstrated no conclusive increase in the maximum DT layer deterioration time. These results are expected to provide useful parameters in the design of cryogenic targets and shroud systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Permeability Calculation in a Fracture Network - 12197

Laminar flow of a viscous fluid in the pore space of a saturated fractured rock medium is considered to calculate the effective permeability of the medium. The effective permeability is determined from the flow field which is calculated numerically by using the finite element method. The computation of permeability components is carried out with a few different discretizations for a number of fracture arrangements. Various features such as flow field in the fracture channels, the convergence of permeability, and the variation of permeability among different fracture networks are discussed. The longitudinal permeability in general appears greater than the transverse ones. The former shows minor variations with fracture arrangement whereas the latter appears to be more sensitive to the arrangement. From the calculations of the permeability in a rock medium with a fracture network (two parallel fractures aligned in the direction of 45-deg counterclockwise from the horizontal and two connecting fractures(narrowing, parallel and widening) the following conclusions are drawn. 1. The permeability of fractured medium not only depends on the primary orientation of the main fractures but also is noticeably influenced by the connecting fractures in the medium. 2. The transverse permeability (the permeability in the direction normal to the direction of the externally imposed macro-scale pressure gradient) is only a fraction of the longitudinal one, but is sensitive to the arrangement of the connecting fractures. 3. It is important to figure out the pattern of the fractures that connect (or cross) the main fractures for reliable calculation of the transverse permeability. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗

Effect of interfacial transition zone on the Young's modulus of carbon nanofiber reinforced cement concrete

Identifying the properties of the region where the cement paste meets the aggregate surface (interfacial transition zone, ITZ) is critical to understanding the strength and fracture behavior of carbon nanofiber (CNF) reinforced cement concrete. In this study, the finite element method is employed to investigate the effect of the ITZ on the Young's modulus of cement concrete made with CNF. The numerical models for cement concrete with and without CNF are constructed based on the digital image processing technique. To consider the interface effect, the concept of ''effective aggregate'' is put forward, namely, the Young's modulus and Poisson's ratio for each aggregate particle are replaced by an effective Young's modulus and effective Poisson's ratio, in which the effect of thickness and the Young's modulus of ITZ is taken into account in an averaged manner. Then, the quantitative nanomechanical mapping (based on atomic force microscopy) technique is adopted to measure the Young's modulus and thickness of ITZ with and without CNFs, which are further used as the input parameters in the numerical model. The numerical simulation results are verified by experimental testing, which indicates that the ITZ effect should be considered when implementing the numerical simulation. In addition, this analysis shows that compared with the plain cement concrete, CNFs can greatly enhance the mechanical properties of ITZ, which will in turn improve the Young's modulus of cement concrete significantly. Finally, the effect of thickness and Young's modulus of ITZ on the Young's modulus of CNF reinforced cement concrete is also discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Fast Reactor Core Analysis Using AZNHEX Neutron Diffusion Code

The AZTLAN platform project emerged as an initiative to rank Mexico at a competitive level in software aspects for the analysis of nuclear reactors, in an effort to achieve autonomy in this area. The AZTLAN platform is conformed by a thermo-hydraulic code (AZTHECA), a neutron transport code (AZTRAN) and two neutron diffusion codes (AZKIND, AZNHEX). AZNHEX is a 3D diffusion code that solves numerically the time dependent neutron diffusion equations in hexagonal z-geometry using the RTN0 nodal finite element method. With this code, it is possible to simulate and analyze light water reactors, such as the Russian VVER, and variety of fast reactors. In this work, a sodium-cooled fast reactor core in steady state, with control rods fully withdrawn, is simulated. The comparison and validation of the values obtained with AZNHEX was done with the results obtained with MCNP6 Monte Carlo code.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validating corrosion models: Influence of physical properties

In chloride containing environments, two metals in physical contact can undergo galvanic corrosion limiting the lifetime of components. Being able to accurately predict galvanic corrosion damage distributions over time has not been widely presented in literature. Therefore, Finite Element Method (FEM) corrosion models were experimentally validated for two galvanic couples as a function of governing equations, environment, anode material, anode:cathode ratio, and time. Carbon steel/stainless steel (CS/SS) and zinc/stainless steel (Zn/SS) galvanic couples were exposed to NaCl solutions at room temperatures for up to 14 days. For the galvanic couples and environments, the Laplace equation with variable conductivity and reactions was sufficient to model the experimental corrosion damage. Increasing the chloride concentration for the CS/SS galvanic couple, regardless of the anode:cathode ratio, decreased the observed and modeled corrosion damage. Decreasing the anode:cathode ratio (i.e., increasing the cathode length), increased the experimentally observed and modeled corrosion damage. For small anode:cathode ratios, the governing equations deviate over long time periods with the Nernst-Plank equation being conservative. For the Zn/SS galvanic couple, the cathode length controls the dominant cathodic reduction reaction. For small cathode lengths, the hydrogen evolution reaction is dominant. For larger cathode lengths, the oxygen reduction reaction is dominant. The results are discussed with regard to the influence of solution chemistry, ohmic drop, and governing reactions. Overall, validated FEM models were presented, and the resultant models, physics, and mechanisms can be applied to other corrosion scenarios with confidence.

Carbon steel↗

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗

A modified moment-fitted integration scheme for X-FEM applications with history-dependent material data

Here, we present a strategy for the numerical integration of partial elements with the eXtended finite element method (X-FEM). The new strategy is specifically designed for problems with propagating cracks through a bulk material that exhibits inelasticity. Following a standard approach with the X-FEM, as the crack propagates into new regions of the domain, elements are split into several copies that contain pieces of the domain on either side of the crack. We examine quadrature rules that have sufficient accuracy to calculate stiffness matrices regardless of the orientation of the crack with respect to the element. This permits the number of integration points within elements to remain constant as a crack propagates, and for state data to be easily transferred between successive discretizations. In order to maintain weights that are strictly positive, we propose an approach that blends moment-fitted weights with volume-fraction based weights. To demonstrate the efficacy of this simple approach, we present results from numerical tests and examples with both elastic and plastic material response.

42 ENGINEERING↗

Development of a graphite probe calorimeter for absolute clinical dosimetry

The aim of this work is to present the numerical design optimization, construction, and experimental proof of concept of a graphite probe calorimeter (GPC) conceived for dose measurement in the clinical environment (U.S. provisional patent 61/652,540). A finite element method (FEM) based numerical heat transfer study was conducted using a commercial software package to explore the feasibility of the GPC and to optimize the shape, dimensions, and materials used in its design. A functioning prototype was constructed inhouse and used to perform dose to water measurements under a 6 MV photon beam at 400 and 1000 MU/min, in a thermally insulated water phantom. Heat loss correction factors were determined using FEM analysis while the radiation field perturbation and the graphite to water absorbed dose conversion factors were calculated using Monte Carlo simulations. The difference in the average measured dose to water for the 400 and 1000 MU/min runs using the TG-51 protocol and the GPC was 0.2% and 1.2%, respectively. Heat loss correction factors ranged from 1.001 to 1.002, while the product of the perturbation and dose conversion factors was calculated to be 1.130. The combined relative uncertainty was estimated to be 1.4%, with the largest contributors being the specific heat capacity of the graphite (type B, 0.8%) and the reproducibility, defined as the standard deviation of the mean measured dose (type A, 0.6%). By establishing the feasibility of using the GPC as a practical clinical absolute photon dosimeter, this work lays the foundation for further device enhancements, including the development of an isothermal mode of operation and an overall miniaturization, making it potentially suitable for use in small and composite radiation fields. It is anticipated that, through the incorporation of isothermal stabilization provided by temperature controllers, a subpercent overall uncertainty will be achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

U-10Mo/Zr Interface Modeling using a Microstructure-Based FEM Approach

The U-10Mo in low enrichments (LEU) has been identified as the most promising alternative to the current highly enriched uranium (HEU) used in the United States’ fleet of high performance research reactors (USHPRRs). The nominal configuration of the new LEU U-10Mo plate fuel comprises a U-10Mo fuel foil enriched to slightly less than 20% U-235 (0.08” to 0.02” thick), a thin Zr interlayer/diffusion barrier (25 m thick) and a relatively thick outer can of 6061 aluminum. Currently, the Zr interlayer is clad by hot roll bonding. Previous studies and observations revealed a thinning of the zirconium (Zr) layer during this fuel fabrication process, which is not desirable from the fuel performance perspective. Coarse UMo grains, dendritic structures, Mo concentration segregation, carbides, and porosity are present in the as-cast material and can lead to a nonuniform UMo/Zr interface. The purpose of the current work is to investigate the effects of these microstructural parameters on the Zr coating variation. A microstructure-based finite-element method model was used in this work, and a study on the effect of homogenization on the interface between U-10Mo and Zr was conducted. The model uses actual backscattered electron–scanning electron microscopy microstructures, Mo concentrations, and mechanical properties to predict the behavior of a representative volume element under compressive loading during the rolling process. The model successfully predicted the experimentally observed thinning of the Zr layer in the as-cast material. The model also uses results from a homogenization model as an input, and a study on the effect of different levels of homogenization on the interface indicated that homogenization helps decrease this thinning. This model can be considered a predictive tool representing a first step for model integration and an input into a larger fuel fabrication performance model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Feasibility of producing a short, high energy s-band linear accelerator using a klystron power source

Purpose: To use a finite-element method (FEM) model to study the feasibility of producing a short s-band (2.9985 GHz) waveguide capable of producing x-rays energies up to 10 MV, for applications in a linac-MR, as well as conventional radiotherapy. Methods: An existing waveguide FEM model developed by the authors' group is used to simulate replacing the magnetron power source with a klystron. Peak fields within the waveguide are compared with a published experimental threshold for electric breakdown. The RF fields in the first accelerating cavity are scaled, approximating the effect of modifications to the first coupling cavity. Electron trajectories are calculated within the RF fields, and the energy spectrum, beam current, and focal spot of the electron beam are analyzed. One electron spectrum is selected for Monte Carlo simulations and the resulting PDD compared to measurement. Results: When the first cavity fields are scaled by a factor of 0.475, the peak magnitude of the electric fields within the waveguide are calculated to be 223.1 MV/m, 29% lower than the published threshold for breakdown at this operating frequency. Maximum electron energy increased from 6.2 to 10.4 MeV, and beam current increased from 134 to 170 mA. The focal spot FWHM is decreased slightly from 0.07 to 0.05 mm, and the width of the energy spectrum increased slightly from 0.44 to 0.70 MeV. Monte Carlo results show d{sub max} is at 2.15 cm for a 10 Multiplication-Sign 10 cm{sup 2} field, compared with 2.3 cm for a Varian 10 MV linac, while the penumbral widths are 4.8 and 5.6 mm, respectively. Conclusions: The authors' simulation results show that a short, high-energy, s-band accelerator is feasible and electric breakdown is not expected to interfere with operation at these field strengths. With minor modifications to the first coupling cavity, all electron beam parameters are improved.

60 APPLIED LIFE SCIENCES↗

Simulated effects of sample size and grain neighborhood on the modeling of extreme value fatigue response

Assessing the size of representative volume elements (RVEs) for fatigue-related applications is challenging. A RVE relevant to random microstructure requires a volume of material that is sufficiently large to capture the grain/phase heterogeneity that captures all statistical moments of the distribution of the driving force for fatigue crack formation at “hot spot” grains. Consequently, the large size of a microstructure RVE required to study fatigue phenomena is largely computationally intractable and difficult to explore. A more realistic objective in this work is to systematically study, as a function of the size of a statistical sample of microstructure, trends towards convergence of the simulated distribution of driving force for fatigue crack formation. Our present work accordingly leverages the recently developed open-source PRISMS-Fatigue framework to examine the trends in convergence of extreme value distributions (EVD) of Fatigue Indicator Parameters (FIPs) in progressively larger polycrystalline microstructure realizations of FCC Al alloy 7075-T6 using crystal plasticity finite element method simulations. The results are compared to the traditional method in which ensembles of statistical volume elements (SVEs) are simulated to build up statistics intended to approximate those associated with a larger volume of material. The convergence of EVDs with increase of size of a SVE of microstructure is closely related to the extent of grain nearest neighbor (NN) interactions. Accordingly, the sensitivity of the local micromechanical response at hot spot grains is quantitatively investigated by systematically varying the orientations of NN grains. Results indicate that SVEs with cubic crystallographic texture tend towards convergence of the EVD of FIPs with tens of thousands of grains while the random and rolled textures require larger volumes. Simple relationships based on microstructure parameters (e.g., Schmid Factor, grain size, NN misorientation) do not completely correlate to fatigue hot spot grains. Finally, the sensitivity of the extreme value fatigue response at hot spot grains extends to the 3rd NN when a single neighborhood grain orientation is altered.

36 MATERIALS SCIENCE↗

Analysis of breathing cracks using vibrations

The influence of half-ellipse-shaped fatigue cracks and machined through-thickness cuts on the impact response of cantilever beam samples were analysed using the Finite Element Method and compared to experimental data for different crack depths. The effect of sample width was analysed for geometries with fatigue cracks (aspect ratio a/c = 0.45) and through-thickness cuts (rectangular cracks). The differences between both types of damage were analysed numerically and experimentally. The effect of sample width was studied and it was found that the behaviour of the cantilever with a fatigue-crack approaches that of the saw cut sample when the fatigue crack reaches the edges of the test piece. Acceleration data were analysed using the Fast Fourier Transform algorithm and Power Density Spectra (PSD) were obtained. For fatigue-cracked samples, harmonics of the frequencies of natural modes of oscillation were found. Analytical modelling was done using a mass-spring-damper system for studying the bilinear behaviour of the stiffness element. Harmonics found were assigned to the breathing behaviour of the crack. The detection threshold using frequency shifts inspection was twice as high as using spectral analysis. This would allow crack detection at half of the fatigue life of the specimen.

42 ENGINEERING↗