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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Local and long-range atomic/magnetic structure of non-stoichiometric spinel iron oxide nanocrystallites

Spinel iron oxide nanoparticles of different mean sizes in the range 10–25 nm have been prepared by surfactant-free up-scalable near- and super-critical hydrothermal synthesis pathways and characterized using a wide range of advanced structural characterization methods to provide a highly detailed structural description. The atomic structure is examined by combined Rietveld analysis of synchrotron powder X-ray diffraction (PXRD) data and time-of-flight neutron powder-diffraction (NPD) data. The local atomic ordering is further analysed by pair distribution function (PDF) analysis of both X-ray and neutron total-scattering data. It is observed that a non-stoichiometric structural model based on a tetragonal γ-Fe 2 O 3 phase with vacancy ordering in the structure (space group P 4 3 2 1 2) yields the best fit to the PXRD and total-scattering data. Detailed peak-profile analysis reveals a shorter coherence length for the superstructure, which may be attributed to the vacancy-ordered domains being smaller than the size of the crystallites and/or the presence of anti-phase boundaries, faulting or other disorder effects. The intermediate stoichiometry between that of γ-Fe 2 O 3 and Fe 3 O 4 is confirmed by refinement of the Fe/O stoichiometry in the scattering data and quantitative analysis of Mössbauer spectra. The structural characterization is complemented by nano/micro-structural analysis using transmission electron microscopy (TEM), elemental mapping using scanning TEM, energy-dispersive X-ray spectroscopy and the measurement of macroscopic magnetic properties using vibrating sample magnetometry. Notably, no evidence is found of a Fe 3 O 4 /γ-Fe 2 O 3 core-shell nanostructure being present, which had previously been suggested for non-stoichiometric spinel iron oxide nanoparticles. Finally, the study is concluded using the magnetic PDF (mPDF) method to model the neutron total-scattering data and determine the local magnetic ordering and magnetic domain sizes in the iron oxide nanoparticles. The mPDF data analysis reveals ferrimagnetic collinear ordering of the spins in the structure and the magnetic domain sizes to be ∼60–70% of the total nanoparticle sizes. The present study is the first in which mPDF analysis has been applied to magnetic nanoparticles, establishing a successful precedent for future studies of magnetic nanoparticles using this technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Measurement of grain boundary strength of Inconel X-750 superalloy using in-situ micro-tensile testing techniques in FIB/SEM system

Grain boundaries (GBs), known as two-dimensional defects, are omnipresent in polycrystalline metallic alloys and thus influence a wide range of mechanical properties under different environmental conditions like irradiation and corrosion. Therefore, quantifying the strength of individual GBs is critical for understanding the degradation of mechanical properties of materials under different conditions. In this study we developed an efficient approach for the fabrication of micro-tensile specimens with a GB almost perpendicular to the tensile direction, which is expected to advance the development of individual GB tensile testing at micro or nanoscale in a wide scope of materials. Further, an in-situ cantilever micro-tensile testing method was developed and used to quantify the strength of a Σ3 GB in Inconel X-750 with the combination of finite element modeling. The average ultimate tensile strength (UTS) of a non-irradiated Σ3 GB is estimated at around 1.4 GPa, comparable to that of a neutron-irradiated Σ3 GB with a dose of ~1.5 dpa (1.3 GPa). Moreover, the in-situ push-to-pull micro-tensile testing technique developed in this work provides valuable insights into the high-angle GB deformation and fracture behavior. This method generates qualitatively similar ductility behavior before and after neutron irradiation as the bulk material testing. However, the ductility and UTS values obtained from this method are different from bulk measurements due to vastly different specimen dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Stable, Ductile and Strong Ultrafine HT-9 Steels via Large Strain Machining

Beyond the current commercial materials, refining the grain size is among the proposed strategies to manufacture resilient materials for industrial applications demanding high resistance to severe environments. Here, large strain machining (LSM) was used to manufacture nanostructured HT-9 steel with enhanced thermal stability, mechanical properties, and ductility. Nanocrystalline HT-9 steels with different aspect rations are achieved. In-situ transmission electron microscopy annealing experiments demonstrated that the nanocrystalline grains have excellent thermal stability up to 700 °C with no additional elemental segregation on the grain boundaries other than the initial carbides, attributing the thermal stability of the LSM materials to the low dislocation densities and strains in the final microstructure. Nano-indentation and micro-tensile testing performed on the LSM material pre- and post-annealing demonstrated the possibility of tuning the material’s strength and ductility. The results expound on the possibility of manufacturing controlled nanocrystalline materials via a scalable and cost-effective method, albeit with additional fundamental understanding of the resultant morphology dependence on the LSM conditions.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Data Systems Dynamic Simulator

The Data System Dynamic Simulator (DSDS) is a discrete event simulation tool. It was developed for NASA for the specific purpose of evaluating candidate architectures for data systems of the Space Station era. DSDS provides three methods for meeting this requirement. First, the user has access to a library of standard pre-programmed elements. These elements represent tailorable components of NASA data systems and can be connected in any logical manner. Secondly, DSDS supports the development of additional elements. This allows the more sophisticated DSDS user the option of extending the standard element set. Thirdly, DSDS supports the use of data streams simulation. Data streams is the name given to a technique that ignores packet boundaries, but is sensitive to rate changes. Because rate changes are rare compared to packet arrivals in a typical NASA data system, data stream simulations require a fraction of the CPU run time. Additionally, the data stream technique is considerably more accurate than another commonly-used optimization technique.

Rouff, Christopher↗
NASA NTRS1994

Non-Reflecting Regions for Finite Difference Methods in Modeling of Elastic Wave Propagation in Plates

Solution of the wave equation using techniques such as finite difference or finite element methods can model elastic wave propagation in solids. This requires mapping the physical geometry into a computational domain whose size is governed by the size of the physical domain of interest and by the required resolution. This computational domain, in turn, dictates the computer memory requirements as well as the calculation time. Quite often, the physical region of interest is only a part of the whole physical body, and does not necessarily include all the physical boundaries. Reduction of the calculation domain requires positioning an artificial boundary or region where a physical boundary does not exist. It is important however that such a boundary, or region, will not affect the internal domain, i.e., it should not cause reflections that propagate back into the material. This paper concentrates on the issue of constructing such a boundary region.

Kishoni, Doron↗
NASA NTRS1959

Source Distribution Method for Unsteady One-Dimensional Flows With Small Mass, Momentum, and Heat Addition and Small Area Variation

A source distribution method is presented for obtaining flow perturbations due to small unsteady area variations, mass, momentum, and heat additions in a basic uniform (or piecewise uniform) one-dimensional flow. First, the perturbations due to an elemental area variation, mass, momentum, and heat addition are found. The general solution is then represented by a spatial and temporal distribution of these elemental (source) solutions. Emphasis is placed on discussing the physical nature of the flow phenomena. The method is illustrated by several examples. These include the determination of perturbations in basic flows consisting of (1) a shock propagating through a nonuniform tube, (2) a constant-velocity piston driving a shock, (3) ideal shock-tube flows, and (4) deflagrations initiated at a closed end. The method is particularly applicable for finding the perturbations due to relatively thin wall boundary layers.

Mirels, Harold↗
NASA NTRS2007

One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure

The development and verification of a one-dimensional material thermal response code with ablation is presented. The implicit time integrator, control volume finite element spatial discretization, and Newton's method for nonlinear iteration on the entire system of residual equations have been implemented and verified for the thermochemical ablation of internally decomposing materials. This study is a continuation of the work presented in "One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure" (AIAA-2006-2910), which described the derivation, implementation, and verification of the constant density solid energy equation terms and boundary conditions. The present study extends the model to decomposing materials including decomposition kinetics, pyrolysis gas flow through the porous char layer, and a mixture (solid and gas) energy equation. Verification results are presented for the thermochemical ablation of a carbon-phenolic ablator which involves the solution of the entire system of governing equations.

Amar, Adam J.↗
NASA NTRS1974

Integral representation of field variables for the finite element solution of viscous flow problems

The kinematics aspects of compressible and incompressible viscous flow problems are recasted into an integral representation for the velocity vector. For the incompressible flow problem, the kinetics aspect is similarly recasted into an integral representation. The result is a system of integral equations ideally suited for the finite element method. The integral representations are shown to permit the confinement of the computation field to the region of nonnegligible vorticity and dilatation which, for the incompressible flow, is identified with the viscous region. This drastic reduction in computation field results in a superior solution speed which is further improved by the use of a flowfield segmentation technique. The new approach is further shown to remove the difficulties associated with previous methods in specifying the far field and extraneous boundary conditions.

Wu, J. C.↗
NASA NTRS1992

Hybrid analysis of deformation using interferometric measurements

Two hybrid methods of deformation analysis combining moire interferometry and finite element techniques are introduced. The hybrid methods are applied to determine the strain fields in the test section of a composite cruciform specimen in order to evaluate the uniformity of the shear strain field. Both methods show that there are significant shear strain gradients in the specimen test section. One of the hybrid methods - the full-field analysis using the random-node mesh generator - detects the presence of material property nonuniformity which is assumed uniform in the boundary-node method.

Lin, S.-Y.↗
NASA NTRS1983

Propeller noise prediction

Analytic propeller noise prediction involves a sequence of computations culminating in the application of acoustic equations. The prediction sequence currently used by NASA in its ANOPP (aircraft noise prediction) program is described. The elements of the sequence are called program modules. The first group of modules analyzes the propeller geometry, the aerodynamics, including both potential and boundary layer flow, the propeller performance, and the surface loading distribution. This group of modules is based entirely on aerodynamic strip theory. The next group of modules deals with the actual noise prediction, based on data from the first group. Deterministic predictions of periodic thickness and loading noise are made using Farassat's time-domain methods. Broadband noise is predicted by the semi-empirical Schlinker-Amiet method. Near-field predictions of fuselage surface pressures include the effects of boundary layer refraction and (for a cylinder) scattering. Far-field predictions include atmospheric and ground effects. Experimental data from subsonic and transonic propellers are compared and NASA's future direction is propeller noise technology development are indicated.

Zorumski, W. E.↗
DOE OSTI2024

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method was developed that could predict the pressure distribution of a ducted source from far field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) was solved using explicit marching techniques. The Von Neumann method was used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach was successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1978

Numerical spatial marching techniques for estimating duct attenuation and source pressure profiles

A numerical method is developed that could predict the pressure distribution of a ducted source from far-field pressure inputs. Using an initial value formulation, the two-dimensional homogeneous Helmholtz wave equation (no steady flow) is solved using explicit marching techniques. The Von Neumann method is used to develop relationships which describe how sound frequency and grid spacing effect numerical stability. At the present time, stability considerations limit the approach to high frequency sound. Sample calculations for both hard and soft wall ducts compare favorably to known boundary value solutions. In addition, assuming that reflections in the duct are small, this initial value approach is successfully used to determine the attenuation of a straight soft wall duct. Compared to conventional finite difference or finite element boundary value approaches, the numerical marching technique is orders of magnitude shorter in computation time and required computer storage and can be easily employed in problems involving high frequency sound.

Baumeister, K. J.↗
NASA NTRS1979

Optimization of multi-element airfoils for maximum lift

Two theoretical methods are presented for optimizing multi-element airfoils to obtain maximum lift. The analyses assume that the shapes of the various high lift elements are fixed. The objective of the design procedures is then to determine the optimum location and/or deflection of the leading and trailing edge devices. The first analysis determines the optimum horizontal and vertical location and the deflection of a leading edge slat. The structure of the flow field is calculated by iteratively coupling potential flow and boundary layer analysis. This design procedure does not require that flow separation effects be modeled. The second analysis determines the slat and flap deflection required to maximize the lift of a three element airfoil. This approach requires that the effects of flow separation from one or more of the airfoil elements be taken into account. The theoretical results are in good agreement with results of a wind tunnel test used to corroborate the predicted optimum slat and flap positions.

Lawrence E Olson↗
NASA NTRS1980

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1)the elimination of large matrix storage requirements, and (2)shorter solution times under most conditions.

Baumeister, K. J.↗
NASA NTRS1980

Time dependent difference theory for sound propagation in axisymmetric ducts with plug flow

The time dependent governing/acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard wall duct with a plug mean flow and spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solutions of the acoustic field. The time dependent finite difference analysis has two advantages over the steady state finite difference and finite element techniques: (1) the elimination of large matrix storage requirements, and (2) shorter solution times under most conditions.

Baumeister, K. J.↗
NASA NTRS2010

The ICARE Method

The ICARE method is a flexible, widely applicable method for systems engineers to solve problems and resolve issues in a complete and comprehensive manner. The method can be tailored by diverse users for direct application to their function (e.g. system integrators, design engineers, technical discipline leads, analysts, etc.). The clever acronym, ICARE, instills the attitude of accountability, safety, technical rigor and engagement in the problem resolution: Identify, Communicate, Assess, Report, Execute (ICARE). This method was developed through observation of Space Shuttle Propulsion Systems Engineering and Integration (PSE&I) office personnel approach in an attempt to succinctly describe the actions of an effective systems engineer. Additionally it evolved from an effort to make a broadly-defined checklist for a PSE&I worker to perform their responsibilities in an iterative and recursive manner. The National Aeronautics and Space Administration (NASA) Systems Engineering Handbook states, engineering of NASA systems requires a systematic and disciplined set of processes that are applied recursively and iteratively for the design, development, operation, maintenance, and closeout of systems throughout the life cycle of the programs and projects. ICARE is a method that can be applied within the boundaries and requirements of NASA s systems engineering set of processes to provide an elevated sense of duty and responsibility to crew and vehicle safety. The importance of a disciplined set of processes and a safety-conscious mindset increases with the complexity of the system. Moreover, the larger the system and the larger the workforce, the more important it is to encourage the usage of the ICARE method as widely as possible. According to the NASA Systems Engineering Handbook, elements of a system can include people, hardware, software, facilities, policies and documents; all things required to produce system-level results, qualities, properties, characteristics, functions, behavior and performance. The ICARE method can be used to improve all elements of a system and, consequently, the system-level functional, physical and operational performance. Even though ICARE was specifically designed for a systems engineer, any person whose job is to examine another person, product, or process can use the ICARE method to improve effectiveness, implementation, usefulness, value, capability, efficiency, integration, design, and/or marketability. This paper provides the details of the ICARE method, emphasizing the method s application to systems engineering. In addition, a sample of other, non-systems engineering applications are briefly discussed to demonstrate how ICARE can be tailored to a variety of diverse jobs (from project management to parenting).

Henke, Luke↗
NASA NTRS2007

A Semi-Analytical Method for Determining the Energy Release Rate of Cracks in Adhesively-Bonded Single-Lap Composite Joints

A semi-analytical method for determining the strain energy release rate due to a prescribed interface crack in an adhesively-bonded, single-lap composite joint subjected to axial tension is presented. The field equations in terms of displacements within the joint are formulated by using first-order shear deformable, laminated plate theory together with kinematic relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. Based on the adhesive stress distributions, the forces at the crack tip are obtained and the strain energy release rate of the crack is determined by using the virtual crack closure technique (VCCT). Additionally, the test specimen geometry from both the ASTM D3165 and D1002 test standards are utilized during the derivation of the field equations in order to correlate analytical models with future test results. The system of second-order differential field equations is solved to provide the adherend and adhesive stress response using the symbolic computation tool, Maple 9. Finite element analyses using J-integral as well as VCCT were performed to verify the developed analytical model. The finite element analyses were conducted using the commercial finite element analysis software ABAQUS. The results determined using the analytical method correlated well with the results from the finite element analyses.

Yang, Charles↗