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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

An energy and potential enstrophy conserving scheme for the shallow water equations

A three-dimensional finite difference scheme for the solution of the shallow water momentum equations which accounts for the conservation of potential enstrophy in the flow of a homogeneous incompressible shallow atmosphere over steep topography as well as for total energy conservation is presented. The scheme is derived to be consistent with a reasonable scheme for potential vorticity advection in a long-term integration for a general flow with divergent mass flux. Numerical comparisons of the characteristics of the present potential enstrophy-conserving scheme with those of a scheme that conserves potential enstrophy only for purely horizontal nondivergent flow are presented which demonstrate the reduction of computational noise in the wind field with the enstrophy-conserving scheme and its convergence even in relatively coarse grids.

Arakawa, A.↗
NASA NTRS1988

Efficient assembly of finite-element subsystems with large relative rotations

A finite element approach is presented for the modeling of rotorcraft undergoing elastic deformation in addition to large rigid body motion with respect to inertial space, with particular attention given to the coupling of the rotor and fuselage subsystems subject to large relative rotations. The component synthesis technique used here allows the coupling of rotors to the fuselage for different rotorcraft configurations. The formulation is general and applicable to any rotorcraft vibration, aeroelasticity, and dynamics problem.

Fuh, Jon-Shen↗
NASA NTRS2004

Reduced-Order Blade Mistuning Analysis Techniques Developed for the Robust Design of Engine Rotors

The primary objective of this research program is to develop vibration analysis tools, design tools, and design strategies to significantly improve the safety and robustness of turbine engine rotors. Bladed disks in turbine engines always feature small, random blade-to-blade differences, or mistuning. Mistuning can lead to a dramatic increase in blade forced-response amplitudes and stresses. Ultimately, this results in high-cycle fatigue, which is a major safety and cost concern. In this research program, the necessary steps will be taken to transform a state-of-the-art vibration analysis tool, the Turbo-Reduce forced-response prediction code, into an effective design tool by enhancing and extending the underlying modeling and analysis methods. Furthermore, novel techniques will be developed to assess the safety of a given design. In particular, a procedure will be established for using eigenfrequency curve veerings to identify "danger zones" in the operating conditions--ranges of rotational speeds and engine orders in which there is a great risk that the rotor blades will suffer high stresses. This work also will aid statistical studies of the forced response by reducing the necessary number of simulations. Finally, new strategies for improving the design of rotors will be pursued. Several methods will be investigated, including the use of intentional mistuning patterns to mitigate the harmful effects of random mistuning, and the modification of disk stiffness to avoid reaching critical values of interblade coupling in the desired operating range. Recent research progress is summarized in the following paragraphs. First, significant progress was made in the development of the component mode mistuning (CMM) and static mode compensation (SMC) methods for reduced-order modeling of mistuned bladed disks (see the following figure). The CMM method has been formalized and extended to allow a general treatment of mistuning. In addition, CMM allows individual mode mistuning, which accounts for the realistic effects of local variations in blade properties that lead to different mistuning values for different mode types (e.g., mistuning of the first torsion mode versus the second flexural mode). The accuracy and efficiency of the CMM method and the corresponding Turbo-Reduce code were validated for an example finite element model of a bladed disk.

Min, James B.↗
NASA NTRS2015

Testing Dissipative Magnetosphere Model Light Curves and Spectra with Fermi Pulsars

We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. comparing its high-energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light cylinder (LC), is assumed to be negligible, while outside the LC it rescales with a finite conductivity (sigma). In our approach we calculate the corresponding high-energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First, we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young gamma-ray pulsars, including features of the phase-resolved spectra. Second, we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The sigma values that best describe each of the pulsars in this group show an increase with the spin-down rate (E˙ ) and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.

acceleration of particles↗
NASA NTRS1987

Finite-element analysis of corner cracks in rectangular bars

Stress intensity factors for a wide range of quarter-elliptical corner cracks in rectangular bars are presented. Cracked configurations were subjected to remote tension, in-plane bending, or out-of-plane bending. Results were compared with other numerical results for a quarter-circular corner crack configuration. The present results along the interior of the crack generally agreed within 3 percent with those from the literature. Some larger differences (3 to 13 percent) were observed near the intersection of the crack front and the free surfaces in the boundary layer region. Analyses were performed to study the effect of mesh refinement in the boundary layer region and the influence of Poisson's ratio on the distribution of stress-intensity factors.

Raju, L. S.↗
DOE OSTI2025

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D code for its forward and adjoint flux solution. Previous work identified the PERSENT features required to be verified to support commercial design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying PERSENT’s ability to correctly calculate four key quantities: perturbation worth distributions, kinetics parameters, sensitivity coefficients, and cross section uncertainty quantification. This manuscript provides the verification tasks and their results with respect to these quantities needed for commercial design activities. For the perturbation worth distributions, hand calculations are deployed to verify the PERSENT calculated results. Similarly, hand calculation of the PERSENT computed kinetics parameters is also used to verify the PERSENT results. In both of these, the input to PERSENT is manipulated to ensure the hand calculation exactly matches the equations PERSENT is calculating. The sensitivity coefficients involve calculating the derivatives of a parameter (such as reactivity worth), with respect to the cross section data. Direct finite difference calculations with DIF3D are used to verify the PERSENT calculated results. For the uncertainty quantification, manufactured input to PERSENT is used to allow an exact hand calculation to reproduce the PERSENT calculated results. The work detailed in this report verified that significant issues were identified for earlier versions of PERSENT for sensitivity coefficients which were corrected in this work and thus version 12.1.0 of PERSENT must be used to reproduce all of the verified work in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Development of Thermoacoustic Sensors for Sodium-cooled Fast Reactors

This Final Report refers to the project “Development of Thermoacoustic Sensors for Sodium-cooled Fast Reactor Systems”, which was led by Westinghouse Electric Company (Westinghouse) and carried out in collaboration with Argonne National Laboratory (ANL) and University of Pittsburgh. Thermo-acoustic Power Sensors (TAPS) are self-powered, wireless sensors envisioned for measuring key parameters, such as local temperature and neutron flux, in a nuclear reactor core. This project was intended to specifically investigate their applicability to Sodium-cooled Fast Reactors (SFR). TAPS are non-invasive (wireless) and passive (self-powered) devices. The passivity derives from their ability to use conditions that “naturally” exist in a nuclear reactor, such as gamma and neutron flux, as power sources. They generate oscillating pressure waves (i.e., sound waves) which, with a frequency and amplitude dependent upon these conditions, can travel through the core and associated structures, and reach the outside of the reactor vessel where a properly designed network of receivers can detect and interpret them. These receivers require a very small amount of power which, during loss of power events, can be provided for example by harvesting gamma radiation energy, thus resulting in a monitoring system that can function both during normal operation and during loss of power events. The project aimed at TAPS development through a series of tasks which are listed and briefly discussed as follows. TASK 1 – Sensor hardware design Subtask 1a: Assessment of sensor applications to SFRs Subtask 1b: Development of sensor functional requirements Subtask 1c: Definition of sensor hardware design specifications Task description: TAPS design was informed by considerations on their application (Subtask 1a), both the ultimate one in an SFR and the actual one in the ANL testing facilities that was intended to be used in support of the project. Considerations were made to identify optimum sensor design features that optimize the sensor size, materials, and output signal, for installation inside an SFR core. These considerations led to the development of Functional Requirements (Subtask 1b) and Design Requirements (Subtask 1c). TASK 2 – Sensor Hardware Manufacture Subtask 2a: Sensor hardware construction drawing development Subtask 2b: Sensor manufacture and assembly Task description: TAPS technical drawings were developed (Subtask 2a) using the Design Requirements established under Task 1. Subsequently, in spite of some problems which ultimately caused the program to be delayed, TAPS manufacturing was completed based on drawings (Subtask 2b). TASK 3 – Development of TAPS Signal Measurement System and TAPS Testing in Water Subtask 3a: Design, assembly and testing of signal measurement system, and TAPS testing in water Subtask 3b: Signal prediction-correction methodology development Task description: An assessment was performed on the techniques that can potentially be used to detect the signals emitted by the TAPS, e.g. a fiber-optic based acoustic signal measurement system, a laser vibrometer system, or an accelerometer-based system. The most suited technology, i.e. the accelerometer-based system, was developed further, and tested in water (Subtask 3a). Moreover, efforts were made to develop the methodology required to determine the actual system temperature and neutron flux distribution using differences between the measured and predicted TAPS responses (Subtask 3b). TASK 4 – Sensor System Testing in Sodium Subtask 4a: Test plan development Subtask 4b: Design, assembly and testing in small-scale sodium facility Subtask 4c: Design, assembly and testing in large-scale sodium and structures facility Task description: Upon proper test plan development (Subtask 4a), the fabricated TAPS was planned to be tested in sodium, by using two sodium facilities at ANL having different size and different purpose. The Under Sodium Viewing (USV) small-scale facility was intended to be used to investigate the effect of sodium on the sensor and its performance (Subtask 4b). The Mechanism Engineering Test Loop (METL) large-scale facility was instead intended to be used to assess the additional effect of prototypical SFR structures, such as fuel assembly mockup or parts of the core restrain structure, on sensor performance (Subtask 4c). As discussed in Section 3.2.2.7, unexpected issues during the TAPS manufacturing process resulted in some activities being delayed, with the TAPS and USV facility developed to the point to be ready for testing in sodium, however without the possibility to actually perform such testing (including the testing in METL) due to the end of the program’s performance period. Overall, through the development and testing (in water only) of two TAPS devices (a First-Generation TAPS followed by an optimized Second-Generation TAPS), the project confirmed the capability of this technology to generate acoustic signals proportional to temperature, which can be detected through a network of accelerometers identified as the best-suited type of receivers for acoustic signal detection. Moreover, the project also developed a computational model to predict the characteristics of the acoustic signals being generated, which combines thermal analysis of the TAPS with Finite Element Modeling (FEM)-aided acoustic characterization of the system. This model was benchmarked against experimental data collected during the project and, although general agreement was obtained, some limitations of the modeling methods were identified, which require additional development. Additional testing is needed in order to assess the effect, on TAPS operation and performance, of environmental changes resulting from the transition from water to liquid sodium. Such testing, which is suggested to be performed in the future, should look specifically at 1) both the effect resulting from the different thermoacoustic behavior of sodium (relative to water) and the effects of higher temperature on TAPS performance, and 2) the performance of the sensor-receiver system when multiple TAPS are used simultaneously and prototypical reactor structures are positioned in the testing environment. The latter testing is needed to assess the effects that potential signal attenuation/ distortion phenomena, as well as potential interference between signals emitted simultaneously, have on the performance of the technology for ultimate application in a nuclear reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1974

Ground-state energies of simple metals

A structural expansion for the static ground-state energy of a simple metal is derived. Two methods are presented, one an approach based on single-particle band structure which treats the electron gas as a nonlinear dielectric, the other a more general many-particle analysis using finite-temperature perturbation theory. The two methods are compared, and it is shown in detail how band-structure effects, Fermi-surface distortions, and chemical-potential shifts affect the total energy. These are of special interest in corrections to the total energy beyond third order in the electron-ion interaction and hence to systems where differences in energies for various crystal structures are exceptionally small. Preliminary calculations using these methods for the zero-temperature thermodynamic functions of atomic hydrogen are reported.

Hammerberg, J.↗
NASA NTRS1984

Static and dynamic structural-sensitivity derivative calculations in the finite-element-based Engineering Analysis Language (EAL) system

The implementation of static and dynamic structural-sensitivity derivative calculations in a general purpose, finite-element computer program denoted the Engineering Analysis Language (EAL) System is described. Derivatives are calculated with respect to structural parameters, specifically, member sectional properties including thicknesses, cross-sectional areas, and moments of inertia. Derivatives are obtained for displacements, stresses, vibration frequencies and mode shapes, and buckling loads and mode shapes. Three methods for calculating derivatives are implemented (analytical, semianalytical, and finite differences), and comparisons of computer time and accuracy are made. Results are presented for four examples: a swept wing, a box beam, a stiffened cylinder with a cutout, and a space radiometer-antenna truss.

Camarda, C. J.↗
NASA NTRS2006

Porous Media Approach for Modeling Closed Cell Foam

In order to minimize boil off of the liquid oxygen and liquid hydrogen and to prevent the formation of ice on its exterior surface, the Space Shuttle External Tank (ET) is insulated using various low-density, closed-cell polymeric foams. Improved analysis methods for these foam materials are needed to predict the foam structural response and to help identify the foam fracture behavior in order to help minimize foam shedding occurrences. This presentation describes a continuum based approach to modeling the foam thermo-mechanical behavior that accounts for the cellular nature of the material and explicitly addresses the effect of the internal cell gas pressure. A porous media approach is implemented in a finite element frame work to model the mechanical behavior of the closed cell foam. The ABAQUS general purpose finite element program is used to simulate the continuum behavior of the foam. The soil mechanics element is implemented to account for the cell internal pressure and its effect on the stress and strain fields. The pressure variation inside the closed cells is calculated using the ideal gas laws. The soil mechanics element is compatible with an orthotropic materials model to capture the different behavior between the rise and in-plane directions of the foam. The porous media approach is applied to model the foam thermal strain and calculate the foam effective coefficient of thermal expansion. The calculated foam coefficients of thermal expansion were able to simulate the measured thermal strain during heat up from cryogenic temperature to room temperature in vacuum. The porous media approach was applied to an insulated substrate with one inch foam and compared to a simple elastic solution without pore pressure. The porous media approach is also applied to model the foam mechanical behavior during subscale laboratory experiments. In this test, a foam layer sprayed on a metal substrate is subjected to a temperature variation while the metal substrate is stretched to simulate the structural response of the tank during operation. The thermal expansion mismatch between the foam and the metal substrate and the thermal gradient in the foam layer causes high tensile stresses near the metal/foam interface that can lead to delamination.

Ghosn, Louis J.↗
NASA NTRS2010

Hypersonic Navier-Stokes Comparisons to Orbiter Flight Data

During the STS-119 flight of Space Shuttle Discovery, two sets of surface temperature measurements were made. Under the HYTHIRM program3 quantitative thermal images of the windward side of the Orbiter with a were taken. In addition, the Boundary Layer Transition Flight Experiment 4 made thermocouple measurements at discrete locations on the Orbiter wind side. Most of these measurements were made downstream of a surface protuberance designed to trip the boundary layer to turbulent flow. In this paper, we use the US3D computational fluid dynamics code to simulate the Orbiter flow field at conditions corresponding to the STS-119 re-entry. We employ a standard two-temperature, five-species finite-rate model for high-temperature air, and the surface catalysis model of Stewart.1 This work is similar to the analysis of Wood et al . 2 except that we use a different approach for modeling turbulent flow. We use the one-equation Spalart-Allmaras turbulence model8 with compressibility corrections 9 and an approach for tripping the boundary layer at discrete locations. In general, the comparison between the simulations and flight data is remarkably good

Candler, Graham V.↗
NASA NTRS2016

Revisiting of Multiscale Static Analysis of Notched Laminates Using the Generalized Method of Cells

Composite material systems generally exhibit a range of behavior on different length scales (from constituent level to macro); therefore, a multiscale framework is beneficial for the design and engineering of these material systems. The complex nature of the observed composite failure during experiments suggests the need for a three-dimensional (3D) multiscale model to attain a reliable prediction. However, the size of a multiscale three-dimensional finite element model can become prohibitively large and computationally costly. Two-dimensional (2D) models are preferred due to computational efficiency, especially if many different configurations have to be analyzed for an in-depth damage tolerance and durability design study. In this study, various 2D and 3D multiscale analyses will be employed to conduct a detailed investigation into the tensile failure of a given multidirectional, notched carbon fiber reinforced polymer laminate. Threedimensional finite element analysis is typically considered more accurate than a 2D finite element model, as compared with experiments. Nevertheless, in the absence of adequate mesh refinement, large differences may be observed between a 2D and 3D analysis, especially for a shear-dominated layup. This observed difference has not been widely addressed in previous literature and is the main focus of this paper.

Materials Engineering↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗
DOE OSTI2020

Understanding the deformation mechanisms in Ni-based superalloys with using crystal plasticity finite element method

Nl-based superalloy is considered as a good candidate due to its excellent resistance to elevated temperature deformation for long term period application. Understanding the deformation and failure mechanisms of Ni-Based superalloys is very helpful for providing design guidelines for processing Ni-based superalloys. Experimental characterization indicates that the deformation mechanisms of Ni based superalloy is strongly microstructure dependent. Besides, damage transform from the void nucleation to the macro cracks by voids growth leading to the failure of the Ni-based superalloys are also showing strong microstructure sensitivity. Therefore, this work focuses on the prediction and comprehension of the deformation and void growth behavior in Ni based superalloy at different working conditions via crystal plasticity finite element modeling and simulation. Physically based crystal plasticity frameworks were developed for newly Ni- based superalloy Haynes 282. It was found that dislocation shearing through the precipitates were acting as the main contributor to the strength of Haynes 282 at room temperature and 815°C. Our analysis of the creeping behavior of Haynes 282 exhibited that resistance of general climb replaced by the resistance induced by the deposited climb dislocation density. In addition, in the study of void growth behavior, our simulation results demonstrated that as the main loading axis perpendicular to the grain boundary (GB), voids grow more slowly on tilt GBs in bicrystals than those in single and bicrystals with twist GBs. And tilt GBs would promote the void grow into irregular shape.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High-speed nanoscale characterization of dewetting via dynamic transmission electron microscopy

The dewetting of thin films can occur in either the solid or the liquid state for which different mass transport mechanisms are expected to control morphological changes. Traditionally, dewetting dynamics have been examined on time scales between several seconds to hours, and length scales ranging between nanometers and millimeters. The determination of mass transport mechanisms on the nanoscale, however, requires nanoscale spatial resolution and much shorter time scales. This study reports the high-speed observation of dewetting phenomena for kinetically constrained Ni thin films on crystalline SrTiO{sub 3} substrates. Movie-mode Dynamic Transmission Electron Microscopy (DTEM) was used for high-speed image acquisition during thin film dewetting at different temperatures. DTEM imaging confirmed that the initial stages of film agglomeration include edge retraction, hole formation, and growth. Finite element modeling was used to simulate temperature distributions within the DTEM samples after laser irradiation with different energies. For pulsed laser irradiation at 18 μJ, experimentally observed hole growth suggests that Marangoni flow dominates hole formation in the liquid nickel film. After irradiation with 13.8 μJ, however, the observations suggest that dewetting was initiated by nucleation of voids followed by hole growth through solid-state surface diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1973

Numerical solution of viscous reacting blunt body flows of a multicomponent mixture.

Hypersonic flow of a chemically nonequilibrium multicomponent mixture past an axisymmetric blunt body is considered on the basis of the exact Navier-Stokes and chemical rate equations. The flowfield, characterized by a sharp shock and a boundary layer, is solved numerically in the region between the shock and the body. The governing equations are first formulated in generalized orthogonal coordinates, and then recast in computational coordinates in which both the shock and the body are two straight lines and one coordinate is stretched to allow higher space resolution near the body. The thermodynamic and transport properties are obtained from the most rigorous theories and reliable data for each component and for the mixture itself. A time-marching second-order finite-difference method is used to solve the equations.

Li, C. P.↗
NASA NTRS1988

The origin of pristine KREEP - Effects of mixing between UrKREEP and the magmas parental to the Mg-rich cumulates

The fact that pristine KREEP basalts generally have bulk-rock molar MgO/(MgO+FeO) ratios lower than average for the lunar crust despite their extraordinarily high incompatible element contents is discussed. Pristine KREEP basalts also have primitive Ni contents compared to mare basalts with comparable REE contents. It is suggested that these facts can be explained through mixing between ancient KREEP precursor materials and primitive Mg-rich melts. Finite-difference models of this mixing, followed by anorthosite assimilation and fractional crystallization provide satisfactory fits to the composition of these basalts. It is suggested that this is a further confirmation of the magmasphere hypothesis.

Warren, Paul H.↗
NASA NTRS1992

Conical detonation waves - A comparison of theoretical and numerical results

A predictive capability for conical oblique detonation waves (ODWs) is discussed with particular attention given to the special case of an axisymmetric conical projectile, which yields necessary conditions for attached waves. A multidimensional numerical model was used to study the influence of a finite Damkoehler number (Da) and viscosity on the formation of a conical ODW. The ODW is predicted for Da on the order of 10, and is found to be important in the ODW formation. For a Da of 1, shock induced combustion is predicted. Simulations of the experimental geometry showed general agreement with experimental flow features, but underpredicted the shock angle. It is suggested that this difference may be caused by inadequate boundary layer resolution, lack of grid orthogonality at the cone surface, and an inappropriate chemical mechanism.

Fort, James A.↗