Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 1,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Radiation forces and torque on a rigid elliptical cylinder in acoustical plane progressive and (quasi)standing waves with arbitrary incidence

This paper presents two key contributions; the first concerns the development of analytical expressions for the axial and transverse acoustic radiation forces exerted on a 2D rigid elliptical cylinder placed in the field of plane progressive, quasi-standing, or standing waves with arbitrary incidence. The second emphasis is on the acoustic radiation torque per length. The rigid elliptical cylinder case is important to be considered as a first-order approximation of the behavior of a cylindrical fluid column trapped in air because of the significant acoustic impedance mismatch at the particle boundary. Based on the rigorous partial-wave series expansion method in cylindrical coordinates, non-dimensional acoustic radiation force and torque functions are derived and defined in terms of the scattering coefficients of the elliptic cylinder. A coupled system of linear equations is obtained after applying the Neumann boundary condition for an immovable surface in a non-viscous fluid and solved numerically by matrix inversion after performing a single numerical integration procedure. Computational results for the non-dimensional force components and torque, showing the transition from the progressive to the (equi-amplitude) standing wave behavior, are performed with particular emphasis on the aspect ratio a/b, where a and b are the semi-axes of the ellipse, the dimensionless size parameter, as well as the angle of incidence ranging from end-on to broadside incidence. The results show that the elliptical geometry has a direct influence on the radiation force and torque, so that the standard theory for circular cylinders (at normal incidence) leads to significant miscalculations when the cylinder cross section becomes non-circular. Moreover, the elliptical cylinder experiences, in addition to the acoustic radiation force, a radiation torque that vanishes for the circular cylinder case. The application of the formalism presented here may be extended to other 2D surfaces of arbitrary shape, such as Chebyshev cylindrical particles with a small deformation, stadiums (with oval shape), or other non-circular geometries.

42 ENGINEERING↗
DOE OSTI2016

Thermodynamic efficiency limits of classical and bifacial multi-junction tandem solar cells: An analytical approach

Bifacial tandem cells promise to reduce three fundamental losses (i.e., above-bandgap, below bandgap, and the uncollected light between panels) inherent in classical single junction photovoltaic (PV) systems. The successive filtering of light through the bandgap cascade and the requirement of current continuity make optimization of tandem cells difficult and accessible only to numerical solution through computer modeling. The challenge is even more complicated for bifacial design. In this paper, we use an elegantly simple analytical approach to show that the essential physics of optimization is intuitively obvious, and deeply insightful results can be obtained with a few lines of algebra. This powerful approach reproduces, as special cases, all of the known results of conventional and bifacial tandem cells and highlights the asymptotic efficiency gain of these technologies.

Alam, Muhammad Ashraful (ORCID:0000000187756043)↗
DOE OSTI2017

On- and off-axis spectral emission features from laser-produced gas breakdown plasmas

Laser-heated gas breakdown plasmas or sparks emit profoundly in the ultraviolet and visible region of the electromagnetic spectrum with contributions from ionic, atomic, and molecular species. Laser created kernels expand into a cold ambient with high velocities during its early lifetime followed by confinement of the plasma kernel and eventually collapse. However, the plasma kernels produced during laser breakdown of gases are also capable of exciting and ionizing the surrounding ambient medium. Two mechanisms can be responsible for excitation and ionization of surrounding ambient: viz. photoexcitation and ionization by intense ultraviolet emission from the sparks produced during the early times of its creation and/or heating by strong shocks generated by the kernel during its expansion into the ambient. In this study, an investigation is made on the spectral features of on- and off-axis emission features of laser-induced plasma breakdown kernels generated in atmospheric pressure conditions with an aim to elucidate the mechanisms leading to ambient excitation and emission. Pulses from an Nd:YAG laser emitting at 1064 nm with 6 ns pulse duration are used to generate plasma kernels. Laser sparks were generated in air, argon, and helium gases to provide different physical properties of expansion dynamics and plasma chemistry considering the differences in laser absorption properties, mass density and speciation. Point shadowgraphy and time-resolved imaging were used to evaluate the shock wave and spark self-emission morphology at early and late times while space and time resolved spectroscopy is used for evaluating the emission features as well as for inferring plasma fundaments at on- and off-axis. Structure and dynamics of the plasma kernel obtained using imaging techniques are also compared to numerical simulations using computational fluid dynamics code. The emission from the kernel showed that spectral features from ions, atoms and molecules are separated in time with an early time temperatures and densities in excess of 35000 K and 4×10 18 /cm 3 with an existence of thermal equilibrium. However, the emission from the off-kernel positions from the breakdown plasmas showed enhanced ultraviolet radiation with the presence of N 2 bands and represented by non-LTE conditions. Finally, our results also highlight that the ultraviolet radiation emitted during early time of spark evolution is the predominant source of the photo-excitation of the surrounding medium.

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

Parameter-Space Survey of Linear G-mode and Interchange in Extended Magnetohydrodynamics

The extended magnetohydrodynamic stability of interchange modes is studied in two configurations. In slab geometry, a local dispersion relation for the gravitational interchange mode (g-mode) with three different extensions of the MHD model [P. Zhu, et al., Phys. Rev. Lett. 101, 085005 (2008)] is analyzed. Our results delineate where drifts stablize the g-mode with gyroviscosity alone and with a two-fluid Ohm’s law alone. Including the two-fluid Ohm’s law produces an ion drift wave that interacts with the g-mode. This interaction then gives rise to a second instability at finite k y . A second instability is also observed in numerical extended MHD computations of linear interchange in cylindrical screw-pinch equilibria, the second configuration. Particularly with incomplete models, this mode limits the regions of stability for physically realistic conditions. But, applying a consistent two-temperature extended MHD model that includes the diamagnetic heat flux density ($\vec{q}$ * ) makes the onset of the second mode occur at larger Hall parameter. For conditions relevant to the SSPX experiment [E.B. Hooper, Plasma Phys. Controlled Fusion 54, 113001 (2012)], significant stabilization is observed for Suydam parameters as large as unity (D s ≲1).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The effect of entrance flow development on vortex formation and wall shear stress in a curved artery model

In this work, we numerically investigate the effect of entrance condition on the spatial and temporal evolution of multiple three-dimensional vortex pairs and the wall shear stress distribution in a curved artery model. We perform this study using a Newtonian blood-analog fluid subjected to a pulsatile flow with two inflow conditions. The first flow condition is fully developed while the second condition is undeveloped (i.e., uniform). We discuss the connection along the axial direction between regions of organized vorticity observed at various cross sections of the model and compare results between the different entrance conditions. We model a human artery with a simple, rigid 180° curved pipe with a circular cross section and constant curvature, neglecting the effects of taper, torsion, and elasticity. Numerical results are computed from a discontinuous high-order spectral element flow solver. The flow rate used in this study is physiological. We observe differences in secondary flow patterns, especially during the deceleration phase of the physiological waveform where multiple vortical structures of both Dean-type and Lyne-type coexist. The results indicate that decreased axial velocities under an undeveloped condition produce smaller secondary flows that ultimately inhibit growth of any interior flow vortices. We highlight the effect of the entrance condition on the formation of these structures and subsequent appearance of abnormal inner wall shear stresses, which suggest there may be a lower prevalence of cardiovascular disease in curved arteries where the flow is rather undeveloped—a potentially physiologically significant result to help understand the influence of blood flow development on disease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

THE ABUNDANCES OF LIGHT NEUTRON-CAPTURE ELEMENTS IN PLANETARY NEBULAE. III. THE IMPACT OF NEW ATOMIC DATA ON NEBULAR SELENIUM AND KRYPTON ABUNDANCE DETERMINATIONS

The detection of neutron(n)-capture elements in several planetary nebulae (PNe) has provided a new means of investigating s-process nucleosynthesis in low-mass stars. However, a lack of atomic data has inhibited accurate trans-iron element abundance determinations in astrophysical nebulae. Recently, photoionization (PI) and recombination data were determined for Se and Kr, the two most widely detected n-capture elements in nebular spectra. We have incorporated these new data into the photoionization code Cloudy. To test the atomic data, numerical models were computed for 15 PNe that exhibit emission lines from multiple Kr ions. We found systematic discrepancies between the predicted and observed emission lines that are most likely caused by inaccurate PI and recombination data. These discrepancies were removed by adjusting the Kr{sup +}–Kr{sup 3+} PI cross sections within their cited uncertainties and the dielectronic recombination rate coefficients by slightly larger amounts. From grids of models spanning the physical conditions encountered in PNe, we derive new, broadly applicable ionization correction factor (ICF) formulae for calculating Se and Kr elemental abundances. The ICFs were applied to our previous survey of near-infrared [Kr iii] and [Se iv] emission lines in 120 PNe. The revised Se and Kr abundances are 0.1–0.3 dex lower than former estimates, with average values of [Se/(O, Ar)] = 0.12 ± 0.27 and [Kr/(O, Ar)] = 0.82 ± 0.29, but correlations previously found between their abundances and other nebular and stellar properties are unaffected. We also find a tendency for high-velocity PNe that can be associated with the Galactic thick disk to exhibit larger s-process enrichments than low-velocity PNe belonging to the thin-disk population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Challenges in nuclear structure theory

The goal of nuclear structure theory is to build a comprehensive microscopic framework in which properties of nuclei and extended nuclear matter, and nuclear reactions and decays can all be consistently described. Due to novel theoretical concepts, breakthroughs in the experimentation with rare isotopes, increased exchange of ideas across different research areas, and the progress in computer technologies and numerical algorithms, nuclear theorists have been quite successful in solving various bits and pieces of the nuclear many-body puzzle and the prospects are exciting. This article contains a brief, personal perspective on the status of the field.

Physics↗
DOE OSTI2025

Including the vacuum energy in stellarator coil design

Being three-dimensional, stellarators have the advantage that plasma currents are not essential for creating rotational-transform; however, the external current-carrying coils in stellarators can have strong geometrical shaping, which can complicate the construction. Reducing the inter-coil electromagnetic forces acting on strongly shaped 3D coils and the stress on the support structure while preserving the favorable properties of the magnetic field is a design challenge. In this work, we recognize that the inter-coil forces are the gradient of the vacuum magnetic energy. We introduce an objective functional built on the usual quadratic flux on a prescribed target surface together with a weighed penalty on the vacuum energy. The Euler–Lagrange equation for stationary states is derived, and numerical illustrations are computed using a modern stellarator optimization framework. A study of the effect of the energy functional on the inter-coil forces is conducted and the energy is shown to be a promising quantity in producing coils with low forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Zinc-blende to rocksalt transition in SiC in a laser-heated diamond-anvil cell

We explore the stability of the ambient pressure zinc-blende polymorph (B3) structure of silicon carbide (SiC) at high pressures and temperatures where it transforms to the rocksalt (B1) structure. We find that the transition occurs ~40 GPa lower than previously measured when heated to moderately high temperatures. A lower transition pressure is consistent with the transition pressures predicted in numerous ab initio computations. We find a large volume decrease across the transition of ~17%, with the volume drop increasing at higher formation pressures, suggesting this transition is volume driven yielding a nearly pressure-independent Clapeyron slope. Such a dramatic density increase occurring at pressure is important to consider in applications where SiC is exposed to extreme conditions, such as in industrial applications or planetary interiors.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Embedding hard physical constraints in neural network coarse-graining of three-dimensional turbulence

In recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of partial differential equations is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into convolutional neural networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. Here, we demonstrate results on three-dimensional fully developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Universal Tripartite Entanglement in One-Dimensional Many-Body Systems

Motivated by conjectures in holography relating the entanglement of purification and reflected entropy to the entanglement wedge cross section, we introduce two related non-negative measures of tripartite entanglement g and h. We prove structure theorems which show that states with nonzero g or h have nontrivial tripartite entanglement. We then establish that in one dimension these tripartite entanglement measures are universal quantities that depend only on the emergent low-energy theory. For a gapped system, we argue that either g≠0 and h=0 or g=h=0, depending on whether the ground state has long-range order. For a critical system, we develop a numerical algorithm for computing g and h from a lattice model. We compute g and h for various CFTs and show that h depends only on the central charge whereas g depends on the whole operator content.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Geometrically necessary dislocation fingerprints of dislocation loop absorption at grain boundaries

Here we present a numerical methodology to compute the Nye-tensor fingerprints of dislocation loop absorption at grain boundaries (GBs) for comparison with TEM observations of irradiated polycrystals. Our approach links atomistic simulations of self-interstitial atom (SIA) prismatic loops gliding toward and interacting with GBs in body-centered cubic iron with experimentally extracted geometrically necessary dislocation (GND) maps to facilitate the interpretation of damage processes. The Nye-tensor analysis is strongly mesh-size dependent—corresponding to resolution-dependent TEM observations. The method computes GND fingerprints from discretized dislocation line segments extracted from molecular dynamics simulations of dislocation loops being absorbed at a GB. Specifically, we perform MD simulation of prismatic loops of two diameters and monitor the three stages of the absorption process: loop glide, the partial, and full absorption of the loops at a [1 0 0] symmetric tilt GB. These methods provide a framework for future investigations of the nature of defect absorption by grain boundaries under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Inferring effective electrostatic interaction of charge-stabilized colloids from scattering using deep learning

In this article, an innovative strategy is presented that incorporates deep auto-encoder networks into a least-squares fitting framework to address the potential inversion problem in small-angle scattering. To evaluate the performance of the proposed approach, a detailed case study focusing on charged colloidal suspensions was carried out. The results clearly indicate that a deep learning solution offers a reliable and quantitative method for studying molecular interactions. The approach surpasses existing deterministic approaches with respect to both numerical accuracy and computational efficiency. Overall, this work demonstrates the potential of deep learning techniques in tackling complex problems in soft-matter structures and beyond.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Development of a Converter Based Microgrid Test Platform

To overcome the limitations of digital simulation (numerical oscillation, limited computing capability of processors, etc.), a converter-based hardware test-bed was developed at CURENT for real-time power grid emulation. However, distribution systems especially microgrids cannot be emulated and tested on the existing hardware test-bed. This paper develops a microgrid test-bed based on the existing hardware test-bed to enable controller testing for microgrids with dynamic boundary. The design and realization of the microgrid hardware test-bed are introduced. The experimental results of the microgrid controller tests are also provided.

Li, Dingrui↗
DOE OSTI2015

Compiler-Directed Transformation for Higher-Order Stencils

As the cost of data movement increasingly dominates performance, developers of finite-volume and finite-difference solutions for partial differential equations (PDEs) are exploring novel higher-order stencils that increase numerical accuracy and computational intensity. This paper describes a new compiler reordering transformation applied to stencil operators that performs partial sums in buffers, and reuses the partial sums in computing multiple results. This optimization has multiple effect son improving stencil performance that are particularly important to higher-order stencils: exploits data reuse, reduces floating-point operations, and exposes efficient SIMD parallelism to backend compilers. We study the benefit of this optimization in the context of Geometric Multigrid (GMG), a widely used method to solvePDEs, using four different Jacobi smoothers built from 7-, 13-, 27-and 125-point stencils. We quantify performance, speedup, andnumerical accuracy, and use the Roofline model to qualify our results. Ultimately, we obtain over 4× speedup on the smoothers themselves and up to a 3× speedup on the multigrid solver. Finally, we demonstrate that high-order multigrid solvers have the potential of reducing total data movement and energy by several orders of magnitude.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Role of Azimuthal Prestress in the Longitudinal Degradation of Nb 3 Sn Superconducting Magnets

Superconducting magnet coils are subject to enormous electro-magnetic forces which push the cables away from the winding pole, and against the surrounding structure. This structure is usually optimized trying to limit the overall motion and the strains experienced by the superconducting elements. To achieve this, preload forces are applied both in the coil cross-section and along its length. If the e.m. forces overcome these preloads, separation between the coil and the pole occurs, resulting in an overall loss of rigidity. During the magnet design process, it is often tempting to treat the optimization of the azimuthal and longitudinal preload systems separately. However, the two are inextricably related: as the cross-section preload increases, friction can prevent any motion in the longitudinal plane, and decreasing it can instead allow dangerous motions in the conductor ends. The latter can result in very high strains and, in Nb 3 Sn conductors, damage that can prevent the magnet to reach the desired performances. In an attempt to define design guidelines, in this paper we use simplified numerical models to compute, as a function of the in-plane prestress, the variation of the peak strains in the end region of the coils. Finally, we investigate the impact of the azimuthal prestress on a real magnet case, the High-Luminosity Nb 3 Sn Quadrupole MQXF.

Nb3Sn↗
DOE OSTI2018

Duality-Free Decomposition Based Data-Driven Stochastic Security-Constrained Unit Commitment

To incorporate the superiority of both stochastic and robust approaches, a data-driven stochastic optimization is employed to solve the security-constrained unit commitment model. This approach makes the most use of the historical data to generate a set of possible probability distributions for wind power outputs and then it optimizes the unit commitment under the worst-case probability distribution. However, this model suffers from huge computational burden, as a large number of scenarios are considered. To tackle this issue, a duality-free decomposition method is proposed in this paper. This approach does not require doing duality, which can save a large set of dual variables and constraints, and therefore reduces the computational burden. In addition, the inner max-min problem has a special mathematical structure, where the scenarios have the similar constraint. Thus, the max-min problem can be decomposed into independent sub-problems to be solved in parallel, which further improves the computational efficiency. A numerical study on an IEEE 118-bus system with practical data of a wind power system has demonstrated the effectiveness of the proposal.

42 ENGINEERING↗