Search NASA⌕ Search

SEARCH · Search NASA

Results for “consistent discretization”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2019

An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables

We derive a formulation of the nonhydrostatic equations in spherical geometry with a Lorenz staggered vertical discretization. The combination conserves a discrete energy in exact time integration when coupled with a mimetic horizontal discretization. The formulation is a version of Dubos and Tort (2014, https://doi.org/10.1175/MWR-D-14-00069.1) rewritten in terms of primitive variables. It is valid for terrain following mass or height coordinates and for both Eulerian or vertically Lagrangian discretizations. The discretization relies on an extension to Simmons and Burridge (1981, https://doi.org/10.1175/1520-0493(1981)109<0758:AEAAMC>2.0.CO;2) vertical differencing, which we show obeys a discrete derivative product rule. This product rule allows us to simplify the treatment of the vertical transport terms. Energy conservation is obtained via a term-by-term balance in the kinetic, internal, and potential energy budgets, ensuring an energy-consistent discretization up to time truncation error with no spurious sources of energy. We demonstrate convergence with respect to time truncation error in a spectral element code with a horizontal explicit vertically implicit implicit-explicit time stepping algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗
DOE OSTI2022

Numerical Study on the Effect of Methane Doping in Hydrogen-Air Rotating Detonation Engines for Various Temperatures and Pressures

Rotating detonation engines (RDEs) have gained attention as a promising technology for future aviation engines. However, the numerical studies of these systems pose severe challenges due to the broad range of spatial and temporal scales. In this study, we use an adaptive mesh refinement based compressible, reactive solver PeleC to resolve the broad range of scales and accurately capture shock and detonation waves using high-resolution numerical schemes. Multi- species transport along with compressible Navier-Stokes equations are solved in the model along with a finite-rate based chemistry model. Embedded boundary method is used to model the complex geometry consisting of discrete fuel nozzles and the combustion chamber. The fuel consisting predominantly of hydrogen is doped with varying levels of methane while air is used as the oxidizer. For a specified total pressure and temperature, the number of stable detonation waves is found to decrease with increasing methane concentration in the fuel mixture. Additionally, no stable detonation solutions are observed for methane composition higher than 20% by volume for the range of operating conditions studied (300-900K, 10 Atm). The increased presence of high temperature zones is also indicative of higher thermal NOx emissions at low methane concentrations. The effect of fuel-air mixture composition and temperature on the detonability, detonation wave structure, mode transitions and their stability are analyzed in this study in addition to their implication on NOx emission.

automatic mesh refinement↗
DOE OSTI2025

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

Second moment methods (SMMs) are developed that are consistent with the discontinuous Galerkin spatial discretization of the discrete ordinates (or S\(_N\)) transport equations. The low-order (LO) diffusion system of equations is discretized with fully consistent P\(_1\), local discontinuous Galerkin (LDG), and interior penalty (IP) methods. A discrete residual approach is used to derive SMM correction terms that make each of the LO systems consistent with the high-order discretization. We show that the consistent methods are more accurate and have better solution quality than independently discretized LO systems, that they preserve the diffusion limit, and that the LDG and IP consistent SMMs can be scalably solved in parallel on a challenging, multimaterial benchmark problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Gauge constrained algorithm of variational discrete action theory at N = 3 for the multiorbital Hubbard model

The recently developed variational discrete action theory (VDAT) provides a systematic variational approach to the ground state of the quantum many-body problem, where the quality of the solution is controlled by an integer N, and increasing N monotonically approaches the exact solution. VDAT can be exactly evaluated in the d = ∞ multiorbital Hubbard model using the self-consistent canonical discrete action theory (SCDA), which requires a self-consistency condition for the integer time Green's functions. Previous work demonstrates that N = 3 accurately captures multiorbital Mott/Hund physics at a cost similar to the Gutzwiller approximation. Here we employ a gauge constraint to automatically satisfy the self-consistency condition of the SCDA at N = 3, yielding an even more efficient algorithm with enhanced numerical stability. We derive closed form expressions of the gauge constrained algorithm for the multiorbital Hubbard model with general density-density interactions, allowing VDAT at N = 3 to be straightforwardly applied to the seven-orbital Hubbard model. We present results and a performance analysis using N = 2 and N = 3 for the SU⁡(2⁢N orb ) Hubbard model in d = ∞ with N orb = 2–8, and compare to numerically exact dynamical mean-field theory solutions where available. Finally, the developments in this work will greatly facilitate the application of VDAT at N = 3 to strongly correlated electron materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Characterizing the microstructures of mammalian enamel by synchrotron phase contrast microCT

The enamel of mammalian teeth is a highly mineralized tissue that must endure a lifetime of cyclic contact and is inspiring the development of next-generation engineering materials. Attempts to implement enamel-inspired structures in synthetic materials have had limited success, largely due to the absence of a detailed understanding of its microstructure. Here, the present work used synchrotron phase-contrast microCT imaging to evaluate the three-dimensional microstructure of enamel from four mammals including Lion, Gray Wolf, Snow Leopard, and Black Bear. Quantitative results of image analysis revealed that the decussation pattern of enamel consists of discrete diazone (D) and parazone (P) bands of rods organized with stacking arrangement of D+/P/D-/P in all mammals evaluated; the D+ and D- refer to distinct diazone bands with juxtaposed rod orientations from the reference plane. Furthermore, the rod orientations in the bands can be described in terms of two principal angles, defined here as the pitch and yaw. While the pitch angle increases from the outer enamel to a maximum (up to ≈ 40°) near the dentin enamel junction, minimal spatial variations are observed in yaw across the enamel thickness. There are clear differences in the decussation parameters of enamel across species that are interpreted here with respect to the structural demands placed on their teeth. The rod pitch and band width of enamel are identified as important design parameters and appear to be correlated with the bite force quotient of the four mammals evaluated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydrogen bonding to the formate ligand: Structural characterization of 2,2′:6′,2″-terpyridine zinc formate·H 2 O

Here, two polymorphs of the 2,2’:6’,2”-terpyridine zinc formate compound (terpy)Zn(O 2 CH) 2 •H 2 O have been structurally characterized by X-ray diffraction. In each case, the formate ligands coordinate via a K 1 -monodentate coordination mode and participate in hydrogen bonding interactions with water molecules. Despite this similarity, however, the polymorphs exhibit different hydrogen bonding motifs. Specifically, while one polymorph forms chains in which water molecules link adjacent (terpy)Zn(O 2 CH) 2 molecules, the other polymorph consists of discrete tetranuclear [(terpy)Zn(O 2 CH) 2 •H 2 O] 4 moieties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Thermodynamic Analysis of Nano-Electrospray Induced Gas Jets

Nano-electrospray (nES) produces a plume of charged liquid droplets which have drag interactions with the surrounding gas after they are emitted from a source capillary. The resulting induced gas flow has often been considered unimportant, but we have demonstrated that a gas microjet of significant velocity up to 10s of m/s can be produced. In this work we present a thermodynamic framework that enables analysis of gas jet generation from electrosprays and introduce the important metrics for such analysis, effectiveness (a measure of momentum transfer from the electrosprayed aerosol to gas) and efficiency (a measure of energy conversion from electrical energy generating the electrospray to gas kinetic energy). This analytical framework is applicable to any two-phase flows consisting of discrete conservative-force-driven particles which exchange momentum with an inert, otherwise quiescent fluid medium to yield a co-flowing two-phase jet. Here, we apply this framework to sprays of water from nano-electrospray emitters to demonstrate that increasing the applied electrical potential difference, increasing liquid mass flowrate, and decreasing droplet size all can increase electrospray induced gas jet strength, but only the latter two do so while also increasing the momentum transfer effectiveness and energy conversion efficiency.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The role of cellular structure, non-equilibrium eutectic phases and precipitates on quasi-static strengthening mechanisms of as-built AlSi10Mg parts 3D printed via laser powder bed fusion

The quasi-static loading strengthening mechanism in the as-built state of an AlSi10Mg alloy 3D-printed via Laser Powder Bed Fusion (LPBF) was thoroughly identified and quantified using state-of-the-art electron microscopy and synchrotron X-ray diffraction techniques. The yield strength was comprehensively modelled through an in-depth characterization and quantification of the microstructural features as well as their effective volumes for strengthening. In particular, the non-equilibrium eutectic network was characterized with a two-fold structure: cell boundary particles as well as intracellular lamellar/fibrous networks, each consisting of discrete phases exhibiting a thru-thickness compositional gradient, a semi-coherent interface with the matrix and an abundance of crystal defects such as nano-sized sub-grains, microstrains and stacking faults. Further, these altogether made the eutectic phase the most potent contributor to the yield strength accounting for ~30–40% of the estimated value (i.e., cell boundary and eutectic network strengthening combined). Precipitates strengthening was identified as the second most potent mechanism via a shearing process only. Moreover, the presented methodology was able to capture the effect of LPBF processing variables on the individual strengthening contributions, e.g., the effect of a lower laser scanning speed on increasing the cell size as well as the mean precipitate size, which are shown to exhibit opposing impacts on strengthening.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Crystal structure of haloxon, C 14 H 14 C l3 O 6 P

The crystal structure of haloxon has been solved and refined using synchrotron X-ray powder diffraction data, and optimized using density functional theory techniques. Haloxon crystallizes in space group P2 1 /n (#14) with a = 19.60382(6), b = 10.05473(3), c = 8.73591(2) Å, β = 92.6617(2)°, V = 1720.088(11) Å 3 , and Z = 4. The structure consists of discrete molecules. The mean planes of the fused ring systems are approximately 0–11 and 011. The rings form staggered stacks perpendicular to these planes. There are no traditional hydrogen bonds in the structure, but several C–H∙∙∙O and C–H∙∙∙Cl hydrogen bonds contribute to the crystal energy. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®).

36 MATERIALS SCIENCE↗
DOE OSTI2024

Crystal structure of valbenazine, C 24 H 38 N 2 O 4

The crystal structure of valbenazine has been solved and refined using synchrotron X-ray powder diffraction data and optimized using density functional theory techniques. Valbenazine crystallizes in space groupP2 1 2 1 2 1 (#19) witha= 5.260267(17),b= 17.77028(7),c= 26.16427(9) Å,V= 2445.742(11) Å 3 , andZ= 4 at 295 K. The crystal structure consists of discrete molecules and the mean plane of the molecules is approximately (8,−2,15). There are no obvious strong intermolecular interactions. There is only one weak classical hydrogen bond in the structure, from the amino group to the ether oxygen atom. Two intramolecular and one intermolecular C–H⋯O hydrogen bonds also contribute to the lattice energy. The powder pattern has been submitted to ICDD for inclusion in the Powder Diffraction File™ (PDF®)

Materials Science↗
DOE OSTI2020

Cas9 interrogates DNA in discrete steps modulated by mismatches and supercoiling

The CRISPR-Cas9 nuclease has been widely repurposed as a molecular and cell biology tool for its ability to programmably target and cleave DNA. Cas9 recognizes its target site by unwinding the DNA double helix and hybridizing a 20-nucleotide section of its associated guide RNA to one DNA strand, forming an R-loop structure. A dynamic and mechanical description of R-loop formation is needed to understand the biophysics of target searching and develop rational approaches for mitigating off-target activity while accounting for the influence of torsional strain in the genome. Here we investigate the dynamics of Cas9 R-loop formation and collapse using rotor bead tracking (RBT), a single-molecule technique that can simultaneously monitor DNA unwinding with base-pair resolution and binding of fluorescently labeled macromolecules in real time. By measuring changes in torque upon unwinding of the double helix, we find that R-loop formation and collapse proceed via a transient discrete intermediate, consistent with DNA:RNA hybridization within an initial seed region. Using systematic measurements of target and off-target sequences under controlled mechanical perturbations, we characterize position-dependent effects of sequence mismatches and show how DNA supercoiling modulates the energy landscape of R-loop formation and dictates access to states competent for stable binding and cleavage. Consistent with this energy landscape model, in bulk experiments we observe promiscuous cleavage under physiological negative supercoiling. The detailed description of DNA interrogation presented here suggests strategies for improving the specificity and kinetics of Cas9 as a genome engineering tool and may inspire expanded applications that exploit sensitivity to DNA supercoiling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Low power implementation of high frequency SiPM readout for Cherenkov and scintillation detectors in TOF-PET

State-of-the-art (SoA) electronic readout for silicon photomultiplier (SiPM)-based scintillation detectors that demonstrate experimental limits in achievable coincidence time resolution (CTR) leverage low noise, high frequency signal processing to facilitate a single photon time response that is near the limit of the SiPMs architecture. This readout strategy can optimally exploit fast luminescence and prompt photon populations, and promising measurements show detector concepts employing this readout can greatly advance PET detector CTR, relative to SoA in clinical systems. However, the technique employs power hungry components which make the electronics chain impractical for channel-dense time-of-flight (TOF)-PET detectors. We have developed and tested a low noise and high frequency readout circuit which is performant at low power and consists of discrete elements with small footprints, making it feasible for integration into TOF-PET detector prototypes. A 3 × 3 mm 2 Broadcom SiPM with this readout chain exhibited sub-100 ps single photon time resolution at 10 mW of power consumption, with a relatively minor performance degradation to 120 ± 2 ps FWHM at 5 mW. CTR measurements with 3 × 3 × 20 mm 3 LYSO and fast LGSO scintillators demonstrated 127 ± 3 ps and 113 ± 2 ps FWHM at optimal power operation and 133 ± 2 ps and 121 ± 3 ps CTR at 5 mW. BGO crystals 3 × 3 × 20 mm 3 in size show 271 ± 5 ps FWHM CTR (1174 ± 14 ps full-width-at-tenth-maximum (FWTM)) at optimal power dissipation and 289 ± 8 ps (1296 ± 33 ps FWTM) at 5 mW. The compact and low power readout topology that achieves this performance thereby offers a platform to greatly advance PET system CTR and also opportunities to provide high performance TOF-PET at reduced material cost.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

3D polycatenated architected materials

Architected materials derive their properties from the geometric arrangement of their internal structural elements. Their designs rely on continuous networks of members to control the global mechanical behavior of the bulk. Here, in this study, we introduce a class of materials that consist of discrete concatenated rings or cage particles interlocked in three-dimensional networks, forming polycatenated architected materials (PAMs). We propose a general design framework that translates arbitrary crystalline networks into particle concatenations and geometries. In response to small external loads, PAMs behave like non-Newtonian fluids, showing both shear-thinning and shear-thickening responses, which can be controlled by their catenation topologies. At larger strains, PAMs behave like lattices and foams, with a nonlinear stress-strain relation. At microscale, we demonstrate that PAMs can change their shapes in response to applied electrostatic charges. The distinctive properties of PAMs pave the path for developing stimuli-responsive materials, energy-absorbing systems, and morphing architectures.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Experimental and Computational Analysis of NEAMS Pebble Bed Reactors

The Advanced High-Temperature Pebble Bed Reactor concept leverages a particle-based fuel format consisting of discrete spherical graphite pebbles arrayed in a packed bed architecture. Thermal regulation achieved via flow of gas (e.g., helium) or liquid (e.g., molten salt) coolants through the void spaces between pebbles in the bed (characteristic pebble diameters: ~ 6 cm, gas cooled; ~ 3 cm, liquid cooled). This design is of interest owing to inherent advantages including passive safety features and highly efficient heat transfer characteristics that enable high power densities to be achieved. Fully leveraging these advantages requires a fundamental understanding of the complex coolant flow structure within the pebble-shaped fuel bed. Unfortunately, complexities associated with non-invasively probing the micro-scale flow phenomena within the complex network of void spaces between pebbles have hindered efforts to characterize the underlying transport phenomena. Without this fundamental foundation, it becomes challenging to rationally construct macroscopic transport models based upon volume averages of micro-scale parameters that accurately capture the observed flow behavior. Consequently, this lack of fundamental knowledge has limited previous work to coarse comparisons between theory and extremely limited experimental data. Moreover, these comparisons have been made only in terms of integral or macroscopic results. Therefore, multipoint measurements with a high level of spatial and temporal resolution are critically needed to fully map the complex flow patterns and to provide data at high spatial density to permit accurate volume averaging in the pebble bed. For this project, we performed a coordinated experimental and computational effort to quantitatively map the full-field velocity and temperature fields in the interstitial spaces within the pebble bed. The major accomplishments of this project included the important knowledge acquired for the flow mixing and heat transfer in the interstitial pore scales, the effects of flow mixing in the near-wall and far-wall regions, and the by-pass flows of pebble bed systems. The produced high-resolution experimental data (not broadly available in existing literature) were then used for the numerical studies enabling the refinement and validation of standard simulation tools (Nek5000).

42 ENGINEERING↗
DOE OSTI2022

Microscopic Visualization of Cell-Cell Adhesion Complexes at Micro and Nanoscale

Considerable progress has been made in our knowledge of the morphological and functional varieties of anchoring junctions. Cell-cell adhesion contacts consist of discrete junctional structures responsible for the mechanical coupling of cytoskeletons and allow the transmission of mechanical signals across the cell collective. The three main adhesion complexes are adherens junctions, tight junctions, and desmosomes. Microscopy has played a fundamental role in understanding these adhesion complexes on different levels in both physiological and pathological conditions. In this review, we discuss the main light and electron microscopy techniques used to unravel the structure and composition of the three cell-cell contacts in epithelial and endothelial cells. It functions as a guide to pick the appropriate imaging technique(s) for the adhesion complexes of interest. We also point out the latest techniques that have emerged. At the end, we discuss the problems investigators encounter during their cell-cell adhesion research using microscopic techniques.

59 BASIC BIOLOGICAL SCIENCES↗