Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffuse layer model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
NASA NTRS2017

Computational Assessment of a 3-Stage Axial Compressor Which Provides Airflow to the NASA 11- by 11-Foot Transonic Wind Tunnel, Including Design Changes for Increased Performance

A 24 foot diameter 3-stage axial compressor powered by variable-speed induction motors provides the airflow in the closed-return 11- by 11-Foot Transonic Wind Tunnel (11-Foot TWT) Facility at NASA Ames Research Center at Moffett Field, California. The facility is part of the Unitary Plan Wind Tunnel, which was completed in 1955. Since then, upgrades made to the 11-Foot TWT such as flow conditioning devices and instrumentation have increased blockage and pressure loss in the tunnel, somewhat reducing the peak Mach number capability of the test section. Due to erosion effects on the existing aluminum alloy rotor blades, fabrication of new steel rotor blades is planned. This presents an opportunity to increase the Mach number capability of the tunnel by redesigning the compressor for increased pressure ratio. Challenging design constraints exist for any proposed design, demanding the use of the existing driveline, rotor disks, stator vanes, and hub and casing flow paths, so as to minimize cost and installation time. The current effort was undertaken to characterize the performance of the existing compressor design using available design tools and computational fluid dynamics (CFD) codes and subsequently recommend a new compressor design to achieve higher pressure ratio, which directly correlates with increased test section Mach number. The constant cross-sectional area of the compressor leads to highly diffusion factors, which presents a challenge in simulating the existing design. The CFD code APNASA was used to simulate the aerodynamic performance of the existing compressor. The simulations were compared to performance predictions from the HT0300 turbomachinery design and analysis code, and to compressor performance data taken during a 1997 facility test. It was found that the CFD simulations were sensitive to endwall leakages associated with stator buttons, and to a lesser degree, under-stator-platform flow recirculation at the hub. When stator button leakages were modeled, pumping capability increased by over 20 of pressure rise at design point due to a large reduction in aerodynamic blockage at the hub. Incorporating the stator button leakages was crucial to matching test data. Under-stator-platform flow recirculation was thought to be large due to a lack of seals. The effect of this recirculation was assessed with APNASA simulations recirculating 0.5, 1, and 2 of inlet flow about stators 1 and 2, modeled as axisymmetric mass flux boundary conditions on the hub before and after the vanes. The injection of flow ahead of the stators tended to re-energize the boundary layer and reduce hub separations, resulting in about 3 increased stall margin per 1 of inlet flow recirculated. In order to assess the value of the flow recirculation, a mixing plane simulation of the compressor which gridded the under-stator cavities was generated using the ADPAC CFD code. This simulation indicated that about 0.65 of the inlet flow is recirculated around each shrouded stator. This collective information was applied during the redesign of the compressor. A potential design was identified using HT0300 which improved overall pressure ratio by removing pre-swirl into rotor 1, replacing existing NASA 65 series rotors with double circular arc sections, and re-staggering rotors and the existing stators. The performance of the new design predicted by APNASA and HT0300 is compared to the existing design.

Turbomachinery↗
NASA NTRS2002

Numerical Simulations of Vortex Generator Vanes and Jets on a Flat Plate

Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Allan, Brian G.↗
DOE OSTI2023

Origin of Rapid Delithiation In Secondary Particles Of LiNi 0.8 Co 0.15 Al 0.05 O 2 and LiNi y Mn z Co 1− y − z O 2 Cathodes

Abstract Most research on the electrochemical dynamics in materials for high‐energy Li‐ion batteries has focused on the global behavior of the electrode. This approach is susceptible to misleading analyses resulting from idiosyncratic kinetic conditions, such as surface impurities inducing an apparent two‐phase transformation within LiNi 0.8 Co 0.15 Al 0.05 O 2 . Here, nano‐focused X‐ray probes are used to measure delithiation operando at the scale of secondary particle agglomerates in layered cathode materials during charge. After an initial latent phase, individual secondary particles undergo rapid, stochastic, and largely uniform delithiation, which is in contrast with the gradual increase in cell potential. This behavior reproduces across several layered oxides. Operando X‐ray microdiffraction (‐XRD) leverages the relationship between Li content and lattice parameter to further reveal that rate acceleration occurs between Li‐site fraction ( x Li ) ≈0.9 and ≈0.5 for LiNi 0.8 Co 0.15 Al 0.05 O 2 . Physics‐based modeling shows that, to reproduce the experimental results, the exchange current density ( i 0 ) must depend on x Li , and that i 0 should increase rapidly over three orders of magnitude at the transition point. The specifics and implications of this jump in i 0 are crucial to understanding the charge‐storage reaction of Li‐ion battery cathodes.

25 ENERGY STORAGE↗
DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗
NASA NTRS2014

Behavior of Particle Depots in Molten Silicon During Float-Zone Growth in Strong Static Magnetic Fields

In the photovoltaics industry, the largest market share is represented by solar cells made from multicrystalline silicon, which is grown by directional solidification. During the growth process, the silicon melt is in contact with the silicon nitride coated crucible walls and the furnace atmosphere which contains carbon monoxide. The dissolution of the crucible coating, the carbon bearing gas, and the carbon already present in the feedstock, lead to the precipitation of silicon carbide, and silicon nitride, at later stages of the growth process. The precipitation of Si3N4 and SiC particles of up to several hundred micrometers in diameter leads to severe problems during the wire sawing process for wafering the ingots. Furthermore the growth of the silicon grains can be negatively influenced by the presence of particles, which act as nucleation sources and lead to a grit structure of small grains and are sources for dislocations. If doped with Nitrogen from the dissolved crucible coating, SiC is a semi conductive material, and can act as a shunt, short circuiting parts of the solar cell. For these reasons, the incorporation of such particles needs to be avoided. In this contribution we performed model experiments in which the transport of intentionally added SiC particles and their interaction with the solid-liquid interface during float zone growth of silicon in strong steady magnetic fields was investigated. SiC particles of 7μm and 60μm size are placed in single crystal silicon [100] and [111] rods of 8mm diameter. This is achieved by drilling a hole of 2mm diameter, filling in the particles and closing the hole by melting the surface of the rod until a film of silicon covers the hole. The samples are processed under a vacuum of 1x10(exp -5) mbar or better, to prevent gas inclusions. An oxide layer to suppress Marangoni convection is applied by wet oxidation. Experiments without and with static magnetic field are carried out to investigate the influence of melt convection on the distribution of particles and their incorporation into the crystal. The field strengths applied by a superconducting magnet are 1T, 3T, 4.5T, and 5T. The increase in field strength dampens the melt flow, and so this study provides comparative data to the crystal growth experiment to be carried out onboard the sounding rocket mission TEXUS 51, where purely diffusive growth condition will be achieved under microgravity conditions.

Jauss, T.↗
NASA NTRS1992

Lidar multiple scattering: Dependence on atmospheric parameters

Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).

Wiegner, Matthias↗
NASA NTRS1993

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗
NASA NTRS1993

Numerical investigations in three-dimensional internal flows

In the present reporting period, the 3D version of the OVERFLOW code was used to solve the flow within the internal portion of the supersonic inlet. The internal portion of this inlet is bounded by an inflow plane containing the leading edge of the sidewalls, the sidewalls, the ramp and cowl surfaces and an outflow plane just downstream of the minimum geometric area of the inlet. Boundary layer bleed was used in the two-dimensional calculations discussed in the previous progress report and that same bleed was applied in the present study. For reference, this bleed corresponds to locations designated as R2 and R3 in the Mach 5 inlet model test. Using the GRIDGEN code, a three dimensional grid was generated that accounted for the viscous effects expected to occur on the sidewall, as well as those known to occur on the ramp and cowl surfaces. The internal flow grid size was 141 streamwise by 101 cross stream by 71 in the lateral direction between sidewalls. Since the flow entering the inlet was not symmetrical, the inlet was solved from sidewall to sidewall (without using a symmetry plane). In addition to the short sidewalls proposed in the Langley geometry database, a set of shorter sidewalls was also investigated in the present study and was shown to have beneficial effects with respect to the flow distortion exiting the supersonic inlet. In addition to these calculations, additional 3D solutions using the OVERFLOW code were obtained for the flow downstream of the throat of the supersonic inlet, including a terminal shock wave system produced by a backpressured subsonic diffuser.

Rose, William C.↗
DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling the electrochemical behavior and interfacial junction profiles of bipolar membranes at solar flux relevant operating current densities

A 1-dimensional, multi-physics model that accounts for the migration and diffusion of solution species, electrostatics, and chemical reactions, in particular water dissociation (WD), at bipolar membrane (BPM) interfaces was developed to study the electrochemical behavior of bipolar membranes (BPMs) at solar flux relevant operating current densities (tens of mA cm -2 ). Significant partial current densities for WD were observed at BPM voltages much less than the equilibrium voltage, e.g., 59 mV × ΔpH from both experiments and modeling. The co-ion leakage across the BPM at pH differentials accounted for the early presence of the partial current density for WD. Two distinctive electric field dependent WD pathways, the un-catalyzed pathway and the catalyzed pathway, were quantitatively and parametrically studied to improve the turn-on potential of the BPM. The catalyzed pathway accounted for the majority of the partial current density for WD at low voltages, while the un-catalyzed pathway dominated the WD at relatively high voltages. Significant WD was observed only within the interfacial CL (<5 nm), in which a large electric field was present. To improve the electrochemical behavior and the turn-on potential of BPMs, the impacts of the pKa of the immobilized WD catalysts, the electric-field dependent rate constant, the thickness of the catalyst layer and fixed charge density in BPMs on the partial current densities for WD were studied systematically. In addition, the electrochemical behavior and concentration profiles of BPMs in a buffered electrolyte were studied and contrasted with those in an un-buffered electrolyte from both modeling and experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Nonequilibrium radiation and chemistry models for aerocapture vehicle flowfields

The primary accomplishments of the project were as follows: (1) From an overall standpoint, the primary accomplishment of this research was the development of a complete gasdynamic-radiatively coupled nonequilibrium viscous shock layer solution method for axisymmetric blunt bodies. This method can be used for rapid engineering modeling of nonequilibrium re-entry flowfields over a wide range of conditions. (2) Another significant accomplishment was the development of an air radiation model that included local thermodynamic nonequilibrium (LTNE) phenomena. (3) As part of this research, three electron-electronic energy models were developed. The first was a quasi-equilibrium electron (QEE) model which determined an effective free electron temperature and assumed that the electronic states were in equilibrium with the free electrons. The second was a quasi-equilibrium electron-electronic (QEEE) model which computed an effective electron-electronic temperature. The third model was a full electron-electronic (FEE) differential equation model which included convective, collisional, viscous, conductive, vibrational coupling, and chemical effects on electron-electronic energy. (4) Since vibration-dissociation coupling phenomena as well as vibrational thermal nonequilibrium phenomena are important in the nonequilibrium zone behind a shock front, a vibrational energy and vibration-dissociation coupling model was developed and included in the flowfield model. This model was a modified coupled vibrational dissociation vibrational (MCVDV) model and also included electron-vibrational coupling. (5) Another accomplishment of the project was the usage of the developed models to investigate radiative heating. (6) A multi-component diffusion model which properly models the multi-component nature of diffusion in complex gas mixtures such as air, was developed and incorporated into the blunt body model. (7) A model was developed to predict the magnitude and characteristics of the shock wave precursor ahead of vehicles entering the Earth's atmosphere. (8) Since considerable data exists for radiating nonequilibrium flow behind normal shock waves, a normal shock wave version of the blunt body code was developed. (9) By comparing predictions from the models and codes with available normal shock data and the flight data of Fire II, it is believed that the developed flowfield and nonequilibrium radiation models have been essentially validated for engineering applications.

Carlson, Leland A.↗
NASA NTRS1995

Velocity Fields of Axisymmetric Hydrogen-Air Counterflow Diffusion Flames from LDV, PIV, and Numerical Computation

Laminar fuel-air counterflow diffusion flames (CFDFs) were studied using axisymmetric convergent-nozzle and straight-tube opposed jet burners (OJBs). The subject diagnostics were used to probe a systematic set of H2/N2-air CFDFs over wide ranges of fuel input (22 to 100% Ha), and input axial strain rate (130 to 1700 Us) just upstream of the airside edge, for both plug-flow and parabolic input velocity profiles. Laser Doppler Velocimetry (LDV) was applied along the centerline of seeded air flows from a convergent nozzle OJB (7.2 mm i.d.), and Particle Imaging Velocimetry (PIV) was applied on the entire airside of both nozzle and tube OJBs (7 and 5 mm i.d.) to characterize global velocity structure. Data are compared to numerical results from a one-dimensional (1-D) CFDF code based on a stream function solution for a potential flow input boundary condition. Axial strain rate inputs at the airside edge of nozzle-OJB flows, using LDV and PIV, were consistent with 1-D impingement theory, and supported earlier diagnostic studies. The LDV results also characterized a heat-release hump. Radial strain rates in the flame substantially exceeded 1-D numerical predictions. Whereas the 1-D model closely predicted the max I min axial velocity ratio in the hot layer, it overpredicted its thickness. The results also support previously measured effects of plug-flow and parabolic input strain rates on CFDF extinction limits. Finally, the submillimeter-scale LDV and PIV diagnostics were tested under severe conditions, which reinforced their use with subcentimeter OJB tools to assess effects of aerodynamic strain, and fueVair composition, on laminar CFDF properties, including extinction.

Pellett, Gerald L.↗
DOE OSTI2025

Multi-X-point radiation and its dynamics in up/down asymmetry in plasma detachment regimes in Wendelstein 7-X

In Wendelstein 7-X (W7-X), stable plasma detachment in ECRH hydrogen plasmas is routinely achieved with the standard magnetic configuration, containing five magnetic islands in the scrape-off layer (SOL). Plasma detachment induced by intrinsic carbon impurities is characterized by a high radiation fraction (ƒ rad ) with significantly reduced divertor heat load and particle flux. In this work, we present the radiation dynamics during the detached plasma (DP) phase, such as the shift of radiation zones in radial and poloidal directions towards the X-points as ƒ rad increases (∼0.6–0.9), as well as the penetration of the radiation layer into the confinement region in the deep DP phase with ƒ rad approaching unity. In particular, a structure of multi-X-point radiation (multi-XPR) with an up/down asymmetry in the DP phase is highlighted, which is revealed by bolometer tomography and is further confirmed by video diagnostics. The multi-XPR structure forms helical 3D bands aligned with W7-X field periodicity. Field reversal experiments show that the brightest XPR displaces between the upper and lower SOL regions, suggesting the E× B drift effects. This paper presents the multi-XPR structure observed in the W7-X plasma in detail for the first time. A simplified model considering the influence of the poloidal E× B drift ( V d ) on the impurity flow in the SOL shows that the downstream drift toward the target or target shadow region ( V d < 0) decreases the impurity content, while the upstream drift toward the LCFS ( V d > 0) increases the impurity content. On this basis, the poloidal drift potentially leads to an up/down asymmetry of impurity density in the SOL despite the symmetry magnetic topology. The dynamics of the up/down asymmetry in the multi-XPR structure is also related to the magnitude V d /D (normalized to the impurity diffusivity), with an additional effect owing to the radial inward shift of the emission zone. These results provide new insights into impurity-induced detachment dynamics and provide a basis for improving 3D modeling of impurity transport, which typically does not consider drifts.

E × B drift↗
DOE OSTI2025

Continuum shock mixture models for Ni+Al multilayers: Individual layers and bulk equations of state

Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni x Al y ⁠, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni x Al y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Second-order closure PBL model with new third-order moments: Comparison with LES data

This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.

Canuto, V. M.↗
NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗
NASA NTRS1995

Time-dependent simulation of oblique MHD cosmic-ray shocks using the two-fluid model

Using a new, second-order accurate numerical method we present dynamical simulations of oblique MHD cosmic-ray (CR)-modified plane shock evolution. Most of the calculations are done with a two-fluid model for diffusive shock acceleration, but we provide also comparisons between a typical shock computed that way against calculations carried out using the more complete, momentum-dependent, diffusion-advection equation. We also illustrate a test showing that these simulations evolve to dynamical equilibria consistent with previously published steady state analytic calculations for such shocks. In order to improve understanding of the dynamical role of magnetic fields in shocks modified by CR pressure we have explored for time asymptotic states the parameter space of upstream fast mode Mach number, M(sub f), and plasma beta. We compile the results into maps of dynamical steady state CR acceleration efficiency, epsilon(sub c). We have run simulations using constant, and nonisotropic, obliquity (and hence spatially) dependent forms of the diffusion coefficient kappa. Comparison of the results shows that while the final steady states achieved are the same in each case, the history of CR-MHD shocks can be strongly modified by variations in kappa and, therefore, in the acceleration timescale. Also, the coupling of CR and MHD in low beta, oblique shocks substantially influences the transient density spike that forms in strongly CR-modified shocks. We find that inside the density spike a MHD slow mode wave can be generated that eventually steepens into a shock. A strong layer develops within the density spike, driven by MHD stresses. We conjecture that currents in the shear layer could, in nonplanar flows, results in enhanced particle accretion through drift acceleration.

Frank, Adam↗