SEARCH · Engineering Papers

Results for “Geometric discontinuities”

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 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Multi-Dimensional Spectral Difference Method for Unstructured Grids

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. It combines the best features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. Universal reconstructions are obtained by distributing unknowns in a geometrically similar manner for all unstructured cells. Placements of the unknown and flux points with various order of accuracy are given for the line, triangular and tetrahedral elements. The data structure of the new method permits an optimum use of cache memory, resulting in further computational efficiency on modern computers. A new pointer system is developed that reduces memory requirements and simplifies programming for any order of accuracy. Numerical solutions are presented and compared with the exact solutions for wave propagation problems in both two and three dimensions to demonstrate the capability of the method. Excellent agreement has been found. The method is simpler and more efficient than previous discontinuous Galerkin and spectral volume methods for unstructured grids.

Liu, Yen↗

NASA NTRS1996

Systematic Underestimation of Earthquake Magnitudes from Large Intracontinental Reverse Faults: Historical Ruptures Break Across Segment Boundaries

Because most large-magnitude earthquakes along reverse faults have such irregular and complicated rupture patterns, reverse-fault segments defined on the basis of geometry alone may not be very useful for estimating sizes of future seismic sources. Most modern large ruptures of historical earthquakes generated by intracontinental reverse faults have involved geometrically complex rupture patterns. Ruptures across surficial discontinuities and complexities such as stepovers and cross-faults are common. Specifically, segment boundaries defined on the basis of discontinuities in surficial fault traces, pronounced changes in the geomorphology along strike, or the intersection of active faults commonly have not proven to be major impediments to rupture. Assuming that the seismic rupture will initiate and terminate at adjacent major geometric irregularities will commonly lead to underestimation of magnitudes of future large earthquakes.

Rubin, C. M.↗

NASA NTRS1993

(abstract) The SIR-C/X-SAR Mission

The Shuttle Imaging Radar-C and X-Band Synthetic Aperture Radar (SIR-C/X-SAR) is a cooperative experiment between the National Aeronautics and Space Administration (NASA), the German Space Agency (DARA), and the Italian Space Agency (ASI). The experiment is the next evolutionary step in NASA's Spaceborne Imaging Radar (SIR) program that began with the Seasat SAR in 1978, and continued with SAR -A in 1981 and SAR-B in 1984. The program will eventually lead to the Earth Observing System (EOS) SAR later in this decade. SIR-C will provide increased capability over Seasat, SIR-A, and SIR-B by acquiring polarimetric images simultaneously at two microwave frequencies: L-band (wavelength 24 cm) and C-band (wavelength 6 cm). X-SAR will operate at X-band (wavelength 3 cm) with VV polarization, resulting in a three-frequency capability. Because radar backscatter is most strongly influenced by objects comparable in size to the radar wavelength, this multifrequency capability will provide information about the Earth's surface over a wide range of scales not discernable with single-frequency experiments. The polarimetric data will yield quantitative information about the surface geometric structure, vegetation dielectric properties, and surface discontinuities. The first SIR-C/X-SAR flight is planned for early in 1994, with two subsequent flights in different seasons to enable data aquisition under different environmental conditions. SIR-C/X-SAR will collect data for a period of 6 to 9 days during each flight. The instrument will be flown in a 57(deg) inclination orbit with a nomimal orbit altitude of 215 km. This altitude will result in a slightly drifting orbit. The latest status of the SIR-C/X-SAR mission and instrument design, science activities, and expected results will be presented.

SIR-C/XSAR radar backscatter surface geometric str↗

NASA NTRS1992

A curvilinear, anisotropic, p-version, brick finite element based on geometric entities

A 'brick' solid finite element is presently developed on the basis of the p-version analysis, and used to demonstrate the FEM concept of 'geometric entities'. This method eliminates interelement discontinuities between low- and high-order elements, allowing very fine control over the shape-function order in various parts of the model. Attention is given to the illustrative cases of a one-element model of an elliptic pipe, and a square cross-section cantilevered beam.

Hinnant, Howard E.↗

DOE OSTI2022

Surface stress tensor and junction conditions on a rotating null horizon

The general form of the surface stress tensor of an infinitesimally thin shell located on a rotating null horizon is derived, when different interior and exterior geometries are joined there. Even though the induced metric on the surface must be the same approached from either side, the first derivatives of the metric need not be. Such discontinuities lead to a Dirac δ-distribution in the Einstein tensor localized on the horizon. For a general stationary axisymmetric geometry the surface stress tensor can be expressed in terms of two geometric invariants that characterize the surface, namely the discontinuities [κ] and [J] of the surface gravity κ and angular momentum density J. The Komar energy and angular momentum are given in coordinates adapted to the Killing symmetries, and the surface contributions to each determined in terms of [κ] and [J] . Guided by these, a simple modification of the original Israel junction conditions is verified directly from the Einstein tensor density to give the correct finite result for the surface stress, when the normal n to the surface is allowed to tend continuously to a null vector. The relation to Israel’s original junction conditions, which fail on null surfaces, is given. The modified junction conditions are suitable to the matching of a rotating “black hole” exterior to any interior geometry joined at the Kerr null horizon surface, even when the surface normal is itself discontinuous and the Barrabès-Israel formalism is also inapplicable. This joining on a rotating null horizon is purely of the matter shell type and does not contain a propagating gravitational shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2005

Numerical Evaluation of P-Multigrid Method for the Solution of Discontinuous Galerkin Discretizations of Diffusive Equations

This paper describes numerical experiments with P-multigrid to corroborate analysis, validate the present implementation, and to examine issues that arise in the implementations of the various combinations of relaxation schemes, discretizations and P-multigrid methods. The two approaches to implement P-multigrid presented here are equivalent for most high-order discretization methods such as spectral element, SUPG, and discontinuous Galerkin applied to advection; however it is discovered that the approach that mimics the common geometric multigrid implementation is less robust, and frequently unstable when applied to discontinuous Galerkin discretizations of di usion. Gauss-Seidel relaxation converges 40% faster than block Jacobi, as predicted by analysis; however, the implementation of Gauss-Seidel is considerably more expensive that one would expect because gradients in most neighboring elements must be updated. A compromise quasi Gauss-Seidel relaxation method that evaluates the gradient in each element twice per iteration converges at rates similar to those predicted for true Gauss-Seidel.

Atkins, H. L.↗

DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Analysis for delamination initiation in postbuckled dropped-ply laminates

The compression strength of dropped-ply, graphite-epoxy laminated plates for the delamination mode of failure is studied by analysis and corroborated with experiments. The nonlinear response of the test specimens is modeled by a geometrically nonlinear finite element analysis. The methodology for predicting delamination is based on a quadratic interlaminar stress criterion evaluated at a characteristic distance from the ply drop-off. The compression strength of specimens exhibiting a linear response is greater than the compression strength of specimens with the same layup exhibiting a geometrically nonlinear response. The analyses for both linear and nonlinear response show that severe interlaminar stress gradients occur in the interfaces at the drop-off because of the thickness/stiffness discontinuity. However, these interlaminar stress distributions are altered in the geometrically nonlinear response such that, with increasing load, their growth at the center of the laminate is retarded while their growth near the unloaded supported edge is increased.

Davila, Carlos G.↗

NASA NTRS2000

Experiment Evaluation of Bifurcation in Sands

The basic principles of bifurcation analysis have been established by several investigators, however several issues remain unresolved, specifically how do stress level, grain size distribution, and boundary conditions affect general bifurcation phenomena in pressure sensitive and dilatant materials. General geometrical and kinematics conditions for moving surfaces of discontinuity was derived and applied to problems of instability of solids. In 1962, the theoretical framework of bifurcation by studying the acceleration waves in elasto-plastic (J2) solids were presented. Bifurcation analysis for more specific forms of constitutive behavior was examined by studying localization in pressure-sensitive, dilatant materials, however, analyses were restricted to plane deformation states only. Bifurcation analyses were presented and applied to predict shear band formations in sand under plane strain condition. The properties of discontinuous bifurcation solutions for elastic-plastic solids under axisymmetric and plane strain loading conditions were studied. The study focused on theory, but also references and comparisons to experiments were made. The current paper includes a presentation of a summary of bifurcation analyses for biaxial and triaxial (axisymmetric) loading conditions. The Coulomb model is implemented using incremental piecewise scheme to predict the constitutive relations and shear band inclination angles. Then, a comprehensive evaluation of bifurcation phenomena is presented based on data from triaxial experiments performed under microgravity conditions aboard the Space Shuttle under very low effective confining pressure (0.05 to 1.30 kPa), in which very high peak friction angles (47 to 75 degrees) and dilatancy angles (30 to 31 degrees) were measured. The evaluation will be extended to include biaxial experiments performed on the same material under low (10 kPa) and moderate (100 kPa) confining pressures. A comparison between the behavior under biaxial and triaxial loading conditions will be presented, and related issues concerning influence of confining pressure will be discussed.

Alshibi, Khalid A.↗

DOE OSTI2022

Quantum simulation of quantum phase transitions using the convex geometry of reduced density matrices

Transitions of many-particle quantum systems between distinct phases at absolute-zero temperature, known as quantum phase transitions, require an exacting treatment of particle correlations. Here in this work, we present a general quantum-computing approach to quantum phase transitions that exploits the geometric structure of reduced density matrices. While typical approaches to quantum phase transitions examine discontinuities in the order parameters, the origin of phase transitions—their order parameters and symmetry breaking—can be understood geometrically in terms of the set of two-particle reduced density matrices (2-RDMs). The convex set of 2-RDMs provides a comprehensive map of the quantum system including its distinct phases as well as the transitions connecting these phases. Because 2-RDMs can potentially be computed on quantum computers at nonexponential cost, even when the quantum system is strongly correlated, they are ideally suited for a quantum-computing approach to quantum phase transitions. We compute the convex set of 2-RDMs for a Lipkin-Meshkov-Glick spin model on IBM superconducting-qubit quantum processors. Even though computations are limited to few-particle models due to device noise, comparisons with a classically solvable 1000-particle model reveal that the finite-particle quantum solutions capture the key features of the phase transitions including the strong correlation and the symmetry breaking.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1985

Geometric-optical Modeling of a Conifer Forest Canopy

The objective of this research is to explore how the geometry of trees in forest stands influences the reflectance of the forest as imaged from space. Most plant canopy modeling has viewed the canopy as an assemblage of plane-parallel layers on top of a soil surface. For these models, leaf angle distribution, leaf area index, and the angular transmittance and reflectance of leaves are the primary optical and geometric parameters. Such models are now sufficiently well developed to explain most of the variance in angular reflectance measurements observed from homogeneous plant canopies. However, forest canopies as imaged by airborne and spaceborne scanners exhibit considerable variance at quite a different scale. Brightness values vary strongly from one pixel to the next primarily as a function of the number of trees they contain. At this scale, the forest canopy is nonuniform and discontinuous. This research focuses on a discrete-element, geometric-optical view of the forest canopy.

Strahler, A. H.↗

NASA NTRS1991

Improved finite-element methods for rotorcraft structures

An overview of the research directed at improving finite-element methods for rotorcraft airframes is presented. The development of a modification to the finite element method which eliminates interelement discontinuities is covered. The following subject areas are discussed: geometric entities, interelement continuity, dependent rotational degrees of freedom, and adaptive numerical integration. This new methodology is being implemented as an anisotropic, curvilinear, p-version, beam, shell, and brick finite element program.

Hinnant, Howard E.↗

NASA NTRS2019

Numerical Predictions of Dust-Induced Heat Flux Augmentation in Hypersonic Blunt-Body Flows Using a Discontinuous Galerkin Multiphase Flow Solver

Recent interest in human-scale missions to Mars has motivated the need for high-fidelity simulations of reentry flows. During a dust storm, there can be high levels of suspended dust in the Martian atmosphere, which cannot only enhance erosion of thermal protection systems but also transfer energy and momentum to the shock layer, thereby significantly augmenting the surface heat flux. Second-order finite-volume schemes are typically employed for hypersonic flow simulations, but such schemes suffer from a number of disadvantages. An attractive alternative is discontinuous Galerkin methods, which benefit from arbitrarily high spatial order of accuracy, geometric flexibility, and other properties. To enable accurate computations of high-speed particle-laden flows, an Euler-Lagrange methodology was developed in which the Eulerian field of the carrier gas is calculated using a discontinuous Galerkin scheme while the disperse phase is treated with Lagrangian particle tracking. We discuss challenges associated with coupling these two formulations and how to handle them. Momentum and energy transfer between the carrier gas and the particle phase is considered, and the importance of accounting for interparticle collisions is assessed. In addition, we describe the physical model of the particle phase and examine effects of its uncertainties on the numerical solution. We demonstrate the performance of the Euler-Lagrange method in representative testcases, with focus on the accurate prediction of particle trajectories and heating augmentation. Quantitative comparisons with experiments are provided.

Ching, Eric J.↗

NASA NTRS1992

Coupled 2D-3D finite element method for analysis of a skin panel with a discontinuous stiffener

This paper describes a computationally efficient analysis method which was used to predict detailed stress states in a typical composite compression panel with a discontinuous hat stiffener. A global-local approach was used. The global model incorporated both 2D shell and 3D brick elements connected by newly developed transition elements. Most of the panel was modeled with 2D elements, while 3D elements were employed to model the stiffener flange and the adjacent skin. Both linear and geometrically nonlinear analyses were performed on the global model. The effect of geometric nonlinearity induced by the eccentric load path due to the discontinuous hat stiffener was significant. The local model used a fine mesh of 3D brick elements to model the region at the end of the stiffener. Boundary conditions of the local 3D model were obtained by spline interpolation of the nodal displacements from the global analysis. Detailed in-plane and through-the-thickness stresses were calculated in the flange-skin interface near the end of the stiffener.

Wang, J. T.↗

NASA NTRS1993

Compression strength failure mechanisms in unidirectional composite laminates containing a hole

Experiments on graphite-epoxy laminated plates containing unloaded small holes show that these laminates are notch insensitive. That is, the uniaxial strength of these laminates with small holes exceeds the strength predicted by a point stress criterion using the stress concentration factor for the in-plane stress field. Laminates containing large holes exhibit notch sensitive behavior and consequently their strength is reasonably well predicted by the stress concentration effect. This hole size effect is manifested both in tension and in compression. Apparently, some mechanism must cause in-plane stress relief for laminates containing small holes. The purpose of this research was to study the influence of geometric nonlinearity on the micromechanical response of a filamentary composite material in the presence of a strain gradient caused by a discontinuity such as a hole. A mathematical model was developed at the micromechanical level to investigate this geometrically nonlinear effect.

Johnson, Eric R.↗

NASA NTRS1992

Adaptive finite element strategies for shell structures

The present paper extends existing finite element adaptive refinement (AR) techniques to shell structures, which have heretofore been neglected in the AR literature. Specific challenges in applying AR to shell structures include: (1) physical discontinuities (e.g., stiffener intersections); (2) boundary layers; (3) sensitivity to geometric imperfections; (4) the sensitivity of most shell elements to mesh distortion, constraint definition and/or thinness; and (5) intrinsic geometric nonlinearity. All of these challenges but (5) are addressed here.

Stanley, G.↗

DOE OSTI2025

An innovative radial gradient material design using hot isostatic pressing for applications in extreme environments

Functionally graded materials (FGMs) are highly advanced continuous or discontinuous structures whose structural and material properties vary along a singular geometric dimension either in the axial or radial direction. Here, the radial gradient FGM design makes for an optimal structural design to incorporate a bi-metallic structure with a copper-based high entropy alloy (Cu-HEA) with good mechanical properties and high irradiation resistance, and Chromium (Cr) with great corrosion resistance. This study focuses on the experimental design of a metal powder loading mechanism to fabricate a bi-metallic radial gradient structure using Cu-HEA and 99.9 % pure Cr metal powders. The powder loading strategy uses custom-designed concentric cylindrical dividers to separate the individual compositions. Two benchtop trial runs were performed for design optimization. The optimized design was then implemented to eventually load the HEA and Cr powders for consolidation via powder metallurgy hot isostatic pressing (PM-HIP). The electron microscopy analysis reveals the successful fabrication of the radial gradient structure with the chemical mapping analysis, demonstrating the gradual composition shift from the HEA at the center to the pure-Cr at the periphery via a three-step gradient.

High Entropy Alloys (HEAs)↗

NASA NTRS1994

Spectral Methods for Computational Fluid Dynamics

As a tool for large-scale computations in fluid dynamics, spectral methods were prophesized in 1944, born in 1954, virtually buried in the mid-1960's, resurrected in 1969, evangalized in the 1970's, and catholicized in the 1980's. The use of spectral methods for meteorological problems was proposed by Blinova in 1944 and the first numerical computations were conducted by Silberman (1954). By the early 1960's computers had achieved sufficient power to permit calculations with hundreds of degrees of freedom. For problems of this size the traditional way of computing the nonlinear terms in spectral methods was expensive compared with finite-difference methods. Consequently, spectral methods fell out of favor. The expense of computing nonlinear terms remained a severe drawback until Orszag (1969) and Eliasen, Machenauer, and Rasmussen (1970) developed the transform methods that still form the backbone of many large-scale spectral computations. The original proselytes of spectral methods were meteorologists involved in global weather modeling and fluid dynamicists investigating isotropic turbulence. The converts who were inspired by the successes of these pioneers remained, for the most part, confined to these and closely related fields throughout the 1970's. During that decade spectral methods appeared to be well-suited only for problems governed by ordinary diSerential eqllations or by partial differential equations with periodic boundary conditions. And, of course, the solution itself needed to be smooth. Some of the obstacles to wider application of spectral methods were: (1) poor resolution of discontinuous solutions; (2) inefficient implementation of implicit methods; and (3) drastic geometric constraints. All of these barriers have undergone some erosion during the 1980's, particularly the latter two. As a result, the applicability and appeal of spectral methods for computational fluid dynamics has broadened considerably. The motivation for the use of spectral methods in numerical calculations stems from the attractive approximation properties of orthogonal polynomial expansions.

Zang, T. A.↗