Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DRAG EFFECT”

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.

225 records · Page 13

Bipolar Membranes for Ion Management in (Photo)Electrochemical Energy Conversion

Conspectus (Photo)electrochemical energy conversion is important in the development of a carbon-neutral energy economy because it can provide a pathway for mitigating the intermittency of renewable energy sources such as wind and solar. In order to operate efficiently, these technologies, which include photoelectrochemical cells, water and CO 2 electrolyzers, fuel cells, and redox flow batteries, require fast charge transfer kinetics at the electrode/electrolyte interface as well as robust ion management in the electrolyte. In conventional electrolyzers and fuel cells, the electrolyte is strongly acidic or basic and ionic current is carried by H + or OH – ions. In contrast, photoelectrodes and electrocatalysts for water splitting are often studied in buffered solutions. The question of ion balance in these systems led us to analyze the polarization losses due to ion concentration gradients in cells that employed various buffer–membrane combinations. Continuously driving the buffer ions across an ionomer membrane not only lowers the buffer capacity of an aqueous electrolyte but also introduces pH gradients that result in significant energy losses. To address the problem, we and other groups have studied the use of reverse-biased bipolar membranes (BPMs) in (photo)electrolytic cells. BPMs consist of an anion exchange layer (AEL) laminated with a cation exchange layer (CEL) and are usually equipped with a catalytic layer in between to accelerate the water dissociation reaction. At the AEL/CEL interface, water dissociates into protons and hydroxide ions, which replenish those consumed at the cathode and anode. Compared to conventional water electrolyzers with proton/anion exchange membranes (PEM/AEM), BPM electrolyzers provide the unique advantage of continuously operating the cathode and anode under different pH conditions, which is desirable when the two electrode reactions have different pH requirements. BPMs also enable the use of buffered electrolytes at pH values that are optimized for electrode stability and product selectivity in applications such as CO 2 electrolysis. Product crossover losses and CO 2 pumping can be dramatically reduced in BPM-based CO 2 electrolyzers, relative to conventional alkaline membranes, by electrostatic repulsion (of anionic products) and electroosmotic drag (for neutral products). BPM-based gas fed CO 2 electrolyzers can achieve high current density, but they suffer from low Faradaic efficiency (FE) due to the acidic local environment of the CEL. This problem can be mitigated by adding an aqueous buffering layer or by creating a weak acid cation exchange film on the CEL face of the membrane. The use of BPMs in fuel cells and redox flow batteries offers some interesting advantages. Configurations with both reverse and forward bias have been studied, but forward bias has been favored due to material compatibility, reaction kinetics, and thermodynamic considerations. The net reaction at the AEL/CEL interface is the acid–base neutralization reaction, which has a high inherent reaction rate constant, but in the BPM is limited to a nanometric space-charge layer and requires efficient catalysis to achieve high current density. Understanding the mechanism of the acid–base neutralization and the opposite process, the water dissociation reaction, will be essential for improving the performance of forward-biased BPMs. In conclusion, this Account reviews our current understanding of the working mechanisms of BPMs and discusses how we can use them to effectively manage ions for various (photo)electrochemical applications.

Yan, Zhifei↗

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗

Droplet-Gas Interactions in Nanoelectrospray Multiphase Flow

Electrospray (ES) is production of a charged droplet plume that transits from an electrically biased supply capillary to a counter electrode under the influence of an electric field. Drag interactions with the surrounding gas decelerate the droplets while transferring the momentum from droplets to accelerate the surrounding gas, resulting in an induced gas flow. This work seeks to characterize the structure of the resulting gas jet and identify the key mechanisms defining the flow structure. The phenomenon of gas jets produced by momentum transfer from nano-electrospray (nES) plumes is explored with schlieren visualization, thermal anemometry, and numerical simulations. Schlieren visualization experiments provide information on the flow structure in support of simulation predictions, and the hot thermistor anemometry measurements of gas velocities outside the spray demonstrate quantitatively validated simulation results. The study reveals the formation of a moderately high velocity coaxial gas jet within the nES plume and provides insight into the overall flow structure of the induced flow. The multiphase electrohydrodynamic simulations enable numerical experimentation to explore the fundamental physics of coupled droplet-gas transport and the resulting flow structure. The simulations, confirmed by the experiments, reveal gas jetting induced by nES with a narrow (sub hundred micrometers in diameter) core originating from the nES liquid-jet breakup region and a surrounding larger-in-extent zone (several hundred micrometers in diameter) of lower velocity gas flow within the electrospray plume. This behavior is due to nES ejecting a stream of droplets from the tip of a narrow liquid cone-jet, where the combined effect of many small droplets transferring momentum to a confined region of gas yields a narrow but high-velocity gas stream. Furthermore, the micro-jets produced by nES have practical utility for mass spectrometry, 3D printing and fabrication, and thermal management.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Airfoil Computational Fluid Dynamics - 9k shapes, 2 AoA's

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 8,996 airfoil shapes, computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this data, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We fixed the linear deformations to be the mean over the database and sampled new shapes over a four-dimensional parameter space of higher-order perturbation. This sampling approaches allows for isolated analysis of non-linear airfoil shape deformations while holding other aspects (e.g., airfoil thickness) approximately constant. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, Reynolds number of 9M, and at two angles of attack, 4 deg. and 12 deg. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_9k_data.ipynb.

9k↗

The dynamics of suspensions of prolate spheroidal particles—Effects of volume fraction

Here, the effect of volume fraction on the dynamics of a gravity driven suspension of prolate spheroidal solid particles in a fully periodic domain is examined by fully resolved numerical simulations, for relatively modest Reynolds numbers (around 20, depending on the volume fraction). Three systems are examined, for volume fractions of 2.2% (20 particles), 5.5% (50 particles), and 9.9% (90 particles). The results show a transition from flow at low volume fractions dominated by the hydrodynamic interactions between the particles and the fluid, although modified by collisions, to flow dominated by collisions at higher volume fractions. The distribution of the particles with respect to each other, as measured by the probability distribution of nearest distances, is nearly random at low volume fractions but at high volume fractions, there is more clustering than for a random distribution. At lower volume fractions most of the particles fall broadside on whereas at the highest volume fraction, their orientation is essentially random.

42 ENGINEERING↗

SARS-CoV-2 M pro Protease Variants of Concern Display Altered Viral Substrate and Cell Host Target Galectin-8 Processing but Retain Sensitivity toward Antivirals

The main protease of SARS-CoV-2 (M pro ) is the most promising drug target against coronaviruses due to its essential role in virus replication. With newly emerging variants there is a concern that mutations in M pro may alter the structural and functional properties of protease and subsequently the potency of existing and potential antivirals. We explored the effect of 31 mutations belonging to 5 variants of concern (VOCs) on catalytic parameters and substrate specificity, which revealed changes in substrate binding and the rate of cleavage of a viral peptide. Crystal structures of 11 M pro mutants provided structural insight into their altered functionality. Additionally, we show M pro mutations influence proteolysis of an immunomodulatory host protein Galectin-8 (Gal-8) and a subsequent significant decrease in cytokine secretion, providing evidence for alterations in the escape of host-antiviral mechanisms. Accordingly, mutations associated with the Gamma VOC and highly virulent Delta VOC resulted in a significant increase in Gal-8 cleavage. Importantly, IC50s of nirmatrelvir (Pfizer) and our irreversible inhibitor AVI-8053 demonstrated no changes in potency for both drugs for all mutants, suggesting M pro will remain a high-priority antiviral drug candidate as SARS-CoV-2 evolves.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reverse phonon thermal flux from an applied electric field

In the phenomenon known as electron drag, a phonon thermal flux is established in a conducting crystal by an electric field applied under isothermal conditions through directed transfer of quasimomentum from the electronic charge current to the phonon subsystem. Prior understanding of this phenomenon involves a thermal current composed of low-frequency acoustic phonons that is in the same direction as the charge current. This results in an increase of the Peltier thermopower. Here, we show that it is also possible to establish a phonon thermal current that is in the opposite direction to the charge current, which we refer to as a reverse phonon thermal flux. We demonstrate this behavior both through a simple qualitative model and through first-principles calculations performed for three materials: p-type 𝜃-TaN, n-type BAs, and n-type Si. The reverse phonon flux is shown to arise through intervalley electron-phonon scattering processes involving high-frequency acoustic phonons. Unlike n-BAs and n-Si, the combined features in the band structure and phonon dispersions of p-type 𝜃-TaN promote a particularly large reverse phonon thermal flux comparable in magnitude to the oppositely directed thermal flux from low-frequency phonons. Finally, this work highlights a previously unrecognized behavior in the coupled electron-phonon system and advances our understanding of the rich physics of transport in solids.

lattice dynamics↗