Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sim”

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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Exploration of Sn 70 Li 30 alloy as possible material for flowing liquid metal plasma facing components

As an advanced alternative to solid materials, Liquid Metals (LM) may offer more resilient and feasible Plasma Facing Components (PFCs). Particularly, regarding the unavoidable material erosion/degradation produced by particle/heat fluxes in future fusion devices where much longer duty cycles are expected. Furthermore, configurations that propose a flowing LM surface can add the advantage of a continuously fresh and clean layer facing the plasma. Although lithium is the most widely tested option, tin-lithium (SnLi) alloys have been proposed to attempt to combine the positive characteristics of both pure elements and ameliorate the specific issues of lithium. In this work, the potential use of Sn70Li30 alloy in such flowing concepts has been explored by addressing several preliminary and mandatory aspects for its utilization. Key issues such as wettability and compatibility of the alloy with relevant substrates have been studied in a multidisciplinary approach. The data obtained from deposited liquid tin-lithium droplets indicates approximate wetting temperatures of 360 °C, 390 °C and 405 °C for the fresh alloy on 316 stainless steel, molybdenum, and tungsten, respectively. However, the alloy contamination appeared to strongly affect the wetting characteristics of materials, increasing their wetting temperature by ~130 °C in the worst observed cases. Interestingly, in some instances, the instability of the liquid alloy surface was observed in the form of sudden gaseous ejection. The deposited droplets were posteriorly characterized in terms of absolute composition and depth profile by Inductively Coupled PlasmaOptical Emission Spectroscopy (ICP-OES) and Secondary Ion Mass Spectrometry (SIMS-ToF). Additionally, the nature and composition of the boundaries between the substrates and alloy microparticles was investigated by Scanning Electron Microscopy (SEM), Energy Dispersive X-Ray Spectroscopy (EDS), and 3D Laser microscopy. The overall results of this post-mortem characterization revealed that first signs of corrosion induced by both alloy elements (lithium-chromium association and iron-tin intermetallic mixing) were present on 316 stainless steel after short exposures (≤3 h) at temperatures lower than 550 °C. Conversely, molybdenum and tungsten showed good compatibility with the alloy in equivalent conditions. The global implications of these results are finally addressed, focusing on the future perspectives and the more viable scenarios for the eventual utilization of these alloys in flowing liquid metal configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Testing of spark plasma sintered porous tungsten under neon glow discharge cleaning conditions in LTX-β

This paper demonstrates that tungsten (W) based powder reconstituted Plasma Facing Components (PFCs) can be treated in situ in a fusion reactor to remove W oxide and carbon (C) contamination. Doing so should ease the challenge of using these materials in a Capillary Porous System (CPS) with lithium (Li) by enabling better wetting and less contamination of the Li by the underlying CPS. Most powder reconstituted materials including 3D printed and sintered PFCs suffer from a high surface contamination from oxides and surface C, which complicates their use with liquid Li, a primary PFC candidate. Spark plasma sintered porous W samples were fabricated to be used as a CPS with liquid Li. The samples were characterized in terms of morphology and surface chemistry. Analysis confirms a high C and oxygen (O) contamination. We present the results of exposing this type of CPS to Glow Discharge Cleaning (GDC) cycles in the Lithium Tokamak Experiment-β (LTX-β). The sample was exposed to neon (Ne) GDC in the midplane of the low-field side of LTX-β and analyzed in vacuo with Temperature Programmed Desorption (TPD) and Secondary Ion Mass Spectrometry (SIMS) to investigate the effects the Ne GDC had on the chemical composition of the sample. The combination of Ne GDC with rapid heating as done in TPD was successful in reducing the W oxides and removing the C contamination.

Capillary porous system↗
DOE OSTI2025

Natural top quark condensation (a redux)

The Nambu--Jona-Lasinio (NJL) model involves a pointlike 4-fermion interaction. While it gives a useful description of chiral dynamics (mainly in QCD), it nonetheless omits the crucially important internal wave-function of a two-body bound state, $\phi(r)$. This becomes significant near critical coupling where $\phi(r)$ extends to large distance, leading to dilution and suppression of induced couplings $\propto \phi(0)$, such as the Yukawa and quartic couplings, as well as reduced fine-tuning of a hierarchy. In top quark condensation, where the BEH boson is a $\bar{t}t$ bound state and we have a UV completion such as topcolor, we must go beyond the NJL model and include effects of $\phi(r)$. We provide a formulation of this for the Brout-Englert-Higgs boson, and find that it leads to an extended $\phi(r)$, a significantly reduced and natural composite scale of $M_0 \sim 6$ TeV, a successful prediction for the quartic coupling, $\lambda$, and fine tuning that is reduced to a few percent, providing a compelling candidate solution to the naturalness problem of the BEH boson. The theory is testable and new physics should begin to emerge on the multi-TeV mass scales and possibly accessible to the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Assessing the fidelity of marine vertebrate microfossil δ18O signatures and their potential for palaeo-ecological and -climatic reconstructions

Conodont biogenic apatite has become a preferred analytical target for oxygen isotope studies investigating ocean temperature and palaeoclimate change in the Palaeozoic. Despite the growing application in geochemical based palaeoenvironmental reconstructions, the paucity or absence of conodont fossils in certain facies necessitates greater flexibility in selection of robust oxygen bearing compounds for analysis. Microvertebrates offer a potential substitute for conodonts from the middle Palaeozoic. Microvertebrate bioapatite is particularly advantageous given a fossil record extending to the present with representatives across freshwater to fully marine environments, thus widening the scope of oxygen isotope studies on bioapatite. However, significant tissue heterogeneity within vertebrates and differential susceptibility of these tissues to diagenetic alteration have been raised as potential problems affecting the reliability of the oxygen isotope ratios as palaeoclimate proxies. Pristine microvertebrate and co-occurring conodont fossils from the Late Devonian and Early Carboniferous of the Lennard Shelf, Canning Basin, Western Australia, were analysed using bulk (gas isotope ratio mass spectrometry) and in-situ (secondary ion mass spectrometry) methodologies, with the latter technique allowing investigation of specific tissues within vertebrate elements. The δ18Oconodont results may be interpreted in terms of palaeolatitudinally and environmentally sensible palaeotemperatures and provide a baseline standard for comparison against δ18Omicrovertebrate values. Despite an absence of obvious diagenetic influences, GIRMS of microvertebrate denticles yielded δ18O values depleted by 2-4 ‰ relative to co-occurring conodonts. SIMS analysis of hypermineralised tissues in both scales and teeth produced δ18O values comparable with those of associated conodonts. The susceptibility of porous phosphatic fossil tissues to microbial activity, fluid interaction and introduction of mineral precipitates post-formation is demonstrated in microvertebrate dentine, which showed significant heterogeneity and consistent depletion in 18O. The hypermineralised tissues present in both teeth and scales appear resistant to many diagenetic processes and indicate potential for palaeoclimatic reconstructions and palaeoecological investigations.

conodont↗
DOE OSTI2022

Hydrodynamics and Load Shedding Behavior of a Variable-Geometry Oscillating Surge Wave Energy Converter (OSWEC)

In order to improve their long-term viability, wave energy converters (WECs) need to be able to shed loads when a threshold wave condition is exceeded. As shown by Tom et al. (2016) [1], provision of adjustable flaps within the body of an oscillating surge wave energy converter (OSWEC) allows wave energy to pass through the device. A control system may then be able to open and close the flaps when waves approaching the device exceed preset thresholds. The variable-geometry OSWEC (VG-OSWEC) concept studied in this paper is a bottom-hinged, rectangular wave paddle with five flaps of elliptical cross-section embedded into the face of the paddle. System ID tests were conducted on this VG-OSWEC device at a 1:14 scale in a wave basin. Free decay tests showed that the damping was distinctly nonlinear when the flaps were fully open, and the natural frequency increased by 40% when compared with the flaps in a fully closed configuration. Tests with regular wave conditions were used to develop the response amplitude operator for the rotational motion about the hinge. These response amplitude operator results when compared with numerical simulations run using WEC-Sim/WAMIT and ANSYS AQWA, show strong agreement with the flap open and closed conditions. The regular-wave condition measurements also show that the wave excitation moment about the hinge was reduced by up to 60% when the flaps were fully open. The experiments serve to demonstrate the potential of the variable geometry design to shed loads and survive harsh ocean environments.

load shedding↗
DOE OSTI2025

On the optimal sizing of power take-off systems for wave energy converters

The power take-off (PTO) system in a wave energy converter (WEC) is the means by which the energy in an ocean wave is converted into useful energy. There has been a general lack of technology convergence in the wave energy industry and PTO design has been largely device-specific with the design process determined by the WEC developer. An optimized PTO design can significantly affect the efficiency, reliability, performance, and overall cost and viability of a device. There remains a need in the marine energy industry for best practices and efficient design processes to be identified as it could greatly benefit developers and researchers in designing an optimal WEC PTO. This research utilizes the open-source Matlab-Simulink based software, WEC-Sim, which was developed by the National Renewable Energy Laboratory, and the Eagle, a high-performance computing system, to model and optimize two WECs of different archetypes and then investigates if any trends exist that can be exploited for greater efficiency in optimal PTO design. The results from this research indicate that there may be an optimal sizing for WEC PTOs for both power and PTO force rating and that substantial reductions can be made in the cost of a device without significant loss in the energy produced. Furthermore this optimal sizing may be independent of WEC type and deployment location. With an optimal, or near optimal, approach to WEC PTO sizing demonstrated, a methodology is proposed to address the challenge of nameplate ratings in the wave energy industry.

13 HYDRO ENERGY↗
DOE OSTI2016

Low-cost plasma immersion ion implantation doping for Interdigitated back passivated contact (IBPC) solar cells

Here, we present progress to develop low-cost interdigitated back contact solar cells with pc-Si/SiO 2 /c-Si passivated contacts formed by plasma immersion ion implantation (PIII). PIII is a lower-cost implantation technique than traditional beam line implantation due to its simpler design, lower operating costs, and ability to run high doses (1E14-1E18 cm -2 ) at low ion energies (20 eV-10 keV). These benefits make PIII ideal for high throughput production of patterned passivated contacts, where high-dose, low-energy implantations are made into thin (20-200 nm) a-Si layers instead of into the wafer itself. For this work symmetric passivated contact test structures (~100 nm thick) grown on n-Cz wafers with pH3 PIII doping gave implied open circuit voltage (iV oc ) values of 730 mV with J o values of 2 fA/cm 2 . Samples doped with B 2 H 6 gave iV oc values of 690 mV and J o values of 24 fA/cm 2 , outperforming BF 3 doping, which gave iV oc values in the 660-680 mV range. Samples were further characterized by SIMS, photoluminescence, TEM, EELS, and post-metallization TLM to reveal micro- and macro-scopic structural, chemical and electrical information.

14 SOLAR ENERGY↗
DOE OSTI2023

Screen printable copper pastes for silicon solar cells

The current work demonstrates the successful metallization of a PERC silicon solar cell with screen-printable copper (Cu) paste that is sintered at elevated temperature in air atmosphere. The existing state of the art in Silicon (Si) solar cell metallization is silver (Ag) paste; Cu cost is 1/100th the cost of Ag and has a comparable conductivity. However, Cu undergoes rapid oxidation at elevated temperatures and the high diffusion of Cu into Si restricts its usage in the metallization of silicon solar cells. In this paper, a Cu paste containing a proprietary mixture of antioxidant additives and diffusion inhibitors was used to make front gridlines on PERC cells. The Cu-printed cells were fired in an IR belt furnace with a peak temperature of 576°C in air atmosphere. The solar cells discussed here reached an efficiency of 19%, which is 88% of the efficiency measured for a commercial Ag-printed cell. The working devices were further characterized using STEM, EDX, ToF-SIMS and long-term suns-V oc studies demonstrating an absence of Cu diffusion. Furthermore, thermally stressing the Cu-printed solar cells yielded a reduction of less than 1.9% of the initial pseudo fill factor (pFF) after 1000 hours at 200°C. Here, the shunt resistance of cells monitored over 3 years remained within 2.5% of its initial value, which demonstrates the effectiveness of the Cu diffusion barrier.

14 SOLAR ENERGY↗
DOE OSTI2025

Thermodynamics, local structure, and transport of protons in triple-conducing oxide, BaCo 0.4 Fe 0.4 Zr 0.1 Y 0.1 O 3-δ (BCFZY4411)

Triple-conducting oxides (TCOs) are an emerging class of mixed ionic and electronically conducting materials that show great promise for oxygen reduction/oxygen evolution (ORR/OER) electrocatalysis—primarily in high-temperature ceramic electrochemical cells— but also in aqueous alkaline environments. Their high activity is attributed, at least in part, to their ability to incorporate and transport three mobile charge carriers: protons, oxygen vacancies, and electron-holes Despite their promise, fundamental studies of TCOs are challenging, as transport dynamics from three charge carriers cannot be fully disentangled via traditional electrical measurement techniques. Characterizing proton dynamics in TCOs is particularly difficult as protons are generally the minority carrier, and their conduction response is typically obscured by the oxygen vacancies and electron holes. Here, we demonstrate successful isolation of the proton behavior in an archetypal TCO, BaCo 0.4 Fe 0.4 Zr 0.1 Y 0.1 O 3-δ (BCFZY4411), using a combination of non-electrical techniques. We determine proton uptake and oxygen non-stoichiometry (δ) using thermogravimetric analysis (TGA). X-ray absorption near edge structure (XANES) and neutron diffraction (ND) are used to validate the oxidation state of Co and the δ values obtained through TGA. We apply 1H solid-state magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) to provide insights into local structure, dynamics, and proton kinetics. Finally, the proton transport properties are further quantified using tracer isotope exchange with time-of-flight secondary ion mass spectrometry (ToF-SIMS). Despite the very low proton concentrations in BCFZY4411 (<0.2% under most conditions), our analysis suggests that the oxygen Manuscript File Click here to view linked References 2 reduction and evolution reactions are nevertheless limited by the oxygen ion kinetics (e.g., oxygen surface exchange) rather than the proton kinetics at the reduced operating temperatures (<500 °C) that are targeted for electrochemical cell applications. These findings provide a comprehensive understanding of proton behavior in BCFZY4411 and pave the way for advancing the fundamental study of TCOs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Mitochondrial-targeted aryl hydrocarbon receptor and the impact of 2,3,7,8-tetrachlorodibenzo-p-dioxin on cellular respiration and the mitochondrial proteome

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor within the Per-Arnt-Sim (PAS) domain superfamily. Exposure to the most potent AHR ligand, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), is associated with various pathological effects including metabolic syndrome. While research over the last several years has demonstrated a role for oxidative stress and metabolic dysfunction in AHR-dependent TCDD-induced toxicity, the role of the mitochondria in this process has not been fully explored. Our previous research suggested that a portion of the cellular pool of AHR could be found in the mitochondria (mitoAHR). Using a protease protection assay with digitonin extraction, we have now shown that this mitoAHR is localized to the inter-membrane space (IMS) of the organelle. TCDD exposure induced a degradation of mitoAHR similar to that of cytosolic AHR. Furthermore, siRNA-mediated knockdown revealed that translocase of outer-mitochondrial membrane 20 (TOMM20) was involved in the import of AHR into the mitochondria. In addition, TCDD altered cellular respiration in an AHR-dependent manner to maintain respiratory efficiency as measured by oxygen consumption rate (OCR). Stable isotope labeling by amino acids in cell culture (SILAC) identified a battery of proteins within the mitochondrial proteome influenced by TCDD in an AHR-dependent manner. Among these, 17 proteins with fold changes ≥ 2 are associated with various metabolic pathways, suggesting a role of mitochondrial retrograde signaling in TCDD-mediated pathologies. Collectively, these studies suggest that mitoAHR is localized to the IMS and AHR-dependent TCDD-induced toxicity, including metabolic dysfunction, wasting syndrome, and hepatic steatosis, involves mitochondrial dysfunction. - Highlights: • The mitoAHR is localized in the mitochondrial intermembrane space. • TOMM20 participates in mitoAHR translocation. • AHR contributes to the maintenance of respiratory control ratio following TCDD exposure. • TCDD-induced AHR-dependent changes in the mitochondrial proteome are identified.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

Scaling diode-pumped, high energy picosecond lasers to kilowatt average powers

Recent results in the development of diode-driven high energy, high repetition rate, picosecond lasers, including the demonstration of a cryogenic Yb:YAG active mirror amplifier that produces 1.5 J pulses at 500 Hz repetition rate (0.75 kW average power) are reviewed. These pulses are compressed resulting in the generation of ${\sim}5~\text{ps}$ duration, 1 J pulses with 0.5 kW average power. A full characterization of this high power cryogenic amplifier, including at-wavelength interferometry of the active region under ${>}1~\text{kW}$ average power pump conditions, is presented. An initial demonstration of operation at 1 kW average power (1 J, 1 kHz) is reported.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Experimental platform for the investigation of magnetized-reverse-shock dynamics in the context of POLAR

The influence of a strong external magnetic field on the collimation of a high Mach number plasma flow and its collision with a solid obstacle is investigated experimentally and numerically. The laser irradiation ( $I\sim 2\times 10^{14}~\text{W}\cdot \text{cm}^{-2}$ ) of a multilayer target generates a shock wave that produces a rear side plasma expanding flow. Immersed in a homogeneous 10 T external magnetic field, this plasma flow propagates in vacuum and impacts an obstacle located a few mm from the main target. A reverse shock is then formed with typical velocities of the order of 15–20 $\pm$ 5 km/s. The experimental results are compared with 2D radiative magnetohydrodynamic simulations using the FLASH code. This platform allows investigating the dynamics of reverse shock, mimicking the processes occurring in a cataclysmic variable of polar type.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Frontogenesis and frontal arrest of a dense filament in the oceanic surface boundary layer

The evolution of upper ocean currents involves a set of complex, poorly understood interactions between submesoscale turbulence (e.g. density fronts and filaments and coherent vortices) and smaller-scale boundary-layer turbulence. Here we simulate the lifecycle of a cold (dense) filament undergoing frontogenesis in the presence of turbulence generated by surface stress and/or buoyancy loss. This phenomenon is examined in large-eddy simulations with resolved turbulent motions in large horizontal domains using ${\sim}10^{10}$ grid points. Steady winds are oriented in directions perpendicular or parallel to the filament axis. Due to turbulent vertical momentum mixing, cold filaments generate a potent two-celled secondary circulation in the cross-filament plane that is frontogenetic, sharpens the cross-filament buoyancy and horizontal velocity gradients and blocks Ekman buoyancy flux across the cold filament core towards the warm filament edge. Within less than a day, the frontogenesis is arrested at a small width, ${\approx}100~\text{m}$ , primarily by an enhancement of the turbulence through a small submesoscale, horizontal shear instability of the sharpened filament, followed by a subsequent slow decay of the filament by further turbulent mixing. The boundary-layer turbulence is inhomogeneous and non-stationary in relation to the evolving submesoscale currents and density stratification. The occurrence of frontogenesis and arrest are qualitatively similar with varying stress direction or with convective cooling, but the detailed evolution and flow structure differ among the cases. Thus submesoscale filament frontogenesis caused by boundary-layer turbulence, frontal arrest by frontal instability and frontal decay by forward energy cascade, and turbulent mixing are generic processes in the upper ocean.

Mechanics↗
DOE OSTI2018

Inclined porous medium convection at large Rayleigh number

High-Rayleigh-number ( $Ra$ ) convection in an inclined two-dimensional porous layer is investigated using direct numerical simulations (DNS) and stability and variational upper-bound analyses. When the inclination angle $\unicode[STIX]{x1D719}$ of the layer satisfies $0^{\circ }<\unicode[STIX]{x1D719}\lesssim 25^{\circ }$ , DNS confirm that the flow exhibits a three-region wall-normal asymptotic structure in accord with the strictly horizontal ( $\unicode[STIX]{x1D719}=0^{\circ }$ ) case, except that as $\unicode[STIX]{x1D719}$ is increased the time-mean spacing between neighbouring interior plumes also increases substantially. Both DNS and upper-bound analysis indicate that the heat transport enhancement factor (i.e. the Nusselt number) $Nu\sim CRa$ with a $\unicode[STIX]{x1D719}$ -dependent prefactor $C$ . When $\unicode[STIX]{x1D719}>\unicode[STIX]{x1D719}_{t}$ , however, where $30^{\circ }<\unicode[STIX]{x1D719}_{t}<32^{\circ }$ independently of $Ra$ , the columnar flow structure is completely broken down: the flow transitions to a large-scale travelling-wave convective roll state, and the heat transport is significantly reduced. To better understand the physics of inclined porous medium convection at large $Ra$ and modest inclination angles, a spatial Floquet analysis is performed, yielding predictions of the linear stability of numerically computed, fully nonlinear steady convective states. The results show that there exist two types of instability when $\unicode[STIX]{x1D719}\neq 0^{\circ }$ : a bulk-mode instability and a wall-mode instability, consistent with previous findings for $\unicode[STIX]{x1D719}=0^{\circ }$ (Wen et al. , J. Fluid Mech. , vol. 772, 2015, pp. 197–224). The background flow induced by the inclination of the layer intensifies the bulk-mode instability during its subsequent nonlinear evolution, thereby favouring increased spacing between the interior plumes relative to that observed in convection in a horizontal porous layer.

Mechanics↗
DOE OSTI2018

Simultaneous direct measurements of concentration and velocity in the Richtmyer–Meshkov instability

The Richtmyer–Meshkov instability (RMI) is experimentally investigated in a vertical shock tube using a broadband initial condition imposed on an interface between a helium–acetone mixture and argon (Atwood number $A\approx 0.7$ ). In the present work, a shear layer is introduced at the interface to serve as a statistically repeatable, broadband initial condition to the RMI, and the density interface is accelerated by either an $M=1.6$ or $M=2.2$ planar shock wave. The development of the ensuing mixing layer is investigated using simultaneous planar laser-induced fluorescence (PLIF) and particle image velocimetry (PIV). PLIF images are processed to reveal the light-gas mole fraction, while PIV particle image pairs yield corresponding two-component planar velocity results. Field structure and distribution are explored through probability density functions (PDFs), and a decomposition is performed on concentration and velocity results to obtain a mean flow field and define fluctuations. Simultaneous concentration and velocity field measurements allow – for the first time in this regime – experimentally determined turbulence quantities such as Reynolds stresses, turbulent mass-flux velocities and turbulent kinetic energy to be obtained. We show that by the latest times the mixing layer has passed the turbulent threshold, and there is evidence of turbulent mixing occurring sooner for the higher Mach number case. Interface measurements show nonlinear growth with a power-law fit to the thickness data, and that integral measurements of mixing layer thickness are proportional to threshold measurements. Spectral analysis demonstrates the emergence of an inertial range with a slope ${\sim}k^{-5/3}$ when considering both density and velocity effects in planar turbulent kinetic energy (TKE) measurements.

Mechanics↗
DOE OSTI2019

Nanoparticle diffusion in sheared cellular blood flow

Using a multiscale blood flow solver, the complete diffusion tensor of nanoparticles (NPs) in sheared cellular blood flow is calculated over a wide range of shear rate and haematocrit. In the short-time regime, NPs exhibit anomalous dispersive behaviors under high shear and high haematocrit due to the transient elongation and alignment of the red blood cells (RBCs). In the long-time regime, the NP diffusion tensor features high anisotropy. Particularly, there exists a critical shear rate ( ${\sim}100~\text{s}^{-1}$ ) around which the shear-rate dependence of the diffusivity tensor changes from linear to nonlinear scale. Above the critical shear rate, the cross-stream diffusivity terms vary sublinearly with shear rate, while the longitudinal term varies superlinearly. The dependence on haematocrit is linear in general except at high shear rates, where a sublinear scale is found for the vorticity term and a quadratic scale for the longitudinal term. Through analysis of the suspension microstructure and numerical experiments, the nonlinear haemorheological dependence of the NP diffusion tensor is attributed to the streamwise elongation and cross-stream contraction of RBCs under high shear, quantified by a capillary number. The RBC size is shown to be the characteristic length scale affecting the RBC-enhanced shear-induced diffusion (RESID), while the NP submicrometre size exhibits negligible influence on the RESID. Based on the observed scaling behaviours, empirical correlations are proposed to bridge the NP diffusion tensor to specific shear rate and haematocrit. As a result, the characterized NP diffusion tensor provides a constitutive relation that can lead to more effective continuum models to tackle large-scale NP biotransport applications.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Langmuir turbulence and filament frontogenesis in the oceanic surface boundary layer

Submesoscale currents, small-scale turbulence and surface gravity waves co-exist in the upper ocean and interact in complex ways. To expose the couplings, the frontogenetic life cycle of an idealized cold dense submesoscale filament interacting with upper ocean Langmuir turbulence is investigated in large-eddy simulations (LESs) based on the incompressible wave-averaged equations. The simulations utilize large domains and fine meshes with $6.4\times 10^{9}$ grid points. Case studies are made with surface winds or surface cooling with waves oriented in across-filament (perpendicular) or down-filament (parallel) directions relative to the two-dimensional filament axis. The currents $u$ , $v$ and $w$ are aligned with the across-filament, down-filament and vertical directions, respectively. Frontogenesis is induced by across-filament Lagrangian secondary circulations in the boundary layer, and it is shown to be strongly impacted by surface waves, in particular the propagation direction relative to the filament axis. In a horizontally heterogeneous boundary layer, surface waves induce both mean and fluctuating Stokes-drift vortex forces that modify a linear, hydrostatic turbulent thermal wind (TTW) approximation for momentum. Down-filament winds and waves are found to be especially impactful, they significantly reduce the peak level of frontogenesis by fragmenting the filament into primary and secondary down-welling sites in a broad frontal zone over a width ${\sim}500~\text{m}$ . At peak frontogenesis, opposing down-filament jets $\langle v\rangle$ overlie each other resulting in a vigorous vertical shear layer $\unicode[STIX]{x2202}_{z}\langle v\rangle$ with large vertical momentum flux $\langle v^{\prime }w^{\prime }\rangle$ . Filament arrest is induced by a lateral shear instability that generates horizontal momentum flux $\langle u^{\prime }v^{\prime }\rangle$ at low wavenumbers. The turbulent vertical velocity patterns, indicative of coherent Langmuir cells, change markedly across the horizontal domain with both across-filament and down-filament winds under the action of submesoscale currents.

Mechanics↗
DOE OSTI2022

Self-similarity of scalar isosurface area density in a temporal mixing layer

Understanding the behaviour of scalar isosurfaces in a turbulent flow is of particular interest for many problems in turbulent mixing that contain sharp interfaces between regions of the flow. Common examples include combustion, where the chemical reactions occur in thin regions within the flow, and the turbulent/non-turbulent interface in shear flows, where a thin region separates the rotational, turbulent motions from the irrotational, non-turbulent background. Recent advances in computing technology allow for in-depth analysis of these interface problems that are difficult to quantify in a laboratory setting. In this paper, the results of a direct numerical simulation of a passive scalar $\varPhi$ evolving on a turbulent, temporally developing mixing layer are described. A novel approach has been taken to calculate the surface area of individual scalar isosurfaces, $A_{iso}$ , throughout the simulation, as well as the mean isosurface area density, $\varSigma$ , as a function of the cross-stream direction and time. A notable finding is that the profiles of $\varSigma$ develop in a self-similar manner when scaled by the Taylor scale of the scalar field, $\lambda _\phi$ . Remarkably, the scaling appears to hold for a wide range of isovalues. A rough scaling argument based on the formal definition of $\varSigma$ and properties of a temporal mixing layer is presented which also exposes a dependence on $\lambda _\phi$ . Based on these results, a possible scaling for the isosurface area is presented as $A_{iso}/A_0 \sim (Re \, Sc)^{1/2}$ , where $Re$ and $Sc$ are local Reynolds and Schmidt numbers, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗