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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Development and verification of local/global analysis techniques for laminated composites

Analysis and design methods for laminated composite materials have been the subject of considerable research over the past 20 years, and are currently well developed. In performing the detailed three-dimensional analyses which are often required in proximity to discontinuities, however, analysts often encounter difficulties due to large models. Even with the current availability of powerful computers, models which are too large to run, either from a resource or time standpoint, are often required. There are several approaches which can permit such analyses, including substructuring, use of superelements or transition elements, and the global/local approach. This effort is based on the so-called zoom technique to global/local analysis, where a global analysis is run, with the results of that analysis applied to a smaller region as boundary conditions, in as many iterations as is required to attain an analysis of the desired region. Before beginning the global/local analyses, it was necessary to evaluate the accuracy of the three-dimensional elements currently implemented in the Computational Structural Mechanics (CSM) Testbed. It was also desired to install, using the Experimental Element Capability, a number of displacement formulation elements which have well known behavior when used for analysis of laminated composites.

Griffin, O. Hayden, Jr.↗

NASA NTRS1989

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

NASA NTRS1990

On estimating scale invariance in stratocumulus cloud fields

Examination of cloud radiance fields derived from satellite observations sometimes indicates the existence of a range of scales over which the statistics of the field are scale invariant. Many methods were developed to quantify this scaling behavior in geophysics. The usefulness of such techniques depends both on the physics of the process being robust over a wide range of scales and on the availability of high resolution, low noise observations over these scales. These techniques (area perimeter relation, distribution of areas, estimation of the capacity, d0, through box counting, correlation exponent) are applied to the high resolution satellite data taken during the FIRE experiment and provides initial estimates of the quality of data required by analyzing simple sets. The results of the observed fields are contrasted with those of images of objects with known characteristics (e.g., dimension) where the details of the constructed image simulate current observational limits. Throughout when cloud elements and cloud boundaries are mentioned; it should be clearly understood that by this structures in the radiance field are meant: all the boundaries considered are defined by simple threshold arguments.

Seze, Genevieve↗

NASA NTRS1992

Inverse design of a proper number, shapes, sizes, and locations of coolant flow passages

During the past several years we have developed an inverse method that allows a thermal cooling system designer to determine proper sizes, shapes, and locations of coolant passages (holes) in, say, an internally cooled turbine blade, a scram jet strut, a rocket chamber wall, etc. Using this method the designer can enforce a desired heat flux distribution on the hot outer surface of the object, while simultaneously enforcing desired temperature distributions on the same hot outer surface as well as on the cooled interior surfaces of each of the coolant passages. This constitutes an over-specified problem which is solved by allowing the number, sizes, locations and shapes of the holes to adjust iteratively until the final internally cooled configuration satisfies the over-specified surface thermal conditions and the governing equation for the steady temperature field. The problem is solved by minimizing an error function expressing the difference between the specified and the computed hot surface heat fluxes. The temperature field analysis was performed using our highly accurate boundary integral element code with linearly varying temperature along straight surface panels. Examples of the inverse design applied to internally cooled turbine blades and scram jet struts (coated and non-coated) having circular and non-circular coolant flow passages will be shown.

George S. Dulikravich↗

DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗

NASA NTRS2023

BOLT II Forebody Surface Measurement Analysis

Surface impression molds of the nosetip/forebody interface, discrete tripping elements, and other surface imperfections were formed on the BOLT II flight vehicle prior to its launch in March 2022. The silicon-rubber molds cast on the interface between the nickel nosetip and the stainless steel forebody were measured using a 3D optical measurement system and a stylus profilometer and were evaluated for step direction and height. The step around the nosetip/forebody interface was found to vary in both magnitude and direction, but values were generally modest with most samples below 60 μm across the span and only two samples, both cast in the gutter, with step values larger than 100 μm. Impressions of the secondary-side tripping elements were strictly scanned with the optical device and similar methods were used to calculate trip height as step height. The measured centerline trip height was 2623 ±12 μm, while the leading-edge, quadrant 2 and 3, trip heights were 897 ±5 μm and 891 ±18 μm, respectively. Surface imperfections on the nose and forebody were also evaluated. Surface roughness values and spectra were estimated, and profiles were extracted across the imperfections to quantify the magnitude of those features.

Trip height↗

DOE OSTI2025

Simulated effect of defect volume and location on very high cycle fatigue of laser beam powder bed fused AlSi10Mg

This study quantifies the interaction between volumetric defect location and size on the very high cycle fatigue (VHCF) of laser beam powder bed fused (LB-PBF) AlSi10Mg. Crystal plasticity finite element method (CPFEM) simulations were used to investigate the effects of defect location and size on the driving force for crack initiation. The CPFEM model was calibrated against uniaxial and cyclic experimental data of LB-PBF AlSi10Mg. Defect characteristics were informed by experimental data from the specimens produced in various geometries to create realistic representative volume elements (RVEs) with equivalent volume fractions of defects. By embedding defects of varying sizes and locations within the RVEs, fatigue indicator parameters (FIPs) were calculated to analyze the impact of defects’ characteristics on fatigue performance. Different combinations of defect volume and locations were generated for various microstructure instantiations, providing insight into extreme value fatigue responses. Larger defect volumes located on free surfaces consistently generated the highest FIPs, suggesting defect size and boundary proximity intensify stress concentration effects. RVEs with multiple smaller defects produced lower FIPs than those with single large critical defects. These findings underscore the critical role of defect characteristics on fatigue life, providing a foundation for future predictive modeling in fatigue-sensitive AM applications.

AlSi10Mg↗

DOE OSTI2021

PRISMS-Fatigue computational framework for fatigue analysis in polycrystalline metals and alloys

Abstract The PRISMS-Fatigue open-source framework for simulation-based analysis of microstructural influences on fatigue resistance for polycrystalline metals and alloys is presented here. The framework uses the crystal plasticity finite element method as its microstructure analysis tool and provides a highly efficient, scalable, flexible, and easy-to-use ICME community platform. The PRISMS-Fatigue framework is linked to different open-source software to instantiate microstructures, compute the material response, and assess fatigue indicator parameters. The performance of PRISMS-Fatigue is benchmarked against a similar framework implemented using ABAQUS. Results indicate that the multilevel parallelism scheme of PRISMS-Fatigue is more efficient and scalable than ABAQUS for large-scale fatigue simulations. The performance and flexibility of this framework is demonstrated with various examples that assess the driving force for fatigue crack formation of microstructures with different crystallographic textures, grain morphologies, and grain numbers, and under different multiaxial strain states, strain magnitudes, and boundary conditions.

Chemistry↗

NASA NTRS1989

Multigrid solution of the Navier-Stokes equations on triangular meshes

A Navier-Stokes algorithm for use on unstructured triangular meshes is presented. Spatial discretization of the governing equations is achieved using a finite element Galerkin approximation, which can be shown to be equivalent to a finite volume approximation for regular equilateral triangular meshes. Integration steady-state is performed using a multistage time-stepping scheme, and convergence is accelerated by means of implicit residual smoothing and an unstructured multigrid algorithm. Directional scaling of the artificial dissipation and the implicit residual smoothing operator is achieved for unstructured meshes by considering local mesh stretching vectors at each point. The accuracy of the scheme for highly stretched triangular meshes is validated by comparing computed flat-plate laminar boundary layer results with the well known similarity solution, and by comparing laminar airfoil results with those obtained from various well-established structured quadrilateral-mesh codes. The convergence efficiency of the present method is also shown to be competitive with those demonstrated by structured quadrilateral-mesh algorithms.

Mavriplis, Dimitri J.↗

NASA NTRS1990

Multigrid solution of the Navier-Stokes equations on triangular meshes

Mavriplis, D. J.↗

NASA NTRS2018

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

DOE OSTI2021

Role of Pairwise Reactions on the Synthesis of Li 0.3 La 0.57 TiO 3 and the Resulting Structure–Property Correlations

The performance of single-ion conductors is highly sensitive to the material’s defect chemistry. Tuning these defects is limited for solid-state reactions as they occur at particle–particle interfaces, which provide a complex evolving energy landscape for atomic rearrangement and product formation. In this report, we investigate the (1) order of addition and (2) lithium precursor decomposition temperature and their effect on the synthesis and grain boundary conductivity of the perovskite lithium lanthanum titanium oxide (LLTO). We use an intimately mixed sol–gel, a solid-state reaction of Li precursor + La 2 O 3 + TiO 2 , and Li precursor + amorphous La 0.57 TiO x as different chemical routes to change the way in which the elements are brought together. The results show that the perovskite can accommodate a wide range of Li deficiencies (upward of 50%) while maintaining the tetragonal LLTO structure, indicating that X-ray diffraction (XRD) is insufficient to fully characterize the chemical nature of the product (i.e., Li-deficient LLTO may behave differently than stoichiometric LLTO). Variations in the relative intensities of different reflections in XRD suggest variations in the La ordering within the crystal structure between synthesis methods. Furthermore, the choice of the precursor and the order of addition of the reactants lower the time required to form a pure phase. Density functional theory calculations of the formation energy of possible reaction intermediates support the hypothesis that a greater thermodynamic driving force to form LLTO leads to a greater LLTO yield. Here, the retention of lithium is correlated with the thermal decomposition temperature of the Li precursor and the starting material mixing strategy. Taking the results together suggests that cations that share a site with Li should be mixed early to avoid ordering. Such cation ordering inhibits Li motion, leading to higher Li ion resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2002

Aerodynamic Design Study of Advanced Multistage Axial Compressor

As a direct response to the need for further performance gains from current multistage axial compressors, an investigation of advanced aerodynamic design concepts that will lead to compact, high-efficiency, and wide-operability configurations is being pursued. Part I of this report describes the projected level of technical advancement relative to the state of the art and quantifies it in terms of basic aerodynamic technology elements of current design systems. A rational enhancement of these elements is shown to lead to a substantial expansion of the design and operability space. Aerodynamic design considerations for a four-stage core compressor intended to serve as a vehicle to develop, integrate, and demonstrate aerotechnology advancements are discussed. This design is biased toward high efficiency at high loading. Three-dimensional blading and spanwise tailoring of vector diagrams guided by computational fluid dynamics (CFD) are used to manage the aerodynamics of the high-loaded endwall regions. Certain deleterious flow features, such as leakage-vortex-dominated endwall flow and strong shock-boundary-layer interactions, were identified and targeted for improvement. However, the preliminary results were encouraging and the front two stages were extracted for further aerodynamic trimming using a three-dimensional inverse design method described in part II of this report. The benefits of the inverse design method are illustrated by developing an appropriate pressure-loading strategy for transonic blading and applying it to reblade the rotors in the front two stages of the four-stage configuration. Multistage CFD simulations based on the average passage formulation indicated an overall efficiency potential far exceeding current practice for the front two stages. Results of the CFD simulation at the aerodynamic design point are interrogated to identify areas requiring additional development. In spite of the significantly higher aerodynamic loadings, advanced CFD-based tools were able to effectively guide the design of a very efficient axial compressor under state-of-the-art aeromechanical constraints.

Larosiliere, Louis M.↗

NASA NTRS2004

Computation of Gust-Cascade Interaction Using the CE/SE Method

The problem 2 in Category 3 of the 4th Computational Aeroacoustic(CAA) Workshop is solved using the space-time conservation element and solution element (CE/SE) method. This problem models rotor-stator interaction in a 2D cascade. It involves complex geometries and flow physics including vortex shedding and acoustic radiation. The parallel version of the 2D nonlinear Euler solver is used with an unstructured triangular mesh to solve this problem. The Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Validation on the Giles approach is performed using Problem 3.1 in the 2nd CAA Workshop. The space-time CE/SE method is a finite volume method with second-order accuracy in both space and time. The flux conservation is enforced in both space and time instead of space only. It has low numerical dissipation and dispersion errors. It uses simple non-reflecting boundary conditions and is compatible with unstructured meshes. It is simple, flexible, and generate reasonably accurate solutions. The CE/SE method has been successfully applied to solve numerous practical problems, especially aeroacoustic problems. Some preliminary numerical results of the benchmark problem 3.2 of the 4th CAA Workshop are shown. The steady-state pressure contour is plotted. The mean pressure distribution on the blade surface is compared with Turbo solution showing a good agreement. The sound pressure level versus the rotor harmonic n at the six designated positions on the blade surface, three locations at inlet plane, and three locations at the outlet plane are plotted. It can be seen that the acoustic response exists only at the excitation frequencies (n = 1,2,3). On the blade surface, the acoustic wave at n = 1 is dominant, while at the inlet and outlet planes, the sound pressure level at n = 2 becomes the largest, which is similar to the results presented. The distribution of sound pressure level at different spatial modes along the z- direction is plotted for n = 1,2,3, respectively. It shows that the spatial modes m = -32 and 22 at n = 1 exponentially decay, and the spatial modes m = 10 at n = 2, m = -42 and 12 at n = 3 propagate both upstream and downstream, which agrees with the prediction based on the linearized theory. Some oscillations are observed, which needs to be investigated further. In the final paper, the numerical results will be compared with a frequency-domain solver LINFLUX solution if it is available.

Wang, X.-Y.↗

DOE OSTI2021

Synthesis and thermodynamics of uranium-incorporated α-Fe 2 O 3 nanoparticles

Hematite nanoparticles were synthesized with U(VI) in circumneutral water through a coprecipitation and hydrothermal treatment process. XRD, TEM, and EXAFS analyses reveal that uranium may aggregate along grain boundaries and occupy Fe sites within hematite. The described synthesis method produces crystalline, single-phase iron oxide nanoparticles absent of surface-bound uranyl complexes. EXAFS data were comparable to spectra from existing studies whose syntheses were more representative of naturally occurring, extended aging processes. Herein this work provides and validates an accelerated method of synthesizing uranium-immobilized iron oxide nanoparticles for further mechanistic studies. High temperature oxide melt solution calorimetry measurements were performed to calculate the thermodynamic stability of uranium-incorporated iron oxide nanoparticles. Increasing uranium content within hematite resulted in more positive formation enthalpies. Standard formation enthalpies of U x Fe 2–2x O 3 were as high as 76.88 ± 2.83 kJ/mol relative to their binary oxides, or -764.04 ± 3.74 kJ/mol relative to their constituent elements, at x = 0.037. Data on the thermodynamic stability of uranium retention pathways may assist in predicting waste uranyl remobilization, as well as in developing more effective methods to retain uranium captured from aqueous environments.

36 MATERIALS SCIENCE↗

DOE OSTI2020

High-Energy X-Ray Diffraction Microscopy in Materials Science

High-energy diffraction microscopy (HEDM) is an implementation of three-dimensional X-ray diffraction microscopy. HEDM yields maps of internal crystal orientation fields, strain states, grain shapes and locations as well as intragranular orientation distributions, and grain boundary character. Because it is nondestructive in hard materials, notably metals and ceramics, HEDM has been used to study responses of these materials to external fields including high temperature and mechanical loading. Currently available sources and detectors lead to a spatial resolution of ~1 μm and an orientation resolution of <0.1°. With the penetration characteristic of high energies ( E ≥ 50 keV), sample cross-section dimensions of ~1 mm can be studied in materials containing elements across much of the Periodic Table. This review describes hardware and software associated with HEDM as well as examples of applications. Overall, these applications include studies of grain growth, recrystallization, texture development, orientation gradients, deformation twinning, annealing twinning, plastic deformation, and additive manufacturing. We also describe relationships to other X-ray-based methods as well as prospects for further development.

36 MATERIALS SCIENCE↗

NASA NTRS2018

Comparing Free-Free and Shaker Table Model Correlation Methods Using Jim Beam

Finite element model correlation as part of a spacecraft program has always been a challenge. For any NASA mission, the coupled system response of the spacecraft and launch vehicle can be determined analytically through a Coupled Loads Analysis (CLA), as it is not possible to test the spacecraft and launch vehicle coupled system before launch. The value of the CLA is highly dependent on the accuracy of the frequencies and mode shapes extracted from the spacecraft model. NASA standards require the spacecraft model used in the final Verification Loads Cycle to be correlated by either a modal test or by comparison of the model with Frequency Response Functions (FRFs) obtained during the environmental qualification test. Due to budgetary and time constraints, most programs opt to correlate the spacecraft dynamic model during the environmental qualification test, conducted on a large shaker table. For any model correlation effort, the key has always been finding a proper definition of the boundary conditions. This paper is a correlation case study to investigate the difference in responses of a simple structure using a free-free boundary, a fixed boundary on the shaker table, and a base-drive vibration test, all using identical instrumentation. The NAVCON Jim Beam test structure, featured in the IMAC round robin modal test of 2009, was selected as a simple, well recognized and well characterized structure to conduct this investigation. First, a free-free impact modal test of the Jim Beam was done as an experimental control. Second, the Jim Beam was mounted to a large 20,000 lbf shaker, and an impact modal test in this fixed configuration was conducted. Lastly, a vibration test of the Jim Beam was conducted on the shaker table. The free-free impact test, the fixed impact test, and the base-drive test were used to assess the effect of the shaker modes, evaluate the validity of fixed-base modeling assumptions, and compare final model correlation results between these boundary conditions.

Modal Test↗