Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2019

High-Frequency Acoustic Mode Identification of Unstable Combustors

Abstract High frequency thermoacoustic instabilities are becoming increasingly problematic in modern combustion systems. Understanding which acoustic mode is being excited is important for understanding potential mechanisms and control approaches — for example, influence of a helical shear layer mode on the flame has profoundly different effects on the first tangential acoustic mode, than a higher order axial mode of similar frequency. Nonetheless, the modal density increases with frequency and it becomes increasingly difficult to determine which acoustic mode is self-excited, based upon frequency calculations alone. Moreover, access issues and cost usually limit the number of pressure probes that can be distributed axially and azimuthally in the combustor. This paper presents a methodology for identifying the acoustic mode by using high temperature pressure transducers flush mounted in a combustion chamber. Modal identification is demonstrated with a siren under non-reacting conditions. The siren is mounted on the chamber to excite longitudinal and azimuthal waves. Five acoustic sensors at different axial and azimuthal locations measure the pressure fluctuations simultaneously. Given the forcing frequency and the speed of sound, the pressure distribution in the combustor is reconstructed in the time domain from the measured data by using a least squares method to determine its mode shapes. In addition, the finite element method (FEM) solver is used to provide the eigenfrequencies and corresponding mode shapes. The test results demonstrate that the mode shapes from the reconstructed data and corresponding frequencies are consistent with those predicted from the FEM, which validates the methodology in this study. In addition, the methodology is extended to practical reacting cases without the siren to determine the acoustic mode shapes of naturally occurring instabilities. In these cases, the modal features have strong stochastic features, such as what appear to be stochastic variations in overall amplitude and relative amplitudes of clockwise and counterclockwise waves.

Kim, J.↗
DOE OSTI2018

Performance of b Jet Identification at $\sqrt{s}$ = 13 TeV with the ATLAS Detector at CERN [Slides]

The efficient identification of jets from bottom quarks (b-jets) is one of the most important techniques for many physics analyses at the Large Hadron Collider, including studies of the Higgs boson, the top quark, and searches beyond the Standard Model. The performance is characterized by b-tagging efficiency (probability to identify a b-jet as such) and the mistag rate (probability to mistakenly accept a non-b-jet).The mistag occurs as a result of finite detector resolution, presence of long-lived particles, and material interactions. As these effects can be different between the experimental data and Monte Carlo (MC) simulation, it is important to measure the b-tagging performance in data and derive the MC correction factors. I will describe various methods to measure the performance of b jet identification in pp collisions at √s = 13 TeV with the ATLAS detector and will report on Direct Tag method which is recently developed in the ATLAS collaboration for Light jet calibration.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Compact finite difference schemes for the time fractional diffusion equation with nonlocal boundary conditions

Fourth-order compact finite difference schemes for solving the time fractional diffusion equation with nonlocal boundary conditions are considered in this paper. The second-order derivative in one space variable is approximated by a fourth-order compact finite difference. Two different time derivatives, the Caputo derivative and Riemann–Liouville derivative are considered in the equations that we discuss. The Caputo time derivative is discretized by the L1-approximation and a second-order approximation is used for the Riemann–Liouville derivative. Stability analysis using the discrete energy method for Scheme I is given, and error estimation of this scheme is derived by conducting the local truncation error analysis, using the stability result. Numerical results are provided to show that both two fully discrete schemes have fourth-order accuracy in space, and 2−γ order of accuracy in time for Scheme I with γ being the order of the fractional derivative, while Scheme II has second-order accuracy in time when the theoretical solutions are smooth with zero initial values.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Inverse aqueous transport modeling for emergency response

ALGE is a three-dimensional, finite-difference aqueous transport model that simulates pollutant fate and transport in lakes, rivers, bays, and estuaries by solving the prognostic equations of mass, momentum, and energy. Its current modeling capabilities include transport of dissolved tracer for a series of predefined basins across the continental United States. Recently, an inverse method (also known as backtracking) has been added to ALGE to provide a possible source of a pollutant should one be detected by a sensor in a body of water and a source is not known. This inverse method is a three step process that uses an algorithm to inverse the flow. We demonstrate the new model’s capabilities through simulating the 2021 Piney Point spill in Tampa Bay, Florida (USA). This involves moving tracer backwards from its detection points, encompassing a potential source area, and applying Bayes’ Theorem and $\frac{𝜒}{𝑄}$ to reduce the area within which the true source could be located.

hydrological modeling↗
DOE OSTI2019

Stochastic multiscale flux basis for Stokes-Darcy flows

Three algorithms are developed for uncertainty quantification in modeling coupled Stokes and Darcy flows. The porous media may consist of multiple regions with different properties. The permeability is modeled as a non-stationary stochastic variable, with its log represented as a sum of local Karhunen-Loève (KL) expansions. The problem is approximated by stochastic collocation on either tensor-product or sparse grids, coupled with a multiscale mortar mixed finite element method for the spatial discretization. A non-overlapping domain decomposition algorithm reduces the global problem to a coarse scale mortar interface problem, which is solved by an iterative solver, for each stochastic realization. In the traditional domain decomposition implementation, each subdomain solves a local Dirichlet or Neumann problem in every interface iteration. To reduce this cost, two additional algorithms based on deterministic or stochastic multiscale flux basis are introduced. The basis consists of the local flux (or velocity trace) responses from each mortar degree of freedom. It is computed for each subdomain independently before the interface iteration begins. The use of the multiscale flux basis avoids the need for subdomain solves on each iteration. The deterministic basis is computed at each stochastic collocation and used only at this realization. The stochastic basis is formed by further looping over all local realizations of a subdomain's KL region before the stochastic collocation begins. It is used over multiple realizations. Numerical tests are presented to illustrate the performance of the three algorithms, with the stochastic multiscale flux basis showing significant savings in computational cost compared to the other two algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗
DOE OSTI2022

Reduced-Order Modeling Method for Phase-Change Thermal Energy Storage Heat Exchangers

Thermal energy storage can facilitate the effective utilization of renewable energy. To speed up the design process of thermal energy storage devices, it is critical to develop fast and accurate modeling methods for phase change material embedded heat exchangers (PCM HXs). This study developed and compared two approximation-assisted reduced-order PCM HX models for the simulation of thermal storage components and systems, which were verified against a validated finite-volume model. They are a pure black-box model, and a grey-box model based on the Number of Transfer Units (or effectiveness-NTU) approach. We used the reduced-order models to predict the performance of standalone PCM HX for 1000 cases with different fluid inlet conditions and HX designs. We also integrated the models with a vapor compression system to predict the compressor energy consumption and total charging time under various conditions. The results show that overall, the black-box model gave more accurate results than the grey-box model. On average, the mean absolute deviation in the PCM HX fluid outlet temperature was 0.05 K and 0.1 K for the black-box model and grey-box model, respectively. The grey-box model gave larger temperature deviation toward the end of the phase change process, due to the simplified two-node PCM representation. On the system level, the mean absolute deviation in compressor energy consumption was 0.2% and 0.3% for the black-box model and grey-box model, respectively. Their mean absolute deviation in total charging time was 1.1% and 2.6%, respectively. In terms of computation efficiency, the system simulation speedup ratio gained by using the reduced-order models was 11 to 57. System simulation time was decreased from an average of 1465 s to an average of 59 s. This shows that the proposed reduced-order modeling methods can be used to predict system performance of PCM HXs with less than 3% accuracy penalty, and 25 times less computational time than finite-volume models, enabling faster design and evaluation of PCM thermal storage devices.

25 ENERGY STORAGE↗
DOE OSTI2025

Tailored Frequency Conversion in Nonlinear Subwavelength Grating Metaslabs

The use of metaslab waveguides consisting of subwavelength gratings (SWGs) is proposed as a highly flexible approach to obtain phase‐matched frequency conversion in nonlinear integrated devices. Control over the phase matching (PM) condition is achieved through modification of the linear modal propagation constants by the SWG. The method is experimentally demonstrated by fabrication of 0.4 mm‐long GaAs metaslabs designed for second harmonic generation of near‐infrared light. A controlled spectral shift of up to 100 nm in the PM condition, relative to a uniform slab is observed. Excellent agreement is found between the experimental results and finite‐difference eigenmode simulations, which predict a tunable frequency conversion bandwidth of ∼300 nm for the design parameters. The different regimes of PM control supported by the structure are discussed. This work highlights the potential of SWGs as building blocks for tailored frequency conversion in nonlinear integrated devices.

Form birefringence↗
DOE OSTI2019

Influence of Single-Point Incremental Force Process Parameters on Forming Characteristics and Microstructure Evolution of AA-6061 Alloy Sheet

Single-point incremental forming (SPIF) process comprises a set of process variables, such as forming tool diameter, vertical pitch, spindle rotation speed, and tool velocity, which may affect the forming behavior of the sheet to be deformed. The objective of this work is to study the effect of SPIF process parameters on forming characteristics and microstructure development for the AA-6061 (T6) aluminum alloy sheet. The SPIF experiments and finite element (FE) simulations were performed at different process parameters to achieve conical shapes from the AA-6061 blanks. The effect of process parameters on forming forces, thickness uniformity in wall region of the cone and surface roughness of the blank was analyzed. A detailed microstructure study was performed to analyze the effect of process variables on microstructure and texture evolution during the SPIF process. This study reveals that the process parameters are likely to influence the texture development especially at high tool diameter and vertical pitch values. Therefore, suitability and consequences of using different combinations of tool diameter and vertical pitch values are discussed for AA-6061 alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A finite element procedure for radio-frequency sheath–plasma interactions based on a sheath impedance model

A finite element code that solves self-consistent radio-frequency (RF) sheath-plasma interaction problems is improved by incorporating a generalized sheath boundary condition in the macroscopic solution scheme. This sheath boundary condition makes use of a complex sheath impedance including both the sheath capacitance and resistance, which enables evaluation of not only the RF voltage across the sheath but also the power dissipation in the sheath. The newly developed finite element procedure is applied to cases where the background magnetic field is perpendicular to the sheath surface in one- and two-dimensional domains filled by uniform low- and high-density plasmas. The numerical results are compared with those obtained by employing the previous capacitive sheath model at a typical frequency for ion cyclotron heating used in fusion experiments. It is shown that for sheaths on the order of 100 V in a high-density plasma, localized RF power deposition can reach a level which causes material damage. It is also shown that the sheath-plasma wave resonances predicted by the capacitive sheath model do not occur when parameters are such that the generalized sheath impedance model substantially modifies the capacitive character of the sheath. Here, possible explanations for the difference in the maximum RF sheath voltage depending on the plasma density are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Reversible Switch in Charge Storage Enabled by Selective Ion Transport in Solid Electrolyte Interphase

Solid-electrolyte interphases (SEIs) in advanced rechargeable batteries ensure reversible electrode reactions at extreme potentials beyond the thermodynamic stability limits of electrolytes by insulating electrons while allowing working ions to transport. Such selective ion transport occurs naturally in biological cell membranes as a ubiquitous prerequisite of many life processes and a foundation of biodiversity. In addition, cell membranes can selectively open and close the ion channels in response to external stimuli (e.g., electrical, chemical, mechanical, thermal), giving rise to “gating” mechanisms that help manage intracellular reactions. We wondered whether the chemistry and structure of SEIs can mimic cell membranes, such that ion gating can be replicated. That is, can SEIs realize a reversible switching between two electrochemical behaviors, i.e., the ion intercalation chemistry of batteries and the ion adsorption of capacitors? Herein, we report such SEIs that result in thermally activated selective ion transport. The function of open/close gate switches is governed by the chemical and structural dynamics of SEIs under different thermal conditions, with precise behaviors as conducting and insulating interphases that enable battery and capacitive processes within a finite temperature window. Such an ion gating function is synergistically contributed by Arrhenius-activated ion transport and SEI dissolution/regrowth. Following the understanding of this new mechanism, we then develop an electrochemical method to heal the SEI layer in situ. As a result, the knowledge acquired in this work reveals the possibility of hitherto unknown biomimetic properties of SEIs, which will guide us to leverage such complexities to design better SEIs for future battery chemistries.

25 ENERGY STORAGE↗
DOE OSTI2017

Improved spatial resolution for spot sampling in thermal desorption atomic force microscopy – mass spectrometry via rapid heating functions

The key to advancing materials is to understand and control their structure and chemistry. However, thorough chemical characterization is challenging since existing techniques characterize only a few properties of the specimen, thereby necessitating multiple measurement platforms to acquire the necessary information. The multimodal combination of atomic force microscopy (AFM) and mass spectrometry (MS) transcends existing analytical capabilities for nanometer scale spatially resolved correlation of the chemical and physical properties of a sample surface. One such hybrid system employs heated AFM cantilevers for thermal desorption (TD) sampling of molecules from a surface and subsequent gas phase ionization and detection of the liberated species by MS. Here in this paper, we report on the use of voltage pulse trains to tailor cantilever heating such that spot sampling size was reduced and mass spectral signal was improved compared to constant voltage, static heating of the cantilever. Desorption efficiency (DE), defined as the quotient of the mass spectral signal intensity and the volume of the desorption crater, was used to judge the effectiveness of a particular tailored heating function. To guide the development and optimization of the heating functions and aid in interpreting experimental results, a 1D finite element model was developed that predicted the cantilever response to different heating functions. Three tailored heating functions that used different combinations, magnitudes, and durations of rectangular voltage pulses, were used for surface spot sampling. The resultant sampling spot size and DE were compared to the same metrics obtained with the conventional method that uses a single voltage pulse. Using a model system composed of a thin film of ink containing pigment yellow 74 as a model system, desorption craters shrunk from 2 μm, using the conventional approach, to 310 nm using the optimum tailored heating function. This same pulsed heating function produced a 381× improvement in the DE and an 8× improvement in spatial resolution compared to the conventional heating approach showing that signal/amount of material sampled was improved significantly by this new cantilever heating strategy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Planetary normal mode computation: Parallel algorithms, performance and reproducibility

This report is an extension of work entitled “Computing planetary interior normal modes with a highly parallel polynomial filtering eigensolver.” by Shi et al., [1] originally presented at the SC18 conference. A highly parallel polynomial filtered eigensolver was developed and exploited to calculate the planetary normal modes. The proposed method is ideally suited for computing interior eigenpairs for large-scale eigenvalue problems as it greatly enhances memory and computational efficiency. In this article, the second-order finite element method is used to further improve the accuracy as only the first-order finite element method was deployed in the previous work. The parallel algorithm, its parallel performance up to 20k processors, and the great computational accuracy are illustrated. The reproducibility of the previous work was successfully performed on the Student Cluster Competition at the SC19 conference by several participant teams using a completely different Mars-model dataset on different clusters. Both weak and strong scaling performances of the reproducibility by the participant teams were impressive and encouraging. The analysis and reflection of their results are demonstrated and future direction is discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

THERMAL SAFETY MARGIN CALCULATION OF THE MP-2 EXPERIMENT IN THE ADVANCED TEST REACTOR

The Mini-Plate 2 (MP-2) irradiation test is a fueled experiment designed for irradiation in multiple test locations in the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). The experiment is a drop-in test where small aluminum-clad fuel plate samples (mini plates) are cooled directly by the ATR Primary Coolant System (PCS) water. The MP-2 fuel plate experiment will be irradiated in several different irradiation locations of the ATR. This fueled experiment contains aluminum-clad fuel mini plates consisting of monolithic U-Mo. Four different types of fuel plates with fuel meat thickness and cladding are part of the MP-2 test. A thermal analysis has been performed on the MP-2 experiment. A method for calculating Departure from Nucleate Boiling Ratio (DNBR) and Flow Instability Ratio (FIR) during a reactivity transient using the commercial finite element and heat transfer code ABAQUS is discussed. At the start of an ATR cycle the heat generation rate of the fueled experiment is high and the heat rate multiplier from the outer shim control cylinders is low, while the reverse is true at the end of the ATR cycle. Thermal analyses at 10 day increments during the cycle calculate the DNBR and FIR during a reactivity transient. This technique calculates DNBR for the fuel plate surfaces and FIR for all water components for each finite element surface and node at various times during the ATR cycle. Heat rates vary with time during the transient calculation that are provided by a detailed physics analysis. Oxide growth on the fuel plates is also incorporated Results from the transient calculations are displayed with the ABAQUS post processor. By calculating these parameters at each location in the finite element model, conservatism is replaced with accuracy. This allows for a greater margin for the thermal hydraulic safety parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

PYRANDA

Pyranda Is fortran 90 code which solves compact finite difference equations on a three dimensional fortran arrays, domain decomposed using the Message Passing Interface (MPI). The fortran source solves the linear system of equations which forms a pentadtagona1 matrix. Derivatives are used to construct hyperbolic partlal differential equation (PDE) operators, namely gradient, divergence, Laplacian, and curl operators. 5. Solution Method: The library provides a python Interface which facilitates problem setup and description. However, the core solver Is written in the fortran 90 language. The solver couples MPI communication and linear solves together so that a direct (non-iterative) solution to the linear system of equations can be computed with minimal communication. Domain decomposition is limited to a uniform data size but provides perfect load balancing. 6. Computer for which the code Is written:

Olson, BrittonJ.↗
DOE OSTI2020

Thermal Model Heat Rate Predictions of the AGR-5/6/7 Experiment

Several fuel and material irradiation experiments have been planned for the U.S. Department of Energy Advanced Gas Reactor Fuel Development and Qualification Program, which supports the development and qualification of tristructural isotropic (TRISO) coated particle fuel for use in high-temperature gas-cooled reactors. The goals of these experiments are to provide irradiation performance data to support fuel process development, qualify fuel for normal operating conditions, support development of fuel performance models and codes, and provide irradiated fuel and materials for post-irradiation examination and safety testing. Originally planned and named as separate fuel experiments, but subsequently combined into a single test train. AGR-5/6/7 will test low-enriched uranium oxycarbide (UCO – a heterogeneous mixture of uranium oxide and uranium carbide) TRISO fuel. The AGR-5/6 portion of the experiment will provide data to support qualification of the selected reference fuel design, while the AGR-7 portion will serve as a margin test, irradiating the fuel beyond normal operating conditions. During irradiation, the temperature in each capsule is controlled by varying the helium / neon gas mixture flowing through the capsule until the desired thermocouple temperature is reached. Further adjustments can be made by changing the lobe power and the neutron filter. Three neutron filters are available to use as the experiment burns up the fissionable fuel. A thermal finite element model has been created for the five capsules comprising the AGR-5/6/7 experiment. Heat rates from a physics analysis are imported into the model along with fast neutron fluence. Graphite shrinkage due to the fast neutron fluence is incorporated into the model. Gas gaps change as a function of fast neutron fluence. Radiation heat transfer is a major contributor in this model. This is a large model with more than 1 million finite element brick elements. More than 150 different parts are modeled in the finite element model and communicate with each other from a heat transfer sense. More than 50 thermocouples are used in the experiment and are compared to actual measurements. The experiment is composed of five separate stainless steel capsules all welded together. This paper presents a method used to predict the gas mixture in the event that all of the thermocouples fail in a particular capsule. Curve fitting of heat rates from previous cycles and projecting them into the future is the basis of the method.

36 MATERIALS SCIENCE↗
DOE OSTI2018

A hybrid finite element-spectral boundary integral approach: Applications to dynamic rupture modeling in unbounded domains

The finite element method (FEM) and the spectral boundary integral method (SBI) have both been widely used in the study of dynamic rupture simulations along a weak interface. Here, we present a hybrid method that combines FEM and SBI through the consistent exchange of displacement and traction boundary conditions, thereby benefiting from the flexibility of FEM in handling problems with nonlinearities or small-scale heterogeneities and from the superior performance and accuracy of SBI. We validate the hybrid method using a benchmark problem from the Southern California Earthquake Center's dynamic rupture simulation validation exercises.We further demonstrate the capability and computational efficiency of the hybrid scheme for resolving off-fault heterogeneities by studying a 2D in-plane shear crack in two different settings: one where the crack is embedded in a high-velocity zone and another where it is embedded in a low-velocity zone. Finally, we discuss the potential of the hybrid method for addressing a wide range of problems in geophysics and engineering.

58 GEOSCIENCES↗