Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Analysi”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Flashback characterization of additively manufactured swirl-stabilized fuel injector with varying surface roughness

This work investigates the effects of as-built surface roughness on the flashback propensity in additively manufactured (AM) swirl-stabilized lean premixed (LPM) fuel injectors. Adoption of AM for rapid prototyping and fabrication of complex fuel-flexible injector designs requires investigating surface roughness effects on the flow and flame stability characteristics of the combustor. Wall roughness increases the near-wall shear, which could alter the boundary layer structure and change the propensity for flame flashback. Accounting for the realistic as-built surface roughness is crucial in carrying out computational modeling and experimental analysis to establish a feedback loop for the precision designing of fuel-flexible injectors. The presented numerical analysis of as-built AM injector is an essential consideration to optimize injector design for aerodynamics and additive manufacturability.

Mohanty, Pratikshya↗
DOE OSTI2021

Excitation spectra of quantum matter without quasiparticles. I. Sachdev-Ye-Kitaev models

Here, we study the low-frequency spectra of complex Sachdev-Ye-Kitaev models at general densities. The analysis applies also to SU(M) magnets with random exchange at large M. The spectral densities are computed by numerical analysis of the saddle-point equations on the real frequency ω axis at zero temperature T. The asymptotic low-ω behaviors are found to be in excellent agreement with the scaling dimensions of irrelevant operators which perturb the conformally invariant critical states. Of possible experimental interest is our computation of the universal spin spectral weight of the SU(M) magnets at low ω and T; this includes a contribution from the time reparametrization mode, which is the boundary graviton of the holographic dual. This analysis is extended to a random t-J model in the following paper [M. Tikhanovskaya et al., following paper, Phys. Rev. B 103, 075142 (2021)].

36 MATERIALS SCIENCE↗
DOE OSTI2026

DEMO-FTES: Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (CRADA Final Report)

The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.

25 ENERGY STORAGE↗
DOE OSTI2023

Dynamic Mooring Field Experiment and Design of a Wave Energy Converter Platform Test System

Abstract This study evaluated a multi-catenary spread mooring system design of a mobile ocean test berth for wave energy converters (WECs), the Ocean Sentinel (OS) instrumentation buoy, through dynamic simulation, numerical analysis, and comparison with measured field motion data of the OS off the OR coast. First, the accuracy of the numerical employed model based on a fully-coupled method of the OS and its mooring lines was validated by comparing predicted mooring tensions to the field measurements. Then, the anchor movability, fatigue damage, and extreme mooring tension of the OS mooring system were analyzed to assess survivability. Field test results show that the numerical model provided accurate predictions of mooring tensions even under environmental conditions of strong wind, current, and waves. Factors affecting the accuracy are discussed. One mooring anchor was shown to have moved significantly during the ocean field test. Mooring fatigue damage was calculated for different levels of sea states. Design strengths of the mooring lines were calculated and analyzed.

Engineering↗
DOE OSTI2023

Crack-Free Joining of High-Strength AA7055 Sheets by Friction Based Self-Piercing Riveting with the Aid of Numerical Design

Unique friction-based self-piercing riveting (F-SPR) was employed to join high-strength, low-ductility aluminum alloy 7055 for lightweight vehicle applications. This study aimed to maximize the joint strength of the AA7055 F-SPR joint while avoiding cracking issues due to low ductility at room temperature. A fully coupled Eulerian–Lagrangian (CEL) model was employed to predict the process temperature during F-SPR, and the temperature field was then mapped onto a 2D axisymmetric equivalent model for accelerated numerical analysis. The geometry, dimensions, and material strength of the rivet, as well as the depth of the die cavity and plunging depth, were investigated to enhance joint formation. Also, a static finite-element analysis model was developed to predict and analyze the stress distribution in the rivet under different mechanical testing loading conditions. Overall, the numerical model showed good agreement with the experiment results, such as joint formation and mechanical joint strength. With the aid of virtual fabrication through numerical modeling, the joint design iterations and process development time of F-SPR were greatly reduced regarding the goal of lightweight, high-strength aluminum joining.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Numerical and experimental analysis of precooled desiccant system

This paper presents a numerical and experimental study of four novel, desiccant air-conditioning systems based on pre-cooling of the ambient air. The purpose of the study is to: 1) establish which configuration is most profitable for the future studies through numerical and experimental methods; 2) verify the concept of precooled system and confirm that it is working in experimental conditions (for this purpose two prototypes of evaporative cooling units were developed and compared for different operation arrangements with desiccant wheel); 3) Analyze practical aspects of the systems, which includes water consumption analysis, economic analysis and establishing operational modes during a cooling season in temperate climate conditions (the proposed solution has a significant advantage during operation in a temperate climate, since it can utilize its first evaporative cooler the only source of cooling energy when the humidity of outdoor air allows to to meet the latent cooling load). Four arrangements of the system with different dew point evaporative air cooler prototypes were selected and developed for comparative analysis. The system was simulated as a whole, also, each component was individually tested on measurement station. develop the optimal arrangement of the system through numerical and experimental analysis. Four different performance indicators were considered: thermal COP (Coefficient Of Performance), the ERR (Energy Efficiency Ratio), coefficient a describing the potential in compensating indoor cooling loads, water consumption ratio. Also costs of the systems were compared through economic analysis. It was established that arrangement which is characterized by best-averaged performance in terms of all considered factors is a unit combined with two different types of evaporative coolers (marked as System B)- this is caused by the favorable arrangement of two units, where one is better suited to the pre-cooler role, while the other is better suited as the post-cooler.

42 ENGINEERING↗
DOE OSTI2023

Effects of additively manufactured surface roughness on small diameter plenums with multiple side discharges for high temperature gas turbine blade cooling applications

Combined heat and power (CHP) applications have significant environmental and economic benefits that are consistent with the goals of the U.S. Department of Energy (DOE). One area that is currently being studied includes the potential benefits of CHP turbine operation at higher turbine inlet temperatures. Internal cooling concepts enabled by additive manufacturing (AM) are of primary interest. Here, the effectiveness of internal cooling is hindered by many factors such as velocity distribution of the cooling air to the hot surface considering impingement cooling. To simulate cooling air exiting from a series of orifices for internal cooling in an airfoil, a straight smooth wall tubing with multiple side discharging orifices is used and compared to additively manufactured tubing (Ti6Al4V Grade 23) with orifice size and spacing as well as inner and outer diameters identical to the smooth wall tubing. Similar to flow discharging form perforated pipes, the flow discharged from individual orifices along the tubes in this study is found to be nonuniformly distributed, and the horizontal (axial direction) momentum can be observed from the experimental data. Discharge velocities have been measured with two-dimensional (2D) particle imaging velocimetry (PIV) and single element hot wire anemometry. Numerical analysis has also been conducted to predict the velocity distributions along the orifices in the smooth wall and additive manufacturing (AM) tubing, which are inherited with surface roughness. Numerical and measured results in this study are compared, presented, and discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Quantitative analysis of Raman spectral parameters for carbon fibers: practical considerations and connection to mechanical properties

Although the literature on the Raman spectra of carbon fibers is vast, no consistent, robust predictive relationship between mechanical properties of carbon fibers and spectral parameters exists. This shortcoming is due to the use of numerous fitting functions to evaluate Raman spectra of carbon fibers and the inconsistencies in establishing the best fitting models in a statistically robust fashion. To address this gap, we present a comprehensive work on the Raman spectra of carbon fibers that combines a vast library of experimental data with a robust numerical analysis and a statistical evaluation of a wide range of suggested fitting models. This manuscript begins with a brief review of the commonly applied fitting models. Then, the Raman spectra of 32 commercially available polyacrylonitrile-based carbon fibers collected at excitation wavelengths 532, 633, and 785 nm are presented and the best fit for all fibers is evaluated based on several statistical criteria in conjunction with numerical calculations and physical arguments. The results suggest that high-performance fibers must be fit with at least five peaks, whereas high-tensile modulus fibers are best fit with at least six distinct peaks. Finally, we employ simultaneous fitting of the Raman spectra of specific fibers and wavelengths and demonstrate that strong correlations exist between mechanical properties and the D1 peak position and shape across the range of evaluated mechanical properties. Finally, we suggest straightforward improvements in fitting analysis procedures that can be implemented to increase coherency in the understanding of the underlying carbon fiber microstructure intuited from Raman spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Exploring orbital angular momentum and spin-orbit correlations for gluons at the Electron-Ion Collider

In our previous work [S. Bhattacharya , ], we introduced a pioneering observable aimed at experimentally detecting the orbital angular momentum (OAM) of gluons. Our focus was on the longitudinal double spin asymmetry observed in exclusive dijet production during electron-proton scattering. We demonstrated the sensitivity of the cos ϕ angular correlation between the scattered electron and proton as a probe for gluon OAM at small x and its intricate interplay with gluon helicity. This current work provides a comprehensive exposition, diving further into the aforementioned calculation with added elaboration and in-depth analysis. We reveal that, in addition to the gluon OAM, one also gains access to the spin-orbit correlation of gluons. We supplement our work with a detailed numerical analysis of our observables for the kinematics of the Electron-Ion Collider. In addition to dijet production, we also consider the recently proposed semi-inclusive diffractive deep inelastic scattering process, which potentially offers experimental advantages over dijet measurements. Finally, we investigate quark-channel contributions to these processes and find an unexpected breakdown of collinear factorization. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Solving the homogeneous Bethe-Salpeter equation with a quantum annealer

The homogeneous Bethe-Salpeter equation (hBSE), describing a bound system in a genuinely relativistic quantum-field theory framework, was solved for the first time by using a D-Wave quantum annealer. After applying standard techniques of discretization, the hBSE, in ladder approximation, can be formally transformed in a generalized eigenvalue problem (GEVP), with two square matrices: one symmetric and the other nonsymmetric. The latter matrix poses the challenge of obtaining a suitable formal approach for investigating the GEVP by means of a quantum annealer, i.e., to recast it as a quadratic unconstrained binary optimization problem. A broad numerical analysis of the proposed algorithms, applied to matrices of dimension up to 64, was carried out by using both the simulated-annealing package and the D-Wave . The numerical results very nicely compare with those obtained with standard classical algorithms, and also show interesting scalability features. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Quantum Control Landscapes Beyond the Dipole Approximation: Controllability, Singular Controls, and Resources

We investigate the control landscapes of closed n-level quantum systems beyond the dipole approximation by including a polarizability term in the Hamiltonian. The latter term is quadratic in the control field. Theoretical analysis of singular controls is presented, which are candidates for producing landscape traps. The results for considering the presence of singular controls are compared to their counterparts in the dipole approximation (i.e., without polarizability). A numerical analysis of the existence of traps in control landscapes for generating unitary transformations beyond the dipole approximation is made upon including the polarizability term. An extensive exploration of these control landscapes is achieved by creating many random Hamiltonians which include terms linear and quadratic in a single control field. The discovered singular controls are all found not to be local optima. This result extends a great body of recent work on typical landscapes of quantum systems where the dipole approximation is made. We further investigate the relationship between the magnitude of the polarizability and the fluence of the control resulting from optimization. It is also shown that including a polarizability term in an otherwise uncontrollable dipole coupled system removes traps from the corresponding control landscape by restoring controllability. We numerically assess the effect of a polarizability term on a known example of a particular three-level Λ-system with a second order trap in its control landscape. It is found that the addition of the polarizability removes the trap from the landscape. The general practical control implications of these simulations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Can Error Mitigation Improve Trainability of Noisy Variational Quantum Algorithms?

Variational Quantum Algorithms (VQAs) are often viewed as the best hope for near-term quantum advantage. However, recent studies have shown that noise can severely limit the trainability of VQAs, e.g., by exponentially flattening the cost landscape and suppressing the magnitudes of cost gradients. Error Mitigation (EM) shows promise in reducing the impact of noise on near-term devices. Thus, it is natural to ask whether EM can improve the trainability of VQAs. In this work, we first show that, for a broad class of EM strategies, exponential cost concentration cannot be resolved without committing exponential resources elsewhere. This class of strategies includes as special cases Zero Noise Extrapolation, Virtual Distillation, Probabilistic Error Cancellation, and Clifford Data Regression. Second, we perform analytical and numerical analysis of these EM protocols, and we find that some of them (e.g., Virtual Distillation) can make it harder to resolve cost function values compared to running no EM at all. As a positive result, we do find numerical evidence that Clifford Data Regression (CDR) can aid the training process in certain settings where cost concentration is not too severe. Our results show that care should be taken in applying EM protocols as they can either worsen or not improve trainability. On the other hand, our positive results for CDR highlight the possibility of engineering error mitigation methods to improve trainability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Pointer States and Quantum Darwinism with Two-Body Interactions

Quantum Darwinism explains the emergence of classical objectivity within a quantum universe. However, to date, most research on quantum Darwinism has focused on specific models and their stationary properties. To further our understanding of the quantum-to-classical transition, it appears desirable to identify the general criteria a Hamiltonian has to fulfill to support classical reality. To this end, we categorize all N-qubit models with two-body interactions, and show that only those with separable interaction of the system and environment can support a pointer basis. We further demonstrate that “perfect” quantum Darwinism can only emerge if there are no intra-environmental interactions. Our analysis is complemented by solving the ensuing dynamics. We find that in systems exhibiting information scrambling, the dynamical emergence of classical objectivity directly competes with the non-local spread of quantum correlations. Our rigorous findings are illustrated through the numerical analysis of four representative models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Computational capacity in hydrodynamic real-time hybrid simulation applied to simulate the dynamic response of floating offshore wind turbines

Real-time hybrid simulation (RTHS) mitigates similitude distortions in model-scale tests of floating offshore wind turbines (FOWTs) by coupling physical experiments with numerical models in real time. The coupling requires faster-than-real-time numerical computations to satisfy temporal similitude with the physical experiment, presenting a bottleneck for using more complex numerical models in RTHS. This paper presents a hydrodynamic-RTHS (hydro-RTHS) framework for FOWTs that simulates the hydrodynamics physically and the aerodynamics numerically with sensor feedback from the physical testing. The framework adapts the three-loop hardware architecture to leverage greater computational resources and mitigate strict temporal requirements, enabling more computationally demanding numerical analyses in hydro-RTHS. The three-loop hardware architecture integrates multiple machines, each dedicated to either numerical analysis or RTHS controls, with a rate-transition algorithm to synchronize the tasks executed across the different machine processors. Virtual and physical tests verified and validated the hydro-RTHS framework, respectively. The ”virtual” tests, which approximates the physical domain numerically, verified the RTHS framework with respect to a numerical full-scale complete FOWT model simulated in the open-source software, OpenFAST. The virtual tests were able to maintain comparable control signals while enabling greater computational resources for the numerical calculations. Real-world physical tests demonstrated that the hydro-RTHS framework computes aerodynamic forces similar to the complete OpenFAST model, validating the hydro-RTHS framework using the three-loop hardware architecture. Findings show that the hydro-RTHS framework with the three-loop hardware architecture is computationally efficient, with reserve capacity to simulate more complex problems due to the customized software, hardware, and rate-transition algorithm.

17 WIND ENERGY↗
DOE OSTI2026

Discrete Superconvergence Analysis for Quantum Magnus Algorithms of Unbounded Hamiltonian Simulation

Motivated by various applications, unbounded Hamiltonian simulation has recently garnered great attention. Quantum Magnus algorithms, designed to achieve commutator scaling for time-dependent Hamiltonian simulation, have been found to be particularly efficient for such applications. When applied to unbounded Hamiltonian simulation in the interaction picture, they exhibit an unexpected superconvergence phenomenon. However, existing proofs are limited to the spatially continuous setting and do not extend to discrete spatial discretizations. Here, in this work, we provide the first superconvergence estimate in the fully discrete setting with a finite number of spatial discretization points N, and show that it holds with an error constant uniform in N. The proof is based on the two-parameter symbol class, which, to our knowledge, is applied for the first time in algorithm analysis. The key idea is to establish a semiclassical framework by identifying two parameters through the discretization number and the time step size rescaled by the operator norm, such that the semiclassical uniformity guarantees the uniformity of both. This approach may have broader applications in numerical analysis beyond the specific context of this work.

Borns-Weil, Yonah [University of California, Berke↗
DOE OSTI2022

A novel approach for tetrahedral-element-based finite element simulations of anisotropic hyperelastic intervertebral disc behavior

Intervertebral discs are microstructurally complex spinal tissues that add greatly to the flexibility and mechanical strength of the human spine. Attempting to provide an adjustable basis for capturing a wide range of mechanical characteristics and to better address known challenges of numerical modeling of the disc, we present a robust finite-element-based model formulation for spinal segments in a hyperelastic framework using tetrahedral elements. We evaluate the model stability and accuracy using numerical simulations, with particular attention to the degenerated intervertebral discs and their likely skewed and narrowed geometry. To this end, 1) annulus fibrosus is modeled as a fiber-reinforced Mooney-Rivlin type solid for numerical analysis. 2) An adaptive state-variable dependent explicit time step is proposed and utilized here as a computationally efficient alternative to theoretical estimates. 3) Tetrahedral-element-based FE models for spinal segments under various loading conditions are evaluated for their use in robust numerical simulations. For flexion, extension, lateral bending, and axial rotation load cases, numerical simulations reveal that a suitable framework based on tetrahedral elements can provide greater stability and flexibility concerning geometrical meshing over commonly employed hexahedral-element-based ones for representation and study of spinal segments in various stages of degeneration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Water desalination through the dewpoint evaporative system

This paper presents the first numerical analysis of a novel Dew Point Desalination unit (DPD), which utilizes the phenomenon of the evaporative cooling process to minimize energy consumption. The proposed solution is based on a modular structure, which allows it to obtain any required water production capacity. The analyses for various operational parameters and various climatic conditions are based on the numerical simulations conducted using a validated mathematical model. The system performance was described with different factors, including specific energy consumption (SEC) and daily water production rate respected to 1 m^3*s^-1 of passing air stream (DWP). It is shown that the DPD system can operate in almost every climate zone, and the desalinated water production consumes less than 1.5 kWh*m^-3 of electric energy It is noted the best performance of the DPD system is achieved in the semi-arid, desert and subtropical climate zones. In these climate zones, the 2-stage DPD system allows the production of desalinated water with an electric energy consumption SEC lower than 0.5 kWh*m^-3. The additional benefit of the DPD system is its modularity. The presented DPD system stands as a module, that can be multiplied as needed to achieve the desired size of the desalination plant.

Desalination↗
DOE OSTI2022

Measuring sub-surface spatially varying thermal conductivity of silicon implanted with krypton

The thermal properties of semiconductors following exposure to ion irradiation are of great interest for the cooling of electronic devices; however, gradients in composition and structure due to irradiation often make the measurement difficult. Furthermore, the nature of spatial variations in thermal resistances due to spatially varying ion irradiation damage is not well understood. In this work, we develop an advancement in the analysis of time-domain thermoreflectance to account for spatially varying thermal conductivity in a material resulting from a spatial distribution of defects. We then use this method to measure the near-surface (≲1 μm) thermal conductivity of silicon wafers irradiated with Kr + ions, which has an approximate Gaussian distribution centered 260 nm into the sample. Our numerical analysis presented here allows for the spatial gradient of thermal conductivity to be extracted via what is fundamentally a volumetric measurement technique. We validate our findings via transmission electron microscopy, which is able to confirm the spatial variation of the sub-surface silicon structure, and provide additional insight into the local structure resulting from the effects of ion bombardment. Thermal measurements found the ion stopping region to have a nearly 50x reduction in thermal conductivity as compared to pristine silicon, while TEM showed the region was not fully amorphized. Our results suggest this drastic reduction in silicon thermal conductivity is primarily driven by structural defects in crystalline regions along with boundary scattering between amorphous and crystalline regions, with a negligible contribution being due to implanted krypton ions themselves.

36 MATERIALS SCIENCE↗