Engineering PapersSearch

SEARCH · Engineering Papers

Results for “hydrodynamic coefficients”

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 37 records · Page 2

Theory of the lattice Boltzmann Method: Dispersion, Dissipation, Isotropy, Galilean Invariance, and Stability

The generalized hydrodynamics (the wave vector dependence of the transport coefficients) of a generalized lattice Boltzmann equation (LBE) is studied in detail. The generalized lattice Boltzmann equation is constructed in moment space rather than in discrete velocity space. The generalized hydrodynamics of the model is obtained by solving the dispersion equation of the linearized LBE either analytically by using perturbation technique or numerically. The proposed LBE model has a maximum number of adjustable parameters for the given set of discrete velocities. Generalized hydrodynamics characterizes dispersion, dissipation (hyper-viscosities), anisotropy, and lack of Galilean invariance of the model, and can be applied to select the values of the adjustable parameters which optimize the properties of the model. The proposed generalized hydrodynamic analysis also provides some insights into stability and proper initial conditions for LBE simulations. The stability properties of some 2D LBE models are analyzed and compared with each other in the parameter space of the mean streaming velocity and the viscous relaxation time. The procedure described in this work can be applied to analyze other LBE models. As examples, LBE models with various interpolation schemes are analyzed. Numerical results on shear flow with an initially discontinuous velocity profile (shock) with or without a constant streaming velocity are shown to demonstrate the dispersion effects in the LBE model; the results compare favorably with our theoretical analysis. We also show that whereas linear analysis of the LBE evolution operator is equivalent to Chapman-Enskog analysis in the long wave-length limit (wave vector k = 0), it can also provide results for large values of k. Such results are important for the stability and other hydrodynamic properties of the LBE method and cannot be obtained through Chapman-Enskog analysis.

Lallemand, Pierre

Mechanistic mass transfer in hollow fiber membrane solvent extraction for bio-based isobutanol

Membrane solvent extraction (MSE) has emerged as a promising method for selectively recovering bioproducts from complex aqueous streams. Bio-isobutanol, a next-generation feedstock for biofuel, remains challenging to recover because of its low concentration and the presence of inhibitory substances. This study explores the potential of hollow fiber (HF) MSE for bio-isobutanol recovery and systematically examines the coupled effects of fiber packing, shell-side flow dynamics, and aqueous chemistry on performance. A resistance-in-series model is applied to understand mass transfer in the HF MSE modules, quantify local resistances, and validate overall performance. The results show that increasing the fiber packing provides a larger interfacial area but induces poor flow distribution and channeling, hindering effective isobutanol transport. Meanwhile, increasing the shell-side velocity improves isobutanol recovery due to reductions in the boundary layer thickness. The presence of salts, added to mimic fermentation broth, increases the partition coefficient through salting-out effects, further improving isobutanol flux. A modified correlation for the shell-side mass transfer coefficient (k s,ϕ+v ), integrating geometric and hydrodynamic effects, was developed and validated. The proposed model achieves highly predictive accuracy (r 2 = 0.9808) across a wide range of conditions, outperforming previous models. The findings provide mechanistic insight into the interaction of geometric packing, hydrodynamics, and chemistry in governing mass transfer in HF MSE. Overall, this work demonstrates the potential of HF MSE for efficient bio-isobutanol recovery and also provides practical guidelines on critical factors (packing fraction, partition coefficient, and shell-side velocity), aiding in the design and scaling of MSE systems for resource recovery.

Aqueous chemistry

Consistent transport coefficients in astrophysics

A consistent theory for dealing with transport phenomena in stellar atmospheres starting with the kinetic equations and introducing three cases (LTE, partial LTE, and non-LTE) was developed. The consistent hydrodynamical equations were presented for partial-LTE, the transport coefficients defined, and a method shown to calculate them. The method is based on the numerical solution of kinetic equations considering Landau, Boltzmann, and Focker-Planck collision terms. Finally a set of results for the transport coefficients derived for a partially ionized hydrogen gas with radiation was shown, considering ionization and recombination as well as elastic collisions. The results obtained imply major changes is some types of theoretical model calculations and can resolve some important current problems concerning energy and mass balance in the solar atmosphere. It is shown that energy balance in the lower solar transition region can be fully explained by means of radiation losses and conductive flux.

Fontenla, Juan M.

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY

Critical fluid dynamics in two and three dimensions

We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS

The Friction of Piston Rings

The coefficient of friction between piston ring and cylinder liner was measured in relation to gliding acceleration, pressure, temperature, quantity of oil and quality of oil. Comparing former lubrication-technical tests, conclusions were drawn as to the state of friction. The coefficients of friction as figured out according to the hydrodynamic theory were compared with those measured by tests. Special tests were made on "oiliness." The highest permissible pressure was measured and the ratio of pressure discussed.

LUBRICATING OILS - OILINESS

Charged particle transport coefficient challenges in high energy density plasmas

High energy density physics (HEDP) and inertial confinement fusion (ICF) research typically relies on computational modeling using radiation-hydrodynamics codes in order to design experiments and understand their results. These tools, in turn, rely on numerous charged particle transport and relaxation coefficients to account for laser energy absorption, viscous dissipation, mass transport, thermal conduction, electrical conduction, non-local ion (including charged fusion product) transport, non-local electron transport, magnetohydrodynamics, multi-ion-species thermalization, and electron-ion equilibration. In many situations, these coefficients couple to other physics, such as imposed or self-generated magnetic fields. Furthermore, how these coefficients combine are sensitive to plasma conditions as well as how materials are distributed within a computational cell. Uncertainties in these coefficients and how they couple to other physics could explain many of the discrepancies between simulation predictions and experimental results that persist in even the most detailed calculations. This paper reviews the challenges faced by radiation-hydrodynamics in predicting the results of HEDP and ICF experiments with regard to these and other physics models typically included in simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Simulation of a High-Confinement Tidal Energy Array (Abstract)

When the projected area of a current turbine array occupies a substantial fraction of a tidal or river channel cross-sectional area, the confinement associated with the channel boundaries substantially increases turbine efficiency. This effect can be exploited to substantially reduce the levelized cost of energy for current turbine arrays. However, the presence of such an array will affect inflow velocity, water depth, and circulation in the surrounding area. This has implications for array performance and environmental impacts. This project combines existing numerical modeling capability from PNNL with an understanding of confined array performance from UW. The work will 1) Simulate tidal currents and water surface elevation in a regional circulation model that includes the Knik Arm site at appropriate resolution; 2) Incorporate a varying thrust coefficient into this model that reflects the outcomes of tidal hydrodynamic simulations of confined turbine arrays; 3) Use the model to determine the effects of array operation on water levels and inflow conditions as a function of different blockage ratios and control strategies; and 4) Identify the “far field” environmental implications of operating an array of this type. UW will use the outputs from this model to evaluate the implications of “feedback effects” and variable inflow conditions on levelized cost of energy estimates for a confinement exploiting array of cross-flow turbines.

16 TIDAL AND WAVE POWER

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity

Nonlinearity Effects of Lateral Density Diffusion Coefficient on Gain-Guided VCSEL Performance

Electron and hole diffusions in the plane of semiconductor quantum wells play an important part in the static and dynamic operations of semiconductor lasers. In this paper, we apply a hydrodynamic model developed from the semiconductor Bloch equations to numerically study the effects of nonlinearity in the diffusion coefficient on single mode operation and direct modulation of a gain-guided InGaAs/GaAs multiple quantum well laser, operating not too far from threshold. We found that a small diffusion coefficient is advantageous for lowering the threshold current and increasing the modulation bandwidth. Most importantly, the effects of nonlinearity in the coefficient can be approximately reproduced by replacing the coefficient with an effective constant diffusion coefficient, which corresponds roughly to the half height density of the density distribution.

Li, Jian-Zhong

Steady-state and transient analysis of a squeeze film damper bearing for rotor stability

A study of the steady-state and transient response of the squeeze film damper bearing is presented. Both the steady-state and transient equations for the hydrodynamic bearing forces are derived. The bearing equivalent stiffness and damping coefficients are determined by steady-state equations. These coefficients are used to find the bearing configuration which will provide the optimum support characteristics based on a stability analysis of the rotor-bearing system. The transient analysis of rotor-bearing systems is performed by coupling the bearing and journal equations and integrating forward in time. The effects of unbalance, cavitation, and retainer springs are included in the analysis. Methods of determining the stability of a rotor-bearing system under the influence of aerodynamic forces and internal shaft friction are discussed with emphasis on solving the system characteristic frequency equation and on producing stability maps. It is shown that for optimum stability and low force transmissability the squeeze bearing should operate at an eccentricity ratio epsilon 0.4.

Barrett, L. E.

Cellular automaton formulation of passive scalar dynamics

Cellular automata modeling of the advection of a passive scalar in a two-dimensional flow is examined in the context of discrete lattice kinetic theory. It is shown that if the passive scalar is represented by tagging or 'coloring' automation particles a passive advection-diffusion equation emerges without use of perturbation expansions. For the specific case of the hydrodynamic lattice gas model of Frisch et al. (1986), the diffusion coefficient is calculated by perturbation.

Chen, Hudong

Bifurcations from stationary to periodic solutions in a low-order model of forced, dissipative barotropic flow

The considered investigation is concerned with periodic solutions in the context of a forced, dissipative, barotropic spectral model truncated to three complex coefficients with constant forcing on only the intermediate scale. It is found that determining a periodic solution of this three-coefficient model also reduces to finding the algebraic roots of a real polynomial. In the derivation of this polynomial, a class of hydrodynamic spectral systems is described for which a periodic solution might be similarly specified. The existence of periodic solutions of the three-coefficient model is controlled by the roots of the stability polynomial of the basic stationary solution, which represents the simplest response to the constant forcing. When the forcing exceeds a critical value, the basic solution becomes unstable. Owing to the nature of the roots of the stability polynomial at critical forcing, bifurcation theory guarantees the existence of a periodic solution.

Mitchell, K. E.

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics

A Hydrodynamic Theory for Spatially Inhomogeneous Semiconductor Lasers: Numerical Results - 2

We present numerical results of the diffusion coefficients (DCs) in the coupled diffusion model derived in the preceding paper for a semiconductor quantum well. These include self and mutual DCs in the general two-component case, as well as density- and temperature-related DCs under the single-component approximation. The results are analyzed from the viewpoint of free Fermi gas theory with many-body effects incorporated. We discuss in detail the dependence of these DCs on densities and temperatures in order to identify different roles played by the free carrier contributions including carrier statistics and carrier-LO phonon scattering, and many-body corrections including bandgap renormalization and electron-hole (e-h) scattering. In the general two-component case, it is found that the self- and mutual- diffusion coefficients are determined mainly by the free carrier contributions, but with significant many-body corrections near the critical density. Carrier-LO phonon scattering is dominant at low density, but e-h scattering becomes important in determining their density dependence above the critical electron density. In the single-component case, it is found that many-body effects suppress the density coefficients but enhance the temperature coefficients. The modification is of the order of 10% and reaches a maximum of over 20% for the density coefficients. Overall, temperature elevation enhances the diffusive capability or DCs of carriers linearly, and such an enhancement grows with density. Finally, the complete dataset of various DCs as functions of carrier densities and temperatures provides necessary ingredients for future applications of the model to various spatially inhomogeneous optoelectronic devices.

Li, Jianzhong

Real time estimation of the heaving and pitching motions of a ship, using a Kalman filter

In the present study the estimation of the heave and pitch motion of a ship is considered, using Kalman filtering techniques. A significant part of the study is devoted to constructing appropriate models for the sea and the ship. The governing equations are obtained from hydrodynamic considerations in the form of linear differential equations with frequency dependent coefficients. In addition, nonminimum phase characteristics are obtained due to the spatial integration of the water wave forces. The resulting transfer matrix function is irrational and nonminimum phase. The conditions for a finite-dimensional approximation are considered and the impact of the various parameters is assessed. A numerical application is considered for a DD-963 destroyer.

Triantafyllou, M.

Real time estimation of ship motions using Kalman filtering techniques

The estimation of the heave, pitch, roll, sway, and yaw motions of a DD-963 destroyer is studied, using Kalman filtering techniques, for application in VTOL aircraft landing. The governing equations are obtained from hydrodynamic considerations in the form of linear differential equations with frequency dependent coefficients. In addition, nonminimum phase characteristics are obtained due to the spatial integration of the water wave forces. The resulting transfer matrix function is irrational and nonminimum phase. The conditions for a finite-dimensional approximation are considered and the impact of the various parameters is assessed. A detailed numerical application for a DD-963 destroyer is presented and simulations of the estimations obtained from Kalman filters are discussed.

Triantafyllou, M. S.

Rough-water Impact-load Investigation of a Chine-immersed V-bottom Model Having a Dead-rise Angle of 10 Degrees

A hydrodynamic rough-water impact-loads investigation of a fixed-trim V-bottom float with a beam-loading coefficient of 5.78 and dead-rise angle of 10 degrees was made at the Langley impact basin. The size of the waves varied from approximately 10 to 60 feet in length and 1 to 2 feet in height. Time histories were obtained showing the position of the model relative to the wave throughout the impact and typical examples are presented. The load coefficient was found to vary primarily with the slope of the impacting wave.

HYDRODYNAMIC CONFIGURATIONS - GENERAL STUDIES