Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shear Layers”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Low Mach number fluctuating hydrodynamics model for ionic liquids

We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

42 ENGINEERING↗
DOE OSTI2020

Evaluating turbulence models at high energy densities

The RESHOCK campaign at Lawrence Livermore National Laboratory has been working on increasing our understanding of the evolution of turbulent, unstable plasma interfaces which are applicable in various HED applications including ICF implosions. A common approach to model such interfaces involves the use of Renolds-Averaged-Navier-Stokes (RANS) models whose parameters have been constrained by theory and experiments. These parameter sets are not unique and the applicability of one set of tuned model parameters to systems with differing shock strengths or material densities has not been well tested. Our new NIF experiments, along with related studies of turbulent shear flows, provide data and model validation for plasma interfaces in the HED regime. Our experiments measure mixing-layer width at an unstable interface in the plasma regime and are specifically designed to challenge mix models of the Reynolds-Averaged-Navier-Stokes (RANS) type where turbulence is assumed to be fully developed, i.e. the flow exhibits a broad spectrum of length scales without memory of the initial condition.4 We utilize precise control of the initial interface conditions (densities on both sides of the interface, an initial rippled perturbation pattern consisting of equally weighted wavelengths between 10 and 20 microns, and initial material compositions) along with a repeatable drive history

36 MATERIALS SCIENCE↗
DOE OSTI2025

Direct ink writing of shear exfoliated two-dimensional nanomaterial- elastomeric multifunctional nanocomposite

Direct ink writing (DIW) of polymer nanocomposites with high loadings of two-dimensional (2D) nanofillers (graphene and hexagonal boron nitride (hBN)) is challenging because of potential clogging, use of hazardous solvents, and agglomeration. Here, in this work, a shear exfoliation and sieving method to prepare DIW ink with high loading of nanofillers produced from low-cost bulk layered materials such as graphite and bulk hBN powder for successful DIW printing without the use of any solvents, binders, or plasticizers. The single-step exfoliation technique resulted in a composite with substantial layer reduction along the c-axis, as confirmed by SEM, TEM, XRD, and Raman analysis. Incorporating exfoliated graphene (40 wt%) increased viscosity by ∼6 orders of magnitude due to enhanced particle–matrix interactions, leading to pronounced yield stress behavior and a yield stress of approximately 1598 Pa, which enabled excellent shape retention during extrusion. Using the DIW technique, porous structures such as desalination membranes, self-sensing bone scaffolds, thermal management coating, and serpentine strain sensors were fabricated. When tested in a direct contact membrane distillation setup, the fabricated membrane demonstrated a promising permeate flux of 21.85 Lm −2 h −1 and a salt rejection of 74.3 %. The fabricated serpentine sensor exhibited stable signal variations under cyclic tensile loading, with a working range of 0–200 % strain and a maximum gauge factor of 43,735. A cell culture test using the printed bone scaffold demonstrated promising cell attachment and proliferation. The DIW printed hBN nanocomposite exhibited reversible shape change under heat, demonstrating potential 4D printing capability and efficient thermal management when exposed to high heat or flame.

Desalination↗
DOE OSTI2021

Fiber/matrix debonding evaluation of SiCf/SiC composites using micropillar compression technique

Third-generation silicon carbide (SiC) composites reinforced by SiC fibers (Hi-Nicalon S [HNS] and Tyranno SA3 [SA3]) are attractive for use in next generation reactors owing to their high strength and chemical inertness at high temperatures, as well as enhanced radiation tolerance under neutron irradiation environments. To optimize composite performance, the interfacial mechanical properties of chemical vapor–infiltrated (CVI) SiC f /SiC composites are investigated in this effort by using a slant interface micropillar compression testing procedure. The micropillar test specimens, containing an inclined pyrolytic carbon (PyC) interphase, are prepared using a focused ion beam. The novel microcompression testing successfully quantifies the debond shear strength and internal friction coefficient of micropillar test samples by using the Mohr-Coulomb formulation. Furthermore, according to four types of SiCf/SiC composite microcompression test results, interfacial properties and debond mechanisms are significantly affected by the PyC layer thickness, the local bonding mechanism of the PyC interphase on the SiC fiber surface, and the surface roughness of fibers. Regardless of PyC thicknesses, SA3-reinforced CVI SiC f /SiC composites are found to have much higher debond shear strengths than HNS-reinforced SiC f /SiC CVI composites. By using this micropillar compression technique alongside analytical methods, we uncover new understandings of PyC interface properties. Additionally, the micropillar test results obtained are correlated with macroscopic mechanical properties of neutron-irradiated CVI SiCf/SiC composites.

42 ENGINEERING↗
DOE OSTI2023

Dynamics of multilayer Rayleigh–Taylor instability at moderately high Atwood numbers

This paper investigates the multilayer Rayleigh–Taylor instability (RTI) using statistically stationary experiments conducted in a gas tunnel. Employing diagnostics such as particle image velocimetry (PIV) and planar laser induced fluorescence (PLIF), we make simultaneous velocity–density measurements to study how dynamics and mixing are linked in this variable density flow. Experiments are conducted in a newly built, blow-down three-layer gas tunnel facility. Mixing between three gas streams is studied, where the top and bottom streams are comprised of air, and the middle stream is an air–helium mixture. Shear is minimized between these streams by matching their inlet velocities. The four experimental conditions investigated here consist of two different density ratios (Atwood numbers 0.3 and 0.6), each investigated at two instability development times (or equivalently, two streamwise locations), and all experiments are with the same middle stream thickness of 3 cm. The growth of the middle layer is measured using laser-based planar Mie scattering visualization. The mixing width is found to grow linearly with time at late times. Various quantitative measures of molecular mixing indicate a very high degree of molecular mixing at late times in the multilayer RTI flow. The vertical turbulent mass flux $a_y$ is calculated. In addition to mostly negative values of $a_y$ , typical of buoyancy-dominated flows due to negative correlation between velocity and density fluctuations, positive regions are also observed in profiles of $a_y$ due to entrainment and erosion at the lower edge of the mixing region. Global energy budgets are calculated for the multilayer RTI flow at late times and it is found that the majority of potential energy released has been dissipated due to viscous effects, and a large value of mixing efficiency ( $\sim$ 60 %) is observed.

Mechanics↗
DOE OSTI2020

Soliton-Dependent Electronic Transport across Bilayer Graphene Domain Wall

Layer-stacking domain wall in bilayer graphene is one type of topological defects that can greatly affect the electronic properties of bilayer graphene and therefore lead to nontrivial transport behaviors. An outstanding question on the layer stacking domain wall is how the electrons hop between two adjacent stacking domains. Here we report the first experimental observation of electronic transport across bilayer graphene domain walls by combining near-field infrared nanoscopy and scanning voltage microscopy techniques. In this work, we observe markedly different electron transport behaviors across the tensile- and shear-type domain walls. The tensile-type domain wall is highly reflective of low-energy incident electrons, but becomes more transparent when the electron density and the Fermi energy are increased by electrostatic gating. In contrast, the shear-type domain wall is always highly transparent at different gate voltages. Such soliton-dependent electronic transport can open up new routes to engineer novel nanoelectronic devices based on layer-stacking domain walls in bilayer graphene.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Ice-shelf ocean boundary layer dynamics from large-eddy simulations

Abstract. Small-scale turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small-scale processes, which regulate heat and salt transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Adaptive Determination of the Optimal Exchange Location in Wall-Modeled Large-Eddy Simulation

Wall-modeled large-eddy simulation introduces a modeling interface (or exchange location) separating the wall-modeled layer from the rest of the domain. The current state-of-the-art is to rely on user expertise when choosing where to place this modeling interface, whether this choice is tied to the grid or not. This paper presents a postprocessing algorithm that determines the exchange location systematically. The algorithm is based on a model for the error in the predicted wall shear stress and a model for the computational cost, and then finds the exchange location that minimizes a combination of the two. Here, the algorithm is tested both a priori and a posteriori using an equilibrium wall model for the flow over a wall-mounted hump, a boundary layer in an adverse pressure gradient, and a shock/boundary-layer interaction. The algorithm produces exchange locations that mostly agree with what an experienced user would suggest, with thinner wall-modeled layers in nonequilibrium flow regions and thicker wall-modeled layers where the boundary layer is closer to equilibrium. This suggests that the algorithm should be useful in simulations of realistic and highly complex geometries.

42 ENGINEERING↗
DOE OSTI2025

Probing Nanorod Assembly and Dynamics in Polymer Nanocomposites in Equilibrium and Shear

Coarse-grained molecular dynamics simulations are used to examine the structure and dynamics of nanorod assemblies in polymer melts under equilibrium and simple shear. We show that as the concentration of nanorods increases, there is a transition from an isotropic phase to a two-phase region in which the nanorods phase separate into a dilute phase and dense bundles of hexagonally packed nanorods. The onset of the two-phase region is below that predicted by Onsager theory, which we attribute to an effective increase in the diameter of nanorods due to a layer of polymer bound to the rod surfaces. Equilibrium simulations show that increasing polymer chain length N enhances nanorod bundling at fixed nanorod concentration. Simulations of systems undergoing simple shear show that flow enhances nanorod alignment and bundling relative to those of equilibrium systems with similar properties. Finally, simulations reveal that increasing nanorod length enhances nanorod alignment under shear at equivalent shear rates but not at equivalent Péclet numbers. Overall, our simulations highlight that polymer matrix-nanorod attraction (i.e., chemistry), shear rate, and matrix chain length are desirable design variables to control the structure of nanorod-containing soft materials under simple shear.

diffusion↗
DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗
DOE OSTI2025

pyTCR: A tropical cyclone rainfall model for python

pyTCR is a climatology software package developed in the Python programming language. It integrates the capabilities of several legacy physical models and increases computational efficiency to allow rapid estimation of tropical cyclone (TC) rainfall consistent with the large-scale environment. Specifically, pyTCR implements a horizontally distributed and vertically integrated model [Zhu et al., 2013] for simulating rainfall driven by TCs. Along storm tracks, rainfall is estimated by computing the cross-boundary-layer, upward water vapor transport caused by different mechanisms including frictional convergence, vortex stretching, large-scale baroclinic effect (i.e., wind shear), topographic forcing, and radiative cooling [Lu et al., 2018]. The package provides essential functionalities for modeling and interpreting spatio-temporal TC rainfall data. pyTCR requires a limited number of model input parameters, making it a convenient and useful tool for analyzing rainfall mechanisms driven by TCs. To sample rare (most intense) rainfall events that are often of great societal interest, pyTCR adapts and leverages outputs from a statistical-dynamical TC downscaling model [Lin et al., 2023] capable of rapidly generating a large number of synthetic TCs given a certain climate. As a result, pyTCR significantly reduces computational effort and improves the efficiency in capturing extreme TC rainfall events at the tail of the distributions from limited datasets. Furthermore, the TC downscaling model is forced entirely by large-scale environmental conditions from reanalysis data or coupled General Circulation Models (GCMs), simplifying the projection of TC-induced rainfall and wind speed under future climate using pyTCR. Finally, pyTCR can be coupled with hydrological and wind models to assess risks associated with independent and compound events (e.g., storm surges and freshwater flooding).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2022

Evolving Interior Mixing Regimes in a Tidal River Plume

Microstructure profiling was utilized to estimate vertical mixing (via vertical turbulent buoyancy flux) during a tidal pulse in the interior Merrimack River plume in calm winds. Multiple stratified shear mixing regimes appear and evolve with time. Initially the plume acts as a nearfield jet, with mixing in the plume (plume layer mixing) and over the plume-ambient interface (nearfield interfacial mixing). As the plume grows, interfacial mixing is suppressed offshore of the nearfield as currents slow, diminishing turbulent exchange between plume and shelf. At the end of ebb, ambient tidal currents reverse direction below plume, initiating another mode of internal, interfacial mixing (coined here as tidal interfacial mixing), allowing exchange between plume and ambient waters offshore. This work highlights previously unreported tidally modulated mixing within the near and midfield of a river plume.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Multi-material additive manufacturing of aluminum 6061-T6 alloy with stainless steel 304: Suppression of intermetallic compounds and interface growth mechanism

Fabricating integrated aluminum-steel structures is difficult because of metallurgical incompatibilities that promote intermetallic compound (IMC) formation at their interfaces. This work introduces a practical route for additively manufacturing 6061-T6 aluminum alloy (AA6061-T6) directly onto 304 stainless steel (SS304) using a high shear strain rate-assisted interfacial bonding mechanism. Through the friction extrusion deposition-based additive manufacturing method, both single- and multi-layer deposits of AA6061-T6 were successfully fabricated on SS304, yielding a robust hybrid multi-material structure. Comprehensive analyses of deposition quality, interface porosity, and bonding performance showed that the interface achieved a tensile strength exceeding 154 MPa under quasi-static loading. Microstructural observations revealed that the deposited aluminum layers experienced severe plastic deformation, resulting in pronounced grain refinement. Importantly, the interfacial zone was found to be free of brittle IMCs, instead containing a thin amorphous layer that evolved through a non-linear reaction-zone growth law. This solid-state additive strategy establishes a promising pathway for lightweight structural systems, nuclear energy component cladding, and multifunctional engineering components, redefining how dissimilar metals can be integrated for advanced applications.

Aluminum alloy↗
DOE OSTI2024

Methodology to determine printability criteria of highly concentrated pastes through rheological characterization

Material extrusion is an additive manufacturing technique that enables the creation of reproducible and complex hardware by depositing a viscous, shear-thinning ink onto a substrate in a custom-pattern via extrusion through a syringe. Here, the ability of an ink to be extruded onto a substrate in many layers, and maintain the desired shape is what defines the printability. Printability is often investigated by formulating, printing, and postmortem analysis of final parts in an iterative manner. Investigations of printability through rheological characterization have often been concerned with inks that straddle the line between printable and too thin, leaving out an entire class of inks that are highly-filled pastes, where extrudability is the limiting factor. Highly-filled pastes continue to pose issues for researchers as the effect of filler morphology, size, loading, and packing fraction on the ink rheology and corresponding printability is not understood. While traditional rheological characterizations may be useful for some inks, we show that protocols utilizing steady-shear, or large-amplitude oscillatory shear are difficult and unreliable for highly-filled pastes. Through transient rheology paired with real-time images we show that each traditional protocol produces inhomogeneous deformations that violate the assumptions that underly common rheological definitions. Instead, we demonstrate metrics measured with small-amplitude oscillatory shear that are correlated to the printability of various ink formulations ranging in loading. The rheological measures that accurately predict the printability of the inks are the axial stress measured at small amplitudes, and the critical stress amplitude above which rheological characterizations become impossible. In addition, we estimate the maximum packing fraction for each filler, based on the exponent common to hard sphere models, and show that the printability of each ink can be predicted by the ratio of the packing fraction to the theoretical maximum. We show how small-amplitude oscillatory shear allows users to develop printability criteria for any ink to enhance the workflow in the development of new inks, increase the performance of material extrusion printing, and improve the stability of printed parts, with less wasted time and materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Swirls and scoops: Ice base melt revealed by multibeam imagery of an Antarctic ice shelf

Knowledge gaps about how the ocean melts Antarctica’s ice shelves, borne from a lack of observations, lead to large uncertainties in sea level predictions. Using high-resolution maps of the underside of Dotson Ice Shelf, West Antarctica, we reveal the imprint that ice shelf basal melting leaves on the ice. Convection and intermittent warm water intrusions form widespread terraced features through slow melting in quiescent areas, while shear-driven turbulence rapidly melts smooth, eroded topographies in outflow areas, as well as enigmatic teardrop-shaped indentations that result from boundary-layer flow rotation. Full-thickness ice fractures, with bases modified by basal melting and convective processes, are observed throughout the area. This new wealth of processes, all active under a single ice shelf, must be considered to accurately predict future Antarctic ice shelf melt.

54 ENVIRONMENTAL SCIENCES↗