Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic efficiency”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Development and Benchmarking of a Hybrid PIC Code For Dense Plasmas and Fast Ignition (Final Report)

Radiation processes play an important role in the study of both fast ignition and other inertial confinement schemes, such as plasma jet driven magneto-inertial fusion, both in their effect on energy balance, and in generating diagnostic signals. In the latter case, warm and hot dense matter may be produced by the convergence of a plasma shell formed by the merging of an assembly of high Mach number plasma jets. This innovative approach has the potential advantage of creating matter of high energy densities in voluminous amount compared with high power lasers or particle beams. An important application of this technology is as a plasma liner for the flux compression of magnetized plasma to create ultra-high magnetic fields and burning plasmas. HyperV Technologies Corp. has been developing plasma jet accelerator technology in both coaxial and linear railgun geometries to produce plasma jets of sufficient mass, density, and velocity to create such imploding plasma liners. An enabling tool for the development of this technology is the ability to model the plasma dynamics, not only in the accelerators themselves, but also in the resulting magnetized target plasma and within the merging/interacting plasma jets during transport to the target. Welch pioneered numerical modeling of such plasmas (including for fast ignition) using the LSP simulation code. Lsp is an electromagnetic, parallelized, plasma simulation code under development since 1995. It has a number of innovative features making it uniquely suitable for modeling high energy density plasmas including a hybrid fluid model for electrons that allows electrons in dense plasmas to be modeled with a kinetic or fluid treatment as appropriate. In addition to in-house use at Voss Scientific, several groups carrying out research in Fast Ignition (LLNL, SNL, UCSD, AWE (UK), and Imperial College (UK)) also use LSP. A collaborative team consisting of HyperV Technologies Corp., Voss Scientific LLC, FAR-TECH, Inc., Prism Computational Sciences, Inc. and Advanced Energy Systems Inc. joined efforts to develop new physics and numerical models for LSP in several key areas to enhance the ability of LSP to model high energy density plasmas (HEDP). This final report details those efforts. Areas addressed in this research effort include: adding radiation transport to LSP, first in 2D and then fully 3D, extending the EMHD model to 3D, implementing more advanced radiation and electrode plasma boundary conditions, and installing more efficient implicit numerical algorithms to speed complex 2-D and 3-D computations. The new capabilities allow modeling of the dominant processes in high energy density plasmas, and further assist the development and optimization of plasma jet accelerators, with particular attention to MHD instabilities and plasma/wall interaction (based on physical models for ion drag friction and ablation/erosion of the electrodes). In the first funding cycle we implemented a solver for the radiation diffusion equation. To solve this equation in 2-D, we used finite-differencing and applied the parallelized sparse-matrix solvers in the PETSc library (Argonne National Laboratory) to the resulting system of equations. A database of the necessary coefficients for materials of interest was assembled using the PROPACEOS and ATBASE codes from Prism. The model was benchmarked against Prism's 1-D radiation hydrodynamics code HELIOS, and against experimental data obtained from HyperV's separately funded plasma jet accelerator development program. Work in the second funding cycle focused on extending the radiation diffusion model to full 3-D, continued development of the EMHD model, optimizing the direct-implicit model to speed up calculations, add in multiply ionized atoms, and improved the way boundary conditions are handled in LSP. These new LSP capabilities were then used, along with analytic calculations and Mach2 runs, to investigate plasma jet merging, plasma detachment and transport, restrike and advanced jet accelerator design. In addition, a strong linkage to diagnostic measurements was made by modeling plasma jet experiments on PLX to support benchmarking of the code. A large number of upgrades and improvements advancing hybrid PIC algorithms were implemented in LSP during the second funding cycle. These include development of fully 3D radiation transport algorithms, new boundary conditions for plasma-electrode interactions, and a charge conserving equation of state that permits multiply ionized high-Z ions. The final funding cycle focused on 1) mitigating the effects of a slow-growing grid instability which is most pronounced in plasma jet frame expansion problems using the two-fluid Eulerian remap algorithm, 2) extension of the Eulerian Smoothing Algorithm to allow EOS/Radiation modeling, 3) simulations of collisionless shocks formed by jet merging, 4) simulations of merging jets using high-Z gases, 5) generation of PROPACEOS EOS/Opacity databases, 6) simulations of plasma jet transport experiments, 7) simulations of plasma jet penetration through transverse magnetic fields, and 8) GPU PIC code development. The tools developed during this project are applicable not only to the study of plasma jets, but also to a wide variety of HEDP plasmas of interest to DOE, including plasmas created in short-pulse laser experiments performed to study fast ignition concepts for inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Influence of turbulence and thermophysical fluid properties on cavitation erosion predictions in channel flow geometries

Cavitation and cavitation-induced erosion have been observed in fuel injectors in regions of high acceleration and low pressure. Although these phenomena can have a large influence on the performance and lifetime of injector hardware, questions still remain on how these physics should be accurately and efficiently represented within a computational fluid dynamics model. While several studies have focused on the validation of cavitation predictions within canonical and realistic injector geometries, it is not well documented what influence the numerical and physical parameters selected to represent turbulence and phase change will have on the predictions for cavitation erosion propensity and severity. In this work, a range of numerical and physical parameters are evaluated within the mixture modeling approach in CONVERGE to understand their influence on predictions of cavitation, condensation and erosion. Particular attention is paid to grid resolution, turbulence model and near-wall treatment, fuel surrogate properties, and non-condensable gas content. Assessment of cavitation predictions are conducted through comparison of measured and predicted mass flow rates and cavitation probability distributions for flow through a channel with a sharp inlet. Predictions for hydrodynamic impact loading and cavitation erosion are compared with the experimentally measured incubation period and critical site for erosion. Based on these findings, recommendations are provided for modeling turbulent cavitating flows, using the single fluid mixture modeling approach, to improve predictions for cavitation-induced erosion. In particular, to capture the fluid dynamic phenomena characterizing cavitation cloud formation, development and shedding, a Large Eddy simulation with grid resolution as fine as 2.50 μm is recommended. The assumed concentration of non-condensable gas content is observed to have a strong influence on the predicted cavitation erosion severity, which motivates the need for dissolved gas concentration measurements for future cavitation erosion experimental studies. Using the best practices established in this work, good agreement is achieved between the measured and predicted cavitation parameters, as well as the critical site for cavitation-induced erosion.

Magnotti, Gina Maureen↗
DOE OSTI2018

A gradient based method for modeling baryons and matter in halos of fast simulations

Fast N-body PM simulations with a small number of time steps such as FastPM or COLA have been remarkably successful in modeling the galaxy statistics, but their lack of small scale force resolution and long time steps cannot give accurate halo matter profiles or matter power spectrum. High resolution N-body simulations can improve on this, but lack baryonic effects, which can only be properly included in hydro simulations. Here we present a scheme to calibrate the fast simulations to mimic the precision of the hydrodynamic simulations or high resolution N-body simulations. The scheme is based on a gradient descent of either effective gravitational potential, which mimics the short range force, or of effective enthalpy, which mimics gas hydrodynamics and feedback. The scheme is fast and differentiable, and can be incorporated as a post-processing step into any simulation. It gives very good results for the matter power spectrum for several of the baryonic feedback and dark matter simulations, and also gives improved dark matter halo profiles. The scheme is even able to find the large subhalos, and increase the correlation coefficient between the fast simulations and the high resolution N-body or hydro simulations. It can also be used to add baryonic effects to the high resolution N-body simulations. While the method has free parameters that can be calibrated on various simulations, they can also be viewed as astrophysical nuisance parameters describing baryonic effects that can be marginalized over during the data analysis. In this view these parameters can be viewed as an efficient parametrization of baryonic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Hydrodynamic Test Requirements Process Improvements

Hydrodynamic testing at Los Alamos National Laboratory would benefit from a process improvement for the requirements process. Cameo was used as a digital solution for requirements management to allow Lead Engineers to track requirements more effectively. This was identified as a process improvement throughout this Capstone project. This report includes a project proposal, business case, literature review, methodology, project plan, data analysis, decision-making report, financial analysis report, discussion, and conclusion. Initially, this project focused on figuring out a solution for the hydrotest requirement process improvements. The scope narrowed to focus on the use of Cameo for requirement capture and management. During this Capstone, four tests had digital models produced for requirements management in Cameo. The initial model was the baseline, with core requirements used across the tests. Commonalities in tests were used and the core requirements allowed for process efficiencies. In the data analysis, it was seen that overall, the implementation of using Cameo for requirements resulted in a decreasing trend for both schedule and normalized cost. Tests have different complexity levels which is also a factor in how long the requirements process will take. Additional data is needed to continue analyzing process improvements. Through decisionmaking and financial analysis, the recommendation was to use Cameo for requirements process improvement. Multiple experts provided feedback for requirements that were then captured within models. Numerous tangible and intangible benefits were identified with this process improvement. For return on investment, the metric of success was schedule reduction, which was overall seen. Four tests were analyzed, so future analysis will be needed. There is also not a great financial risk because the main cost would be purchasing more licenses. Individuals must generate requirements whether using this software or not. Overall, training is needed to help improve the skillsets of Lead Engineers but is already being supported on a regular frequency. A desktop guide was started associated with explaining the process, however, it is a work in progress. The team plans on adding additional information in the digital models to help status when requirements are met using verification methods and artifacts. From working on this project an improved understanding of Cameo and requirements was the result. There are future opportunities to extend the usage for requirements management and progress will continue after this project.

99 GENERAL AND MISCELLANEOUS↗
NASA NTRS1994

Mars atmospheric loss and isotopic fractionation by solar-wind-induced sputtering and photochemical escape

We examine the effects of the loss of Mars atmospheric constituents by solar-wind-induced sputtering and by photochemical escape during the past 3.8 billion years. Sputtering is capable of efficiently removing species from the upper atmosphere, including the light noble gases; nitrogen and oxygen are removed by photochemical processes as well. Due to diffusive separation (by mass) above the homopause, removal from the top of the atmosphere will fractionate the isotopes of each species, with the lighter mass being preferentially lost. For carbon and oxygen, this allows us to determine the size of nonatmospheric reservoirs which mix with the atmosphere; these reservoirs can be CO2 adsorbed in the regolith and H2O in the polar ice caps. We have constructed both simple analytical models and time-dependent models of the loss of volatiles from and supply to the martian atmosphere. Both argon and neon require continued replenishment from outgassing over geologic time. For argon, sputtering loss explains the fractionation of (Ar-36)/(Ar-38) without requiring a distinct epoch of hydrodynamic escape (although fractionation of Xe isotopes still requires very early hydrodynamic loss). For neon, the current (Ne-22)/(Ne-20) ratio represents a balance between loss to space and continued resupply from the interior; the similarity of the ratio to the terrestrial value is coincidental. For nitrogen, the loss by both sputtering and photochemical escape would produce a fractionation of (N-15)/(N-14) larger than observed; an early, thicker carbon dioxide atmosphere could mitigate the nitrogen loss and produce the observed fractionation, as could continued outgassing of juvenile nitorgen. Based on the isotopic constraints, the total amount of carbon dioxide lost over geologic time is probably on the order of tens of millibars rather than a substantial fraction of a bar. The total loss from solar-wind-induced sputtering and photochemical escape, therefore, does not seem able to explain the loss of a putative thick, early atmosphere withput requiring formation of extensive surface carbonate deposits or other nonatmospheric reservoirs for CO2.

Jakosky, Bruce M.↗
DOE OSTI2021

Real-Time Nonlinear Model Predictive Controller for Multiple Degrees of Freedom Wave Energy Converters with Non-Ideal Power Take-Off

An increase in wave energy converter (WEC) efficiency requires not only consideration of the nonlinear effects in the WEC dynamics and the power take-off (PTO) mechanisms, but also more integrated treatment of the whole system, i.e., the buoy dynamics, the PTO system, and the control strategy. It results in an optimization formulation that has a nonquadratic and nonstandard cost functional. This article presents the application of real-time nonlinear model predictive controller (NMPC) to two degrees of freedom point absorber type WEC with highly nonlinear PTO characteristics. The nonlinear effects, such as the fluid viscous drag, are also included in the plant dynamics. The controller is implemented on a real-time target machine, and the WEC device is emulated in real-time using the WECSIM toolbox. The results for the successful performance of the design are presented for irregular waves under linear and nonlinear hydrodynamic conditions.

Haider, Ali S. (ORCID:0000000246513157)↗
NASA NTRS2013

The 0.5-2.22 micrometer Scattered Light Spectrum of the Disk around TW Hya: Detection of a Partially Filled Disk Gap at 80 AU*

We present a 0.5-2.2 micrometer scattered light spectrum of the circumstellar disk around TW Hya from a combination of spatially resolved Hubble Space Telescope STIS spectroscopy and NICMOS coronagraphic images of the disk. We investigate the morphology of the disk at distances greater than 40 AU over this wide range of wavelengths, and identify the presence of a depression in surface brightness at approximately 80 AU that could be caused by a gap in the disk. Additionally, we quantify the surface brightness, azimuthal symmetry, and spectral character of the disk as a function of radius. Our analysis shows that the scattering efficiency of the dust is largely neutral to blue over the observed wavelengths. We model the disk as a steady α-disk with an ad hoc gap structure. The thermal properties of the disk are selfconsistently calculated using a three-dimensional radiative transfer code that uses ray tracing to model the heating of the disk interior and scattered light images. We find a good fit to the data over a wide range of distances from the star if we use a model disk with a partially filled gap of 30% depth at 80 AU and with a self-similar truncation knee at 100 AU. The origin of the gap is unclear, but it could arise from a transition in the nature of the disk's dust composition or the presence of a planetary companion. Based on scalings to previous hydrodynamic simulations of gap-opening criteria for embedded proto-planets, we estimate that a planetary companion forming the gap could have a mass between 6 and 28 solar mass.

radiative transfer↗
DOE OSTI2015

A high-order vertex-based central ENO finite-volume scheme for three-dimensional compressible flows

High-order discretization methods offer the potential to reduce the computational cost associated with modeling compressible flows. However, it is difficult to obtain accurate high-order discretizations of conservation laws that do not produce spurious oscillations near discontinuities, especially on multi-dimensional unstructured meshes. A novel, high-order, central essentially non-oscillatory (CENO) finite-volume method that does not have these difficulties is proposed for tetrahedral meshes. The proposed unstructured method is vertex-based, which differs from existing cell-based CENO formulations, and uses a hybrid reconstruction procedure that switches between two different solution representations. It applies a high-order k-exact reconstruction in smooth regions and a limited linear reconstruction when discontinuities are encountered. Both reconstructions use a single, central stencil for all variables, making the application of CENO to arbitrary unstructured meshes relatively straightforward. The new approach was applied to the conservation equations governing compressible flows and assessed in terms of accuracy and computational cost. For all problems considered, which included various function reconstructions and idealized flows, CENO demonstrated excellent reliability and robustness. Up to fifth-order accuracy was achieved in smooth regions and essentially non-oscillatory solutions were obtained near discontinuities. The high-order schemes were also more computationally efficient for high-accuracy solutions, i.e., they took less wall time than the lower-order schemes to achieve a desired level of error. In one particular case, it took a factor of 24 less wall-time to obtain a given level of error with the fourth-order CENO scheme than to obtain the same error with the second-order scheme.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

THE 0.5-2.22 {mu}m SCATTERED LIGHT SPECTRUM OF THE DISK AROUND TW Hya: DETECTION OF A PARTIALLY FILLED DISK GAP AT 80 AU

We present a 0.5-2.2 {mu}m scattered light spectrum of the circumstellar disk around TW Hya from a combination of spatially resolved Hubble Space Telescope STIS spectroscopy and NICMOS coronagraphic images of the disk. We investigate the morphology of the disk at distances >40 AU over this wide range of wavelengths, and identify the presence of a depression in surface brightness at {approx}80 AU that could be caused by a gap in the disk. Additionally, we quantify the surface brightness, azimuthal symmetry, and spectral character of the disk as a function of radius. Our analysis shows that the scattering efficiency of the dust is largely neutral to blue over the observed wavelengths. We model the disk as a steady {alpha}-disk with an ad hoc gap structure. The thermal properties of the disk are self-consistently calculated using a three-dimensional radiative transfer code that uses ray tracing to model the heating of the disk interior and scattered light images. We find a good fit to the data over a wide range of distances from the star if we use a model disk with a partially filled gap of 30% depth at 80 AU and with a self-similar truncation knee at 100 AU. The origin of the gap is unclear, but it could arise from a transition in the nature of the disk's dust composition or the presence of a planetary companion. Based on scalings to previous hydrodynamic simulations of gap-opening criteria for embedded proto-planets, we estimate that a planetary companion forming the gap could have a mass between 6 and 28 M{sub Circled-Plus }.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Modeling the Effects of Tidal Energy Extraction on Estuarine Hydrodynamics in a Stratified Estuary

A three-dimensional coastal ocean model with a tidal turbine module was used in this paper to study the effects of tidal energy extraction on temperature and salinity stratification and density driven two-layer estuarine circulation. Numerical experiments with various turbine array configurations were carried out to investigate the changes in tidally mean temperature, salinity and velocity profiles in an idealized stratified estuary that connects to coastal water through a narrow tidal channel. The model was driven by tides, river inflow and sea surface heat flux. To represent the realistic size of commercial tidal farms, model simulations were conducted based on a small percentage of the total number of turbines that would generate the maximum extractable energy in the system. Model results indicated that extraction of tidal energy will increase the vertical mixing and decrease the stratification in the estuary. Extraction of tidal energy has stronger impact on the tidally-averaged salinity, temperature and velocity in the surface layer than the bottom. Energy extraction also weakens the two-layer estuarine circulation, especially during neap tides when tidal mixing the weakest and energy extraction is the smallest. Model results also show that energy generation can be much more efficient with higher hub height with relatively small changes in stratification and two-layer estuarine circulation.

Yang, Zhaoqing↗
DOE OSTI2025

Thermal energy transport in laser-driven high x-ray conversion efficiency metallic silver nanowire foams

Here, we report on record brightness from Ag x-ray emission obtained using a novel laser-produced plasma source. The reported K-shell conversion efficiency of nearly 1% with a radiant energy of ≈ 0.6 kJ/sr from Ag ions is the highest presently recorded and is about twofold greater than more conventional metal-lined cavity targets. He-like Ag ions are the dominant radiators at ≈ 22.7 keV, which contrast the x-ray sources from other pulse power facilities that produce K-shell x-rays from Auger processes in near-neutral ions driven by the nonthermal hot electrons produced in the Z-pinch implosion. The reported x-ray emission was produced from a 4-mm-diameter, 4-mm-long underdense silver nanowire target at 10⁢ mg/⁢cm 3 . The National Ignition Facility laser beams deposited ≈ 1000 kJ of 3ω light into the target with a ≈ 400 TW, 2.5 ns square pulse. Finally, we show via comparisons of data with radiation-hydrodynamics simulations that flux-limited classical thermal energy transport models are inadequate in correctly modeling the behavior of these non-local thermodynamic equilibrium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Neutrino Oscillations in Supernovae

A (core-collapse) supernova is the final death call of a massive star during which a gigantic amount of a special kind of elementary particles known as neutrinos are released for a brief period of time. An earlier work by the PI and his collaborators has shown that the dense neutrino gas deep inside the supernova can experience "flavor transformation" or "oscillations" in a collective fashion, a macroscopic quantum phenomenon spanning over distances of 100-1000 kilometers. During this process neutrinos of different energies and propagation directions change their identities simultaneously. This phenomenon may be a key to understanding the mechanism of supernova explosion and the origin of many of the chemical elements found on the earth which were produced in the supernovae exploded in the past. The previous understanding of neutrino oscillations in supernovae was based on a stationary 3-dimensional model with 1 spatial dimension and 2 momentum dimensions. In this "neutrino bulb model" neutrinos are emitted from a spherically symmetric surface and they are decoupled from the motion of matter and the production of chemical elements. In order to accurately model the neutrino oscillations in supernovae and determine its physical consequences, sophisticated 7-dimensional (1 temporal dimension + 3 spatial dimensions + 3 momentum dimensions) modelling of neutrino transport must be performed self-consistently together with the hydrodynamic simulation of the flow of energy and matter inside the supernova and with the calculation of the production and destruction of chemical elements through a large nucleosynthesis network. Although it will be long journey to achieve the lofty goal, our project has made important progress in three different fronts towards this goal. In the first front we explored new computing hardware on which the next generation simulation codes may be run. We ported the neutrino bulb code to the platform of Graphics Processing Unit (GPU) and then later developed a new bulb code that can be run on supercomputers equipped with the Intel Xeon CPUs and/or the Xeon Phi coprocessors. In the second front we developed new algorithms that can be used towards computing neutrino oscillations in multi-dimensional models. We have developed a moment expansion method which is approximately 10 times more efficient than the old algorithm. We have also developed numerical codes to solve neutrino oscillations in toy models with two spatial dimensions. In the last front we used analytic methods to study neutrino oscillations in more sophisticated models, e.g. of multiple spatial dimensions and/or with time dependence. By analyzing the neutrino oscillations in toy models with two spacial dimensions we found that the approximate spacial translation symmetry, which is similar to the spherical symmetry in the neutrino bulb model, can be broken spontaneously as neutrino oscillations begin. By analyzing a time-dependent neutrino bulb model we found that certain modes of neutrino oscillations that are suppressed in the stationary model can become active in the time-dependent model. We have also developed theories which give insights to how ordinary matter may affect neutrino oscillations in supernovae and other astrophysical environments.

Duan, Huaiyu↗
DOE OSTI2015

Laser-ablated loading of solid target through foams of overcritical density

The main objective of the work is to study the conversion of the laser pulse energy into the energy of the hydrodynamic motion of matter in a solid target following the initial absorption of laser radiation in a layer of porous material. Results of experiments on plane massive aluminum targets, coated with a layer of porous plastic with density greater than the critical density of the plasma created, are presented. Experiments were carried out on the laser installation ABC of the Research Center ENEA-Frascati; the targets were irradiated by a beam of the fundamental harmonic of Nd-laser radiation with an energy of about 50 kJ, intensity of 10{sup 13 }W/cm{sup 2}, and 3 ns duration. The experimental method consisted in measuring the volume of the craters created on the aluminum surface behind various thicknesses and densities of the porous absorber of laser radiation. On the basis of these measurements and of an advanced analytical model, quantitative conclusions are made on how the efficiency of laser energy transfer to the solid part of the target (laser-ablated loading) depends on thickness and density of the porous absorber.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hydrodynamics of a multicomponent vesicle under strong confinement

Herein we numerically investigate the hydrodynamics and membrane dynamics of a multicomponent vesicle in two strongly confined geometries. This serves as a simplified model for red blood cells undergoing large deformations while traversing narrow constrictions. We propose a new parameterization for the bending modulus that remains positive for all lipid phase parameter values. For a multicomponent vesicle passing through a stenosis, we establish connections between various properties: lipid phase coarsening, size and flow profile of the lubrication layers, excess pressure, and the tank-treading velocity of the membrane. For a multicomponent vesicle passing through a contracting channel, we find that the lipid always phase separates so that the vesicle is stiffer in the front as it passes through the constriction. For both cases of confinement we find that lipid coarsening is arrested under strong confinement, and resumes at a high rate upon relief from extreme confinement. The results may be useful for efficient sorting lipid domains using microfluidic flows by controlled release of vesicles passing through strong confinement.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2022

After DART: Using the first full-scale test of a kinetic impactor to inform a future planetary defense mission

NASA’s Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos’s response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor β, showing that a particular direction-specific β will be directly determined by the DART results, and that a related direction-specific β is a figure of merit for a kinetic impact mission. The DART β determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos’s near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in-situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction.

Asteroids↗
DOE OSTI2016

Deflection by kinetic impact: Sensitivity to asteroid properties

Impacting an asteroid with a spacecraft traveling at high speed delivers an impulsive change in velocity to the body. In certain circumstances, this strategy could be used to deflect a hazardous asteroid, moving its orbital path off of an Earth-impacting course. However, the efficacy of momentum delivery to asteroids by hypervelocity impact is sensitive to both the impact conditions (particularly velocity) and specific characteristics of the target asteroid. We numerically model asteroid response to kinetic impactors under a wide range of initial conditions, using an Adaptive Smoothed Particle Hydrodynamics code. Impact velocities spanning 1–30 km/s were investigated, yielding, for a particular set of assumptions about the modeled target material, a power-law dependence consistent with a velocity-scaling exponent of μ = 0.44. Target characteristics including equation of state, strength model, porosity, rotational state, and shape were varied, and corresponding changes in asteroid response were documented. Moreover, the kinetic-impact momentum-multiplication factor, β, decreases with increasing asteroid cohesion and increasing porosity. Although increased porosity lowers β, larger porosities result in greater deflection velocities, as a consequence of reduced target masses for asteroids of fixed size. Porosity also lowers disruption risk for kinetic impacts near the threshold of disruption. Including fast (P = 2.5 h) and very fast (P = 100 s) rotation did not significantly alter β but did affect the risk of disruption by the impact event. Asteroid shape is found to influence the efficiency of momentum delivery, as local slope conditions can change the orientation of the crater ejecta momentum vector. Our results emphasize the need for asteroid characterization studies to bracket the range of target conditions expected at near-Earth asteroids while also highlighting some of the principal uncertainties associated with the kinetic-impact deflection strategy.

58 GEOSCIENCES↗
DOE OSTI2014

Fingering convection induced by atomic diffusion in stars: 3D numerical computations and applications to stellar models

Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature, and detailed structure can affect the excitation process of various stellar pulsation modes and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, as well as its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection (which neglect radiative transfer and complex opacity effects) and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (within the scope of the aforementioned approximation) and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. prescription in the Toulouse-Geneva Evolution Code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al., that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases it ultimately triggers dynamical convection. The direct effects of radiative acceleration on the dynamics of fingering convection (especially in the nonlinear regime) remain to be added in the future to improve the quantitative predictions of the model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Heat Integration Optimization and Dynamic Modeling Investigation for Advancing the Coal-Direct Chemical Looping Process

The purpose of the project is to address the optimization and startup operation of a modular coal direct chemical looping (CDCL) combustion system integrated with a steam cycle for power generation to reduce the risks involved in further scale-up of the technology. The modular reactor design of the CDCL process provides flexibility in the fabrication of the reactor and in its operating capacity (i.e. turndown ratio) at the cost of a more complex heat exchange network (HEN) design and integration. To address the technology gaps and advance the efficiency and economic feasibility of the CDCL technology, the project will perform a detailed and comprehensive analysis of the integration of a modular CDCL reactor system and a steam cycle system under both static and transient conditions via HEN process performance simulations and system dynamic modeling, respectively. The scope of work consists of 1) Experimental and computational studies of the CDCL combustor reactor 2) Comprehensive static (i.e. steady-state) system HEN design analysis in CDCL 550 MWe commercial unit for power generation and 3) Dynamic modeling of site specific design of 10MWe CDCL large pilot plant. The project team has successfully developed and validated a kinetic model for the oxidation of oxygen carriers in the combustor using the unreacted shrinking core model (UCSM). The model is capable of capturing the oxidation kinetics of fully or partially reduced oxygen carrier particles. A computational fluid dynamics (CFD) model is developed to simulate the hydrodynamics, heat transfer, and chemical reaction occurring in the CDCL combustor. The model is developed in MFIX and ANSYS Fluent. Key aspects of CDCL combustor operation, including heat transfer, oxygen carrier oxidation, and the transport of oxygen carrier particles, are simulated using this CFD model. The HEN for a commercial scale 550 MWe CDCL power plant is simulated and optimized using ASPEN Plus. Practical design considerations are incorporated based on industrial experiences. The performance and cost for the commercial CDCL plant is updated based on these analyses. A dynamic model for the 10 MWe CDCL pilot plant is developed in ProTRAX simulation software. The model is based on the pilot plant design developed in project DE-FE0027654 “10 MWe CDCL Large Pilot Plang – Pre-FEED Study” and the steam cycle data obtained from Dover Light & Power plant. The transient behaviors during pilot plant load variation are simulated using the dynamic model.

01 COAL, LIGNITE, AND PEAT↗