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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Acoustic wave propagation in heterogeneous structures including experimental validation

A finite element model was developed to solve for the acoustic pressure and energy fields in a heterogeneous suppressor. The derivations from the governing equations assumed that the material properties could vary with position resulting in a heterogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between different materials. For a two-media region consisting of part air and part bulk absorber, a model was used to describe the bulk absorber properties in two directions. Complex metallic structures inside the air duct are simulated by simply changing element properties from air to the structural material in a pattern to describe the desired shapes. To verify the numerical theory, experiments were conducted without flow in a rectangular duct with a single folded cavity mounted above the duct and absorbing material mounted inside a cavity. Changes in a nearly plane wave sound field were measured on the wall opposite the absorbing cavity. Fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber, as a function of distance along the duct. The finite element model provides a convenient method for evaluating the acoustic properties of bulk absorbers.

Baumeister, Kenneth J.↗
DOE OSTI2022

Advanced Characterization of Fuel-cladding Chemical Interaction between U-10Zr Fuel and HT9 Cladding Tested in Fast Flux Test Facility

Fuel cladding chemical interaction (FCCI) can greatly accelerate the cladding failure. However, due to the limited space in a fuel cladding assembly, it is historically challenging to gain an mechanistical understanding of the formation mechanism of FCCI and its influence on fuel and cladding performance. With the imminent need to qualify U-10Zr based metallic fuel cladded by HT-9 for advanced reactors demonstration project, it is of vital importance to use advanced characterization method to study FCCI in a unprecedent detailed manner and gain better mechanism understanding of FCCI. Mechanistic Fuel Failure (MFF) series of prototypic fuel elements irradiated in FFTF [1, 2] provides the best samples to study FCCI since the MFF-series assemblies had an axial fuel height the same as proposed length by industry partners. Jason et al. [3] has performed preliminary post-irradiation examination on a MFF fuel pin sample, which was extracted from the HT9 cladded U-10at.%Zr MFF-3 pin MFF-3 pin (#193045) at an axial location of X/L = 0.98. This sample has a peak burnup of 5.7 at.% and peak inner cladding temperature (PICT) of around 615 °C during in-core testing. Scanning electron microscope examination has identified visible FCCI region on more than half of the HT9 circumference [3]. The most striking feature are grain boundary attacking by apparently lanthanides (Lns) rich phase. However, SEM cannot provide accurate assessment of phase and concentration of grain boundary phases and prevented a better understanding of the formation mechanism of such attach. This study, by pairing transmission electron microscope (TEM) characterization and atom probe tomography (APT) techniques with in-situ micro-tensile testing in scanning electron microscope (SEM), aims at gaining in-depth understanding on the formed FCCI region. The identified FCCI region roughly consists of multilayers as illustrated in Figure 1 (c). The main findings are: (1) layer-B shows observable lanthanides (Lns) infiltration along grain boundaries and mechanical softening due to FCCI- and irradiation-induced microstructural and microchemistry changes, particularly the recovery of martensitic lath structure and dissolution of pre-existing M23C6 together with the formation of coarsened Laves phases, (Fe, Cr)2(Mo, W); (2) layer-C is Fe depleted but Lns significantly enriched, becoming very brittle; (3) layer-D is mainly composed of UFe2 and Lns; (4) three FCCI-induced intermetallic U-Fe-Zr phases, ? (Fe0.5Zr0.32U0.18), e (Fe0.3Zr0.4U0.3), ? (Fe0.06Zr0.23U0.71), were identified near layer-E; (5) the ? (Fe0.5Zr0.32U0.18) phase was characterized to be a face centered cubic (FCC) crystal structure. These results will help to better understanding the governing mechanism of FCCI and facilitating the development of theoretical model for assessing the performance of metallic fuel and cladding integrity.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Finite Element Modeling for Compression Strength Prediction of Honeycomb Cores with Geometric Imperfections Measured using X-ray CT Imaging

Aluminum honeycomb cores have been used extensively in composite sandwich panels due to having high bending rigidity while maintaining low density. During manufacturing of the aluminum honeycomb, the thin metal walls are susceptible to imperfections that deviate from an ideal honeycomb shape. Computational tools to quantify imperfections and investigate their effects on quasi-static compression and impact response can inform design criteria to construct safer and lighter launch vehicle structures. This paper presents the results of a study on finite element modeling of metallic honeycomb cores (HCC) with geometric imperfections measured using X-ray Computer Tomography (CT), to predict the compression response. Finite element (FE) models are constructed with measured imperfections with appropriate boundary conditions. It is shown that the average behavior of larger 10x11 cell models can be predicted by sampling single cell models from the larger domain. Amplification of imperfections for cells at the center of sampled specimens are mode switched from the larger out-of-plane imperfection to the mode shape of the adjacent cells resulting in strengthening of the center specimens.

aluminum honeycomb↗
NASA NTRS1975

Vertical velocity structure and geometry of clear air convective elements

The paper discusses observations of individual convective elements with a high-power narrow-beam scanning radar, an FM-CW radar, and an acoustic sounder, including the determination of the vertical air velocity patterns of convective structures with the FM-CW radar and acoustic sounder. Data are presented which link the observed velocity structure and geometrical patterns to previously proposed models of boundary layer convection. It is shown that the high-power radar provides a clear three-dimensional picture of convective cells and fields over a large area with a resolution of 150 m, where the convective cells are roughly spherical. Analysis of time-height records of the FM-CW radar and acoustic sounder confirms the downdraft-entrainment mechanism of the convective cell. The Doppler return of the acoustic sounder and the insect-trail slopes on FM-CW radar records are independent but redundant methods for obtaining the vertical velocity patterns of convective structures.

Rowland, J. R.↗
NASA NTRS2022

Multi-Scale Modelling of the Bound Metal Deposition Manufacturing of Ti6Al4V

Nonlinear shrinkage of the metal part during manufacturing by bound metal deposition, both on the ground and under microgravity, is considered. A multi-scale physics-based approach is developed to address the problem. It spans timescales from atomistic dynamics on the order of nanoseconds to full-part shrinkage on the order of hours. This approach enables estimation of the key parameters of the problem, including the widths of grain boundaries, the coefficient of surface diffusion, the initial redistribution of particles during the debinding stage, the evolution of the microstructure from round particles to densely-packed grains, the corresponding changes in the total and chemical free energies, and the sintering stress. The method has been used to predict shrinkage at the levels of two particles, of the filament cross-section, of the sub-model, and of the whole green, brown, and metal parts.

Nonlinear shrinkage↗
DOE OSTI2021

The influence of preexisting host rock damage on fault network localization

The transition from stable to unstable fracture propagation occurs when fractures begin to interact and link. Thus, fracture network coalescence controls how rocks and engineered structures fail. Here, to constrain the factors that influence localization in shear zones under brittle conditions, we build discrete element method models with a rough fault embedded in a shear zone. We add varying numbers of diffuse, randomly-placed weaknesses to examine the influence of diffuse damage on fracture network localization. The number of weaknesses controls the localization behavior of the fault network and the final fault geometry. We quantify localization using the Gini coefficient of the fracture volume, which measures the nonuniformity in a population. Each model generally increases in localization toward failure. However, models with more diffuse damage experience delocalization phases that are superimposed on the overall trend of increasing localization. The observed link between delocalization and host rock damage may help explain the varying localization of low magnitude seismicity in southern California. Models with more diffuse damage produce more complex fault geometries comprised of several parallel strands of wing cracks. The propagation of these wing cracks reduces the shear stress acting on the model boundaries, indicating that this fracture development increases the mechanical efficiency of the system.

58 GEOSCIENCES↗
NASA NTRS1996

Coupled Structural, Thermal, Phase-change and Electromagnetic Analysis for Superconductors, Volume 2

Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromag subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermel and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. Volume 1 describes mostly formulation specific problems. Volume 2 describes generalization of those formulations.

Felippa, C. A.↗
DOE OSTI2024

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗
NASA NTRS1996

Coupled Structural, Thermal, Phase-Change and Electromagnetic Analysis for Superconductors

Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermal and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. This volume, Volume 1, describes mostly formulations for specific problems. Volume 2 describes generalization of those formulations.

Felippa, C. A.↗
NASA NTRS2024

X-57 Cruise Motor GVT Using Fixed-Base Correction Technique

The National Aeronautics and Space Administration (NASA) Armstrong Flight Research Center (AFRC) completed a modal survey of the X-57 Maxwell aircraft cruise motor system to help inform cruise motor redesign efforts. X-57 Maxwell was an electric propulsion demonstrator aircraft developed by NASA to inform airworthiness standards for electrified-aircraft. The cruise motor system modal survey was completed in spring of 2023 utilizing the fixed-base correction (FBC) ground vibration test (GVT) technique developed by ATA Engineering, to decouple the motor modes from the aircraft modes. Previously during the full aircraft GVT, a detailed modal assessment of the cruise motors was not performed. Owing to the X-57 project’s phase in the aircraft development cycle when the motor redesign effort occurred, the cruise motor GVT could only be performed with the cruise motor system installed on the aircraft, with most of its installation hardware (wiring, baffling, sensors, etc.) attached. An impact hammer was used to provide excitation input at various locations within the tight confines of the cruise motor installation. To better support motor redesign efforts, the FBC methodology was utilized to fix, separate and de-couple the cruise motor modes from aircraft modal response. During the GVT, this required additional impact tap tests on candidate fixed-boundary points for each degree of freedom (DOF) to be fixed. Additional triaxial accelerometers installed at the candidate points were used to compute frequency response functions (FRFs) in X, Y, and Z directions to enable those DOFs to be numerically fixed. Test data was acquired using Hottinger Brüel & Kjær’s LAN-XI data acquisition hardware and BK Connect software. FBC post-test processing was performed using the Structural Modification Using Frequency Response Functions (SMURF) technique with ATA Engineering’s Interface between MATLAB, Analysis, Test (IMAT) software. Utilizing the FBC technique relieved test engineers from having to instrument the entire aircraft to identify and separate aircraft response from cruise motor modes of interest. The FBC technique also permitted structural analysis engineers to omit secondary components from their finite element model (FEM) of the cruise motor system. This FBC modal survey was successful, and the first time NASA AFRC utilized the FBC method on an aircraft rather than a test fixture, and also using an impact hammer rather than multiple shakers allowing significant project schedule and cost savings

Modal Survey↗
NASA NTRS2021

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗
DOE OSTI2024

Multi deep learning-based stochastic microstructure reconstruction and high-fidelity micromechanics simulation of time-dependent ceramic matrix composite response

A multi deep learning-based framework is developed for efficient, automated microstructure reconstruction and generation of stochastic representative volume elements (SRVEs) with periodic boundary conditions (PBCs) for accurate modeling of ceramic matrix composite (CMC) response. The methodology comprises a convolutional neural network coupled with regression layers to act as a vanilla regression network for semantic segmentation of the microstructure, allowing accurate characterization of the phases and their distributions at the microscale. Scanning electron microscope and confocal microscope are used to obtain C/SiNC and SiC/SiNC CMCs micrographs for vanilla regression testing. Microstructure variability in terms of fiber volume fraction and porosity are quantified through the output regression layer, ensuring accurate representation of material variability in SRVE construction. Generative adversarial network (GAN) and its variants are designed to produce high-fidelity SRVE, spanning CMCs microstructure variability space. A circular padding algorithm is developed to generate SRVEs with PBCs during training of GANs. The accuracy of the generated SRVEs is established through micromechanics simulations, where an efficient formulation of the high-fidelity generalized methods of cells (HFGMC) approach is used to compute the effective mechanical properties. Furthermore, an iterative algorithm is implemented in the HFGMC solver to simulate time-dependent deformation of SiC/SiNC subjected to creep loading conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiphase tin equation of state using density functional theory

In this work, we perform density functional theory (DFT) calculations of five solid phases and the liquid phase of tin. The calculations include cold curves of the five solid phases, phonon calculations in the quasiharmonic approximation over a range of volumes for each solid phase, and DFT-based molecular dynamics (DFT-MD) simulations of the liquid phase, including those of the melt curve using the Z method. Using the DFT results, we construct a tabular multiphase sesame equation of state for tin, referred to as sesame 2162. Comparisons to experimental data are made and show a high level of agreement in isobaric data, isothermal data, shock data, and phase boundary measurements, including measurements of the melt curve. The 2162 EOS will be useful for hydrodynamics simulations and has been designed with an eye toward hydrodynamics simulations that incorporate materials strength models and allow for modeling of the kinetics of phase transitions.

3-dimensional systems↗
NASA NTRS2014

Multiscale Analysis of Structurally-Graded Microstructures Using Molecular Dynamics, Discrete Dislocation Dynamics and Continuum Crystal Plasticity

A multiscale modeling methodology is developed for structurally-graded material microstructures. Molecular dynamic (MD) simulations are performed at the nanoscale to determine fundamental failure mechanisms and quantify material constitutive parameters. These parameters are used to calibrate material processes at the mesoscale using discrete dislocation dynamics (DD). Different grain boundary interactions with dislocations are analyzed using DD to predict grain-size dependent stress-strain behavior. These relationships are mapped into crystal plasticity (CP) parameters to develop a computationally efficient finite element-based DD/CP model for continuum-level simulations and complete the multiscale analysis by predicting the behavior of macroscopic physical specimens. The present analysis is focused on simulating the behavior of a graded microstructure in which grain sizes are on the order of nanometers in the exterior region and transition to larger, multi-micron size in the interior domain. This microstructural configuration has been shown to offer improved mechanical properties over homogeneous coarse-grained materials by increasing yield stress while maintaining ductility. Various mesoscopic polycrystal models of structurally-graded microstructures are generated, analyzed and used as a benchmark for comparison between multiscale DD/CP model and DD predictions. A final series of simulations utilize the DD/CP analysis method exclusively to study macroscopic models that cannot be analyzed by MD or DD methods alone due to the model size.

Saether, Erik↗
NASA NTRS2015

Investigation of Micro-Scale Architectural Effects on Damage of Composites

This paper presents a three-dimensional, energy based, anisotropic, stiffness reduction, progressive damage model for composite materials and composite material constituents. The model has been implemented as a user-defined constitutive model within the Abaqus finite element software package and applied to simulate the nonlinear behavior of a damaging epoxy matrix within a unidirectional composite material. Three different composite microstructures were considered as finite element repeating unit cells, with appropriate periodicity conditions applied at the boundaries. Results representing predicted transverse tensile, longitudinal shear, and transverse shear stress-strain curves are presented, along with plots of the local fields indicating the damage progression within the microstructure. It is demonstrated that the damage model functions appropriately at the matrix scale, enabling localization of the damage to simulate failure of the composite material. The influence of the repeating unit cell geometry and the effect of the directionality of the applied loading are investigated and discussed.

damage↗
DOE OSTI2021

Method for Rapid Modeling of Distortion in Laser Powder Bed Fusion Metal Additive Manufacturing Parts

The simulation and modeling of part-level distortion and residual stress in diverse metal additive manufacturing (AM) geometries has great potential to enable the rapid adoption of this technology in engineering design. Moreover, the use of additive manufacturing component libraries (CLs) offer a computationally efficient means of quantifying these part-level defects resultant from AM processing. We report on how the individual simulations of simple shapes, potential entries in a CL, can be superimposed to provide an indication of distortion and residual stresses in complex geometries. Laser powder bed fusion AM was used to construct test geometries of varied shapes and their combinations in the form of CLs in an effort to characterize location-dependent and feature-dependent distortion distributions. Blue light scanning was used to experimentally measure 3D distortions in order to investigate the interaction between the component shapes and local boundary conditions. Overall, part-level distortions were highly dependent on test component geometry, local boundary conditions, and shape combination. Commercial finite element software was used to verify experimental trends and to make predictions of distortion. The use of CLs resulted in over 20 times savings in computational cost while reproducing overall trends in distortion for test geometry assemblies. Therefore, it is anticipated that the use of CLs for L-PBF AM geometries has demonstrated potential to facilitate efficient simulations of full component AM assemblies, thereby reducing the need for costly trial-and-error-type experimental analysis.

36 MATERIALS SCIENCE↗
NASA NTRS1999

History of Ozone Research: From Schonbein to the Present

In 1840, C.F. Schonbein recognized that the smell generated in several different electrical and chemical processes was a single substance. He named this substance "ozein" from the Greek for "to smell". This substance we know today as ozone. Several periods can be distinguished in the continued development of our understanding of ozone. Throughout the late 19th century, the identity and properties of ozone were established and described. Ozone was recognized to be a constituent of normal air and tests were established to measure its concentration. Its disinfectant properties were recognized. New methods were developed for making ozone in the laboratory. In 1879, ultraviolet spectroscopic techniques were applied to the measurement of the solar spectrum and it was discovered by Comu that the spectrum was cut off at about 300 nm wavelength. Hartley suggested, based on laboratory measurements, that this cutoff was due to ozone in the atmosphere which he correctly asserted was somewhere in the upper atmosphere. This began the period of development of the amount and distribution of ozone throughout the atmosphere. In 1930, Chapman put forward the first theory of the formation and destruction of ozone. By the mid-1960s it was becoming obvious that the description of the chemical loss term was inadequate. By the early 1970s the chemical destruction of ozone by the oxides of hydrogen, nitrogen, chlorine, and bromine was recognized as an essential element in the chemical balance determining the ozone concentration. Today, ozone is a broad research project which crosses the boundaries of traditional disciplines. Stratospheric ozone loss due to chlorofluorocarbons is a newsworthy item. The Antarctic ozone hole opens up every spring. The provisions of the Montreal Protocol were agreed upon by countries around the world and promise to reduce the future levels of ozone-destroying chlorine in the stratosphere. Ozone concentrations in polluted cities are a subject of local and national regulations to limit its potential for causing health problems and corrosive effects on materials. A lot is known about ozone but many new questions are still being developed.

Stolarski, Richard S.↗
NASA NTRS1993

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 1

This research program has dealt with the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements; (2) finite element modeling of the electromagnetic problem; (3) coupling of thermal and mechanical effects; and (4) computer implementation and solution of the superconductivity transition problem. The research was carried out over the period September 1988 through March 1993. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles; (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements; and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects; and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The grant has fully supported the thesis work of one doctoral student (James Schuler, who started on January 1989 and completed on January 1993), and partly supported another thesis (Carmelo Militello, who started graduate work on January 1988 completing on August 1991). Twenty-three publications have acknowledged full or part support from this grant, with 16 having appeared in archival journals and 3 in edited books or proceedings.

Felippa, Carlos A.↗