Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary 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

Solidification of II-VI Compounds in a Rotating Magnetic Field

This project is aimed at using a rotating magnetic field (RMF) to control fluid flow and transport during directional solidification of elemental and compound melts. Microgravity experiments have demonstrated that small amounts of residual acceleration of less than a micro-g can initiate and prolong fluid flow, particularly when there is a static component of the field perpendicular to the liquid solid interface. Thus a true diffusion boundary layer is not formed, and it becomes difficult to verify theories of solidification or to achieve diffusion controlled solidification. The RMF superimposes a stirring effect on an electrically conducting liquid, and with appropriate field strengths and frequencies, controlled transport of material through a liquid column can be obtained. As diffusion conditions are precluded and complete mixing conditions prevail, the technique is appropriate for traveling solvent zone or float zone growth methods in which the overall composition of the liquid can be maintained throughout the growth experiment. Crystals grown by RMF techniques in microgravity in previous, unrelated missions have shown exceptional properties. The objective of the project is two-fold, namely (1) using numerical modeling to simulate the behavior of a solvent zone with applied thermal boundary conditions and demonstrate the effects of decreasing gravity levels, or an increasing applied RMF, or both, and (2) to grow elements and II-VI compounds from traveling solvent zones both with and without applied RMFs, and to determine objectively how well the modeling predicts solidification parameters. Numerical modeling has demonstrated that, in the growth of CdTe from a tellurium solution, a rotating magnetic field can advantageously modify the shape of the liquid solid interface such that the interface is convex as seen from the liquid. Under such circumstances, the defect structure is reduced as any defects which are formed tend to grow out and not propagate. The flow of liquid, however, is complex due to the competing flow induced by the rotating magnetic field and the buoyancy driven convection. When the acceleration forces are reduced to one thousandth of gravity, the flow pattern is much simplified and well controlled material transport through the solvent zone can be readily achieved. Triple axis diffractometry and x-ray synchrotron topography have demonstrated that there is no significant improvement in crystal quality for HgCdTe grown on earth from a tellurium solution when a rotating magnetic field is applied. However, modeling shows that the flow in microgravity with a rotating magnetic field would produce a superior product.

Gillies, D. C.↗

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

High-temperature Majorana fermions in magnet-superconductor hybrid systems

Magnet-superconductor hybrid (MSH) structures represent one of the most promising platforms to realize, control, and manipulate Majorana modes using scanning tunneling methods. By depositing either chains or islands of magnetic atoms on the surface of a conventional, elemental superconductor such as Pb or Re, topological superconducting phases can emerge. They feature either localized Majorana bound states at the chain ends or dispersing chiral Majorana modes at the island's boundary. However, some of these experiments have not reached the spectral resolution to clearly distinguish between topological Majorana and trivial Shiba states due to very small superconducting gap sizes and experiments performed at sub-Kelvin temperatures. Here we consider superconducting substrates with unconventional spin-singlet pairing, including high-temperature d -wave and extended s -wave superconductors. We derive topological phase diagrams and compute edge states for cylinder and island geometries and discuss their properties. Several time-reversal invariant topological superconducting phases of the Zhang-Kane-Mele type are found and discussed. Addiotnally, we review one-dimensional MSH structures and show that parameters to realize topologically nontrivial magnetic chains embedded into a larger, two-dimensional substrate differ from the purely one-dimensional case. Quite generally we find that unconventional superconducting substrates work as well as the conventional s -wave substrates to realize topological phases. In particular, iron-based pnictide and chalcogenide superconductors are the most promising class of substrates for future high-temperature MSH systems.

36 MATERIALS SCIENCE↗

Assessing Nonlinear Polarization in Electrochemical Cells using AC Impedance Spectroscopy

AC impedance spectroscopy is widely used to evaluate performance limitations in energy storage and conversion devices (e.g., batteries, supercapacitors, and fuel cells). This work shows that integrating the resistive elements in an equivalent circuit as functions of steady-state current enables one to recover overpotentials associated with different processes (e.g., ion migration, charge transfer, and diffusion) in nonlinear electrochemical power supplies. Closed form expressions for diffusion overpotentials are derived using this method for transmissive and reflective boundary conditions and three electrode symmetries (planar, cylindrical, and spherical). Discussion is also extended to macroscopically homogenous porous electrodes which are relevant for most real-world devices. Overall, the approach described herein is a powerful tool to identify rate-limiting steps and guide material/component design.

25 ENERGY STORAGE↗

Flow prediction over a transport multi-element high-lift system and comparison with flight measurements

Accurate prediction of surface-pressure distributions, merging boundary-layers, and separated-flow regions over multi-element high-lift airfoils is required to design advanced high-lift systems for efficient subsonic transport aircraft. The availability of detailed measurements of pressure distributions and both averaged and time-dependent boundary-layer flow parameters at flight Reynolds numbers is critical to evaluate computational methods and to model the turbulence structure for closure of the flow equations. Several detailed wind-tunnel measurements at subscale Reynolds numbers were conducted to obtain detailed flow information including the Reynolds-stress component. As part of a subsonic-transport high-lift research program, flight experiments are conducted using the NASA-Langley B737-100 research aircraft to obtain detailed flow characteristics for support of computational and wind-tunnel efforts. Planned flight measurements include pressure distributions at several spanwise locations, boundary-layer transition and separation locations, surface skin friction, as well as boundary-layer profiles and Reynolds stresses in adverse pressure-gradient flow.

Vijgen, P. M. H. W.↗

A unit-cell model of textile composite beams for predicting stiffness properties

Flexural stiffness properties of a textile composite beam are obtained from a finite-element model of the unit cell. Three linearly independent deformations, namely, pure extension, pure bending and pure shear, are applied to the unit cell. The top and bottom surfaces of the beam are assumed to be traction free. Periodic boundary conditions on the lateral boundaries of the unit cell are enforced by multi-point constraint elements. From the forces acting on the unit cell, the flexural stiffness coefficients of the composite beam are obtained. The difficulties in determining the transverse shear stiffness are discussed, and a modified approach is presented. The methods are first verified by applying them to isotropic and bimaterial beams for which the results are known, and then illustrated for a simple plain-weave textile composite.

Sankar, Bhavani V.↗

Elastic interactions of a fatigue crack with a micro-defect by the mixed boundary integral equation method

In this paper, the mixed boundary integral equation method is developed to study the elastic interactions of a fatigue crack and a micro-defect such as a void, a rigid inclusion or a transformation inclusion. The method of pseudo-tractions is employed to study the effect of a transformation inclusion. An enriched element which incorporates the mixed-mode stress intensity factors is applied to characterize the singularity at a moving crack tip. In order to evaluate the accuracy of the numerical procedure, the analysis of a crack emanating from a circular hole in a finite plate is performed and the results are compared with the available numerical solution. The effects of various micro-defects on the crack path and fatigue life are investigated. The results agree with the experimental observations.

Lua, Yuan J.↗

Parallel Aeroelastic Analysis Using ENSAERO and NASTRAN

A high fidelity parallel static structural analysis capability is created and interfaced to the multidisciplinary analysis package ENSAERO-MPI of Ames Research Center. This new module replaced ENSAERO's lower fidelity simple finite element and modal modules. Full aircraft structures may be more accurately modeled using the new finite element capability. Parallel computation is performed by breaking the full structure into multiple substructures. This approach is conceptually similar to ENSAERO's multi-zonal fluid analysis capability. The new substructure code is used to solve the structural finite element equations for each substructure in parallel. NASTRAN/COSMIC is utilized as a front end for this code. Its full library of elements can be used to create an accurate and realistic aircraft mode. It is used to create the stiffness matrices for each sub-structure. The new parallel code then uses an iterative preconditioned conjugate gradient method to solve the global structural equations for the sub-structure boundary nodes. Results are presented for a wing-body configuration.

Eldred, Lloyd B.↗

Convection cells in the early lunar magma ocean - Trace-element evidence

A method was proposed for recognizing suites of samples which are related to the differentiation of distinct and regionally separated magmas in the early lunar surface ocean. The sample suites were identified by three Cl-P2O5 compositional groups. New experimental data continue to add samples to these groups. Each magma was associated with a convection cell the boundaries of which are proposed. Independent geochemical evidence supporting a limited number of early magmas is cited. Two new possible Cl-P2O5 groups are identified. Samples not indigenous to a given site and their possible source regions are identified. It is proposed that many Apollo 17 basalts are not locally derived.

Jovanovic, S.↗

Cadmium Selective Etching in CdTe Solar Cells Produces Detrimental Narrow-Gap Te in Grain Boundaries

Recent advances in design and processing technology have made possible commercialization of polycrystalline (px)-CdTe as a photovoltaic absorber. Grain boundaries (GBs) are the most prominent structural defects in these devices and undergo significant changes during device fabrication. However, the effects of device fabrication processes on these GBs are not entirely understood. Prevailing models of GBs in thin-film photovoltaics consider individual GBs to have homogeneous properties in their area. Here, using an aberration-corrected scanning transmission electron microscope (STEM)-based low-loss and core-loss electron energy-loss spectroscopy (EELS), we show that back-surface etching of CdTe leads to inhomogeneity within individual grain boundaries. We observe that etching the back surface leads to the conversion of a region of GBs from CdTe to an elemental Te, which has an only 0.33 eV band gap, as deep as 1 μm from the back surface. The presence of elemental Te in GBs this deep into the absorber layer will increase recombination in the absorber layer and limit the extractable open-circuit voltage, thus reducing device efficiency. However, additive methods for back contact formation such as deposition of Te, ZnTe, or other materials preserve the CdTe stoichiometry of the GBs. Thus, especially for the next generations of CdTe-based cells having longer minority carrier diffusion length and/or thinner absorber layers, additive back contacting methods are superior.

36 MATERIALS SCIENCE↗

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗

27 Al NMR diffusometry of Al 13 Keggin nanoclusters

Although nanometer-sized aluminum hydroxide clusters (i.e., ϵ-Al 13 , [Al 13 O 4 (OH) 24 (H 2 O) 12 ] 7+ ) command a central role in aluminum ion speciation and transformations between minerals, measurement of their translational diffusion is often limited to indirect methods. Here, in this work, 27 Al pulsed field gradient stimulated echo nuclear magnetic resonance (PFGSTE NMR) spectroscopy has been applied to the AlO 4 core of the ϵ-Al 13 cluster with complementary theoretical simulations of the diffusion coefficient and corresponding hydrodynamic radii from a boundary element-based calculation. The tetrahedral AlO 4 center of the ϵ-Al 13 cluster is symmetric and exhibits only weak quadrupolar coupling, which results in favorable T 1 and T 2 27 Al NMR relaxation coefficients for 27 Al PFGSTE NMR studies. Stokes–Einstein relationship was used to relate the 27 Al diffusion coefficient of the ϵ-Al 13 cluster to the hydrodynamic radius for comparison with theoretical simulations, dynamic light scattering from literature, and previously published 1 H PFGSTE NMR studies of chelated Keggin clusters. This first-of-its-kind observation proves that 27 Al PFGSTE NMR diffusometry can probe symmetric Al environments in polynuclear clusters of greater molecular weight than previously considered.

27Al↗

Assessment of Lost Circulation Material Particle-Size Distribution on Fracture Sealing: A Numerical Study

Summary Lost circulation materials (LCMs) are essential to combat fluid loss while drilling and may put the whole operation at risk if a proper LCM design is not used. The focus of this research is understanding the function of LCMs in sealing fractures to reduce fluid loss. One important consideration in the success of fracture sealing is the particle-size distribution (PSD) of LCMs. Various studies have suggested different guidelines for obtaining the best size distribution of LCMs for effective fracture sealing based on limited laboratory experiments or field observations. Hence, there is a need for sophisticated numerical methods to improve the LCM design by providing some predictive capabilities. In this study, computational fluid dynamics (CFD) and discrete element methods (DEM) numerical simulations are coupled to investigate the influence of PSD of granular LCMs on fracture sealing. Dimensionless variables were introduced to compare cases with different PSDs. We validated the CFD-DEM model in reproducing specific laboratory observations of fracture-sealing experiments within the model boundary parameters. Our simulations suggested that a bimodally distributed blend would be the most effective design in comparison to other PSDs tested here.

Engineering↗

Compositional Homogeneity in Rapidly Solidified 316L stainless steel produced by laser powder-bed fusion processing

The segregation of alloying elements in rapidly solidified 316L stainless steel produced using additive manufacturing technologies is widely reported throughout literature. The microstructure of thin walls manufactured from 316L stainless steel using laser powder bed fusion are examined for the presence of alloy segregation, i.e. the rejection of solutes to grain boundaries. Metallographic methods including mechanical polishing and chemical etching are used to prepare samples for microanalysis m

Jankowski, Alan F. [Sandia National Laboratories (↗

Bio-, Magneto- and event-stratigraphy across the K-T boundary

Determining the time and the time structure of rare events in geology can be accomplished by applying three different and independent stratigraphic methods: Biostratigraphy, magneto-stratigraphy and event-stratigraphy. The optimal time resolution of the two former methods is about 1000 years, while by means of event-stratigraphy a resolution of approximately one year can be achieved. For biostratigraphy across the Cretaceous-Tertiary (K-T) boundary micro- and nannofossils have been found best suited. The qualitative and quantitative analyses of minerals and trace elements across the K-T boundary show anomalies on a millimeter scale and permit conclusions regarding the time structure of the K-T event itself. The results of the analyses find a most consistent explanation by the assumption of an extraterrestrial impact. The main portion of the material rain from the atmosphere evidently was deposited within a short time. The long-time components consist of the finest portion of the material rain from the atmosphere and the transported and redeposited fall-out.

Preisinger, A.↗

Discrete Roughness Effects on Shuttle Orbiter at Mach 6

Discrete roughness boundary layer transition results on a Shuttle Orbiter model in the NASA Langley Research Center 20-Inch Mach 6 Air Tunnel have been reanalyzed with new boundary layer calculations to provide consistency for comparison to other published results. The experimental results were previously obtained utilizing the phosphor thermography system to monitor the status of the boundary layer via global heat transfer images of the Orbiter windward surface. The size and location of discrete roughness elements were systematically varied along the centerline of the 0.0075-scale model at an angle of attack of 40 deg and the boundary layer response recorded. Various correlative approaches were attempted, with the roughness transition correlations based on edge properties providing the most reliable results. When a consistent computational method is used to compute edge conditions, transition datasets for different configurations at several angles of attack have been shown to collapse to a well-behaved correlation.

Berry, Scott A.↗