Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak solutions”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical and boundary condition effects on the prediction of detonation engine behavior using detailed numerical simulations

High-fidelity numerical simulations of an experimental rotating detonation engine with discrete fuel/air injection were conducted. A series of configurations with different feed-plenum pressures but with constant equivalence ratio were studied. Detailed chemical kinetics for the hydrogen/air system is used. A resolution study for the full rotating detonation engine (RDE) system simulation is also conducted. Two kinds of boundary conditions, a total pressure boundary and a constant mass flow rate boundary, are used to assess the effects of the inlet boundary. As mass flow rate is increased, the total pressure boundary causes more error in the axial pressure distribution while the constant mass flow rate gives a better solution for all cases ran. The simulations confirm experimental findings, and reproduce qualitative as well as some of the quantitative trends. These results demonstrate that a) fuel-air mixing is highly non-uniform within the detonation chamber, leading to variations in local equivalence ratio, b) the fuel and oxidizer injectors experience significant backflow as the detonation wave passes over, but recover at different rates which further augments the inefficiencies in mixing, and c) parasitic combustion in the mixing region makes the detonation wave weak by extending the reaction zone across the wave.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Fluctuations analysis of spin susceptibility: Néel ordering revisited in dynamical mean field theory

Here in this paper, we revisit the antiferromagnetic (AF) phase diagram of the single-band three-dimensional half-filled Hubbard model on a simple cubic lattice studied within the dynamical mean field theory (DMFT). Although this problem has been investigated extensively in the literature, a comprehensive understanding of the impact of the different one- and, in particular, two-particle local correlation functions of DMFT on the AF transition temperature is still missing. We have, hence, performed a fluctuation analysis of T N with respect to different local bosonic fluctuations (charge, spin, particle-particle) contained in the two-particle vertex of DMFT. Our results indicate that, beyond weak coupling, the screening of the DMFT vertex by local fluctuations leads to an enhancement of T N with respect to a random phase approximation (RPA) like calculation where this vertex is replaced by the bare interaction. The overall suppression of T N in DMFT with respect to RPA is then solely due to the incoherence introduced by the DMFT self-energy in the one-particle Green's functions. This illustrates the Janus-faced role of the local moment formation in the DMFT solution of the Hubbard model, which leads to completely opposite effects in the one- and two-particle correlation functions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The Monomeric Alanediyl : AlAr i Pr8 (Ar i Pr8 = C 6 H-2,6-(C 6 H 2 -2,4,6-Pr i 3 ) 2 -3,5-Pr i 2 ): An Organoaluminum(I) Compound with a One-Coordinate Aluminum Atom

Reduction of the aluminum iodide AlI 2 Ar i Pr8 ( 1 ; Ar i Pr8 = C 6 H-2,6-(C 6 H 2 -2,4,6-Pr i 3 ) 2 -3,5-Pr i 2 ) with 5% w/w Na/NaCl in hexanes gave a dark red solution from which the monomeric alanediyl : AlAr i Pr8 ( 2 ) was isolated in ca. 28% yield as yellow-orange crystals. Compounds 1 and 2 were characterized by X-ray crystallography, electronic and NMR spectroscopy, and theoretical calculations. The Al atom in 2 is one-coordinate, and the compound displays two absorptions in its electronic spectrum at 354 and 455 nm. It reacts with H 2 under ambient conditions to give the aluminum hydride {AlH(μ-H)Ar i Pr8 } 2 , probably via a weakly bound dimer of 2 as an intermediate.

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

Compare linear-system solver and preconditioner stacks with emphasis on GPU performance and propose phase-2 NGP solver development pathway

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind project are Nalu-Wind, which is an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, and OpenFAST, which is a whole-turbine simulation code. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and a momentum equation for the velocity. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. In this report we evaluated GPU-performance baselines for the linear solvers in the Trilinos and hypre solver stacks using two representative Nalu-Wind simulations: an atmospheric boundary layer precursor simulation on a structured mesh, and a fixed-wing simulation using unstructured overset meshes. Both strong-scaling and weak-scaling experiments were conducted on the OLCF supercomputer Summit and similar proxy clusters. We focused on the performance of multi-threaded Gauss-Seidel and two-stage Gauss-Seidel that are extensions of classical Gauss-Seidel; of one-reduce GMRES, a communication-reducing variant of the Krylov GMRES; and algebraic multigrid methods that incorporate the afore-mentioned methods. The team has established that AMG methods are capable of solving linear systems arising from the fixed-wing overset meshes on CPU, a critical intermediate result for ExaWind FY20 Q3 and Q4 milestones. For the fixed-wing strong-scaling study (model with 3M grid-points), the team identified that Nalu-Wind simulations with the new Trilinos and hypre solvers scale to modest GPU counts, maintaining above 70% efficiency up to 6 GPUs. However, there still remain significant bottlenecks to performance: matrix assembly (hypre), AMG setup (hypre and Trilinos) In the weak-scaling experiments (going from 0.4M to 211M gridpoints), it's shown that the solver apply phases are faster on GPUs, but that Nalu-Wind simulation times grow, primarily due to the multigrid-setup process. Finally, based on the report outcomes, we propose a linear solver path-forward for the remainder of the ExaWind project. Near term, the NREL team will continue their work on GPU-based linear-system assembly. They will also investigate how the use of alternatives to the NVIDIA UVM (unified virtual memory) paradigm affects performance. Longer term, the NREL team will evaluate algorithmic performance on other types of accelerators and merge their improvements back to the main hypre repository branch. Near term, the Trilinos team will address performance bottlenecks identified in this milestone, such as implementing a GPU-based segregated momentum solve and reusing matrix graphs across linear-system assembly phases. Longer term, the Trilinos team will do detailed analysis and optimization of multigrid setup.

17 WIND ENERGY↗
DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗
DOE OSTI2024

Lyapunov exponent as a signature of dissipative many-body quantum chaos

A distinct feature of Hermitian quantum chaotic dynamics is the exponential increase of certain out-of-time-order correlation (OTOC) functions around the Ehrenfest time with a rate given by a Lyapunov exponent. Physically, the OTOCs describe the growth of quantum uncertainty that crucially depends on the nature of the quantum motion. Here, we employ the OTOC in order to provide a precise definition of dissipative quantum chaos. For this purpose, we compute analytically the Lyapunov exponent for the vectorized formulation of the large- q limit of a q -body Sachdev-Ye-Kitaev model coupled to a Markovian bath. These analytic results are confirmed by an explicit numerical calculation of the Lyapunov exponent for several values of q ≥ 4 based on the solutions of the Schwinger-Dyson and Bethe-Salpeter equations. We show that the Lyapunov exponent decreases monotonically as the coupling to the bath increases and eventually becomes negative at a critical value of the coupling signaling a transition to a dynamics which is no longer quantum chaotic. Therefore, a positive Lyapunov exponent is a defining feature of dissipative many-body quantum chaos. The observation of the breaking of the exponential growth for sufficiently strong coupling suggests that dissipative quantum chaos may require in certain cases a sufficiently weak coupling to the environment. Published by the American Physical Society 2024

García-García, Antonio M.↗
DOE OSTI2021

On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compare these results to Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

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

Reassessment of Dose Contributions from the M-Area Glass Special Waste Form Buried in the E-Area Low-Level Waste Facility

During a recent revision to the groundwater (GW) and inadvertent human intruder (IHI) screening analysis report (Aleman and Hamm, 2023), weaknesses within the screening logic were observed resulting in U-235 being brought back into the list of radionuclides requiring inventory limits within the M-Area Glass SWF. In addition, while reviewing the earlier vadose zone (VZ) transport analysis efforts, an input deck error was observed in the U-234G VZ transport simulation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

The Role of Dislocations in the Anelasticity of the Upper Mantle

Dislocation‐based dissipation mechanisms potentially control the viscoelastic response of Earth's upper mantle across a variety of geodynamic contexts, including glacial isostatic adjustment, postseismic creep, and seismic‐wave attenuation. However, there is no consensus on which dislocation‐based, microphysical process controls the viscoelastic behavior of the upper mantle. Although both intergranular (plastic anisotropy) and intragranular (backstress) mechanisms have been proposed, there is currently insufficient laboratory data to discriminate between those mechanisms. Here, we present the results of forced‐oscillation experiments in a deformation‐DIA apparatus at confining pressures of 3–7 GPa and temperatures of 298–1370 K. Our experiments tested the viscoelastic response of polycrystalline olivine—the main constituent of the upper mantle—at stress amplitudes from 70 to 2,800 MPa. Mechanical data are complemented by microstructural analyses of grain size, crystallographic preferred orientation, and dislocation density. We observe amplitude‐ and frequency‐dependent attenuation and modulus relaxation and find that numerical solutions of the backstress model match our results well. Therefore, we argue that interactions among dislocations, rather than intergranular processes (e.g., plastic anisotropy or grain boundary sliding), control the viscoelastic behavior of polycrystalline olivine in our experiments. In addition, we present a linearized version of the constitutive equations of the backstress model and extrapolate it to conditions typical of seismic‐wave propagation in the upper mantle. Our extrapolation demonstrates that the backstress model can explain the magnitude of seismic‐wave attenuation in the upper mantle, although some modification is required to explain the weak frequency dependence of attenuation observed in nature and in previous experimental work.

Hein, Diede [Univ. of Minnesota, Minneapolis, MN (↗
DOE OSTI2023

Facile and dynamic cleavage of every iron–sulfide bond in cuboidal iron–sulfur clusters

Nature employs weak-field metalloclusters to support a wide range of biological processes. The most ubiquitous metalloclusters are the cuboidal Fe–S clusters, which are comprised of Fe sites with locally high-spin electronic configurations. Such configurations enhance rates of ligand exchange and imbue the clusters with a degree of structural plasticity that is increasingly thought to be functionally relevant. Here, we examine this phenomenon using isotope tracing experiments. Specifically, we demonstrate that synthetic [Fe 4 S 4 ] and [MoFe 3 S 4 ] clusters exchange their Fe atoms with Fe 2+ ions dissolved in solution, a process that involves the reversible cleavage and reformation of every Fe–S bond in the cluster core. This exchange is facile—in most cases occurring at room temperature on the timescale of minutes—and documented over a range of cluster core oxidation states and terminal ligation patterns. In addition to suggesting a highly dynamic picture of cluster structure, these results provide a method for isotopically labeling pre-formed clusters with spin-active nuclei, such as 57 Fe. Such a protocol is demonstrated for the radical S -adenosyl- l -methionine enzyme, RlmN.

Science & Technology - Other Topics↗
DOE OSTI2025

Investigating intermolecular interactions among CO 2 , water and PEEK-ionene membrane using cryo ToF-SIMS and isotopic labeling

Cryogenic time-of-flight secondary ion mass spectrometry (cryo ToF-SIMS) has emerged as a powerful tool for investigating molecular interactions, speciation, and dynamics in materials for CO 2 capture. In this study, we apply cryo ToF-SIMS to probe interactions between CO 2 , water, and PEEK-ionene membranes—a promising material for direct CO 2 capture due to its selectivity, durability, and efficiency. Despite this potential, the mechanisms governing CO 2 diffusion and the influence of water vapor on CO 2 behavior remain unclear. To address this, we loaded PEEK-ionene membranes with 13 CO 2 and D 2 O and employed cryo ToF-SIMS to visualize the 3D distribution of CO 2 and water within the membrane. While prior studies suggest that 13 CO 2 is absorbed under ambient conditions, our cryo ToF-SIMS analysis revealed no enhancement of the 13 C/ 12 C ratio, suggesting weak CO 2 -membrane interactions. As a result, CO 2 vaporizes even at low temperatures (−140°C) under vacuum conditions. In contrast, D 2 O displayed a relatively homogeneous distribution in the membrane, suggesting stronger water-membrane interactions via hydrogen bonding (18–20 kJ/mol). Interestingly, CO 2 was not detected in D 2 O-loaded membranes, indicating minimal interference from water vapor on CO 2 diffusion. As a comparison, the cryo ToF-SIMS data show that 13 CO 2 can readily react with a basic Na 2 CO 3 aqueous solution to form NaH 13 CO 3 . These findings demonstrate cryo ToF-SIMS as a critical technique for understanding gas-water-membrane interactions, offering insights for membrane functionalization to improve CO 2 capture efficiency.

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

High‐order symmetric cubature rules for tetrahedra and pyramids

Abstract In this article, we present an algorithm to construct high‐order fully symmetric cubature rules for tetrahedral and pyramidal elements, with positive weights and integration points that are in the interior of the domain. Cubature rules are fully symmetric if they are invariant to affine transformations of the domain. We divide the integration points into symmetry orbits where each orbit contains all the points generated by the permutation stars. These relations are represented by equality constraints. The construction of symmetric cubature rules require the solution of nonlinear polynomial equations with both inequality and equality constraints. For higher orders, we use an algorithm that consists of five sequential phases to produce the cubature rules. In the literature, symmetric numerical integration rules are available for the tetrahedron for orders p = 1 – 10, 14 , and for the pyramid up to p = 10 . We have obtained fully symmetric cubature rules for both of these elements up to order p = 20 . Numerical tests are presented that verify the polynomial‐precision of the cubature rules. Convergence studies are performed for the integration of exponential, weakly singular, and trigonometric test functions over both elements with flat and curved faces. With increase in p , improvements in accuracy is realized, though nonmonotonic convergence is observed.

Jaśkowiec, Jan↗
DOE OSTI2022

Using electron channeling contrast imaging to inform and improve the growth of high-efficiency GaAs solar cells on nanopatterned GaAs substrates

Patterned substrates provide opportunities for reducing the cost of high-efficiency III-V devices by incorporating mechanically weak layers beneficial for substrate reuse (e.g. by spalling). In this work, the functionality of electron channeling contrast imaging (ECCI) as a tool to efficiently understand and mitigate defect formation is exemplified by developing a process in which high-quality III-V material can be grown on nanopatterned GaAs substrates. Reactive ion etching used in the patterning process was found to damage the GaAs substrate surface, leading to the formation of stacking faults in the epitaxial material as observed by ECCI. Etching the patterned substrates in a 1 NH4OH: 1 H2O2: 50 DI H2O solution for 10 s prior to growth removed the substrate surface damage and stacking faults were no longer present. Growth of solar cell device structures initially produced samples with many macroscale flaws creating shunts in the devices, which complicated the assessment of material quality by device measurements. However, ECCI revealed that the epitaxial material surrounding macroscale flaws was free from any crystallographic defects such as stacking faults and threading dislocations. With this knowledge, we focused on refining the patterning process to eliminate the macroscale flaws. Solar cells were then grown on the improved nanopatterned substrates and exhibited device structures with defect densities less than 5 x 105 cm -2 and average conversion efficiency of 24.8%, nearly identical to devices grown on unpatterned epi-ready substrates (25.0%).

14 SOLAR ENERGY↗
DOE OSTI2023

Investigating microstructural evolution in SolidStir™ extruded oxide dispersion strengthened 14YWT alloy fuel cladding tube

Structural components for extreme environments require advanced materials and manufacturing processes. One such advanced material is oxide dispersion strengthened (ODS) Fe–14 wt.%Cr–3W–0.4Ti–0.3Y 2 O 3 (14YWT), which is developed for components to be used in Gen IV and fusion nuclear reactors. However, the conventional manufacturing processes to produce components, such as fuel cladding tubes, either do not retain the desired microstructural attributes or are costly and less efficient. As a result, a more attractive option is a novel manufacturing process developed on the principle of friction stir welding or processing (FSW/P) and commercially referred to as SolidStir™ Extrusion (SSE). Here this study used the SSE technique to study the manufacturing and microstructural evolution of a fuel cladding tube made of ODS 14YWT alloy. The cladding extrusion by the SSE technique involved the use of ball-milled 14YWT powders with Y, Ti, and O in the solid solution and the use of a specially designed W-25Re-Hf tool for consolidation and extrusion of the powders. A scanning electron microscope (SEM), electron backscattered diffraction (EBSD), and transmission electron microscope (TEM) were utilized for microstructural characterization. A Keyence microscope captured macrostructural photographs of the extruded tube. EBSD examination of the extruded tube on both transverse and longitudinal cross-sections showed the presence of dynamically recrystallized grains and revealed that the average grain size of the transverse cross-section was smaller than that of the longitudinal cross-section. The texture was weak and consisted of some amount of shear texture. The presence of nano-oxide clusters or precipitates of Y, Ti, and O (pyrochlores) in the extruded tube was determined by TEM. Aging heat treatment caused the average precipitate size to decrease and the density to increase compared to the as-processed condition. The results indicated that SSE is a viable tool for manufacturing fuel cladding tubes with the desired microstructural attributes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Crystal Growth, Single Crystal Structure, and Biological Activity of Thiazolo-Pyridine Dicarboxylic Acid Derivatives

Four novel TPDCA derivatives were prepared via a supersaturation method combining TPDCA with water, N-methyl-2-pyrrolidone (NMP), Na(PO 2 H 2 ), and ammonia solution: 2(C 9 H 7 NO 5 S)H 2 O (1), (C 9 H 7 NO 5 S)C 5 H 9 NO (2), (C 9 H 7 NO 5 S)Na(PO 2 H 2 ) (3), and (C 9 H 5 NO 5 S)(NH 4 )2(H 2 O) (4). Their crystal structures were determined by single-crystal X-ray diffraction. Compounds (1) and (2) crystallize in the monoclinic space groups P2 1 and P2 1 /c, respectively, whereas compounds (3) and (4) crystallize in the triclinic space group P1 - . Weak and moderate hydrogen bonds were detected in the four compounds. In the biological tests, (1) and (3) exhibited significant antibacterial activity against Escherichia coli and Staphylococcus aureus; in addition, (1) was cytotoxic against leukemia HL-60 cells with the IC 50 value of 158.5 ± 12.5 μM.

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

Relaxor ferroelectric polymer exhibits ultrahigh electromechanical coupling at low electric field

Electromechanical (EM) coupling—the conversion of energy between electric and mechanical forms—in ferroelectrics has been used for a broad range of applications. Ferroelectric polymers have weak EM coupling that severely limits their usefulness for applications. Here, we introduced a small amount of fluorinated alkyne (FA) monomers (<2 mol %) in relaxor ferroelectric poly(vinylidene fluoride-trifluoroethylene-chlorofluoroethylene) (PVDF-TrFE-CFE) terpolymer that markedly enhances the polarization change with strong EM coupling while suppressing other polarization changes that do not contribute to it. Under a low–dc bias field of 40 megavolts per meter, the relaxor tetrapolymer has an EM coupling factor (k 33 ) of 88% and a piezoelectric coefficient (d 33 ) >1000 picometers per volt. These values make this solution-processed polymer competitive with ceramic oxide piezoelectrics, with the potential for use in distinct applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Hydraulic performance and wave transmission through nature-inspired perforated hollow-base piles breakwater

Nature-based Solutions (NbS) are critical to achieving long-term sustainable coastlines in deltaic coastlines, particularly susceptible to severe impacts under climate change. In recent years, coastal erosion has accelerated in the Vietnam Mekong Delta (VMD), particularly in the West Sea (from Ca Mau to Kien Giang). Due to the characteristics of muddy coastlines and eroded mangrove forests in this region, only a few engineering solutions have effectively reduced waves and promoted sediment accumulation to facilitate beach accretion. However, these solutions require significant investment and cannot be relocated or reused. Therefore, this study introduces a coastal mangrove-inspired perforated hollow-base piles breakwater (PHBPB), which uses the principle of biomimicry by studying the ability of mangrove groups to protect the coastline, stimulate sediment deposition, and over time, restore mangrove forests in the VMD. The hydraulic parameters of PHBPB were tested in 260 scenarios using a physical model, which included varying parameters such as crest width (B), pile row numbers (n c ), crest freeboard (R c ), and wave parameters (H s , T p ) in determining the capabilities of wave transmission, reflection, and wave dissipation. The results showed that the process of wave energy dissipation through perforated hollow-base breakwaters without piles (PHBBWP) is governed by three main parameters, including the relative water depth of the crest (R c /H m0,i ), the relative width of the crest (B/H m0,i ), and the wave slope at the structure location (s m = H m0,i /L m ). In contrast, for PHBPB, the primary influencing factors are the relative submergence depth or the length of the submerged pile portion (R c /H m0,i ) and the relative width of the pile system (X b /L m ). Finally, the hollow structure is suitable for the soft, weak mud geology in this region with poor load-bearing capacity, allows for rapid installation, is reasonably cost-effective, and can be relocated and reused once regeneration has been completed in each stretch of coastline.

54 ENVIRONMENTAL SCIENCES↗