Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equations of motion”

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 415 records · Page 23

Phonon-Induced Localization of Excitons in Molecular Crystals from First Principles

The spatial extent of excitons in molecular systems underpins their photophysics and utility for optoelectronic applications. Phonons are reported to lead to both exciton localization and delocalization. However, a microscopic understanding of phonon-induced (de)localization is lacking, in particular, how localized states form, the role of specific vibrations, and the relative importance of quantum and thermal nuclear fluctuations. Here, we present a first-principles study of these phenomena in solid pentacene, a prototypical molecular crystal, capturing the formation of bound excitons, exciton-phonon coupling to all orders, and phonon anharmonicity, using density functional theory, the ab initio GW-Bethe-Salpeter equation approach, finite-difference, and path integral techniques. In this study, we find that for pentacene zero-point nuclear motion causes uniformly strong localization, with thermal motion providing additional localization only for Wannier-Mott-like excitons. Anharmonic effects drive temperature-dependent localization, and, while such effects prevent the emergence of highly delocalized excitons, we explore the conditions under which these might be realized

74 ATOMIC AND MOLECULAR PHYSICS↗

Exact and locally implicit source term solvers for multifluid-Maxwell systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3 S + 3) × (3 S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. In this paper, we present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗

Interplay between Ultrafast Electronic and Librational Dynamics in Liquid Nitrobenzene Probed with Two-Color Four-Wave Mixing

Here, we present an experimental and theoretical study of the interplay between ultrafast electron dynamics and librational dynamics in liquid nitrobenzene. A femtosecond ultraviolet pulse and two femtosecond near-infrared pulses interact with nitrobenzene molecules, generating a four-wave mixing nonlinear signal measured in the Optical Kerr Effect geometry. The signal is measured to be nonzero only at negative time delays, corresponding to the near-infrared pulses arriving before the ultraviolet pulse. We perform time-dependent Quantum Master Equation calculations with classical libration to simulate the experiment. The simulations support the conclusion that the near-infrared pulses launch librational motion while creating electronic coherences resulting in a libration-modulated electronic nonlinear response. The analysis of the phase-matched four-wave mixing signals suggests a nonparametric process leaving the molecules in an excited electronic state, providing new insight into ultrafast nonlinear optical interactions in liquids and advancing toward probing ultrafast electronic coherences in complex molecular liquids.

Shivaram, Niranjan [Lawrence Berkeley National Lab↗

Electronic, Ionic, and Mixed Conduction in Polymeric Systems

Polymers that simultaneously transport electrons and ions are paramount to drive the technological advances necessary for next-generation electrochemical devices, including energy storage devices and bioelectronics. However, efforts to describe the motion of ions or electrons separately within polymeric systems become inaccurate when both species are present. In this paper, we highlight the basic transport equations necessary to rationalize mixed transport and the multiscale materials properties that influence their transport coefficients. Potential figures of merit that enable a suitable performance benchmark in mixed conducting systems independent of end application are discussed. Practical design and implementation of mixed conducting polymers require an understanding of the evolving nature of structure and transport with ionic and electronic carrier density to capture the dynamic disorder inherent in polymeric materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Final-stage densification kinetics of direct current–sinteredZrB2

Final-stage sintering was analyzed for nominally phase pure zirconium diboride synthesized by borothermal reduction of high-purity ZrO2. Analysis was conducted on ZrB2 ceramics with relative densities greater than 90% using the Nabarro–Herring stress–directed vacancy diffusion model. Temperatures of 1900°C or above and an applied uniaxial pressure of 50 MPa were required to fully densify ZrB2 ceramics by direct current sintering. Ram travel data were collected and used to determine the relative density of the specimens during sintering. Specimens sintered between 1900 and 2100°C achieved relative densities greater than 97%, whereas specimens sintered below 1900°C failed to reach the final stage of sintering. The average grain size ranged from 1.0 to 14.7 μm. The activation energy was calculated from the slope of an Arrhenius plot that used the Kalish equation. The activation energy was 162 ± 34 kJ/mol, which is consistent with the activation energy for dislocation movement in ZrB2. The diffusion coefficients for dislocation motion that controls densification were 5.1 × 10−6 cm2/s at 1900°C and 5.1 × 10−5 cm2/s at 2100°C, as calculated from activation energy and average grain sizes. This study provides evidence that the dominant mechanism for final-stage sintering of ZrB2 ceramics is dislocation motion.

36 MATERIALS SCIENCE↗

Validation of Local Structural Loads Computed by OpenFAST Against Measurements From the FOCAL Experimental Campaign

This work presents the validation of the local structural load modeling capability in OpenFAST for floating substructures based on data from the FOCAL experimental campaign. Previously, OpenFAST could only represent the floating substructure as a rigid body, and though this approach can model the global response of the floater in most cases, it is not able to capture the structural loads within the floater's individual members. Consideration of local substructure loads is important for some floating designs, because the pursuit of cost reduction often results in lighter and more flexible structures. To address this limitation, the HydroDyn (hydrodynamics) and SubDyn (substructure dynamics) modules of OpenFAST have been recently extended to account for the flexibility of floating substructures. To validate this new capability, we compare the results obtained by OpenFAST with data measured during the FOCAL experimental campaign, which analyzed a 1:70 scale performance-matched model of the IEA 15-MW reference turbine atop a modified University of Maine VolturnUS-S semisubmersible in a wave basin under the action of both wind and waves. For the purposes of the present work, the most important feature of the experiment is the presence of load cells at the root of each pontoon, and our objective is to assess how well those loads are reproduced by OpenFAST. To model the distributed hydrodynamic and hydrostatic loads along the floating substructure, we adopt a strip-theory approach based on the Morison equation, and we discuss the impact of different hydrodynamic modeling options (wave stretching, MacCamy-Fuchs correction, and second-order wave kinematics) on both motions and loads. For simplicity, we focus on wave-only conditions, both regular and irregular. The results demonstrate good overall agreement for the loads at the root of the pontoons for the waves analyzed in this work, especially given the assumptions and simplifications inherent to a simple strip-theory model.

floating offshore wind turbine↗

Validation of Local Structural Loads Computed by OpenFAST Against Measurements From the Focal Experimental Campaign: Preprint

This work presents the validation of the local structural load modeling capability in OpenFAST for floating substructures based on data from the FOCAL experimental campaign. Previously, OpenFAST could only represent the floating substructure as a rigid body, and though this approach can model the global response of the floater in most cases, it is not able to capture the structural loads within its individual members. Consideration of local substructure loads is important for some floating designs, as the pursuit of cost reduction often results in lighter and more flexible structures. To address this limitation, the HydroDyn (hydrodynamics) and SubDyn (substructure dynamics) modules of OpenFAST have been recently extended to account for the flexibility of floating substructures. To validate this new capability, we compare the results obtained by OpenFAST with data measured during the FOCAL experimental campaign, which analyzed a 1:70 scale performance-matched model of the IEA 15-MW reference turbine atop a modified University of Maine VolturnUS-S semisubmersible in a wave basin under the action of both wind and waves. For the purposes of the present work, the most important feature of the experiment is the presence of load cells at the root of each pontoon, and our objective is to assess how well those loads are reproduced by OpenFAST. To model the distributed hydrodynamic and hydrostatic loads along the floating substructure, we adopt a strip-theory approach based on the Morison equation, and we discuss the impact of different hydrodynamic modeling options (wave stretching, MacCamy-Fuchs correction, and second-order wave kinematics) on both motions and loads. For simplicity, we focus on wave-only conditions, both regular and irregular. The results demonstrate good overall agreement for the loads at the root of the pontoons for the waves analyzed in this work, especially given the assumptions and simplifications inherent to a simple strip-theory model.

floating offshore wind turbine↗

Inviscid damping of an elliptical vortex subject to an external strain flow

Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. Here, the experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using E x B strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Some mathematics for quasi-symmetry

Quasi-symmetry of a steady magnetic field means integrability of first-order guiding-center motion by a spatial symmetry. Here, we derive many restrictions on the possibilities for a quasi-symmetry. We also derive an analog of the Grad–Shafranov equation for the flux function in a quasi-symmetric magnetohydrostatic field.

97 MATHEMATICS AND COMPUTING↗

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control

Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.

42 ENGINEERING↗

N-particle irreducible actions for stochastic fluids

We construct one- and two-particle irreducible (1PI and 2PI) effective actions for the stochastic fluid dynamics of a conserved density undergoing diffusive motion. We compute the 1PI action in one-loop order and the 2PI action in two-loop approximation. We derive a set of Schwinger-Dyson equations and regularize the resulting equations using Pauli-Villars fields. We numerically solve the Schwinger-Dyson equations for a non-critical fluid. We find that higher-loop effects summed by the Schwinger-Dyson renormalize the non-linear coupling. We also find indications of a diffuson-cascade, the appearance of n-loop correction with smaller and smaller exponential suppression.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pronghorn Porous Media Model Validation with Pressure Drop Measurements

The verification and validation (V&V) of Pronghorn is imperative to assert its accuracy when predicting the fluid velocity, temperature, and pressure in high temperature gas-cooled reactors. Pronghorn is a coarse-mesh, intermediate-fidelity, and multidimensional thermal-hydraulics (TH) code developed by the Idaho National Laboratory (INL). New pebble bed experiments are used to observe the details of the fluid motion and pressure drop in the porous bed under the reactor normal operation. This paper focuses on the validation of the Pronghorn compressible and incompressible Navier-Stokes equations using the pressure drop measurements performed at the engineering-scale pebble bed facility at the Texas A&M university (TAMU). Various pressure drop correlations and porosity functions are implemented in both Pronghorn and STAR-CCM+ to compare the pressure drop due to the combined viscous and inertial resistances in the porous bed. The correlations accounting for the near-wall effect are also utilized to observe if the pressure drop estimates can be improved. Pronghorn porous media models predict the pressure drop well relative to the STAR-CCM+ simulation results and 1D correlations, and both the finite element method (FEM) and finite volume method (FVM) perform accurately. Pronghorn models are also validated with the experimental measurements given the different Reynolds number ranges and specific aspect ratios. The likelihood of the statistical significance between the pressure drop measurements and specific correlations or simulations is low provided that the overlap of their confidence intervals is more than the half of a single arm. Several validation metrics are reasonable in regard to the similar studies from other literature. The precise average pebble bed porosity estimation has much impact on the pressure drop, and the Foumeny and Montillet (dense packing) models carry out the accurate pressure drop prediction by considering the near-wall effect.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fluctuation dynamo in a weakly collisional plasma

The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle–particle collisions allows departures from local thermodynamic equilibrium. These departures – pressure anisotropies – exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present an extensive numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart, particularly when the magnetic field is dynamically weak. If instead the parallel viscous stress is left unabated – a situation relevant to recent kinetic simulations of the fluctuation dynamo and, we argue, to the early stages of the dynamo in a magnetized ICM – the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics reminiscent of the saturated state of the large-Prandtl-number ( ${Pm}\gtrsim {1}$ ) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches simulated dynamo growth rates and magnetic-energy spectra. A prediction of this model, confirmed by our numerical simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. Finally, in all cases that do exhibit a viable dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Staircase formation by resonant and non-resonant transport of potential vorticity

The E × B staircase is a quasi-periodic pattern of pressure profile corrugations. In this work, we present a new mechanism for E × B staircase formation that involves resonant transport versus non-resonant transport. We start from a potential vorticity evolution system and use quasi-linear theory, a model dispersion relation, and a bi-Lorentzian spectrum approximation, to construct the relation between the fluxes and the profiles. With these fluxes, we close the profile evolution equations and the extended turbulence intensity evolution equation, which together constitute a turbulence-profile evolution system. In this system, the Doppler effect from the E × B mean flow can cause resonance between trapped ion precession motion and the trapped ion mode, which drives a resonant transport contribution to the fluxes. The profiles will be flattened where the resonant transport is switched on. In contrast, for the regions of non-resonant transport, profiles are steeper. A quasi-periodic pattern of profile corrugation (the E × B staircase) spontaneously emerges in this system, which is the two states mentioned above, arranged as alternating layers in space. The feedback processes during the staircase pattern formation are identified. An estimate of the critical value of the boundary heat flux is obtained, above which the staircase formation will be triggered. An estimate scaling of the step size in the staircase pattern is obtained. The resonant turbulent transport is also a mechanism for collisionless saturation of zonal flow. Furthermore, this work is related to internal transport barrier formation and suggests some new scenarios, such as an enhanced confined L mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Diagnosing the DARHT Electron Beam-Target Interaction and Hydrodynamic Expansion [Dissertation]

Quantitative electron beam-target interaction studies on the Dual-Axis Radiographic Hydrodynamic Test (DARHT) electron linear induction accelerators have only been performed recently. This includes characterization of the temperature, density, pressure, extent, and expansion velocity of the plasma plume. The results presented in this dissertation present a detailed and unified overview of the accelerator systems, target heating physics, beam transport, and diagnostic tools. Additional information includes calibration sources, radiation hydrodynamics, spectroscopic-quality radiation transport modeling, experimental measurements, and analyses of electron beam driven aluminum experiments. The first set of spatially and temporally resolved spectroscopic measurements of electron beam driven aluminum are presented. Contamination quantification analyses are used to understand the origin of the strong Na-I 3p-3s lines that are observed in absorption within the aluminum plasma continuum. These results inform the creation of the first spectroscopic-quality radiation transport model that links several atomic physics codes to interpret the conditions from which the Na-I lines originate. A good agreement is found between the surface analysis results and the model which confirms the concentration of the sodium present within the aluminum alloy foil material. It also demonstrates, for the first time on electron beam-driven target experiments, the ability to interpret plasma conditions from measured absorption lines. In a second experimental campaign which focuses on pure aluminum, the Al-I 3p-4s and 3p3d doublets are both measured in emission. A detailed analysis of the Al-I 3p-4s doublet reveals that the lines undergo moderate self-absorption. A simple model of the self-absorption effect is successfully used to match measured spectra at various temperature/density/plasma scale length combinations. These measurements led to the realization of the minimum density that can be resolved by the spectrometers for the Al-I 3p-4s lines due to the large slit width required to observe a signal on aluminum. These measurements also demonstrate the sensitivity of visible and long wave UV spectroscopy to minor changes in both temperature and density. The simple self-absorption model will be useful for analysis of other beam-target interaction experiments with spectra exhibiting either self-absorption or full self-reversal. Substantial headway has been made on the modeling front by linking together several codes needed to model both the energy deposition, hydrodynamic motion, and atomic kinetics to produce synthetic spectral calculations that are compared with experimental measurements. There exists ample space for improvement, especially with benchmarking the hydrodynamics codes and equation-of-state tables with experimental measurements. The X-ray diagnostics required to make these new measurements along with the simulation capabilities required to interpret the results are under development and will be the subject of future studies.

43 PARTICLE ACCELERATORS↗