Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Method and Apparatus for Measuring Near-Angle Scattering of Mirror Coatings

Disclosed herein is a method of determining the near angle scattering of a sample reflective surface comprising the steps of: a) splitting a beam of light having a coherence length of greater than or equal to about 2 meters into a sample beam and a reference beam; b) frequency shifting both the sample beam and the reference beam to produce a fixed beat frequency between the sample beam and the reference beam; c) directing the sample beam through a focusing lens and onto the sample reflective surface, d) reflecting the sample beam from the sample reflective surface through a detection restriction disposed on a movable stage; e) recombining the sample beam with the reference beam to form a recombined beam, followed by f) directing the recombined beam to a detector and performing heterodyne analysis on the recombined beam to measure the near-angle scattering of the sample reflective surface, wherein the position of the detection restriction relative to the sample beam is varied to occlude at least a portion of the sample beam to measure the near-angle scattering of the sample reflective surface. An apparatus according to the above method is also disclosed.

Chipman, Russell A.↗
DOE OSTI2021

A Continuous-Discontinuous Galerkin Method for Electromagnetic Simulations Based on an All-Frequency Stable Formulation

In this paper, a potential-based partial-differential formulation, called the all-frequency stable formulation, is presented for the accurate and robust simulation of electromagnetic problems at all frequencies. Due to its stability from (near) dc to microwave frequencies, this formulation can be applied to simulate wide-band and multiscale problems without encountering the infamous low-frequency breakdown issue or applying basis function decompositions such as the tree-cotree splitting technique. To provide both efficient and flexible numerical solutions to the electromagnetic formulation, a mixed continuous-discontinuous Galerkin (CDG) method is proposed and implemented. In regions with homogeneous media, the continuous Galerkin method is employed to avoid the introduction of duplicated degrees of freedom (DoFs) on the elemental interfaces, while on the interfaces of two different media, the discontinuous Galerkin method is applied to permit the jump of the normal components of the electromagnetic fields. Numerical examples are provided to validate and demonstrate the proposed numerical solver for problems in a wide electromagnetic spectrum.

Yan, Su↗
DOE OSTI2019

Modulating Carrier Transport via Defect Engineering in Solar Water Splitting Devices

Developing sustainable solar water splitting devices requires efficient separation and transport of photogenerated carriers. In this Perspective, we examine carrier transport in semiconductor photoelectrodes using bismuth vanadate as a primary model system and highlight strategies to significantly improve carrier delivery through defect engineering. To improve electron transport in low-mobility semiconductors, we introduce two distinct bulk doping methods, by homogeneously enhancing the bulk defect level or by forming distributed homojunctions with graded doping. Next, we demonstrate the use of structural boundaries as extrinsic pathways for fast electron transport, thus providing novel insights to engineer materials’ macroscopic conductivity. Third, we describe the importance of interface design in terms of structural, chemical, and electronic matching at the back contact to suppress carrier recombination. Lastly, we highlight the methods for surface defect control via passivation and catalysis and give a brief outlook of research challenges and opportunities for solar water splitting.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Initial Split Disk Mechanical Testing Studies with a Digital Image Correlator

Nuclear graphite plays a pivotal role in the advanced reactor program. Consequently, developing a viable method to measure and evaluate small samples is essential. The Digital Image Correlation (DIC) technique allows for the examination of much smaller components in real-time applications, particularly concerning reactor capsule design and placement within the reactor. These samples become crucial for understanding the behavior of nuclear graphite before and after radiation exposure. They serve as a means to analyze changes and effectively measure the mechanical strength properties, which may or may not alter following irradiation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Effect of Sample Thickness on the Tensile Strength of Small Graphite Discs

This report formally documents the completion of the Advanced Reactor Technologies (ART) Level 3 Milestone (M3TG-24OR0501054), “Continue activities related to Split Disc-DIC - complete analysis of effect of sample thickness on one fine grain graphite,” due May 31, 2024. Details within this report outline the status of activities aimed at elucidating the effect of sample thickness on the tensile strength measurement using the small-disc, split-disc testing standard in support of the US Department of Energy’s ART Graphite R&D Program. The ASTM D8289, Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite, was developed to provide an alternative means for testing tensile strength on smaller specimens, which are compatible with available irradiation capsule volumes, and other compatibility measurements. ASTM D8289 specifies that acceptable specimen diameter can range from 6 to 12.7 mm and that the maximum allowed sample thickness should be half of the diameter. However, information is limited with respect to the effect of thickness on the measured splitting tensile strength. This report documents efforts to understand the effect of sample thickness on splitting tensile strength. The work involved testing Ø12.7 mm samples of fine-grain graphites 2114 and IG-110 of different thicknesses (6.35, 5, 4, and 3 mm). The digital image correlation (DIC) method was applied to the samples, along with the ASTM D8289 Standard, to help interpret the measured results.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A comprehensive modeling method for proton exchange membrane electrolyzer development

Hydrogen attracts significant interests as an effective energy carrier that can be derived from renewable sources. Hydrogen production using a proton-exchange membrane (PEM) electrolyzer can efficiently convert renewable power via water splitting in wide scales—from large, centralized generation to on-site production. Mathematical models with multiple scales and fidelities facilitate the continuing improvements of PEM electrolyzer development to improve performance, cost, and reliability. The model scopes and methods are presented in this paper, which also introduces a comprehensive PEM electrolysis modeling tool based on computational fluid dynamics (CFD) software, ANSYS/Fluent. The modeling tool incorporates electrochemical model of a PEM electrolysis cell to simulate the performance of coupled thermal-fluid, species transport, and electrochemical processes in a product-scale cell or stack by leveraging the powerful meshing generation and CFD solver of ANSYS/Fluent. The thermal-fluid modeling includes liquid water/gas two-phase flow and simulates a PEM electrolysis cell by using Fluent user-defined functions as add-on modules accounting for PEM-specific species transport and electrochemical processes. The modeling outcomes expediate PEM electrolyzer scaling up from basic material development and laboratory testing.

08 HYDROGEN↗
DOE OSTI2026

Investigation of nonlinear mechanical behavior of two superfine-grained graphites with DIC assisted disc splitting test

Using standardized uniaxial tensile test specimens is not practical due to the limited volumes in irradiation capsules, molten salt degradation facilities, or oxidation apparatus. The ASTM International Standard Test Method for Tensile Strength Estimate by Disc Compression of Manufactured Graphite (ASTM D8289) was developed to provide a convenient way for estimating tensile strength. Unlike the traditional uniaxial tensile test in ASTM C749 (International Standard Test Method for Tensile Stress–Strain of Carbon and Graphite), which uses dog-bone shaped specimens larger than 12.95 mm × 120.65 mm, the ASTM D8289 standard uses smaller discs with diameters of 6–12.7 mm. A digital image correlation (DIC) system, which uses a full-field noncontact surface displacement measurement technique, was applied along with the ASTM D8289 disc splitting test on IG-110 samples and compared with previous measurements from Mersen 2114 samples. Results confirmed that the DIC technique can measure the surface displacement/strain on these small (Ø6 mm × 3 mm) graphite specimens with good repeatability. However, Mersen 2114 and IG-110 samples exhibited strain discrepancies when DIC measurements were compared with analytical and finite element simulation values. The loading history and strain results also indicated different mechanical behaviors between Mersen 2114 and IG-110, particularly the nonlinear behavior of the IG-110 samples. Good agreement was observed by comparing the splitting tensile strengths of two superfine-grained grades with results from other work. The specimen size effect is discussed when comparing the splitting tensile strength with corresponding uniaxial tensile strength of these two graphite grades.

Lin, Lianshan [ORNL] (ORCID:0000000203399219)↗
DOE OSTI2019

Ultrastrong coupling of vibrationally dressed organic Frenkel excitons with Bloch surface waves in a one-sided all-dielectric structure

We demonstrate a transition from weak, to strong, to ultrastrong coupling of Frenkel molecular excitons and Bloch surface wave photons at room temperature using a one-sided, all-dielectric optical structure. The all-dielectric structure comprises an organic semiconductor thin film on the surface of a distributed Bragg reflector. We investigate the evolution of multiple vibronic polariton branches and their dominant absorption peaks as a function of coupling and in-plane momentum, which is absent in previous ultrastrong coupling systems. Measurements are interpreted using both the transfer matrix method and a coupled-oscillator model without the rotating wave approximation. The dependence of Rabi splitting on the number of excitons and electrical field amplitude is also modeled showing a transition to ultrastrong coupling at film thicknesses ≥ 50 nm. Furthermore, this low-loss polaritonic structure enables us to study phenomena such as organic exciton-polariton dynamics, ultralong range polariton propagation, and high efficiency energy transport in the ultrastrong coupling regime at room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1977

A computerized traffic control algorithm to determine optimal traffic signal settings

An algorithm was developed to optimally control the traffic signals at each intersection using a discrete time traffic model applicable to heavy or peak traffic. Off line optimization procedures were applied to compute the cycle splits required to minimize the lengths of the vehicle queues and delay at each intersection. The method was applied to an extensive traffic network in Toledo, Ohio. Results obtained with the derived optimal settings are compared with the control settings presently in use.

Seldner, K.↗
NASA NTRS1977

Backscattering from a Gaussian distributed, perfectly conducting, rough surface

The problem of scattering by random surfaces possessing many scales of roughness is analyzed. The approach is applicable to bistatic scattering from dielectric surfaces, however, this specific analysis is restricted to backscattering from a perfectly conducting surface in order to more clearly illustrate the method. The surface is assumed to be Gaussian distributed so that the surface height can be split into large and small scale components, relative to the electromagnetic wavelength. A first order perturbation approach is employed wherein the scattering solution for the large scale structure is perturbed by the small scale diffraction effects. The scattering from the large scale structure is treated via geometrical optics techniques. The effect of the large scale surface structure is shown to be equivalent to a convolution in k-space of the height spectrum with the following: the shadowing function, a polarization and surface slope dependent function, and a Gaussian factor resulting from the unperturbed geometrical optics solution. This solution provides a continuous transition between the near normal incidence geometrical optics and wide angle Bragg scattering results.

Brown, G. S.↗
NASA NTRS1990

Matrix-based, finite-difference algorithms for computational acoustics

A compact numerical algorithm is introduced for simulating multidimensional acoustic waves. The algorithm is expressed in terms of a set of matrix coefficients on a three-point spatial grid that approximates the acoustic wave equation with a discretization error of O(h exp 5). The method is based on tracking a local phase variable and its implementation suggests a convenient coordinate splitting along with natural intermediate boundary conditions. Results are presented for oblique plane waves and compared with other procedures. Preliminary computations of acoustic diffraction are also considered.

Davis, Sanford↗
NASA NTRS2016

The Relaxation Matrix for Symmetric Tops with Inversion Symmetry. I. Effects of Line Coupling on Self-Broadened v (sub 1) and Pure Rotational Bands of NH3

The Robert-Bonamy formalism has been commonly used to calculate half-widths and shifts of spectral lines for decades. This formalism is based on several approximations. Among them, two have not been fully addressed: the isolated line approximation and the neglect of coupling between the translational and internal motions. Recently, we have shown that the isolated line approximation is not necessary in developing semi-classical line shape theories. Based on this progress, we have been able to develop a new formalism that enables not only to reduce uncertainties on calculated half-widths and shifts, but also to model line mixing effects on spectra starting from the knowledge of the intermolecular potential. In our previous studies, the new formalism had been applied to linear and asymmetric-top molecules. In the present study, the method has been extended to symmetric-top molecules with inversion symmetry. As expected, the inversion splitting induces a complete failure of the isolated line approximation. We have calculated the complex relaxation matrices of selfbroadened NH3. The half-widths and shifts in the ν1 and the pure rotational bands are reported in the present paper. When compared with measurements, the calculated half-widths match the experimental data very well, since the inapplicable isolated line approximation has been removed. With respect to the shifts, only qualitative results are obtained and discussed. Calculated off-diagonal elements of the relaxation matrix and a comparison with the observed line mixing effects are reported in the companion paper (Paper II).

Trajectory models↗
DOE OSTI2024

Divide and conquer: Learning chaotic dynamical systems with multistep penalty neural ordinary differential equations

Forecasting high-dimensional dynamical systems is a fundamental challenge in various fields, such as geosciences and engineering. Neural Ordinary Differential Equations (NODEs), which combine the power of neural networks and numerical solvers, have emerged as a promising algorithm for forecasting complex nonlinear dynamical systems. However, classical techniques used for NODE training are ineffective for learning chaotic dynamical systems. In this work, we propose a novel NODE-training approach that allows for robust learning of chaotic dynamical systems. Here, our method addresses the challenges of non-convexity and exploding gradients associated with underlying chaotic dynamics. Training data trajectories from such systems are split into multiple, non-overlapping time windows. In addition to the deviation from the training data, the optimization loss term further penalizes the discontinuities of the predicted trajectory between the time windows. The window size is selected based on the fastest Lyapunov time scale of the system. Multi-step penalty(MP) method is first demonstrated on Lorenz equation, to illustrate how it improves the loss landscape and thereby accelerates the optimization convergence. MP method can optimize chaotic systems in a manner similar to least-squares shadowing with significantly lower computational costs. Our proposed algorithm, denoted the Multistep Penalty NODE, is applied to chaotic systems such as the Kuramoto-Sivashinsky equation, the two-dimensional Kolmogorov flow, and ERA5 reanalysis data for the atmosphere. It is observed that MP-NODE provide viable performance for such chaotic systems, not only for short-term trajectory predictions but also for invariant statistics that are hallmarks of the chaotic nature of these dynamics.

Chaotic dynamical systems↗
DOE OSTI2021

DEM analysis of the effect of interface transition zone on dynamic splitting tensile behavior of high-strength concrete based on multi-phase model

Highlights: • Crushable aggregates with the real geometry shape are generated by using clump-cluster method. • Flat-Joint Model is adopted as the constitutive law. • The tensile strength is almost not affected by ITZ strength when its strength exceeds 0.8 times that of mortar. • A modified formula that can reflect the functional relationship between the ITZ strength and the TDIF is proposed. The tensile strength is one of the main parameters in the design of concrete structures and has an important influence on the cracking performance and durability of concrete materials. As the weak region in concrete, the strength of the interface transition zone (ITZ) directly affects the dynamic mechanical response of concrete. Meanwhile, considering that the tensile performance of concrete is far weaker than its compressive performance, the study of the strength of ITZ is of great significance to the dynamic tensile fracture characteristics and mechanical properties of concrete. In this paper, a method combining 3D scanning technology and Clump-Cluster was proposed to generate a real geometry shape model of crushable aggregate. On this basis, the 3D discrete element model (DEM) of high-strength concrete was constructed using the Flat-Joint Model (FJM). The effects of the ITZ strength on the failure mode, number of microcracks, tensile strength and strain rate effect of high-strength concrete under strain rates from 0.01/s to 100/s were studied by splitting tensile test. The simulation results indicate that the number of microcracks generated decreases with the increase of ITZ strength, and the proportion of ITZ cracks decreases obviously. When the ITZ strength exceeds the critical value of 0.8 times that of mortar, the tensile strength is almost not affected by ITZ strength. In addition, based on the numerical simulation results, a modified formula that can reflect the functional relationship between the ITZ strength and the dynamic increasing factor of the concrete tensile strength (TDIF) is proposed, and the rationality of the TDIF modified formula is verified.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine-Learning-Driven Discovery of Water Splitting BaFe 2 O 4 and Human-in-the-Loop Improvement via Al-Substitution for Increased Thermal Stability

Thermochemical hydrogen (TCH) production offers a promising method for converting thermal energy into hydrogen fuel through heat-driven redox cycles of metal oxides. Here, in this work a defect graph neural network (dGNN) was used to predict oxygen vacancy formation energies ΔH V O combined with Materials Project predictions of oxygen chemical potential stability to screen candidate oxides via high-throughput database analysis. BaFe 2 O 4 was identified as a promising material for experimental validation based on its predicted ΔH V O , oxygen chemical potential stability range, and potential for tunable substitutions to improve thermal properties. Experimental validation using thermogravimetric analysis (TGA), stagnation flow reactor (SFR), X-ray diffraction (XRD), and electron microscopy confirmed positive water-splitting behavior but also revealed limitations in thermal stability under aggressive reduction conditions. To address this, a human-in-the-loop modification strategy was employed introducing Al substitution in BaFe 2–x Al x O 4 ; this modification improves thermal stability, alters the crystal structure and enhances overall performance. These results demonstrate a combined computational and experimental workflow in which machine learning accelerates identification of promising candidates, while targeted experimental design enables optimization of functional performance. This approach advances the development of robust, cost-effective TCH materials and highlights the importance of integrating data-driven discovery with human-guided materials design in paving the way for scalable hydrogen production technologies.

organic↗
NASA NTRS2002

Multigrid Solution of the Navier-Stokes Equations at Low Speeds with Large Temperature Variations

Multigrid methods for the Navier-Stokes equations at low speeds and large temperature variations are investigated. The compressible equations with time-derivative preconditioning and preconditioned flux-difference splitting of the inviscid terms are used. Three implicit smoothers have been incorporated into a common multigrid procedure. Both full coarsening and semi-coarsening with directional fine-grid defect correction have been studied. The resulting methods have been tested on four 2D laminar problems over a range of Reynolds numbers on both uniform and highly stretched grids. Two of the three methods show efficient and robust performance over the entire range of conditions. In addition none of the methods have any difficulty with the large temperature variations.

Sockol, Peter M.↗
DOE OSTI2019

Development of the Local (Site) Symmetry Method in the Supercell Model for a Crystal with an Impurity

Symmetry aspects of the periodic model of a crystal with a point defect (a supercell model or a extended unit cell (EUC) model) are considered; splitting of the Wyckoff positions in a primitive crystal cell with the introduction of supercell and the change of the one-electron states classification over k-vector (Brillouin zone folding) are discussed. When considering a point defect in a crystal in the supercell model, it is necessary to take into account the symmetry of the one-electron states of the original crystal at the top of the valence band and at the bottom of the conduction band. The selected supercell should reproduce these states. Each specific selection of a supercell corresponds to a specific splitting of the Wyckoff positions of the original crystal and, as a result, the possibility to place a defect in positions with different point symmetry and to perform a calculation without taking into account the point symmetry of the crystal with a defect at all (site symmetry method). By the results of the calculation of a crystal with a defect in the supercell model, taking into account site symmetry, it is possible to determine the real symmetry of the crystal with a defect, which is essential for the interpretation of experimental data. A copper impurity in a lithium node in a LiCl crystal retains the cubic symmetry of the atom being replaced; the impurity of an iron atom in a titanium site with cubic local symmetry in a SrTiO{sub 3} crystal lowers the symmetry to tetragonal; and the H site in the CsPbI{sub 3} crystal is characterized by the complete removal of point symmetry with the formation of an I{sub 2}{sup −} ion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗