Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermodynamic Modeling of Heat Engines Including Heat Transfer and Compression–Expansion Irreversibilities

In this paper, a thermodynamic model based on an endoreversible engine approach is developed to analyze the performance of heat engines operating under different thermodynamic cycles. The model considers finite heat transfer rate, variable heat source and sink temperatures, and irreversibilities associated with the expansion and compression. Expressions for the maximum power and efficiency at maximum power output are obtained as a function of hot and cold reservoir temperatures, the equivalent isentropic efficiency of compression and expansion components, and the effective conductance ratio between heat exchangers. In all cases, the Curzon–Ahlborn efficiency is retrieved at constant reservoir temperatures and neglected compression–expansion irreversibilities. The proposed model allows assessing the effect of isentropic efficiencies and heat exchanger design and operation characteristics for different thermodynamic cycles.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

Real-time estimators for scattering observables: A full account of finite-volume errors for quantum simulation

The real-time correlators of quantum field theories can be directly probed through new approaches to simulation, such as quantum computing and tensor networks. This provides a new framework for computing scattering observables in lattice formulations of strongly interacting theories, such as lattice quantum chromodynamics. In this paper, we prove that the proposal of real-time estimators of scattering observables is universally applicable to all scattering observables of gapped quantum field theories. All finite-volume errors are exponentially suppressed, and the rate of this suppression is controlled by the regulator considered, namely, a displacement of the spectrum of the theory into the complex plane. A partial restoration of Lorentz symmetry by averaging over different boosts gives an additional suppression of finite volume errors. Our results also apply to the simulation of wave packet scattering, where a similar averaging is performed to construct the wave packets that regulate the finite volume effects. This result represents a necessary key step toward determining a broad class of scattering observables via quantum computing that are currently inaccessible via classical computing. Such observables are relevant for various applications, including hadron spectroscopy, hadron structure, and precision tests of the Standard Model. We also comment on potential applications of our results to traditional computational schemes.

Burbano, Ivan M. [University of California, Berkel↗
DOE OSTI2023

How thermal fluctuations influence the function of the FeMo cofactor in nitrogenase enzymes

The catalytic mechanism of N 2 fixation by nitrogenase remains unresolved in how the strong N≡N bond is activated and why the reductive elimination of H 2 is required. Here, we use density functional theory and physiologically relevant thermal simulations to elucidate the mechanism of the complete nitrogenase catalytic cycle. Over the accumulation of four reducing equivalents, we find that protons and electrons transfer to the FeMo cofactor to weaken and break its bridge Fe–S bond, leading to temporary H 2 S formation that exposes the Fe sites to weakly bind N 2 . Remarkably, we find that subsequent H 2 formation is responsible for chemical activation to an N=N double bond accompanied by a low barrier for H 2 release. We emphasize that finite temperature effects smooth out mechanistic differences between DFT functionals observed at 0 K, thus leading to a consistent understanding as to why H formation is an obligatory step in N 2 adsorption and activation.

DFT↗
DOE OSTI2025

Impact of representative ground motion level on seismic PSA with the boundary between overestimation and underestimation

One commonly used approach in seismic probabilistic safety assessment (PSA) is the discrete method. This method follows the standard PSA framework and can be applied to various models, such as multi-unit models, while reducing computational costs using standard software. However, due to the inability to subdivide intervals infinitely, the discrete method approximates with a finite number of subintervals. In practice, different numbers of subintervals are applied, and the representative ground motion level is selected based on expert judgment. When employing a smaller number of subintervals, it is important to take caution to prevent underestimation. This study analyzes the impact of the representative ground motion level on seismic risk. It confirms that underestimation can occur with a small number of subintervals depending on the representative ground motion level. This study also proposes a method for determining the boundary of underestimation and overestimation. The method is demonstrated through examples, providing a mathematical foundation for selecting appropriate representative ground motion levels. By avoiding underestimation, this research helps prevent the oversight of significant risk contributors and enhances the understanding of seismic risk.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Beam impedance and heating analysis of the diagnostic stripline

In this work the heat transfer analysis has been performed for the stripline kicker, design of which is used at many accelerator facilities as a part of the transverse bunch-by-bunch feed-back system, including NSLS-II. The experimental data of the electrode temperature have been collected from a special diagnostic stripline designed with six infrared view ports. The temperature analysis is based on the power loss calculations and the loss factor computed by the 3-D electromagnetic simulation code GdfidL. To estimate the power loss, we analyze the real part of the longitudinal impedance and the loss factor for two different regimes. We use the ANSYS finite element code to correlate the temperature with the power loss. Lambertson and Shafer formalisms of the longitudinal and transverse beam impedances are discussed for the simplified stripline geometries. The required characteristic impedances and the geometric factors are determined analytically for the simplified geometries and compared with the results for the designed stripline geometry using the 2-D POISSON code by solving the Laplace’s equation in two dimensions. Lambertson’s equations are compared with numerical results. We introduce the frequency range, within which we can apply the analytical approach.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach

Here we apply conformal mapping to find the evolving shapes of a dissolving cylinder in a potential flow. Similar equations can be used to describe melting in a flowing liquid phase. Results are compared with microfluidic experiments and numerical simulations. Shapes predicted by conformal mapping agree almost perfectly with experimental observations, after a modest (20 %) rescaling of the time. Finite-volume simulations show that the differences with experiment are connected to the underlying assumptions of the analytical model: potential flow and diffusion-limited dissolution. Approximate solutions of the equations describing the evolution of the shape of the undissolved solid can be derived from a Laurent expansion of the mapping function from the unit circle. Asymptotic expressions for the evolution of the area of the disk and the shift in its centre of mass have been derived at low and high Péclet number. Analytic approximations to the leading-order Laurent coefficients provide additional insight into the mechanisms underlying pore-scale dissolution.

58 GEOSCIENCES↗
DOE OSTI2020

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. Here, we also present results for the energy as a function of temperature in the continuum limit.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

$2+1$ dimensional Floquet systems and lattice fermions: Exact bulk spectral equivalence

A connection has recently been proposed between periodically driven systems known as Floquet insulators in continuous time and static fermion theories in discrete time. This connection has been established in a (1+1) ( 1 + 1 ) -dimensional free theory, where an explicit mapping between the spectra of a Floquet insulator and a discrete-time Dirac fermion theory has been formulated. Here we investigate the potential of static discrete-time theories to capture Floquet physics in higher dimensions, where so-called anomalous Floquet topological insulators can emerge that feature chiral edge states despite having bulk bands with zero Chern number. Starting from a particular model of an anomalous Floquet system, we provide an example of a static discrete-time theory whose bulk spectrum is an exact analytic match for the Floquet spectrum. The spectra with open boundary conditions in a particular strip geometry also match up to finite-size corrections. However, the models differ in several important respects. The discrete-time theory is spatially anisotropic, so that the spectra do not agree for all lattice terminations, e.g. other strip geometries or on half spaces. This difference can be attributed to the fact that the static discrete-time model is quasi-one-dimensional in nature and therefore has a different bulk-boundary correspondence than the Floquet model.

Iadecola, Thomas (ORCID:0000000251456441)↗
DOE OSTI2024

Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems

The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.

Alfven waves↗
DOE OSTI2021

Thermal and Mechanical Design of a High-Voltage Power Electronics Package: Preprint

Thermal and mechanical design aspects of a power electronic package with 5-kV GaN devices are presented. Through finite element analyses, the impact of different cooling configurations and device location on the thermal performance and reliability of the package was investigated. It was found that placing the devices closer to the direct-bond-copper substrate as opposed to a centered approach in the proposed double-sided-cooling configuration resulted in improved heat dissipation. This approach also reduced the total number of attachment layers, thereby likely improving the reliability of the package. Furthermore, simulations revealed that the device location had a negligible impact on the thermomechanical behavior of the attachment layers, as they are more prone to the local coefficient of thermal expansion mismatch.

27 ARPA - Advanced Research Projects Agency-Energy↗
DOE OSTI2020

Crushing behaviors of unidirectional carbon fiber reinforced plastic composites under dynamic bending and axial crushing loading

The dynamic crushing characteristics of unidirectional carbon fiber reinforced plastic composites under two loading types, dynamic three-point bending and axial crushing, are investigated by experiment and finite element simulation. Hat section samples with two different layup orientations are tested at various impact velocities to investigate its effects. The experiment results show that delamination plays a critical role in dynamic bending deformation. Different layup orientations lead to similar crushing bending behaviors, but remarkable variations in delamination pattern evolutions. In contrast, layup orientation exerts more significant effects on failure modes and energy absorption in axial crushing. A finite element model with multiple layers of thick shell element and cohesive element is established, and simulations are performed for both dynamic three-point bending and axial crushing. The model can capture the crushing behaviors in the two loading types reasonably well. The damage evolution and energy absorption through various mechanisms at different loading conditions are discussed.

42 ENGINEERING↗
DOE OSTI2020

Generalized hydrodynamic analysis of transport through a finite open nanopore for two-component single-file systems

Single-file diffusion (SFD) in finite open nanopores is characterized by nonzero spatially varying tracer diffusion coefficients within a generalized hydrodynamic description. This contrasts with infinite SFD systems where tracer diffusivity vanishes. In standard tracer counterpermeation (TCP) analysis, two reservoirs, each containing a different species, are connected to opposite ends of a finite pore. Here, we implement an extended TCP analysis to allow the two reservoirs to contain slightly different mixtures of the two species. Then, determination of diffusion fluxes through the pore allows extraction of diffusion coefficients for near-constant partial concentrations of the two species. This analysis is applied for a lattice-gas model describing two-component SFD through a finite linear pore represented by a one-dimensional array of cells. Two types of particles, A and B, can hop only to adjacent empty cells with generally different rates, $h_A$ and $h_B$. Particles are noninteracting other than exclusion of multiple cell occupancy. Results reveal generalized hydrodynamic tracer diffusion coefficients which adopt small values inversely proportional to pore length in the pore center, but which are strongly enhanced near pore openings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Computing the Central Charge of the 3D Ising CFT Using Quantum Finite Elements

The 3D Ising conformal field theory (CFT) describes different physical systems, such as uniaxial magnets or fluids, at their critical points. In absence of an analytical solution for the 3D Ising model, the scaling dimensions and operator product expansion (OPE) coefficients characterizing this CFT must be determined numerically. The currently most-cited values for these quantities have been obtained from the conformal bootstrap, while lattice calculations have so far only produced reliable results for the scaling dimensions involved in calculating the critical exponents. Using Quantum Finite Elements to investigate critical \(\phi^4\)-theory on \(\mathbb{R}\times\mathbb{S}^2\), we show in this work that it is possible to extract scaling dimensions and OPE coefficients of the 3D Ising CFT by fitting the lattice four-point function with expectations from the operator product expansion for the radially quantized CFT and extrapolating to the continuum limit. This way, we have for the first time been able to use Monte Carlo simulations to compute the central charge of the theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Framework for phase transitions between the Maxwell and Gibbs constructions at finite temperature

The characteristics of the hadron-to-quark first-order phase transition differ depending on whether charge neutrality is locally or globally fulfilled. In 𝛽-equilibrated matter, these two possibilities correspond to the Maxwell and Gibbs constructions. Recently, we presented a new framework in which a continuously varying parameter allows one to describe a first-order phase transition in intermediate scenarios to the two extremes of fully local and fully global charge neutrality. In this work, we extend the previous framework to finite temperatures and out-of-𝛽 equilibrium conditions, making it available for simulations of core-collapse supernovae and binary neutron star mergers. We investigate its impact on key thermodynamic quantities across a range of baryon densities, temperatures, and electron fractions. We find that when matter is not in 𝛽 equilibrium, the pressure in the mixed phase is not constant even for the case of fully local charge neutrality. Moreover, we compute the thermal index using three different approaches, demonstrating that the finite-temperature extension of an equation of state using a constant thermal index can be ill defined when applied to the mixed phase.

QCD phase transitions↗
DOE OSTI2025

Shear viscosity from perturbative quantum chromodynamics to the hadron resonance gas at finite baryon, strangeness, and electric charge densities

Through model-to-data comparisons from heavy-ion collisions, it has been shown that the quark gluon plasma has an extremely small shear viscosity at vanishing densities. At large baryon densities, significantly less is known about the nature of the shear viscosity from quantum chromodynamics (QCD). Within heavy-ion collisions, there are three conserved charges: baryon number (B), strangeness (S), and electric charge (Q). Here we calculate the shear viscosity in two limits using perturbative QCD (pQCD) and an excluded-volume hadron resonance gas at finite BSQ densities. We then develop a framework that interpolates between these two limits such that shear viscosity is possible to calculate across a wide range of finite BSQ densities. We find that the pQCD and hadron resonance gas calculations have different BSQ density dependencies such that a rather nontrivial shear viscosity appears at finite densities.

Phenomenology↗
DOE OSTI2024

Advanced multimaterial shape optimization methods as applied to advanced manufacturing of wind turbine generators

Abstract Currently, many utility‐scale wind turbine generator original equipment manufacturers are dependent on imported rare earth permanent magnets, which are susceptible to market risks from cost instability. To lower the production costs of these generators and stay competitive in the market, several small wind manufacturers are pursuing continuous improvements to both generator design and manufacturing. However, traditional design and manufacturing methods have yielded marginal improvements in wind power performance. This work presents novel methods to redesign a baseline 15‐kW wind turbine generator with reduced rare‐earth permanent magnets by leveraging cutting‐edge three‐dimensional (3D) printed polymer‐bonded permanent magnets and steel. Symmetric, asymmetric, and multimaterial‐magnet parametrization methods are introduced for shape optimization. We extend the symmetric and asymmetric methods to the back iron in the stator to further investigate the impact and opportunities for performance improvements with lesser active materials. We employ a design‐of‐experiments approach with parametric computer‐aided design for shape generation and evaluate different designs by magneto‐thermal modeling and finite‐element analysis. We use adaptive sampling technique to identify better performing designs with lesser magnet mass, higher efficiency, and lower cogging torque when compared with the baseline generator. Asymmetric pole designs resulted in a magnet mass in the range of 4.77–5.37 kg, which was 27%–35% lighter than the baseline generator, suggesting that a new design freedom exists that can be enabled by advanced manufacturing, such as 3D printing. Shaping the back iron in the stator resulted in material savings in electrical steel of up to 14.62 kg, which was 20% lighter than the baseline stator. We conducted a structural analysis to evaluate an optimized asymmetric rotor design from the point of view of mechanical integrity and air‐gap stiffness. The magnetically optimal shape profile was shown as having a positive impact on the radial stiffness, and an optimal solution was discovered to reduce the structural mass by nearly 30 kg, which was 29% lighter than the baseline.

17 WIND ENERGY↗
DOE OSTI2025

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗