Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “One-dimensional Model”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Theory and numerics of subspace approximation of eigenvalue problems

Large-scale eigenvalue problems arise in various fields of science and engineering and demand computationally efficient solutions. In this study, we investigate the subspace approximation for parametric linear eigenvalue problems, aiming to mitigate the computational burden associated with high-fidelity systems. Furthermore, we provide general error estimates under non-simple eigenvalue conditions, establishing some theoretical foundations for understanding the convergence behavior of subspace approximations. Numerical examples, including problems with one-dimensional to three-dimensional spatial domain and one-dimensional to two-dimensional parameter domain, are presented to demonstrate the efficacy of reduced basis method in handling parametric variations in boundary conditions and coefficient fields to achieve significant computational savings while maintaining high accuracy, making them promising tools for practical applications in large-scale eigenvalue computations.

Eigenvalue problems↗
DOE OSTI2025

Isochronous and period-doubling diagrams for symplectic maps of the plane

Symplectic mappings of the plane serve as key models for exploring the fundamental nature of complex behavior in nonlinear systems. Central to this exploration is the effective visualization of stability regimes, which enables the interpretation of how systems evolve under varying conditions. While the area-preserving quadratic Hénon map has received significant theoretical attention, a comprehensive description of its mixed parameter-space dynamics remain lacking. This limitation arises from early attempts to reduce the full two-dimensional phase space to a one-dimensional projection, a simplification that resulted in the loss of important dynamical features. Consequently, there is a clear need for a more thorough understanding of the underlying qualitative aspects. This paper aims to address this gap by revisiting the foundational concepts of reversibility and associated symmetries, first explored in the early works of G.D. Birkhoff. We extend the original framework proposed by Hénon by adding a period-doubling diagram to his isochronous diagram, which allows to represents the system’s bifurcations and the groups of symmetric periodic orbits that emerge in typical bifurcations of the fixed point. A qualitative and quantitative explanation of the main features of the region of parameters with bounded motion is provided, along with the application of this technique to other symplectic mappings, including cases of multiple reversibility. Modern chaos indicators, such as the Reversibility Error Method (REM) and the Generalized Alignment Index (GALI), are employed to distinguish between various dynamical regimes in the mixed space of variables and parameters. These tools prove effective in differentiating regular and chaotic dynamics, as well as in identifying twistless orbits and their associated bifurcations. Additionally, we discuss the application of these methods to real-world problems, such as visualizing dynamic aperture in accelerator physics, where our findings have direct relevance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2025

Driven Majorana modes: A route to synthetic $p_x +ip_y$ superconductivity

We propose a protocol to realize synthetic $p_x +ip_y$ superconductors in one-dimensional topological systems that host Majorana fermions. By periodically driving a localized Majorana mode across the system, our protocol realizes a topological pumping of Majorana fermions, analogous to the adiabatic Thouless pumping of electrical charges. Importantly, similar to the realization of a Chern insulator through Thouless pumping, we show that pumping of Majorana zero modes could lead to a $p_x +ip_y$ superconductor in the two dimensions of space and synthetic time. The Floquet theory is employed to map the driven one-dimensional system to a two-dimensional synthetic system by considering frequency as a new dimension. We demonstrate such Floquet $p_x +ip_y$ superconductors using the Kitaev $p$-wave superconductor chain, a prototypical 1D topological system, as well as its more realistic realization in the 1D Kondo lattice model as examples. We further show the appearance of Majorana $π$ mode at the Floquet zone boundary in an intermediate drive frequency region. Our work suggests a driven magnetic spiral coupled to a superconductor as a promising platform for the realization of novel topological superconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Feasibility Study on Implementing a Staggered-Grid Finite Volume Method for System Analysis Code Development Under the MOOSE Framework

Here, this work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key to the test bed is the implementation of high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. The test bed utilized a more flexible code structure to enable the finite volume method implementation and direct interacting with the solver package, instead of using the natively supported finite element method by the framework. Using a suite of selected test problems with different problem sizes and levels of complexity, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. For a complex reactor model, transient simulation was performed using the newly developed finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development.

MOOSE↗
DOE OSTI2026

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗
DOE OSTI2025

Band-Selective Spin-Charge Separation across the Charge Density Wave Transition in Quasi-1D NbSe 3

Non-Fermi liquid (non-FL) phase is a pivotal enigma in understanding intriguing quantum phases in strongly correlated systems, such as high-temperature superconductivity. Tomonaga-Luttinger liquid (TLL) theory, designed for one-dimensional (1D) systems, serves as one of the microscopic frameworks that elucidates non-FL behavior. Despite its theoretical concreteness, comprehensive experimental verification has remained incomplete. In particular, addressing the persistence of the TLL nature within ordered phases, such as charge density wave (CDW), has posed a significant challenge. We report the observation of TLL characteristics across the CDW transitions in a quasi-1D material NbSe 3 using angle-resolved photoemission spectroscopy. Spin-charge separation, a disentanglement of spin and charge degrees of freedom of electrons, is clearly observed within the band, accompanied by the suppression of spectral weight following a power-law behavior with anticipated temperature scaling. Surprisingly, the spin-charge separation persists even below the CDW transition temperatures, indicating the validity of the TLL nature in the CDW phase. In conclusion, our findings offer a unique opportunity to explore the interplay between the TLL and ordered phases, which could be connected to the interplay between non-FL and ordered phases in strongly correlated systems.

1-dimensional systems↗
DOE OSTI2026

Intertwined charge, spin, and orbital degrees of freedom under electronic correlations in the one-dimensional Fe 3+ chalcogenide chain

Motivated by recent developments in the study of quasi-one-dimensional iron systems with Fe 2+ , we comprehensively study an Fe 3+ chalcogenide chain system. Based on first-principles calculations, the Fe 3+ chain has a similar electronic structure to that discussed before for the Fe 2+ chain, because of the similar Fe⁢𝑋 4 (𝑋 = S or Se) tetrahedron-chain geometry. Furthermore, a three-orbital electronic Hubbard model for this chain was constructed using the density matrix renormalization group method. A robust antiferromagnetic coupling was unveiled in the chain direction. In addition, in the intermediate electronic correlation 𝑈/𝑊 region, we found an interesting orbital-selective Mott phase with the coexistence of localized and itinerant electrons (𝑈 is the on-site Hubbard repulsion, while 𝑊 is the electronic bandwidth) based on the orbital-selective behavior observed in the charge fluctuations. Furthermore, we do not observe any obvious pairing tendency in the Fe 3+ chain in the electronic-correlation 𝑈/𝑊 region, where superconducting pairing tendencies were reported before in iron ladders. This suggests that superconductivity is unlikely to emerge in the Fe 3+ systems. Finally, our results clearly establish the similarities and differences between Fe 2+ and Fe 3+ iron chains, as well as iron ladders.

density matrix renormalization group↗
DOE OSTI2025

Quasi-1D Coulomb Drag in the Nonlinear Regime

One-dimensional Coulomb drag has been an essential tool to probe the physics of interacting Tomonaga-Luttinger liquids. To date, most experimental work has focused on the linear regime while the predictions for Luttinger liquids beyond the linear response theory remain largely untested. In this Letter, we report measurements of reciprocal momentum transfer induced Coulomb drag between vertically coupled quasi-one-dimensional quantum wires in the nonlinear regime. Measurements were performed at ultralow temperatures between wires only 15 nm apart. Our results reveal a nonlinear dependence of the drag voltage as a function of the drive current superimposed with an oscillatory contribution, in agreement with theoretical predictions for Coulomb drag between Tomonaga-Luttinger liquids. Additionally, the observed current-voltage characteristics exhibit a nonmonotonic temperature dependence, further corroborating the presence of non-Fermi-liquid behavior in our system. In conclusion, these findings are observed both in the single and in the multiple subband regimes and in the presence of disorder, extending the onset of this behavior beyond the clean single channel Tomonaga-Luttinger regime where the predictions were originally formulated.

1-dimensional systems↗
DOE OSTI2025

Advanced System Thermal Fluids Solver Development for SAM

This work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key is the implementation of a high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. Leveraging the existing capabilities of the SAM code, significant code coverages were established in the finite volume method code. This in turn allows for a suite of test problems with different problem sizes and levels of complexity to be used to quantify the performance improvement of the finite volume method code. As evidently shown in this study, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE for the wide range of selected problems. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. In addition, for a complex reactor model, transient simulation was performed using the finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development. In this work, short-term priority development and testing items were identified, and long-term code adoption and integration plans were made for the eventual deployment of the finite volume method in the SAM code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A North Sea in Situ Evaluation of the Fitch Wind Farm Parameterization Within the Mellor-Yamada-Nakanishi-Niino and 3D Planetary Boundary Layer Schemes

Wind resource assessments and wind power forecasts that account for wind farm wakes are sensitive to the choice of planetary boundary layer (PBL) scheme. This work compares the one-dimensional Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme with a three-dimensional PBL (3DPBL) scheme, evaluating predictions made with both schemes against two sets of North Sea in situ observations of wind farm wakes. The optimal PBL scheme varies based on the observations (FINO1 tower vs. aircraft), the quantity of interest (wind speed vs. turbulence kinetic energy [TKE]), and the error metric (bias, centered root mean square error [cRMSE], R2, and earth mover's distance [EMD]). Whereas 3DPBL wind speeds outperform MYNN wind speeds with respect to the cRMSE at the FINO1 site located at a single point within the turbine rotor layer, 3DPBL TKE bias is larger than MYNN TKE bias when compared to aircraft observations taken 100 m above a wind farm. Wind speeds in the aircraft region are ambiguous with regard to which PBL scheme is optimal. Aircraft MYNN wind speeds outperform 3DPBL wind speeds with respect to R2 and cRMSE but underperform with respect to bias and EMD. Future evaluations across broader temporal and spatial scales may offer further insight into model differences.

17 WIND ENERGY↗
NASA NTRS1971

Transient-forced convection film boiling on an isothermal flat plate

An approach for the solution of transient-forced convection film boiling on an isothermal flat plate using the boundary layer model is developed. The similarity variables are used to convert the governing partial differential equations to ordinary ones. The results of numerical solutions of these ordinary equations indicate that the transient process can be classified as one-dimensional conduction, intermediate, and the steady state regions. The time required for the one-dimensional conduction and the time necessary to attain a steady state condition are obtained. The influence of interfacial shear is seen to be negligible while the Prandtl Number and the ratio (C sub p delta T divided by h sub fg times Pr) have major influence. The use of local similarity approximations for the intermediate regime facilitates prediction of complete boundary layer growth. Using the ratio of time at any instant to the steady state time as abscissa, the curves representing the boundary layer growth can be merged into a single mean curve within 5 percent. Further, the analysis shows that the average rate of heat transfer during transient is 50 to 100 percent higher than those at steady state. The average rate of vapor convected away is 10 to 15 percent lower than at steady state while the average rate of accumulation to form the vapor layer is 1 to 14 times larger. Further, the total heat transferred during transient increases and the evaporation decreases for increasing values of C sub p delta T divided by h sub fg times Pr

Nagendra, H. R.↗