Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized finite differences”

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

On the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers: Selective mode-type excitation and specific mode enhancement

Here, we investigate the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers. Non-uniform polarization is represented as a non-uniform electromechanical coupling coefficient using a polarization function. Governing equations are derived for the forced in-plane vibrations of a thin wafer under the assumption of generalized plane stress. The effect of non-uniform polarization is explicitly obtained in the forcing terms of the governing equations. These equations are then recast into a variational weak form that is then solved using the finite element method to obtain the displacement fields for different modes. The electromechanical response of a non-uniformly polarized piezoelectric wafer is derived in terms of the out-of-plane displacement profile on the surface of the wafer. Using the derived analytical expression, a necessary and sufficient condition for the presence/absence of a vibrational mode in the electromechanical impedance spectrum is obtained. Based on this condition, criteria for selective mode-type excitation and specific mode enhancement of vibrational modes in the electromechanical impedance spectrum are postulated. Selective mode-type excitation of in-plane extensional and shear modes is demonstrated for a square wafer and that of in-plane bending modes is demonstrated for a rectangular wafer. Specific mode enhancement is demonstrated for both square and rectangular wafers. In addition, it is also demonstrated how the criteria can be used to suppress specific vibrational modes in the electromechanical impedance spectrum. The proposed methodology of using non-uniformly polarized piezoelectric wafers finds application in the design of single element transducers with multi-frequency operation, frequency-tuned receivers/sensors, acoustic holograms, designing acoustic beams of prescribed shape/lobes, and other non-traditional applications such as information storage.

36 MATERIALS SCIENCE↗
NASA NTRS1982

Fluid dynamical description of relativistic nuclear collisions

On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Nix, J. R.↗
NASA NTRS1982

Recent applications of the transonic wing analysis computer code, TWING

An evaluation of the transonic-wing-analysis computer code TWING is given. TWING utilizes a fully implicit approximate factorization iteration scheme to solve the full potential equation in conservative form. A numerical elliptic-solver grid-generation scheme is used to generate the required finite-difference mesh. Several wing configurations were analyzed, and the limits of applicability of this code was evaluated. Comparisons of computed results were made with available experimental data. Results indicate that the code is robust, accurate (when significant viscous effects are not present), and efficient. TWING generally produces solutions an order of magnitude faster than other conservative full potential codes using successive-line overrelaxation. The present method is applicable to a wide range of isolated wing configurations including high-aspect-ratio transport wings and low-aspect-ratio, high-sweep, fighter configurations.

Subramanian, N. R.↗
NASA NTRS2014

Modeling Sound Propagation Through Non-Axisymmetric Jets

A method for computing the far-field adjoint Green's function of the generalized acoustic analogy equations under a locally parallel mean flow approximation is presented. The method is based on expanding the mean-flow-dependent coefficients in the governing equation and the scalar Green's function in truncated Fourier series in the azimuthal direction and a finite difference approximation in the radial direction in circular cylindrical coordinates. The combined spectral/finite difference method yields a highly banded system of algebraic equations that can be efficiently solved using a standard sparse system solver. The method is applied to test cases, with mean flow specified by analytical functions, corresponding to two noise reduction concepts of current interest: the offset jet and the fluid shield. Sample results for the Green's function are given for these two test cases and recommendations made as to the use of the method as part of a RANS-based jet noise prediction code.

Leib, Stewart J.↗
DOE OSTI2025

General relativistic hydrodynamic simulations of binary strange star mergers

We perform fully general-relativistic simulations of binary strange star mergers considering two different approaches for thermal effects. The first uses a cold equation of state (EOS) derived from a modified version of the MIT bag model which is then supplemented by a Γ-law correction. The second approach employs a microphysical description of the finite-temperature effects. We describe results obtained with the two treatments, highlighting the influence of thermal effects. We find that the postmerger dynamics differs significantly in the two cases, leading to quantitative differences in the postmerger gravitational-wave spectrum and ejecta mass. The peak frequency of the postmerger gravitational-wave emission is consistent with the established quasi-universal relations for binary neutron star mergers and as a result, our simulations cannot distinguish between mergers of neutron stars and those of strange stars. Our models with realistic treatment of finite-temperature effects produce a significant amount of ejecta ≳0.02 M ⊙ ​. Here, the resulting flux of strangelets near the Earth, computed assuming that all neutron star mergers are in fact strange-stars mergers and that the binary considered here is representative, is in tension with experimental upper limits. As such, our results tentatively disfavor a scenario in which strange-quark matter is the lowest energy state of matter.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermal conductivity modeling of monodispersed microspheres using discrete element method

Particle beds are widely used in various systems and processes, such as particle heat exchangers, granular flow reactors, and additive manufacturing. Accurate modeling of the thermal conductivity of particle beds and understanding of their heat transfer mechanisms are important. However, previous models were based on simple cubic packing of particles, which could not accurately represent the actual heat transfer processes under certain conditions. Here, we examine the effect of the packing structure on the thermal conductivity of particle beds. We use monodispersed silica microspheres with average particle sizes ranging from 23 to 330 μm as a model material. We employ a transient hot-wire technique to measure the thermal conductivity of the particle beds with packing density of 43%–57% within a temperature range of room temperature to 500 °C and under N2 gaseous pressures of 20–760 Torr. We then use a discrete element method (DEM) to obtain the realistic packing structure of the particles, which is then fed into a finite-element model (FEM) to calculate the thermal conductivity, with the consideration of solid conduction, gas conduction, and radiation heat transfer. Our results show that the thermal conductivity model based on the more realistic random packing structure derived from the DEM shows better agreement with the experimental data compared to that based on the simple cubic-packing structure. The combined DEM and FEM methodology can serve as a useful tool to predict the effective thermal conductivity of particle beds and to quantify different heat transfer mechanisms under various conditions.

14 SOLAR ENERGY↗
NASA NTRS2019

Uniaxial Tensile Properties of AS4 3D Woven Composites with Four Different Resin Systems: Experimental Results & Analysis - Property Computations

As a part of the NASA Composite Technology for Exploration project, eight different AS4 3D orthogonal woven composite panels were manufactured and were subjected to mechanical testing including uniaxial tension along the weaves' warp direction. Each set, with four different resin systems (KCR-IR6070, EP2400, RTM6, and RS-50), included weave architectures designed using 12K and 6K AS4 carbon fiber yarns. For the tension testing conducted at Room Temperature Ambient (RTA) conditions, the elastic modulus and strength of these eight panels (as-processed and thermally cycled) were measured and compared while the potential evolution of micro-cracking before and after thermal cycling were monitored via optical microscopy and X-Ray Computed Tomography. The data set also included test results of the as-processed materials at Elevated Temperature Wet (ETW) conditions. In the second part of this study, efforts were made to compute elastic constants for AS4 6K/RTM6 and AS4 12K/RTM6 materials by implementing a finite element approach and the Multiscale Generalized Method of Cells (MSGMC) technique developed at NASA Glenn Research Center. Digimat-FE was used to model the weave architectures, assign properties, calculate yarn properties, create the finite element mesh, and compute the elastic properties by applying periodic boundary conditions to finite element models of each repeating unit cell. The required input data for MSGMC was generated using Matlab® from Digimat exported weave information. Experimental and computational results were compared, and the differences and limitations in correlating to the test data were briefly discussed.

Property Computations↗
NASA NTRS2019

Uniaxial Tensile Properties of AS4 3D Woven Composites with Four Different Resin Systems: Experimental Results and Analysis: Property Computations

As a part of the NASA Composite Technology for Exploration project, eight different AS4 3D orthogonal woven composite panels were manufactured and were subjected to mechanical testing including uniaxial tension along the weaves' warp direction. Each set, with four different resin systems (KCR-IR6070, EP2400, RTM6, and RS-50), included weave architectures designed using 12K and 6K AS4 carbon fiber yarns. For the tension testing conducted at Room Temperature Ambient (RTA) conditions, the elastic modulus and strength of these eight panels (as-processed and thermally-cycled) were measured and compared while the potential evolution of micro-cracking before and after thermal cycling were monitored via optical microscopy and X-Ray Computed Tomography. The data set also included test results of the as-processed materials at Elevated Temperature Wet (ETW) conditions. In the second part of this study, efforts were made to compute elastic constants for AS4 6K/RTM6 and AS4 12K/RTM6 materials by implementing a finite element approach and the Multiscale Generalized Method of Cells (MSGMC) technique developed at NASA Glenn Research Center. Digimat-FE was used to model the weave architectures, assign properties, calculate yarn properties, create the finite element mesh, and compute the elastic properties by applying periodic boundary conditions to finite element models of each repeating unit cell. The required input data for MSGMC was generated using Matlab® from Digimat exported weave information. Experimental and computational results were compared, and the differences and limitations in correlating to the test data were briefly discussed.

Property Computations↗
DOE OSTI2021

Explaining an unusual electromigration behavior—A comprehensive experimental and theoretical analysis using finite element method

In metallic interconnects, it is generally assumed that electromigration (EM) failure location is independent of the applied electrical current and always occurs at the highest-current-density area. Our experiments show otherwise. We designed an Al interconnect that alters its failure location by only varying the applied current density. The failure occurs near the high for a current above 2 × 10 7 A/cm 2 , but at a location with 59% of the maximum for lower current densities. Thermoreflectance thermal imaging is employed to gather time-dependent high-resolution spatial temperature distributions of the Al interconnect during EM. More importantly, we propose a computationally inexpensive 2D finite element method that tracks EM evolution in time and matches well with the observations from different experimental conditions. A detailed analysis covering the major driving forces of EM is carried out to understand the complex physics behind EM. The atomic depletion rate contributed by each force is quantitatively studied. By examining the results from every tested experimental condition, the model reveals that the temperature gradient is the key reason causing atomic depletion near the failure location. Graphical illustrations and qualitative analysis are provided to intuitively show the key findings of our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Generalized Korn’s inequalities for piecewise $H¹$ and $H²$ vector fields

The purpose of this paper is to construct a new class of discrete generalized Korn’s inequalities for piecewise $H^1$ vector fields and piecewise $H^2$ vector fields in three-dimensional space. The resulting Korn’s inequalities are different from the standard Korn’s inequalities, as they involve the trace-free symmetric gradient operator, in place of the usual symmetric gradient operator. Furthermore, it is anticipated that the new generalized Korn’s inequalities will be useful for the analysis of a broad range of finite element methods, including mixed finite element methods and discontinuous Galerkin methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1978

Some calculations of transonic potential flow for the NACA 64A006 airfoil with oscillating flap

A method for calculating the transonic flow over steady and oscillating airfoils was developed by Isogai. It solves the full potential equation with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. In this paper the method is described in general terms, and results are compared with experimental data for both steady flow and for oscillations at several values of reduced frequency. Good agreement for static pressures is shown for subcritical speeds, with increasing deviation as Mach number is increased into the supercritical speed range. Fair agreement with experiment was obtained at high reduced frequencies with larger deviations at low reduced frequencies.

Bennett, R. M.↗
NASA NTRS1979

Documentation of the Fourth Order Band Model

A general circulation model is presented which uses quadratically conservative, fourth order horizontal space differences on an unstaggered grid and second order vertical space differences with a forward-backward or a smooth leap frog time scheme to solve the primitive equations of motion. The dynamic equations for motion, finite difference equations, a discussion of the structure and flow chart of the program code, a program listing, and three relevent papers are given.

Kalnay-Rivas, E.↗
NASA NTRS1981

The challenge of unsteady separating flows

It is noted that the general fluid dynamic problem of unsteady separation at most practical Reynolds numbers remains an unsolved one and that no completely reliable prediction techniques exist at the present time. The modern design engineer must therefore draw from a combination of approximate theories, empirical correlations of data, and finite difference programs based on uncertain physical modeling of turbulence. An attempt is made to describe the basic features of several representative classes of problems for which unsteady effects produce strong or unusual changes in the separation characteristics of the flow. The analysis concerns itself largely with external flow, and emphasis is placed on the physical phenomena involved.

Mccroskey, W. J.↗
NASA NTRS1983

Calculation of unsteady small disturbance transonic flow at arbitrary reduced frequency using an extended integral equation method.

A novel analytical-numerical method for calculating unsteady small disturbance transonic flow over airfoils has been developed. The method uses an extended integral equation technique, based on both the velocity potential and the acceleration potential, to predict unsteady aerodynamic loading on airfoils oscillating in subcritical transonic free stream conditions. The formulation is an extension of the work of Sivaneri and Harris (1980) for steady, non-lifting flows and utilizes the linear theory of Landahl (1961) for decoupling of steady and unsteady components. The analytical-numerical procedure involves several intnegrating schemes and applies to general frequencies of oscillations. The technique is illustrated by computing the transonic flow about parabolic arc airfoils. Specific unsteady results for reduced frequencies based on semi-chord of 0.01, 0.1, 0.3, 0.4 and 0.6 are given. Comparison of results with those obtained by an ADI finite difference scheme is made.

Lamah, C. A.↗
NASA NTRS1989

Computing Viscoplastic Behavior Of A Material

Finite-element implementation developed for Robinson's unified model of viscoplasticity. Implemented via MARC general-purpose finite-element computer program by incorporating all of nonlinearity of material into initial load vector and treating vector as pseudo body force in finite-element equilibrium equations. In technique, each finite increment split into several equal subincrements, constitutive equations integrated by explicit Euler forward-difference method. HYPELA subroutine in MARC includes adaptive integration scheme selecting optimal sizes of subincrements. Spline function smooths discontinuous boundaries in Robinson's model to facilitate finite-element numerical calculations. Method used to calculate elastic, plastic, and creep deformations in solid materials at various temperatures; also used to study thermomechanical behavior of such things as hot pressure vessels and highly stressed components in engines.

Arya, V. K.↗
NASA NTRS2022

A Benchmark Example for Delamination Propagation Predictions Based on the Calibrated End-Loaded Split Specimen

A benchmark example based on the Calibrated End-Loaded Split (C-ELS) specimen is developed and used to assess the performance of recently developed delamination propagation capabilities in the Abaqus/Standard finite element code. The C-ELS specimen has the advantage of a longer region of stable delamination propagation compared to the existing mode II benchmark case. The new benchmark example may therefore provide a better assessment tool by enabling more stable crack growth in regions further away from the boundary conditions or load application. First, a benchmark result is created manually using two-dimensional finite element models of the C-ELS specimen with different delamination lengths. Second, the performance of the delamination propagation capabilities in the Abaqus/Standard finite element code are assessed by comparing the results to the benchmark case. Two examples with different starter delamination lengths are studied. A shorter starter length is chosen to create a scenario with unstable delamination propagation. A longer delamination causes stable delamination propagation. Detailed results from three-dimensional analyses with aligned and mis-aligned meshes and two levels of mesh refinement are provided for several permutations of numerical input parameters. In general, good agreement can be achieved between the results obtained from the quasi-static propagation analysis and the benchmark analysis. However, particular non-default settings are found to be most reliable, accurate, and numerically efficient. Numerical artifacts including anomalous unreleased nodes in the crack wake and zig-zag crack fronts occur, and further development of the Abaqus/Standard VCCT propagation may be required. Use of the benchmark case to assess the continuous improvements in one finite element code illustrates the value of establishing benchmark solutions.

delamination↗
DOE OSTI2017

A Comparison of Grid-based and SPH Binary Mass-transfer and Merger Simulations

There is currently a great amount of interest in the outcomes and astrophysical implications of mergers of double degenerate binaries. In a commonly adopted approximation, the components of such binaries are represented by polytropes with an index of n = 3/2. We present detailed comparisons of stellar mass-transfer and merger simulations of polytropic binaries that have been carried out using two very different numerical algorithms—a finite-volume “grid” code and a smoothed-particle hydrodynamics (SPH) code. We find that there is agreement in both the ultimate outcomes of the evolutions and the intermediate stages if the initial conditions for each code are chosen to match as closely as possible. We find that even with closely matching initial setups, the time it takes to reach a concordant evolution differs between the two codes because the initial depth of contact cannot be matched exactly. There is a general tendency for SPH to yield higher mass transfer rates and faster evolution to the final outcome. We also present comparisons of simulations calculated from two different energy equations: in one series, we assume a polytropic equation of state and in the other series an ideal gas equation of state. In the latter series of simulations, an atmosphere forms around the accretor, which can exchange angular momentum and cause a more rapid loss of orbital angular momentum. In the simulations presented here, the effect of the ideal equation of state is to de-stabilize the binary in both SPH and grid simulations, but the effect is more pronounced in the grid code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Solidification and crystallographic texture modeling of laser powder bed fusion Ti-6Al-4V using finite difference-monte carlo method

Laser powder bed fusion (LPBF) additive manufacturing makes near-net-shaped parts with reduced material cost and time, rising as a promising technology to fabricate Ti-6Al-4V, a widely used titanium alloy in aerospace and medical industries. However, LPBF Ti-6Al-4V parts produced with 67° rotation between layers, a scan strategy commonly used to reduce microstructure and property inhomogeneity, have varying grain morphologies and weak crystallographic textures that change depending on processing parameters. Here, this study predicts LPBF Ti-6Al-4V solidification at three energy levels using a finite difference-Monte Carlo method and validates the simulations with large-area electron backscatter diffraction (EBSD) scans. The developed model accurately shows that a <001> texture forms at low energy and a <111> texture occurs at higher energies parallel to the build direction but with a lower strength than the textures observed from EBSD. A validated and well-established method of combining spatial correlation and general spherical harmonics representation of texture is developed to calculate a difference score between simulations and experiments. The quantitative comparison enables effective fine-tuning of nucleation density (N 0 ) input, which shows a nonlinear relationship with increasing energy level. Future improvements in texture prediction code and a more comprehensive study of N 0 with different energy levels will further advance the optimization of LPBF Ti-6Al-4V components. These developments contribute a novel understanding of crystallographic texture formation in LPBF Ti-6Al-4V, the development of robust model validation and calibration pipeline methodologies, and provide a platform for mechanical property prediction and process parameter optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗