Purge dynamics of fuel cells.
Anisotropic porous media model for fluid motion in rectangular fuel cell cavities, analyzing pressure, velocity, stream functions and purge time
Gidaspow, D.↗
SEARCH · Search NASA
Results for “ANISOTROPIC FLUID”
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.
Ten-moment fluid modeling of the Weibel instability
We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.
Kuldinow, D. A. (ORCID:0000000319730196)↗
Dispersion relation and instability for an anisotropic nonuniform flowing plasma
A generalized formula for wave instability is developed for an anisotropic nonuniform plasma with finite flows and temperatures. Six-moment fluid equations are solved to give the analytic expression for wave instability in arbitrarily nonuniform plasmas. The analytic formula explicitly states the dependence of wave instability on the nonuniformities of number density, flow velocity, and anisotropic or isotropic pressure. The accuracy of the formalism is verified by a numerical calculation of implicit dispersion relations in complex Fourier space. The analysis shows that nonuniformity plays a critical role in plasma instability, while the flow velocity and anisotropic pressures determine the growth rate of the instability. Lastly, the instability diagram and associated instability criterion for anisotropy-driven instability are introduced as applications of the formalism.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Purge dynamics of fuel cells
Fuel cell cavities, analyzing fluid motion purge dynamics based on anisotropic porous media model
Gidaspow, D.↗
Continuation of tailored composite structures of ordered staple thermoplastic material
The search for the cost effective composite structure has motivated the investigation of several approaches to develop composite structure from innovative material forms. Among the promising approaches is the conversion of a planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. A framework was established which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. Predictions for the effective viscosities of a hyper-anisotropic medium consisting of collimated, discontinuous fibers suspended in viscous matrix were extended to capture the characteristics of typical polymers including non-Newtonian behavior and temperature dependence. In addition, the influence of fiber misorientation was also modeled by compliance averaging to determine ensemble properties for a given orientation distribution. A design tool is presented for predicting the effect of material heterogeneity on the performance of curved composite beams such as those used in aircraft fuselage structures. Material heterogeneity can be induced during manufacturing processes such as sheet forming and stretch forming of thermoplastic composites. This heterogeneity can be introduced in the form of fiber realignment and spreading during the manufacturing process causing radial and tangential gradients in material properties. Two analysis procedures are used to solve the beam problems. The first method uses separate two-dimensional elasticity solutions for the stresses in the flange and web sections of the beam. The separate solutions are coupled by requiring that forces and displacements match section boundaries. The second method uses an approximate Rayleigh-Ritz technique to find the solutions for more complex beams. Analyses are performed for curved beams of various cross-sections loaded in pure bending and with a uniform distributed load. Preliminary results show that the geometry of the beam dictates the effect of heterogeneity on performance. The role of heterogeneity is larger in beams with a small average radius-to-depth ration, R/t, where R is the average radius of the beam and t is the difference between the inside and outside radii. Results of the anlysis are in the form of stresses and displacements and are compared to both mechanics of materials and numerical solutions obtained using finite element analysis.
Santare, Michael H.↗
Interference Provides Clarity: Direct Observation of 2D Materials at Fluid–Fluid Interfaces
Monolayer particles of two-dimensional (2D) materials represent a scientifically and technologically interesting class of anisotropic particles with colloidal-scale lateral sizes but sub-nanometer thicknesses. This atomic-scale thickness leads to interesting phenomena that can be exploited in next-generation thin-film technologies, and fluid-fluid interfaces provide a potentially scalable platform to confine, assemble, and deposit functional thin films of 2D materials. However, directly observing how these materials interact and assemble into a given film morphology is experimentally challenging because of their sub-nanometer thicknesses. Here, we demonstrate the ability to directly observe graphene, molybdenum disulfide (MoS2), and hexagonal boron nitride (h-BN) particles at fluid-fluid interfaces using interference reflection microscopy (IRM). Monolayer MoS2 and graphene particles demonstrated >10% optical contrast at an air-water interface, which allowed us to quantitatively analyze in situ images of self-assembled MoS2 particles and to map trajectories of interacting graphene particles. Additionally, the Brownian motion of a graphene particle was tracked and analyzed in the context of passive microrheology theory for 2D particle probes. Our results demonstrate how IRM can be used to obtain quantitative spatiotemporal information regarding the self-assembly and dynamics of 2D materials at fluid-fluid interfaces. It will have a significant impact on our ability to investigate systems of atomically thin particles at fluid-fluid interfaces, an area that has fundamental scientific importance and materials science applications but has suffered from a lack of direct, in situ observation techniques.
2D material↗
Gneiss-charnockite transformation at Kottavattam, Southern Kerala (India)
At Kottavattam, leucocratic granitic garnet-biotite gneisses (age less than 2 Ga) were partially transformed to coarse-grained charnockite along a system of conjugate fractures (N70E and N20W) and the foliation planes (N60 to 80W; dip 80 to 90 SW) about 550 m.y. ago. To examine and quantify changes in fabric, mineralogy, pore fluids and chemical composition associated with this process, large rock specimens showing gneiss-charnockite transition were studied in detail. The results of the present study corroborate the concept that charnockite formation at Kottavattam is an internally-generated phenomenon and was not triggered by the influx of carbonic fluids from a deep-seated source. It is suggested that charnockitization was caused by the following mechanism: (1) near-isothermal decompression during uplift of the gneiss complex led to an increase of the pore fluid pressure (P sub fluid greater than P sub lith) which - in a regime of anisotropic stress - triggered or at least promoted the development of conjugate fractures; (2) the simultaneous release of pore fluids from bursting fluid inclusions and their escape into the developing fracture system resulted in a drop of fluid pressure; and (3) the internal generation and buffering of the fluids and their, probably, limited migration in an entirely granitic rock system explains the absence of any significant metasomatic mass transfer.
Raith, M.↗
Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model
The global dynamics of the ionized and neutral components in the environment of Io plays an important role in the interaction of Jupiter's corotating magnetospheric plasma with Io. The stationary simulation of this problem was done in the MHD and the electrodynamics approaches. One of the main significant results from the simplified two-fluid model simulations was a production of the structure of the double-peak in the magnetic field signature of the I0 flyby that could not be explained by standard MHD models. In this paper, we develop a method of kinetic ion simulation. This method employs the fluid description for electrons and neutrals whereas for ions multilevel, drift-kinetic and particle, approaches are used. We also take into account charge-exchange and photoionization processes. Our model provides much more accurate description for ion dynamics and allows us to take into account the realistic anisotropic ion distribution that cannot be done in fluid simulations. The first results of such simulation of the dynamics of ions in the Io's environment are discussed in this paper.
Lipatov, Alexander S.↗
Ten-moment fluid model for low-temperature magnetized plasmas
In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.
Kuldinow, Derek Amur (ORCID:0000000319730196)↗
Topological structure and dynamics of three-dimensional active nematics
Topological structures are effective descriptors of the nonequilibrium dynamics of diverse many-body systems. For example, motile, point-like topological defects capture the salient features of two-dimensional active liquid crystals composed of energy-consuming anisotropic units. We dispersed force-generating microtubule bundles in a passive colloidal liquid crystal to form a three-dimensional active nematic. Light-sheet microscopy revealed the temporal evolution of the millimeter-scale structure of these active nematics with single-bundle resolution. The primary topological excitations are extended, charge-neutral disclination loops that undergo complex dynamics and recombination events. Our work suggests a framework for analyzing the nonequilibrium dynamics of bulk anisotropic systems as diverse as driven complex fluids, active metamaterials, biological tissues, and collections of robots or organisms.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Interface Morphology During Crystal Growth: Effects of Anisotropy and Fluid Flow
The effect of a parallel shear flow and anisotropic interface kinetics on the onset of instability during growth from a supersaturated solution is analyzed. The model used for anisotropy is based on the microscopic picture of step motion. A shear flow (linear Couette flow or asymptotic suction profile) parallel to the crystal solution interface in the same direction as the step motion decreases interface stability. A shear flow counter to the step motion enhances stability and for sufficiently large shear rates the interface is absolutely morphologically stable. For large wave numbers, the perturbed flow field can be neglected and a simple analytic approximation for the stability-instability demarcation is found.
Coriell, S. R.↗
Turbulent Output-Based Anisotropic Adaptation
Controlling discretization error is a remaining challenge for computational fluid dynamics simulation. Grid adaptation is applied to reduce estimated discretization error in drag or pressure integral output functions. To enable application to high O(10(exp 7)) Reynolds number turbulent flows, a hybrid approach is utilized that freezes the near-wall boundary layer grids and adapts the grid away from the no slip boundaries. The hybrid approach is not applicable to problems with under resolved initial boundary layer grids, but is a powerful technique for problems with important off-body anisotropic features. Supersonic nozzle plume, turbulent flat plate, and shock-boundary layer interaction examples are presented with comparisons to experimental measurements of pressure and velocity. Adapted grids are produced that resolve off-body features in locations that are not known a priori.
Park, Michael A.↗
Effects of Kinetic Processes in Shaping Io's Global Plasma Environment: A 3D Hybrid Model
The global dynamics of the ionized and neutral gases in the environment of Io plays an important role in the interaction of Jupiter s corotating magnetospheric plasma with Io. Stationary simulations of this problem have already been done using the magnetohydrodynamics (MHD) and the electrodynamics approaches. One of the major results of recent simplified two-fluid model simulations [Saur, J., Neubauer, F.M., Strobel, D.F., Summers, M.E., 2002. J. Geophys. Res. 107 (SMP5), 1-18] was the production of the structure of the double-peak in the magnetic field signature of the Io flyby. These could not be explained before by standard MHD models. In this paper, we present a hybrid simulation for Io with kinetic ions and fluid electrons. This method employs a fluid description for electrons and neutrals, whereas for ions a particle approach is used. We also take into account charge-exchange and photoionization processes and solve self-consistently for electric and magnetic fields. Our model may provide a much more accurate description for the ion dynamics than previous approaches and allows us to account for the realistic anisotropic ion velocity distribution that cannot be done in fluid simulations with isotropic temperatures. The first results of such a simulation of the dynamics of ions in Io s environment are discussed in this paper. Comparison with the Galileo IO flyby results shows that this approach provides an accurate physical basis for the interaction and can therefore naturally reproduce all the observed salient features.
Lipatov, Alexander S.↗
A solar-wind model including proton thermal anisotropy.
A relatively simple model is proposed to study the magnetohydrodynamic expansion of the solar wind. The interplanetary space is divided into an 'inner region' and an 'outer region.' In the inner region the solar wind is one-fluid and thermally isotropic. The solar-wind protons are assumed to be completely collisionless in the outer region. Thus the solar wind is two-fluid in nature, and the protons are thermally anisotropic in the outer region. It is assumed that throughout the interplanetary space electrons are thermally isotropic due to collisions or other processes. The electron heat flux is supposed to obey the well-known relationships that the flux is parallel to the magnetic-field vector and is directly proportional to the spatial derivative of the electron temperature.
Whang, Y. C.↗
Long discontinuous fiber composite structure: Forming and structural mechanics
Cost effective composite structure has motivated the investigation of several new approaches to develop composite structure from innovative material forms. Among the promising new approaches is the conversion of planar sheet to components of complex curvature through sheet forming or stretch forming. In both cases, the potential for material stretch in the fiber direction appears to offer a clear advantage in formability over continuous fiber systems. In the present study, the authors have established a framework which allows the simulation of the anisotropic mechanisms of deformation of long discontinuous fiber laminates wherein the matrix phase is a viscous fluid. The initial study focuses upon the establishment of micromechanics models for prediction of the effective anisotropic viscosities of the oriented fiber assembly in a viscous matrix. Next, the developed constitutive relation is employed through an analogy with incompressible elasticity to exercise the finite element technique for determination of local fiber orientation and laminate thickness after forming. Results are presented for the stretch bending of a curved beam from an arbitrary composite laminate and the bulging of a clamped sheet. Structural analyses are conducted to determine the effect of microstructure on the performance of curved beams manufactured from long discontinuous fiber composites. For the purposes of this study, several curved beams with ideal and non-ideal microstructures are compared for response under pure bending. Material parameters are determined from a separate microstructural analysis.
Pipes, R. B.↗
The problem of exact interior solutions for rotating rigid bodies in general relativity
The (3 + 1) dyadic formalism for timelike congruences is applied to derive interior solutions for stationary, axisymmetric, rigidly rotating bodies. In this approach the mathematics is formulated in terms of three-space-covariant, first-order, vector-dyadic, differential equations for a and Omega, the acceleration and angular velocity three-vectors of the rigid body; for T, the stress dyadic of the matter; and for A and B, the 'electric' and 'magnetic' Weyl curvature dyadics which describe the gravitational field. It is shown how an appropriate ansatz for the forms of these dyadics can be used to discover exact rotating interior solutions such as the perfect fluid solution first published in 1968. By incorporating anisotropic stresses, a generalization is found of that previous solution and, in addition, a very simple new solution that can only exist in toroidal configurations.
Wahlquist, H. D.↗
Kinetic structure of rotational discontinuities: Implications for the magnetopause
Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.
Krauss-Varban, D.↗
Extending the single-fluid solvability conditions for more general plasma systems
We extend the single-fluid solvability conditions to include plasma systems at arbitrary $\beta$ with arbitrary flows and external forcing terms. This treatment includes both the isotropic and the anisotropic cases. The generalized conditions that result can be used to generate certain classes of single-fluid equilibria.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗