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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Low-energy electron-atom bremsstrahlung

This paper extends recent work providing an elementary calculation of bremsstrahlung and opacity associated with the scattering of low-energy electrons by neutral atoms and molecules. The method applies when the scattering potential is 'hard' or when the collision time is short, applying the classical soft-photon emission probability formula for arbitrary bremsstrahlung photon energy. However, now, in addition to correcting the probability factor for finite bremsstrahlung photon energy, another factor corrects for the reduced phase space available to the outgoing electron. The bremsstrahlung cross section and opacity are then computed directly from the elastic scattering cross section, determined experimentally or calculated; a small (approximately 10 percent) correction is computed from the anisotropic term in the elastic scattering cross section. The opacity is evaluated for electron scattering by H, He, and H2 using experimentally determined values for the elastic scattering cross section, and is compared with more elaborate calculations. The agreement is good (within 10 percent), indicating an accuracy for the general method comparable to variations among the results of different elaborate theoretical computations. The agreement seems to validate the basic approximation of short collision time even at large bremsstrahlung photon energy for electron energies and temperatures up to a few eV.

Gould, R. J.↗

NASA NTRS1991

Supersonic propulsion simulation by incorporating component models in the large perturbation inlet (LAPIN) computer code

An approach to simulating the internal flows of supersonic propulsion systems is presented. The approach is based on a fairly simple modification of the Large Perturbation Inlet (LAPIN) computer code. LAPIN uses a quasi-one dimensional, inviscid, unsteady formulation of the continuity, momentum, and energy equations. The equations are solved using a shock capturing, finite difference algorithm. The original code, developed for simulating supersonic inlets, includes engineering models of unstart/restart, bleed, bypass, and variable duct geometry, by means of source terms in the equations. The source terms also provide a mechanism for incorporating, with the inlet, propulsion system components such as compressor stages, combustors, and turbine stages. This requires each component to be distributed axially over a number of grid points. Because of the distributed nature of such components, this representation should be more accurate than a lumped parameter model. Components can be modeled by performance map(s), which in turn are used to compute the source terms. The general approach is described. Then, simulation of a compressor/fan stage is discussed to show the approach in detail.

Cole, Gary L.↗

NASA NTRS2003

Modeling and Validation of Damped Plexiglas Windows for Noise Control

Windows are a significant path for structure-borne and air-borne noise transmission in general aviation aircraft. In this paper, numerical and experimental results are used to evaluate damped plexiglas windows for the reduction of structure-borne and air-borne noise transmitted into the interior of an aircraft. In contrast to conventional homogeneous windows, the damped plexiglas windows were fabricated using two or three layers of plexiglas with transparent viscoelastic damping material sandwiched between the layers. Transmission loss and radiated sound power measurements were used to compare different layups of the damped plexiglas windows with uniform windows of the same nominal thickness. This vibro-acoustic test data was also used for the verification and validation of finite element and boundary element models of the damped plexiglas windows. Numerical models are presented for the prediction of radiated sound power for a point force excitation and transmission loss for diffuse acoustic excitation. Radiated sound power and transmission loss predictions are in good agreement with experimental data. Once validated, the numerical models were used to perform a parametric study to determine the optimum configuration of the damped plexiglas windows for reducing the radiated sound power for a point force excitation.

Buehrle, Ralph D.↗

NASA NTRS2024

Finite Element Simulation of the Direct Energy Deposition using Comsol Multiphysics

Direct Energy Deposition (DED) is an emerging technology extensively employed in metal Additive Manufacturing (AM). Despite its widespread industrial application, mathematical modeling in this domain remains challenging. This complexity arises from the intricate nature of the modeling approach and the nonlinear behavior of material parameters across a broad temperature range. Consequently, experimentalists often resort to a trial-and-error method to achieve structures with desired properties. However, this approach can be time-consuming and may not always yield parts with the requisite characteristics, highlighting the necessity for mathematical modeling to enhance manufacturing success. This study focuses on thin-walled manufacturing with a single bead thickness, adjusting laser parameters to produce such structures. The laser cladding speed, typical for DED manufacturing, is considered to be on the order of centimeters per second. The various wall thicknesses and heights are explored to discern the general characteristics of walls manufactured via this DED approach. Our modeling approach diverges from the conventional DED modeling based on activation-deactivation of predefined mesh domains, commonly implemented in many finite element codes. Instead, a method is proposed that effectively models layer cladding and melt pool dynamics, enabling predictions of the microstructure in the resulting structures. The formation mechanisms of cellular, dendritic columnar, and stray (equiaxed) grains, which arise from the interplay between nucleation and growth from the surface is analyzed. The study also examines the feasibility of microstructure formation to demonstrate the various thermal characteristics inherent in wall manufacturing. Our results illustrate how different process parameters influence the temperature gradient and cooling rate of the molten pool, subsequently affecting the primary dendrite arm spacing (PDAS). This modeling technique allows for the investigation of diverse thermal conditions, facilitating the prediction of microstructure and residual stresses in the manufactured parts.

modeling↗

DOE OSTI2020

PRISMS-PF: A general framework for phase-field modeling with a matrix-free finite element method

Abstract A new phase-field modeling framework with an emphasis on performance, flexibility, and ease of use is presented. Foremost among the strategies employed to fulfill these objectives are the use of a matrix-free finite element method and a modular, application-centric code structure. This approach is implemented in the new open-source PRISMS-PF framework. Its performance is enabled by the combination of a matrix-free variant of the finite element method with adaptive mesh refinement, explicit time integration, and multilevel parallelism. Benchmark testing with a particle growth problem shows PRISMS-PF with adaptive mesh refinement and higher-order elements to be up to 12 times faster than a finite difference code employing a second-order-accurate spatial discretization and first-order-accurate explicit time integration. Furthermore, for a two-dimensional solidification benchmark problem, the performance of PRISMS-PF meets or exceeds that of phase-field frameworks that focus on implicit/semi-implicit time stepping, even though the benchmark problem’s small computational size reduces the scalability advantage of explicit time-integration schemes. PRISMS-PF supports an arbitrary number of coupled governing equations. The code structure simplifies the modification of these governing equations by separating their definition from the implementation of the numerical methods used to solve them. As part of its modular design, the framework includes functionality for nucleation and polycrystalline systems available in any application to further broaden the phenomena that can be used to study. The versatility of this approach is demonstrated with examples from several common types of phase-field simulations, including coarsening subsequent to spinodal decomposition, solidification, precipitation, grain growth, and corrosion.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Anti-symmetric plasma moment equations with symbiotic discrete counterparts

We construct a set of fluid moment equations with inherent consistency and numerical stability, conceived by exploiting the anti-symmetric nature of the plasma flow operator. The obtained equations can be interpreted as an alternative to the traditional Eulerian and Lagrangian representations - one in which plasma flows generate in infinitesimal rotations of generalized fluid moments $\sqrt{n}$, $\sqrt{nv}$, $\sqrt{p}$. The resulting model has a discrete analog, with positivity preservation, and exact mass, momentum, and energy conservation. This property is illustrated in simulations of seeded blob propagation, where we verify numerical conservation to machine accuracy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Development of 3D electromagnetic modeling tools for airborne vehicles

The main goal of this project is to develop methodologies for scattering by airborne composite vehicles. Although our primary focus continues to be the development of a general purpose code for analyzing the entire structure as a single unit, a number of other tasks are also pursued in parallel with this effort. These tasks are important in testing the overall approach and in developing suitable models for materials coatings, junctions and, more generally, in assessing the effectiveness of the various parts comprising the final code. Here, we briefly discuss our progress on the five different tasks which were pursued during this period. Our progress on each of these tasks is described in the detailed reports (listed at the end of this report) and the memoranda included. The first task described below is, of course, the core of this project and deals with the development of the overall code. Undoubtedly, it is the outcome of the research which was funded by NASA-Ames and the Navy over the past three years. During this year we developed the first finite element code for scattering by structures of arbitrary shape and composition. The code employs a new absorbing boundary condition which allows termination of the finite element mesh only 0.3 lambda from the outer surface of the target. This leads to a remarkable reduction of the mesh size and is a unique feature of the code. Other unique features of this code include capabilities to model resistive sheets, impedance sheets and anisotropic materials. This last capability is the latest feature of the code and is still under development. The code has been extensively validated for a number of composite geometries and some examples are given. The validation of the code is still in progress for anisotropic and larger non-metallic geometries and cavities. The developed finite element code is based on a Galerkin's formulation and employs edge-based tetrahedral elements for discretizing the dielectric sections and the region between the target and the outer mesh termination boundary (ATB). This boundary is placed in conformity with the target's outer surface, thus resulting in additional reduction of the unknown count.

Volakis, John L.↗

NASA NTRS1995

Nonlinear interaction between a pair of oblique modes in a supersonic mixing layer: Long-wave limit

The nonlinear interaction between a pair of symmetric, oblique, and spatial instability modes is studied in the long-wave limit using asymptotic methods. The base flow is taken to be a supersonic mixing layer whose Mach number is such that the corresponding vortex sheet is marginally stable according to Miles' criterion. It is shown that the amplitude of the mode obeys a nonlinear integro-differential equation. Numerical solutions of this equation show that, when the obliqueness angle is less than pi/4, the effect of the nonlinearity is to enhance the growth rate of the instability. The solution terminates in a singularity at a finite streamwise location. This result is reminiscent of that obtained in the vicinity of the neutral point by other authors in several different types of flows. On the other hand, when the obliqueness angle is more than pi/4, the streamwise development of the amplitude is characterized by a series of modulations. This arises from the fact that the nonlinear term in the amplitude equation may be either stabilizing or destabilizing, depending on the value of the streamwise coordinate. However, even in this case the amplitude of the disturbance increases, though not as rapidly as in the case for which the angle is less than pi/4. Quite generally then, the nonlinear interaction between two oblique modes in a supersonic mixing layer enhances the growth of the disturbance.

Balsa, Thomas F.↗

DOE OSTI2024

Universal fine grained asymptotics of free and weakly coupled quantum field theory

We give a rigorous proof that in any free quantum field theory with a finite group global symmetry G, on a compact spatial manifold, at sufficiently high energy, the density of states ρ α (E) for each irreducible representation α of G obeys a universal formula as conjectured by Harlow and Ooguri. We further prove that this continues to hold in a weakly coupled quantum field theory, given an appropriate scaling of the coupling with temperature. This generalizes similar results that were previously obtained in (1 + 1)-D to higher spacetime dimension. We discuss the role of averaging in the density of states, and we compare and contrast with the case of continuous group G, where we prove a universal, albeit different, behavior.

Discrete Symmetries↗

DOE OSTI2017

An iterative phase-space explicit discontinuous Galerkin method for stellar radiative transfer in extended atmospheres

In this work, a phase-space discontinuous Galerkin (PSDG) method is presented for the solution of stellar radiative transfer problems. It allows for greater adaptivity than competing methods without sacrificing generality. The method is extensively tested on a spherically symmetric, static, inverse-power-law scattering atmosphere. Results for different sizes of atmospheres and intensities of scattering agreed with asymptotic values. The exponentially decaying behavior of the radiative field in the diffusive-transparent transition region, and the forward peaking behavior at the surface of extended atmospheres were accurately captured. The integrodifferential equation of radiation transfer is solved iteratively by alternating between the radiative pressure equation and the original equation with the integral term treated as an energy density source term. In each iteration, the equations are solved via an explicit, flux-conserving, discontinuous Galerkin method. Finite elements are ordered in wave fronts perpendicular to the characteristic curves so that elemental linear algebraic systems are solved quickly by sweeping the phase space element by element. Two implementations of a diffusive boundary condition at the origin are demonstrated wherein the finite discontinuity in the radiation intensity is accurately captured by the proposed method. This allows for a consistent mechanism to preserve photon luminosity. The method was proved to be robust and fast, and a case is made for the adequacy of parallel processing. In addition to classical two-dimensional plots, results of normalized radiation intensity were mapped onto a log-polar surface exhibiting all distinguishing features of the problem studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Nucleation and structural growth of cluster crystals

We study the nucleation of crystalline cluster phases in the generalized exponential model with exponent n = 4. Due to the finite value of this pair potential for zero separation, at high densities the system forms cluster crystals with multiply occupied lattice sites. Here, we investigate the microscopic mechanisms that lead to the formation of cluster crystals from a supercooled liquid in the low-temperature region of the phase diagram. Using molecular dynamics and umbrella sampling, we calculate the free energy as a function of the size of the largest crystalline nucleus in the system, and compare our results with predictions from classical nucleation theory. Employing bond-order parameters based on a Voronoi tessellation to distinguish different crystal structures, we analyze the average composition of crystalline nuclei. We find that even for conditions where a multiply occupied fcc crystal is the thermodynamically stable phase, the nucleation into bcc cluster crystals is strongly preferred. Furthermore, we study the particle mobility in the supercooled liquid and in the cluster crystal. In the cluster crystal, the motion of individual particles is captured by a simple reaction-diffusion model introduced previously to model the kinetics of hydrogen bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

A THz driven split-ring resonator based ultrafast relativistic electron streak camera

The use of sub-wavelength metal structures to locally enhance high frequency electromagnetic fields, generally known as plasmonics, enables breakthrough opportunities across diverse fields of research such as nonlinear optics, biosensing, photovoltaics and others. Here we study the application of sub-wavelength metallic resonators tuned in the THz frequency range for manipulation and diagnostics of relativistic electron beams. In this work, we report on the use of a double-sided split-ring structure driven by a near single cycle THz field generated by optical rectification to impart a time-dependent angular deviation (streak) on a 4.5 MeV electron beam. Electrons passing through the small gap reveal field enhancement factors larger than 10, in good agreement with finite difference time domain simulations. This work paves the way for further application of high frequency metallic structures in compact particle accelerators such as for THz-based relativistic electron streaking at fs and sub-fs temporal resolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Mechanical Loads on Spent Nuclear Fuel in the General 30 cm Package Drop Scenario

The U.S. Department of Energy Office of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) research program is guided by the high-level goal of closing prioritized knowledge gaps related to spent nuclear fuel (SNF) storage and transportation, which are summarized by Saltzstein et al. (2020). One of the high-priority knowledge gaps is the identification and quantification of mechanical loads that are expected to affect SNF during normal conditions of transportation and storage to inform the range of physical SNF test programs. This report uses modeling and analysis methods to estimate the mechanical loads on spent nuclear fuel (SNF) in the general 30 cm package drop scenario. The drop scenario assumes impact limiters are in place in the transportation configuration and the impact surface is perfectly rigid. The goal of this analysis is to consider the universe of potential mechanical loading conditions that can happen to SNF and present the results in a manner that is useful for materials testing, decision making, and regulatory rule making purposes. This study uses validated finite element models and methods to perform a broad parametric study of key variables that can affect the mechanical loads on SNF during a hypothetical package free drop scenario. Physical drop test data from a cask and fuel assembly drop test campaign is the basis for model validation. Additionally, the results of the parametric study are used to inform a damage model, which uses multiple nonlinear regression to estimate the relationships between input variables and output response. The parametric finite element analyses consider thousands of input variable combinations, while the damage model estimates millions of combinations. The breadth of this study provides confidence that the potential range of mechanical loads that SNF might experience during the general 30 cm package drop scenario are characterized well enough to consider this knowledge gap closed. While this report documents the overall peak values calculated in this study, the 95 th percentile values, the histograms, and the observed trends are equally important. This study covered a large range of SNF temperatures, room temperature to 300°C, and burnups, 10 GWd/MTU to 62 GWd/MTU. Each temperature and burnup combination has a different cladding yield strain, so it is more meaningful to summarize the calculated cladding strain response as its factor of safety, which is defined relative to the yield strain. The factor of safety is calculated as the yield strain divided by the peak cladding strain. A factor of safety greater than unity indicates that the cladding remains below yield, whereas a value less than unity is indicative of plastic deformation. In all cases of this study a safety factor over 1.0 was calculated, although in the most limiting case at 300°C the safety factor was only 1.01, which suggests that yielding could occur when additional loads like rod internal pressure are included. When the temperature is restricted to 200°C the limiting safety factor increases to 1.28, which has significant margin to accommodate internal pressure and potential local cladding defects that could cause a local stress concentration. An important trend in the calculated fuel rod mechanical loads is that the 2 nd highest loaded fuel rod in an assembly tends to be significantly lower than the highest loaded rod. The implication is that even if one rod in an assembly experiences a failure the loads would have to be significantly higher to cause two or more rods to fail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1976

Crack propagation and arrest in pressurized containers

The problem of crack propagation and arrest in a finite volume cylindrical container filled with pressurized gas is considered. It is assumed that the cylinder contains a symmetrically located longitudinal part-through crack with a relatively small net ligament. The net ligament suddenly ruptures initiating the process of fracture propagation and depressurization in the cylinder. Thus the problem is a coupled gas dynamics and solid mechanics problem the exact formulation of which does not seem to be possible. The problem is reduced to a proper initial value problem by introducing a dynamic fracture criterion which relates the crack acceleration to the difference between a load factor and the corresponding strength parameter. The results indicate that generally in gas filled cylinders fracture arrest is not possible unless the material behaves in a ductile manner and the container is relatively long.

Erdogan, F.↗

NASA NTRS1982

A harmonic analysis method for unsteady transonic flow and its application to the flutter of airfoils

A finite difference method for solving the unsteady transonic flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The differential equation for the unsteady velocity potential is linear with spatially varying coefficients and with the time variable eliminated by assuming harmonic motion. A study is presented of the shock motion associated with an oscillating airfoil and its representation by the harmonic procedure. The effects of the shock motion and the resulting pressure pulse are shown to be included in the harmonic pressure distributions and the corresponding generalized forces. Analytical and experimental pressure distributions for the NACA 64A010 airfoil are compared for Mach numbers of 0.75, 0.80 and 0.842. A typical section, two-degree-of-freedom flutter analysis of a NACA 64A010 airfoil is performed. The results show a sharp transonic bucket in one case and abrupt changes in instability modes.

Ehlers, F. E.↗

NASA NTRS1983

On one-dimensional stretching functions for finite-difference calculations

The class of one-dimensional stretching functions used in finite-difference calculations is studied. For solutions containing a highly localized region of rapid variation, simple criteria for a stretching function are derived using a truncation error analysis. These criteria are used to investigate two types of stretching functions. One an interior stretching function, for which the location and slope of an interior clustering region are specified. The simplest such function satisfying the criteria is found to be one based on the inverse hyperbolic sine. The other type of function is a two-sided stretching function, for which the arbitrary slopes at the two ends of the one-dimensional interval are specified. The simplest such general function is found to be one based on the inverse tangent. Previously announced in STAR as N80-25055

Vinokur, M.↗

NASA NTRS1984

Finite element thermal-structural analysis of cable-stiffened space structues

Finite element thermal-structural analyses of large, cable-stiffened space structures are presented. A computational scheme for the calculation of prestresses in the cable-stiffened structures is also described. The determination of thermal loads on orbiting space structures due to environment heating is discussed briefly. Three finite element structural analysis techniques are presented for the analysis of prestressed structures. Linear, stress stiffening, and large displacement analysis techniques were investigated. These three techniques were employed for analysis of prestressed cable structures at different prestress levels. The analyses produced similar results at small prestress, but at higher prestress, differences between the results became significant. For the cable-stiffened structures studied, the linear analysis technique may not provide acceptable results. The stress stiffening analysis technique may yield results of acceptable accuracy depending upon the level of prestress. The large displacement analysis technique produced accurate results over a wide range of prestress and is recommended as a general analysis technique for thermal-structural analysis of cable-stiffened space structures.

Thornton, E. A.↗

NASA NTRS1993

High Temperature Composite Analyzer (HITCAN) theoretical manual, version 1.0

This manual outlines some of the theoretical aspects embedded in the computer code, HITCAN (High Temperature Composite ANalyzer). HITCAN is a general purpose computer program for predicting nonlinear global structural and local stress-strain response of arbitrarily oriented, multilayered high temperature metal matrix composite structures. This code combines composite mechanics and laminate theory with an internal data base of the constituents (matrix, fiber and interface) material properties. The thermal and mechanical properties of the constituents are considered to be nonlinearly dependent on several parameters including temperature, stress and stress rate. The computational procedure for the analysis of the composite structure uses the finite element method. HITCAN consists of three modules: COBSTRAN, METCAN and MHOST. COBSTRAN generates the geometry (pre-processor) and defines the layup of the different plies. METCAN computes the material behavior of the composite and of the constituents. Finally, MHOST is a finite element program based on the mixed iterative solution technique. MHOST has a library for 2D and 3D isoparametric elements. HITCAN is written in FORTRAN 77 computer language and has been configured and executed on the NASA Lewis Research Center CRAY XMP and YMP computers.

Lackney, J. J.↗