Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Predictive analytics”

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,045 records · Page 58

Simulation of Taylor flow evaporation for bubble-pump applications

Single-pressure absorption systems incorporate bubble-pump generators (BPGs) for refrigerant separation and passive fluid circulation. In conventional spot-heated BPGs, heat is transferred over a small area, requiring high source temperatures. Distributed-heated BPGs receive thermal input over most of the component surface, enabling low temperature operation. In this investigation, a Volume-of-Fluid phase-change simulation formulation is developed and validated. Here this approach is applied to the evaporating Taylor flow pattern in distributed-heated BPGs. A 2-D axisymmetric simulation is performed, which yields detailed information about the developing heat transfer and two-phase flow phenomena. Results are used to assess predicted trends and sub-models from a 1-D segmented BPG model. Close agreement is obtained between segmented model and simulation results for bubble rise velocity (5–7% deviation), bubble and slug lengths, void fraction (3%), and hydrodynamic pressure drop (18%). Specifying average Taylor bubble lengths from the simulation as an input to the segmented model reduces hydrodynamic pressure drop deviation to 6%. Simulated flow-evaporation heat transfer coefficients are significantly higher than those predicted using analytic models from the literature. A new flow evaporation heat transfer correlation that accounts for developing slug flow effects is proposed, and yields close agreement with simulation results for heat transfer coefficient (AAD = 11%) and overall heat transfer rate (2%). Overall, this investigation provides validation for a distributed-heated BPG modeling approach, which can enable passive refrigeration for diverse applications.

42 ENGINEERING↗

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SPH modeling of biomass granular flow: Theoretical implementation and experimental validation

The commercialization of biomass-derived energy is impeded by flowability challenges arising from the feeding and handling of granular biomass materials in full-scale biorefineries. To overcome these obstacles, a robust and accurate model to simulate the flow of granular biomass is indispensable. However, conventional mesh-based numerical codes are limited by inherent mesh distortion in simulating large deformation that commonly occurs in granular biomass handling. Here, in this study, we propose a graphics processing unit (GPU)-accelerated meshless Smoothed Particle Hydrodynamics (SPH) code to model the flow of granular biomass materials. A modified void ratio-based mass conversation, a hybrid particle-to-particle/surface frictional boundary treatment, and a hypoplastic constitutive model are implemented. Four numerical examples, an elastic block sliding on inclined planes, sand column collapse, Angle of Repose, and axial compression tests for pine chips, were simulated using the developed SPH code. The results demonstrate good agreement between numerical predictions and analytical and experimental data for all four examples, validating the SPH code and increasing confidence that it can be applied to simulate more complex granular biomass handling processes, such as hopper feeding or auger conveyance.

09 BIOMASS FUELS↗

The influence of environmental microseismicity on detection and interpretation of small-magnitude events in a polar glacier setting

Glacial environments exhibit temporally variable microseismicity. To investigate how microseismicity influences event detection, we implement two noise-adaptive digital power detectors to process seismic data from Taylor Glacier, Antarctica. We add scaled icequake waveforms to the original data stream, run detectors on the hybrid data stream to estimate reliable detection magnitudes and compare analytical magnitudes predicted from an ice crack source model. We find that detection capability is influenced by environmental microseismicity for seismic events with source size comparable to thermal penetration depths. When event counts and minimum detectable event sizes change in the same direction (i.e. increase in event counts and minimum detectable event size), we interpret measured seismicity changes as ‘true’ seismicity changes rather than as changes in detection. Generally, one detector (two degree of freedom (2dof)) outperforms the other: it identifies more events, a more prominent summertime diurnal signal and maintains a higher detection capability. We conclude that real physical processes are responsible for the summertime diurnal inter-detector difference. One detector (3dof) identifies this process as environmental microseismicity; the other detector (2dof) identifies it as elevated waveform activity. Our analysis provides an example for minimizing detection biases and estimating source sizes when interpreting temporal seismicity patterns to better infer glacial seismogenic processes.

54 ENVIRONMENTAL SCIENCES↗

Effect of Macromolecular Architecture on the Morphology of Polystyrene–Polyisoprene Block Copolymers

The impact of block connectivity on the morphologies of four block copolymers of varying architecture containing polystyrene (PS) and polyisoprene (PI) has been studied. The volume fraction of PS and molecular weight are held constant while varying the architecture from a linear PS–PI diblock copolymer to three different miktoarm star architectures: PS2PI, PSPI2, and PS2PI2. Morphologies of the PS2PI and PSPI2 miktoarm stars are different from those observed for the linear copolymer and dependent on the connectivity of the copolymer blocks. The change in morphology with connectivity indicates that combining two chains at a junction point leads to chain crowding, where subsequent excluded volume effects drive the change in morphology for each sample. The PS2PI2 miktoarm star exhibits the same morphology as the linear diblock but with a reduction in the size of the domains. The extent of the decrease in domain size indicates that chain stretching impacts the formation of this morphology. Experimentally observed morphologies for different chain architectures are generally consistent with three-dimensional self-consistent-field theory simulations, taking into account conformational asymmetry and experimental uncertainty in the copolymer composition. Furthermore, these results generally agree with analytical theory predictions that account for architectural and conformational asymmetry.

Polymer Science↗

Observing the onset of pressure-driven K-shell delocalization

The gravitational pressure in many astrophysical objects exceeds one gigabar (one billion atmospheres), creating extreme conditions where the distance between nuclei approaches the size of the K shell. This close proximity modifies these tightly bound states and, above a certain pressure, drives them into a delocalized state. Both processes substantially affect the equation of state and radiation transport and, therefore, the structure and evolution of these objects. Still, our understanding of this transition is far from satisfactory and experimental data are sparse. Here, in this work, we report on experiments that create and diagnose matter at pressures exceeding three gigabars at the National Ignition Facility where 184 laser beams imploded a beryllium shell. Bright X-ray flashes enable precision radiography and X-ray Thomson scattering that reveal both the macroscopic conditions and the microscopic states. The data show clear signs of quantum-degenerate electrons in states reaching 30 times compression, and a temperature of around two million kelvins. At the most extreme conditions, we observe strongly reduced elastic scattering, which mainly originates from K-shell electrons. We attribute this reduction to the onset of delocalization of the remaining K-shell electron. With this interpretation, the ion charge inferred from the scattering data agrees well with ab initio simulations, but it is significantly higher than widely used analytical models predict.

79 ASTRONOMY AND ASTROPHYSICS↗

The effect of the H{sub 2} density on the electron energy distribution in H{sup −} ion sources

The electron energy distribution in H{sup −} ion sources is studied. By three-dimensional Monte Carlo simulations (3D-MCC) and cross sections of collisions, the crucial collisions in the JAEA 10A are chosen to develop theoretical models for the low electron temperature and the high electron temperature, respectively. The effects of the H{sub 2} density on the low electron temperature and the high electron temperature in H{sup −} ion sources are predicted with analytic expressions and verified by 3D-MCC. The results show that both the low and high electron temperatures become lower while n{sub H2} increases.

43 PARTICLE ACCELERATORS↗

Photon assisted current in molecular nanojunctions with novel types of contacts

We propose new approaches to coherent control of transport via molecular junctions, which bypasses several of the hurdles to experimental realization of optically manipulated nanoelectronics noted in the previous literature. The first method is based on the application of intrinsic semiconductor contacts and optical frequencies below the semiconductor bandgap. Our analytical theory predicts a new phenomenon, referred to as coherent destruction of induced tunnelling, which extends the phenomenon of coherent destruction of tunnelling frequently discussed in the previous literature. We also propose to use graphene electrodes as a platform for effective photon assisted tunneling through molecular conduction nanojunctions. We predict dramatic increasing currents evaluated at side-band energies ∼nħω (n is a whole number) related to the modification of graphene gapless spectrum under the action of external electromagnetic field of frequency ω. Our results illustrate the potential of semiconductor and graphene contacts in coherent control of photocurrent.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Temperature dependent low-field measurements of the magnetocaloric ΔT with sub-mK resolution in small volume and thin film samples

We present temperature dependent ΔT measurements of the magnetocaloric effect in a thin film sample of Gd, employing magnetomodulation and detection of thermal radiation. A bulk sample of the metamagnetic material LaFe{sub 11.05}Co{sub 0.91}Si{sub 1.04} shows a strong broadening of the ΔT peak for increasing field amplitudes between 4 and 45 mT. Bulk Gd in comparison shows only a weak broadening. All investigated samples exhibit a clear quadratic dependence of ΔT on the external field H{sub ext} at the ΔT peak maximum, contrary to earlier predictions. An analytic expression is derived that interpolates between the H{sub ext}{sup 2}-behavior at low and the well-known H{sub ext}{sup 2/3}-behavior at high fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Broadband resonances in indium-tin-oxide nanorod arrays

There is currently much discussion within the nanophotonics community regarding the origin of wavelength selective absorption/scattering of light by the resonances in nanorod arrays. Here, we report a study of resonances in ordered indium-tin-oxide nanorod arrays resulting from waveguide-like modes. We find that with only a 2.4% geometrical coverage, micron-length nanorod arrays interact strongly with light across a surprisingly wide band from the visible to the mid-infrared, resulting in less than 10% transmission. Simulations show excellent agreement with our experimental observations. The field profile in the vicinity of the rods obtained from simulations shows that the electric field is mainly localized on the surfaces of the nanorods for all resonances. Based on our analysis, the resonances in the visible are different in character from those in the infrared. When light is incident on the array, part of it propagates in the space between the rods and part of it is guided within the rods. The phase difference (interference) at the ends of the rods forms the basis for the resonances in the visible region. The resonances in the infrared are Fabry-Perot-like resonances involving standing surface waves between the opposing ends of the rods. Simple analytical formulae predict the spectral positions of these resonances. It is suggested that these phenomena can be utilized for wavelength-selective photodetectors, modulators, and nanorod-based solar cells.

14 SOLAR ENERGY↗

Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration

Turbulent Richtmyer–Meshkov instability (RMI) is investigated through a series of high resolution three-dimensional simulations of two initial conditions with eight independent codes. The simulations are initialised with a narrowband perturbation such that instability growth is due to non-linear coupling/backscatter from the energetic modes, thus generating the lowest expected growth rate from a pure RMI. By independently assessing the results from each algorithm and computing ensemble averages of multiple algorithms, the results allow a quantification of key flow properties as well as the uncertainty due to differing numerical approaches. A new analytical model predicting the initial layer growth for a multimode narrowband perturbation is presented, along with two models for the linear and non-linear regimes combined. Overall, the growth rate exponent is determined as θ=0.292±0.009, in good agreement with prior studies; however, the exponent is decaying slowly in time. Also, θ is shown to be relatively insensitive to the choice of mixing layer width measurements. The asymptotic integral molecular mixing measures Θ=0.792±0.014, Ξ=0.800±0.014, and Ψ=0.782±0.013 are lower than some experimental measurements but within the range of prior numerical studies. The flow field is shown to be persistently anisotropic for all algorithms, at the latest time having between 49% and 66% higher kinetic energy in the shock parallel direction compared to perpendicular and does not show any return to isotropy. The plane averaged volume fraction profiles at different time instants collapse reasonably well when scaled by the integral width, implying that the layer can be described by a single length scale and thus a single θ. In conclusion, quantitative data given for both ensemble averages and individual algorithms provide useful benchmark results for future research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implosion of auto-magnetizing helical liners on the Z facility

In the first auto-magnetizing liner implosion experiments on the Z Facility, precompressed internal axial fields near 150 T were measured and 7.2-keV radiography indicated a high level of cylindrical uniformity of the imploding liner's inner surface. An auto-magnetizing (AutoMag) liner is made of discrete metallic helical conductors encapsulated in insulating material. Here, the liner generates internal axial magnetic field as a 1–2 MA, 100–200 ns current prepulse flows through the helical conductors. After the prepulse, the fast-rising main current pulse causes the insulating material between the metallic helices to break down ceasing axial field production. After breakdown, the helical liner, nonuniform in both density and electrical conductivity, implodes in 100 ns. In-flight radiography data demonstrate that while the inner wall maintains cylindrical uniformity, multiple new helically oriented structures are self-generated within the outer liner material layers during the implosion; this was not predicted by simulations. Furthermore, liner stagnation was delayed compared to simulation predictions. An analytical implosion model is compared with experimental data and preshot simulations to explore how changes in the premagnetization field strength and drive current affect the liner implosion trajectory. Both the measurement of >100 T internal axial field production and the demonstration of cylindrical uniformity of the imploding liner's inner wall are encouraging for promoting the use of AutoMag liners in future MagLIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-stream instability of a single long electron bunch in a storage ring

A multistream instability is observed experimentally during the longitudinal expansion of an electron beam in a storage ring. The instability is observed when the beam expands in length to several times the circumference of the ring so that portions of the beam overlap. While portions of the beam overlap in physical space, they form multiple streams and remain separate in velocity space. The streams become unstable as their number increases and their separation in velocity decreases. An analytical theory predicts the onset of the instability, in agreement with simulations and measurements, over a wide range of peak line-charge densities (10.3 pC/m–1.8 nC/m) and bunch lengths. Finally, this work extends previous calculations to include the dynamic nonlinear elongation of the bunch, with a given initial length, and defines an onset criterion for the filament velocity separation for the beam to be unstable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analytic investigation of combined effects of anisotropic thermal transport and energetic particles on stability of resistive plasma resistive wall mode

The combined effects of anisotropic thermal transport and trapped energetic particles (EPs) on the stability of the resistive plasma resistive wall mode (RPRWM) are investigated by an energy-principle based analytical model. The results qualitatively confirm that of a recent toroidal modeling study [Bai et al., Phys. Plasmas 27, 072502 (2020)], in which the thermal transport can stabilize the RPPWM depending on the parameters of both the plasma equilibrium and energetic particles. The analytical model predicts a complete stabilization of the RPRWM in highly resistive plasmas, at sufficiently high EPs’ pressure, and finite plasma flow. The stabilizing effect of thermal transport originates from its enhancement of energy dissipations associated with both the resistive layer and the trapped energetic particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

End-exclusion zones in strongly stretched, molten polymer brushes of arbitrary shape

Theories of strongly stretched polymer brushes, particularly the parabolic brush theory, are valuable for providing analytically tractable predictions for the thermodynamic behavior of surface-grafted polymers in a wide range of settings. Furthermore, the parabolic brush limit fails to describe polymers grafted to convex curved substrates, such as the surfaces of spherical nanoparticles or the interfaces of strongly segregated block copolymers. It has previously been shown that strongly stretched curved brushes require a boundary layer devoid of free chain ends, requiring modifications of the theoretical analysis. While this “end-exclusion zone” has been successfully incorporated into the descriptions of brushes grafted onto the outer surfaces of cylinders and spheres, the behavior of brushes on surfaces of arbitrary curvature has not yet been studied. We present a formulation of the strong-stretching theory for molten brushes on the surfaces of arbitrary curvature and identify four distinct regimes of interest for which brushes are predicted to possess end-exclusion zones, notably including regimes of positive mean curvature but negative Gaussian curvature. Through numerical solutions of the strong-stretching brush equations, we report predicted scaling of the size of the end-exclusion zone, the chain end distribution, the chain polarization, and the free energy of stretching with mean and Gaussian surface curvatures. Through these results, we present a comprehensive picture of how the brush geometry influences the end-exclusion zones and exact strong-stretching free energies, which can be applied, for example, to model the full spectrum of brush geometries encountered in block copolymer melt assembly.

36 MATERIALS SCIENCE↗

Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Materials science↗

BINARY-DISK INTERACTION: GAP-OPENING CRITERIA

We study the interaction of an equal mass binary with an isothermal circumbinary disk motivated by the evidence of the formation of massive black hole binaries surrounded by gas, after a major merger of gas-rich galaxies. We focus on the torques that the binary produces on the disk and how the exchange of angular momentum can drive the gap formation on it. We propose that the angular momentum exchange between the binary and the disk is through the gravitational interaction of the binary and a (tidally formed) global non-axisymmetric perturbation in the disk. Using this interaction, we derive an analytic criterion for the gap formation in the disk that can be expressed either via the characteristic velocities of the binary-disk system or in terms of the structural parameters h/a and M(< r)/M{sub bin}. Using numerical simulations we show that the simulations where the binary opens a gap in the disk and the simulations where the disk does not have a gap are distributed into two well separated regions. Our analytic criterion predicts a shape of the threshold between these two regions that is consistent with our simulations and the other ones in the literature. We propose an analogy between the regime without (with) a gap in the disk and the Type I (Type II) migration that is observed in simulations of planet-disk interaction (extreme mass ratio binary), emphasizing that the interaction that drives the formation of a gap on the disk is different in the regime that we analyze (comparable mass binary).

79 ASTRONOMY AND ASTROPHYSICS↗