Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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

CRADA Final Report: CRADA Number NFE-22-09330 with General Fusion

General Fusion is developing a magnetized target fusion (MTF) approach that involves compressing an initial magnetically confined plasma inside a cavity formed in liquid metal. This approach builds from concepts initially developed under the Linus program at the U.S. Naval Research Laboratory and combines it with advances from compact toroid experiment (CTX) and sustained spheromak plasma experiment (SSPX) in compact toroid plasmas and coaxial Marshall gun systems. Modeling the tokamak during compression is central to designing a successful MTF device. The plasma is formed by coaxial helicity injection in the General Fusion device. Immediately after formation, the plasma has a diverted tokamak configuration with a single null. As the wall moves inwards, the plasma is repelled from the conducting surface and driven inwards by currents induced by its magnetic field in the liquid metal wall. As the liquid metal closes (or bridges) the opening of the coaxial plasma injector, the magnetic field topology alters to remove the null. Due to this, the plasma moves from a diverted to a wall-limited configuration. The liquid metal liner continues to close in and change shape, reducing in radius by a factor of ten at the peak of plasma compression. A model of the MTF plasma must be able to handle this continually varying geometry, and to be predictive, it must faithfully include the real imperfections arising in the process. In this project, we pursued a Monte Carlo approach to closures for MHD by computing kinetic electron trajectories in an MHD plasma background from simulations of GF devices. This requires enhancing the capabilities of the KORC-T code for running large ensembles of kinetic trajectories by porting it to GPU architectures and enabling workflows for large ensembles on OLCF machines. With these capabilities, it is possible to produce a large library of kinetic calculations of electron orbits evolving in plasma configurations spanning the magnetic configurations and plasma density profiles, including non-axisymmetry, arising in the General Fusion’s existing PI3 spherical tokamak device. Using ensembles will capture particles passing a single point in space in a given magnetic configuration, and the entire dataset will cover a range of global magnetic field geometries. By sampling around many starting points, this dataset will capture the spatial dependence of the plasma parameters. From this large dataset, it is possible to produce a reduced model for the kinetic effects not captured in MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Comparison between the South Pacific Convergence Zone and South Atlantic Convergence Zone in January 1979

Using one of the FGGE data sets, the eddy potential and eddy kinetic energy budgets during the period of January 10-27, 1979 were derived for the regions containing the South Pacific Convergence Zone (SPCZ) and the South Atlantic Convergence Zone. A complete set of energy equations developed by Brennan and Vincent (1980) was partitioned into zonal and eddy components; the eddy forms of energy were derived from a grid point's departure from the zonal average over a 5-deg latitude belt. The results of the analysis indicate that the eddy kinetic energy content (KE) is greater than the eddy available potential energy content (AE) throughout the period. Both AE and KE show large losses beginning when the SPCZ enters a decaying state: AE shows a gradual but slow decrease until after January 22, when it decreases rapidly; KE oscillates, reaching a maximum on January 22 before a rapid decrease.

Vincent, Dayton G.↗
DOE OSTI2023

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗
NASA NTRS1976

Hydrogen-fueled scramjets: Potential for detailed combustor analysis

Combustion research related to hypersonic scramjet (supersonic combustion ramjet) propulsion is discussed from the analytical point of view. Because the fuel is gaseous hydrogen, mixing is single phase and the chemical kinetics are well known; therefore, the potential for analysis is good relative to hydro-carbon fueled engines. Recent progress in applying two and three dimensional analytical techniques to mixing and reacting flows indicates cause for optimism, and identifies several areas for continuing effort.

Beach, H. L., Jr.↗
NASA NTRS1987

A cure-rate model for the Shuttle filament-wound case

An epoxy and carbon fiber composite has been used to produce a light-weight rocket case for the Space Shuttle. A kinetic model is developed which can predict the extent of epoxy conversion during the winding and curing of the case. The model accounts for both chemical and physical kinetics. In the model, chemical kinetics occur exclusively up to the time the transition temperature equals the reaction temperature. At this point the resin begins to solidify and the rate of this process limits the rate of epoxy conversion. A comparison of predicted and actual epoxy conversion is presented for isothermal and temperature programmed cure schedules.

Cagliostro, D. E.↗
NASA NTRS1994

Numerical investigation of supersonic turbulent boundary layers with high wall temperature

A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Guo, Y.↗
NASA NTRS1992

Computation of Turbulent Recirculating Flow in Channels, and for Impingement Cooling

Fully elliptic forms of the transport equations have been solved numerically for two flow configurations. The first is turbulent flow in a channel with transverse rectangular ribs, and the second is impingement cooling of a plane surface. Both flows are relevant to proposed designs for active cooling of hypersonic vehicles using supercritical hydrogen as the coolant. Flow downstream of an abrupt pipe expansion and of a backward-facing step were also solved with various near-wall turbulence models as benchmark problems. A simple form of periodicity boundary condition was used for the channel flow with transverse rectangular ribs. The effects of various parameters on heat transfer in channel flow with transverse ribs and in impingement cooling were investigated using the Yap modified Jones and Launder low Reynolds number k-epsilon turbulence model. For the channel flow, predictions were in adequate agreement with experiment for constant property flow, with the results for friction superior to those for heat transfer. For impingement cooling, the agreement with experiment was generally good, but the results suggest that improved modelling of the dissipation rate of turbulence kinetic energy is required in order to obtain improved heat transfer prediction, especially near the stagnation point. The k-epsilon turbulence model was used to predict the mean flow and heat transfer for constant and variable property flows. The effect of variable properties for channel flow was investigated using the same turbulence model, but comparison with experiment yielded no clear conclusions. Also, the wall function method was modified for use in the variable properties flow with a non-adiabatic surface, and an empirical model is suggested to correctly account for the behavior of the viscous sublayer with heating.

Chang, Byong Hoon↗
NASA NTRS2006

Observations of Shock Diffusion and Interactions in Supersonic Freestreams with Counterflowing Jets

One of the technical challenges in long-duration space exploration and interplanetary missions is controlled entry and re-entry into planetary and Earth atmospheres, which requires the dissipation of considerable kinetic energy as the spacecraft decelerates and penetrates the atmosphere. Efficient heat load management of stagnation points and acreage heating remains a technological challenge and poses significant risk, particularly for human missions. An innovative approach using active flow control concept is proposed to significantly modify the external flow field about the spacecraft in planetary atmospheric entry and re-entry in order to mitigate the harsh aerothermal environments, and significantly weaken and disperse the shock-wave system to reduce aerothermal loads and wave drag, as well as improving aerodynamic performance. To explore the potential benefits of this approach, we conducted fundamental experiments in a trisonic blow down wind tunnel to investigate the effects of counterflowing sonic and supersonic jets against supersonic freestreams to gain a better understanding of the flow physics of the interactions of the opposing flows and the resulting shock structure.

Daso, Endwell O.↗
NASA NTRS2018

Blunt-Body Paradox and Improved Application of Transient-Growth Framework

The “Reshotko-Tumin transition criterion" based on optimal transient growth successfully correlates laboratory measurements of roughness induced transition over blunt body configurations. Even though transient growth has not been conclusively linked to the measured onset of transition, the above correlation denotes the only available physics-based model for subcritical transition in blunt body flows, since the latter do not support any modal instabilities at typical experimental conditions. Unlike other established models based on empirical curve fits that are valid for a specific subclass of datasets, the optimal-growth-based transition criterion appears to provide a reasonable correlation with measurements in various wind tunnel and ballistic range facilities and for a broad range of surface temperature ratios. This paper is focused on optimal growth calculations that improve upon significant shortcomings of the computations underlying the Reshotko-Tumin correlation. The improved framework is applied to leeward transition over a spherical section forebody that was tested in the Mach 6 Adjustable Contour Expansion wind tunnel at Texas A&M University. The computed results highlight the significance of nonparallel basic state evolution, curvature terms, and an optimization procedure that varies both inflow and outflow locations of the transient growth interval. More important, the results indicate that the modified correlation is very close to its original form, and hence, that the accuracy of the transient-growth-based transition criterion is not compromised by using a more thorough theoretical framework. Yet the results also show that the optimal energy gain up to the predicted transition onset location can be rather small, highlighting the need to further investigate the optimal growth criterion for additional experimental configurations and to also uncover the in-depth physics underlying blunt body transition.

Freestream conditions↗
DOE OSTI2022

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

SCExAO/MEC and CHARIS Discovery of a Low-mass, 6 au Separation Companion to HIP 109427 Using Stochastic Speckle Discrimination and High-contrast Spectroscopy

We report the direct imaging discovery of a low-mass companion to the nearby accelerating A star, HIP 109427, with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument coupled with the Microwave Kinetic Inductance Detector Exoplanet Camera (MEC) and CHARIS integral field spectrograph. CHARIS data reduced with reference star point spread function (PSF) subtraction yield 1.1–2.4 μm spectra. MEC reveals the companion in Y and J band at a comparable signal-to-noise ratio using stochastic speckle discrimination, with no PSF subtraction techniques. Combined with complementary follow-up L {sub p} photometry from Keck/NIRC2, the SCExAO data favors a spectral type, effective temperature, and luminosity of M4–M5.5, 3000–3200 K, and log{sub 10}(L/L{sub ⊙})=−2.28{sub −0.04}{sup +0.04}, respectively. Relative astrometry of HIP 109427 B from SCExAO/CHARIS and Keck/NIRC2, and complementary Gaia–Hipparcos absolute astrometry of the primary favor a semimajor axis of 6.55{sup +3.0} {sub −0.48} au, an eccentricity of 0.54{sub −0.15}{sup +0.28}, an inclination of 66.7{sub −14}{sup +8.5} degrees, and a dynamical mass of 0.280{sub −0.059}{sup +0.18} M {sub ⊙}. This work shows the potential for extreme AO systems to utilize speckle statistics in addition to widely used postprocessing methods to directly image faint companions to nearby stars near the telescope diffraction limit.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS2000

Does Warming a Lysozyme Solution Cook Ones Data?

Chicken egg white lysozyme has a well characterized thermally driven phase transition. Between pH 4.0 and 5.2, the transition temperature, as defined by the point where the tetragonal and orthorhombic solubility are equal, is a function of the pH, salt (precipitant) type and concentration, and most likely of the buffer concentration as well. This phase transition can be carried out with protein solution alone, prior to the initiation of the crystallization process. We have now measured the kinetics of this process and investigated its reversibility. An aliquot of a stock protein solution is held at a given temperature, and at periodic intervals used to set up batch crystallization experiments. The batch solutions were incubated at 20 C until macroscopic crystals were obtained, at which point the number of crystals in each well were counted. The transition effects increased with temperature, slowly falling off at 30 C with a half time (time to approx. 1/2 the t = 0 number of crystals) of approx. 5 hours, and an estimated half time of approx. 0.5 hours at 43 C. Further, the process was not reversible by simple cooling. After holding a lysozyme solution at 37 C (prior to addition of precipitant) for 16 hours, then cooling and holding it at 4 C, no return to the pre-warmed nucleation kinetics are observed after at least 4 weeks. Thus every thermal excursion above the phase transition point results in a further decrease in the nucleation rate of that solution, the extent being a function of the time and temperature. Orthorhombic lysozyme crystals apparently do not undergo the flow-induced growth cessation of tetragonal lysozyme crystals. We have previously shown that putting the protein in the orthorhombic form does not affect the averaged face growth kinetics, only nucleation, for tetragonal crystals. We may be able to use this differential behavior to elucidate how flow affects tile lysozyme crystal growth process.

Pusey, Marc↗
DOE OSTI2023

Saddle point search with dynamic active volume

Sampling potential energy surface (PES) is critical for many problems in materials science, chemistry, physics, and biology and requires highly efficient saddle point searches (SPS). In the study, we introduce the concept of dynamic active volume (DAV) in addition to the active volume in self-evolving atomistic kinetic Monte Carlo (SEAKMC). The DAV method has further reduced the dimensionality of the PES at the elevation stage of a SPS. At the subsequent converging stage, the dynamic boundary in DAV is lifted to allow the system to converge to the right location of a saddle point. Coupled with the dimer method, the DAV method not only significantly reduces the time cost for a given search attempt, but also dramatically increases the probability of finding relevant saddle points for PES sampling. Further, a Python software package within the framework SEAKMC (SEAKMC_py) with the DAV method has been developed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Validation of SOLPS-ITER simulations with kinetic, fluid, and hybrid neutral models for JET-ILW low-confinement mode plasmas

For JET L-mode plasmas in low-recycling conditions (electron temperature at the outer strike point, T e,ot ≳ 30 eV ), SOLPS-ITER simulations agree within the error bars for the experimental profiles at the low-field side (LFS) divertor target. The peak Balmer-α (Dα ) emission in the LFS divertor agrees within the error bars of the KS3 filterscope diagnostic, but is approximately 30% lower than the peak value of the KT1 spectrometer. Simulations have been performed with fluid, kinetic, and hybrid models for the neutrals. The large fluid-kinetic discrepancies of more than a factor 2 are successfully corrected by using a hybrid fluid-kinetic approach, for which kinetic atoms are transferred to the fluid population when the local Knudsen number of the atom becomes smaller than a user-defined transition Knudsen number Kn t . The hybrid-kinetic discrepancies are limited to a few % for Kn t ≤ 100 . When increasing the upstream density to high-recycling conditions, at the onset of detachment ( T e,ot ≈ 5 eV ), the simulations predict more than a factor 2 lower peak ion saturation current to the LFS divertor than the experiments. Also the Dα emission is underpredicted with approximately a factor 2. For these high-recycling conditions, the fluid-kinetic discrepancies are limited to maximum 50%, which are again corrected by using the hybrid approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Investigation of Dynamic Oxygen Adsorption in Molten Solder Jetting Technology

Surface tension forces play a critical role in fluid dynamic phenomena that are important in materials processing. The surface tension of liquid metals has been shown to be very susceptible to small amounts of adsorbed oxygen. Consequently, the kinetics of oxygen adsorption can influence the capillary breakup of liquid-metal jets targeted for use in electronics assembly applications, where low-melting-point metals (such as tin-containing solders) are utilized as an attachment material for mounting of electronic components to substrates. By interpreting values of surface tension measured at various surface ages, adsorption and diffusion rates of oxygen on the surface of the melt can be estimated. This research program investigates the adsorption kinetics of oxygen on the surface of an atomizing molten-metal jet. A novel oscillating capillary jet method has been developed for the measurement of dynamic surface tension of liquids, and in particular, metal melts which are susceptible to rapid surface degradation caused by oxygen adsorption. The experimental technique captures the evolution of jet swells and necks continuously along the jet propagation axis and is used in conjunction with an existing linear, axisymmetric, constant-property model to determine the variation of the instability growth rate, and, in turn, surface tension of the liquid as a function of surface age measured from the exit orifice. The conditions investigated so far focus on a time window of 2-4ms from the jet orifice. The surface properties of the eutectic 63%Sn-37%Pb solder alloy have been investigated in terms of their variation due to O2 adsorption from a N2 atmosphere containing controlled amounts of oxygen (from 8 ppm to 1000 ppm). The method performed well for situations where the oxygen adsorption was low in that time window. The value of surface tension for the 63Sn-37Pb solder in pure nitrogen was found to be 0.49 N/m, in good agreement with previously published work. A characteristic time of O(1ms) or less was determined for the molten-metal surface to be saturated by oxygen at 1000 ppm concentration in N2.

Megaridis, Constantine M.↗
DOE OSTI2023

Kinetic model describing self-limiting CO 2 diffusion in supported amine adsorbents

A reaction–diffusion shrinking core model describing the decay in diffusivity of supported amine sorbents upon CO 2 sorption under both simulated direct air capture and point source capture conditions is described. The decay in CO 2 diffusivity is associated with crosslinking in the aminopolymer samples and general pore blockage in the amino-silane derived samples, which occurs as CO 2 is adsorbed. The model is used to extract four kinetic parameters that govern the CO 2 uptake kinetics and working capacity: an apparent reaction rate constant, an initial effective diffusivity, and two dimensionless decay parameters. Ideally, an initially reaction limited system would allow for direct determination of the intrinsic reaction rate constant; however, sorption experiments suggest mass transfer resistances related to gas mixing, external boundary layers and intraparticle diffusion are present. Reaction rate constants are determined and agree well with theoretical values predicted with the Eyring equation parameterized using density functional theory energies from literature sources. The kinetic performance is expressed as the average effective diffusivity as a function of average conversion, which can be correlated to the dispersion of sorption sites on the support and the morphology of the active sorbent phase. Four supports are impregnated or grafted with amines, SBA-15, single-walled zeolite nanotubes (ZNT), Syloid SiO 2 , and γ-Al 2 O 3 . Due to its pore structure, γ-Al 2 O 3 supported amines can be modeled at the μm scale or at the nm scale, where the shell balance is on the μm-sized macroporous particle aggregate or on the nm-sized amine film on the surface of the Al 2 O 3 nanoparticles, which comprise the spherical particle aggregates. Faster diffusion rates are maintained under 400 ppm rather than 10% CO 2 due to a slower reaction rate giving a slower decay in diffusivity. In conclusion, this work provides a first principles kinetic analysis of CO 2 sorption where previous models are semi-empirical and use arbitrary kinetic parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗