Numerical simulation of a methane-oxygen rotating detonation rocket engine
Not provided.
Energy & Fuels↗
SEARCH · Engineering Papers
Results for “Field stability”
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.
Exploring the Structural, Dynamic, and Functional Properties of Metal–Organic Frameworks through Molecular Modeling
This review spotlights the role of atomic-level modeling in research on metal-organic frameworks (MOFs), especially the key methodologies of density functional theory (DFT), Monte Carlo (MC) simulations, and molecular dynamics (MD) simulations. The discussion focuses on how periodic and cluster-based DFT calculations can provide novel insights into MOF properties, with a focus on predicting structural transformations, understanding thermodynamic properties and catalysis, and providing information or properties that are fed into classical simulations such as force field parameters or partial charges. Classical simulation methods, highlighting force field selection, databases of MOFs for high-throughput screening, and the synergistic nature of MC and MD simulations, are described. By predicting equilibrium thermodynamic and dynamic properties, these methods offer a wide perspective on MOF behavior and mechanisms. Additionally, the incorporation of machine learning (ML) techniques into quantum and classical simulations is discussed. These methods can enhance accuracy, expedite simulation setup, reduce computational costs, as well as predict key parameters, optimize geometries, and estimate MOF stability. By charting the growth and promise of computational research in the MOF field, the aim is to provide insights and recommendations to facilitate the incorporation of computational modeling more broadly into MOF research.
36 MATERIALS SCIENCE↗
Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria
The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton’s equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Anisotropic elastoplastic response of double-porosity media
Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.
42 ENGINEERING↗
Differentiation strategies of soil rare and abundant microbial taxa in response to changing climatic regimes
Despite the important roles of soil microbes, especially the most diverse rare taxa in maintaining community diversity and multifunctionality, how different climate regimes alter the stability and functions of the rare microbial biosphere remains unknown. We reciprocally transplanted field soils across a latitudinal gradient to simulate climate change and sampled the soils annually after harvesting the maize over the following 6 years (from 2005 to 2011). By sequencing microbial 16S ribosomal RNA gene amplicons, we found that changing climate regimes significantly altered the composition and dynamics of soil microbial communities. A continuous succession of the rare and abundant communities was observed. Rare microbial communities were more stable under changing climatic regimes, with lower variations in temporal dynamics, and higher stability and constancy of diversity. More nitrogen cycling genes were detected in the rare members than in the abundant members, including amoA, napA, nifH, nirK, nirS, norB and nrfA. Random forest analysis and receiver operating characteristics analysis showed that rare taxa may act as potential contributors to maize yield under changing climatics. The study indicates that the taxonomically and functionally diverse rare biosphere has the potential to increase functional redundancy and enhance the ability of soil communities to counteract environmental disturbances. With ongoing global climate change, exploring the succession process and functional changes of rare taxa may be important in elucidating the ecosystem stability and multifunctionality that are mediated by microbial communities.
59 BASIC BIOLOGICAL SCIENCES↗
SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies
Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
An Accurate SUPG-stabilized Continuous Galerkin Discretization for Anisotropic Heat Flux in Magnetic Confinement Fusion
We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov–Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Geometries with mismatched branes
We study Randall-Sundrum two brane setups with mismatched brane tensions. For the vacuum solutions, boundary conditions demand that the induced metric on each of the branes is either de Sitter, Anti-de Sitter, or Minkowski. For incompatible boundary conditions, the bulk metric is necessarily time-dependent. This introduces a new class of time-dependent solutions with the potential to address cosmological issues and provide alternatives to conventional inflationary (or contracting) scenarios. We take a first step in this paper toward such solutions. One important finding is that the resulting solutions can be very succinctly described in terms of an effective action involving only the induced metric on either one of the branes and the radion field. But the full geometry cannot necessarily be simply described with a single coordinate patch. We concentrate here on the time- dependent solutions but argue that supplemented with a brane stabilization mechanism one can potentially construct interesting cosmological models this way. This is true both with and without a brane stabilization mechanism.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Rheology and stability of nanoparticle-stabilized CO 2 foam under reservoir conditions
This work investigates the rheology and stability of nanoparticle-stabilized CO 2 foams under reservoir conditions (high temperature and high pressure) for fracturing applications. The effects of different parameters on foam apparent viscosity and foam stability were experimentally investigated including the effects of nanoparticle concentration, salinity, foam quality (Γ), shear rate and temperature. The power law model was applied to calculate foam apparent viscosity owing to its pseudo-plastic behavior, and the changes of foam heights over time were used to evaluate foam stability. Results showed that the CO 2 foam apparent viscosity featured a mountain-shaped curve versus Γ, with the peak apparent viscosity obtained at 70% foam quality. The increase of salinity (up to 11%) in the continuous phase improved both foam stability and rheology. Higher nanoparticle concentration could contribute to better foam stability, but there was a threshold concentration, above which the foam apparent viscosity remained stabilized. The CO 2 foam stabilized by nanoparticles displayed a shear-thickening behavior as the experiment flow rate increased from 6 mL/min to 18 mL/min. Further studies showed that as the total flow rate increased, the CO 2 foam became finer-textured with better stability. Elevated temperatures could undermine foam apparent viscosity and long-term stability. The results of this study could provide guidelines as to the design of foam fracturing systems for field applications.
03 NATURAL GAS↗
Colloidal Synthesis of Palladium Nanocluster‐Decorated Cs 3 Sb 2 Cl 9 Perovskite Heterostructural Nanorods for Enhanced CO 2 Photoreduction
Developing efficient and sustainable photocatalysts for CO 2 reduction remains a significant challenge, particularly with environmentally benign materials. Here, in this study, we report the first one-step synthesis of metal–lead-free perovskite heterostructural nanocrystals by decorating Cs 3 Sb 2 Cl 9 perovskite nanorods with size-controlled Pd nanoclusters via a one-step hot-injection method. The resulting Pd-Cs 3 Sb 2 Cl 9 heteronanorods (HNRs) exhibit strong interfacial electronic coupling, enhanced charge separation, and excellent colloidal stability. Transient absorption spectroscopy and DFT calculations reveal a built-in electric field that drives directional electron transfer from the perovskite host to the Pd domains. Under UV irradiation, the Pd-Cs 3 Sb 2 Cl 9 HNRs demonstrate excellent CO 2 photoreduction activity with high CH 4 selectivity, achieving a record apparent quantum yield (AQY) of 2.62% among halide perovskite nanocrystal-based systems with a large electronic yield of 689.3 ± 12.2 µmol·g cat −1 . In situ spectroscopic monitoring and Gibbs free energy analysis further unveil a Pd-facilitated reaction pathway involving stabilization of key intermediates. This work introduces a new class of lead-free perovskite-based heterostructures through a facile one-step synthesis strategy and offers a new design principle for next-generation photocatalysts for solar fuel production.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Thermal properties of field-assisted-sintered SiCN–Y 2 O 3 composites
Polymer-derived amorphous SiCN has excellent high-temperature stability and properties. To reduce the shrinkage during pyrolysis and to improve the high-temperature oxidation resistance, Y 2 O 3 was added as a filler. In this study, polymer-derived SiCN–Y 2 O 3 composites were fabricated by mixing a polymeric precursor of SiCN with Y 2 O 3 submicron powders in different ratios. The mixtures were cross-linked and pyrolyzed in argon. SiCN–Y 2 O 3 composites were processed using field-assisted sintering technology at 1350°C for 5 min under vacuum. Dense SiCN–Y 2 O 3 composite pellets were successfully made with relative density higher than 98% and homogeneous microstructure. Due to low temperature and short time of the heat-treatment, the grain growth of Y 2 O 3 was substantially inhibited. The Y 2 O 3 grain size was ~1 μm after sintering. The composites’ heat capacity, thermal diffusivity, and thermal expansion coefficients were characterized as a function of temperature. The thermal conductivity of the composites ceramics decreased as the amount of amorphous SiCN increased and the coefficient of thermal expansion (CTE) of the composites increased with Y 2 O 3 content. However, the thermal conductivity and CTE did not follow the rule of mixture. This is likely due to the partial oxidation of SiCN and the resultant impurity phases such as Y 2 SiO 5 , Y 2 Si 2 O 7 , and Y 4.67 (SiO 4 ) 3 O.
36 MATERIALS SCIENCE↗
Controlled Nucleation and Stabilization of Ferroelectric Domain Wall Patterns in Epitaxial (110) Bismuth Ferrite Heterostructures
Ferroelectric domain walls, topological entities separating domains of uniform polarization, are promising candidates as active elements for nanoscale memories. In such applications, controlled nucleation and stabilization of domain walls are critical. In this work, using in situ transmission electron microscopy and phase-field simulations, a controlled nucleation of vertically oriented 109° domain walls in (110)-oriented BiFeO 3 (BFO) thin films is reported. In the switching experiment, reversed domains that are nucleated preferentially at the nanoscale edges of the “crest and sag” pattern-like electrode under external bias subsequently grow into a stable stripe configuration. In addition, when triangular pockets (with an in-plane polarization component) are present, these domain walls are pinned to form stable flux-closure domains. Phase field simulations show that i) field enhancement at the edges of the electrode causes site-specific domain nucleation, and ii) the local electrostatics at the domain walls drives the formation of flux closure domains, thus stabilizing the striped pattern, irrespective of the initial configuration. The results demonstrate how flux closure pinning can be exploited in conjunction with electrode patterning and substrate orientation to achieve a desired topological defect configuration. These insights constitute critical advancements in exploiting domain walls in next generation ferroelectronic devices.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
FIMP dark matter from flavon portals
We investigate the phenomenology of a non-thermal dark matter (DM) candidate in the context of flavor models that explain the hierarchy in the masses and mixings of quarks and leptons via the Froggatt-Nielsen (FN) mechanism. A flavor-dependent U(1) FN symmetry explains the fermion mass and mixing hierarchy, and also provides a mechanism for suppressed interactions of the DM, assumed to be a Majorana fermion, with the Standard Model (SM) particles, resulting in its FIMP (feebly interacting massive particle) character. Such feeble interactions are mediated by a flavon field through higher dimensional operators governed by the U(1) FN charges. We point out a natural stabilizing mechanism for the DM within this framework with the choice of half-integer U(1) FN charge n for the DM fermion, along with integer charges for the SM fermions and the flavon field. In this flavon portal scenario, the DM is non-thermally produced from the decay of the flavon in the early universe which becomes a relic through the freeze-in mechanism. We explore the allowed parameter space for this DM candidate from relic abundance by solving the relevant Boltzmann equations. We find that reproducing the correct relic density requires the DM mass to be in the range (100 – 300) keV for n = 7.5 and (3 – 10) MeV for n = 8.5 where n is the U(1) FN charge of the DM fermion.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Plasma Facing Components with Capillary Porous System and Liquid Metal Coolant Flow
Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [Khodak and Maingi (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile, as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a Fusion Nuclear Science Facility [Kessel et al. (2019)]: 10T toroidal field and 10 MW/m2 peak incident heat flux. Computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.
High Heat Flux↗
Large-eddy simulation of non-vaporizing sprays using the spectral-element method
Predictive simulations of high-pressure sprays require accurate representation of the turbulent gaseous flow field generated by liquid jet. Typically, the accuracy that can be obtained with low-order numerical methods (e.g. finite volume, finite element) is limited by stability issues in fine grids and the order of convergence of the method. In this work, we resolve the turbulent flow field in an Eulerian manner using the high-order spectral element method, coupled with a Lagrangian parcels approach to model the atomizing liquid jet. Large eddy simulations of single-hole sprays under non-evaporative conditions were conducted and compared against experimental data from Margot et al. (2008) and Spray A data from the Engine Combustion Network. The sensitivity of liquid penetration and droplet sizes to different breakup model parameters was studied. The effect of different numerical parameters, such as polynomial order of the solution (grid resolution), on liquid penetration was also analyzed. The method achieved grid-independent results using p-refinement, achieving finer resolution (by a factor of x 1.7 - x 3.5) in the gas-phase solution than in state-of-the-art simulations using the finite-volume method. Here the results showed good agreement with experimental data, demonstrating the ability of the current method to accurately capture liquid penetration and the shape of the spray.
97 MATHEMATICS AND COMPUTING↗
Boundary conditions manipulation of polar vortex domains in BiFeO 3 membranes via phase-field simulations
Polar vortex domains have recently become an emergent research field due to the abundant physical phenomena and potential applications in high-density memories. Here, we explore the mechanisms of creating polar vortex domains in the BiFeO3 (BFO) membranes subjected to different boundary conditions using phase-field simulations. A major difference is that the vortex in membrane can be stabilized even under short-circuit electrical boundary conditions compared to vortex in other systems, such as thin film or superlattice. We found that (a) the formation of polar vortex domains at the membrane interior under bending is mainly driven by the reduction of elastic energy under short-circuit boundary condition, and the vortex chirality (namely, clockwise and counterclockwise) could be identified by n-shape and u-shape bending; (b) in the unbent open-circuit BFO membrane case, exotic trapezoid-shaped vortex nanodomains form at the terminations of 109 degrees domain walls (DWs) and partially charged 71 degrees DWs, which is driven by the local depolarization field and the interplay among electrostatic, elastic, and gradient and Landau energies. We also examine Kittel's law by establishing the dependence of vortex periods on the membrane thickness. Furthermore, these results give further understanding of the effect of boundary conditions on the formation of polar vortex domains, guiding experimental designs of vortex-based high-density memories.
36 MATERIALS SCIENCE↗
Linear instabilities in the hot-ion regime in a high-field spherical tokamak
Abstract The present operating high-field compact spherical tokamak ST40 is an important step toward an ST-based fusion reactor (Gryaznevich and Asunta 2017 Fusion Eng. Des. 123 177; McNamara 2023 Nucl. Fusion 63 054002). Temperature and density profiles and their uncertainties in recent hot ion ST40 plasmas with central ion temperatures exceeding 8.6 keV, have been inferred using an integrated analysis of several diagnostics including line-of-sight integration and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. A linear gyrokinetic stability analysis has been carried out to identify the most unstable micro-instabilities in these hot ion plasmas. In one of these plasmas, it is found that linear growth rates of both ion- and electron-scale ( k θ ρ s ⩾ 0.2 where k θ is the poloidal wavenumber and ρ s is the ion gyro-radius with sound speed) modes decrease from the edge toward the core of the plasma (i.e. from ρ = 0.8 to 0.3 where ρ is the square root of the normalized magnetic flux), while the change of linear growth rates at k θ ρ s < 0.2 is non-monotonic. In particular, at ρ = 0.3 no unstable mode was found at k θ ρ s > 5 , and about an order of magnitude reduction in the maximum linear growth rate in the ion-scale wavenumber range of 0.2 ⩽ k θ ρ s < 1 is seen at ρ = 0.3 compared with the ρ = 0.8 location. It is found through parametric scans that while the ion temperature gradient (ITG) mode/trapped electron mode (TEM) and kinetic ballooning mode (KBM)/kinetic shear Alfvén (KSA) mode are more important in the plasma core, ubiquitous mode (UM) is found to be the important ion-scale instability at ρ ⩾ 0.6 . At ρ = 0.8, two instabilities are found at the electron scale, with one being the typical electrostatic electron temperature gradient (ETG) mode at ρ e scale and beyond and the other being an unidentified low-real-frequency instability (at the intermediate scale between electron and ion gyroradii) which can propagate in the ion direction and is driven by ETG. The existence of UM at the ion scale and the low-real-frequency instability and ETG mode at the electron scale at ρ = 0.8 is supported by an extensive linear gyrokinetic stability analysis of another ST40 hot ion plasma in a wider parametric range, taking into account the uncertainties in the experimental profiles.
ubiquitous mode↗
Electrical Potential Investigation of Reversible Metastability and Irreversible Degradation of CdTe Solar Cells
In this paper, we report on the stability of CdTe devices with a structure of TCO/MZO/CdSeTe/CdTe/back-contact. The device showed reversible transitions between the light-soak state (LSS) with the best device efficiency and the dark-soak state (DSS) with an inferior efficiency. However, it showed an irreversible degradation state (DgS) driven by long-hour light soaking at an elevated temperature. We have investigated transitions between these three states from the perspective of the electric field by nm-resolution potential imaging across the devices using Kelvin probe force microscopy (KPFM). The results exhibit different anomalous electric field profiles. At the LSS, the electric field exhibits a main peak inside the CdSeTe layer instead of the MZO/CdSeTe heterointerface, illustrating a buried homojunction (BHJ) of the device. At the DSS, a large electric field peak at the MZO/CdSeTe interface was measured, which probably resulted in the inferior fill factor at the DSS. At the DgS, the electric field peak at the MZO/CdSeTe interface increased further and a third electric field was measured at the back contact of the device. Device modeling using COMSOL software, in alignment with both the electric field and device current-voltage curves, elucidates that a low n-doped CdSeTe in the region near the MZO/CdSeTe interface caused the BHJ in the LSS and a loss of MZO doping and/or increase of the conduction band offset spike due to long-term stress caused the increased electric field near the MZO/CdSeTe interface at the DgS.
CdTe solar cell↗