Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 703 records · Page 39

Absolute flux measurements for swift atoms

While a torsion balance in vacuum can easily measure the momentum transfer from a gas beam impinging on a surface attached to the balance, this measurement depends on the accommodation coefficients of the atoms with the surface and the distribution of the recoil. A torsion balance is described for making absolute flux measurements independent of recoil effects. The torsion balance is a conventional taut suspension wire design and the Young modulus of the wire determines the relationship between the displacement and the applied torque. A compensating magnetic field is applied to maintain zero displacement and provide critical damping. The unique feature is to couple the impinging gas beam to the torsion balance via a Wood's horn, i.e., a thin wall tube with a gradual 90 deg bend. Just as light is trapped in a Wood's horn by specular reflection from the curved surfaces, the gas beam diffuses through the tube. Instead of trapping the beam, the end of the tube is open so that the atoms exit the tube at 90 deg to their original direction. Therefore, all of the forward momentum of the gas beam is transferred to the torsion balance independent of the angle of reflection from the surfaces inside the tube.

Fink, M.↗

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Effect of f -element complexation on the radiolysis of 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP])

A systematic study of the impact on the chemical reactivity of the oxidising n-dodecane radical cation (RH˙ + ) with f-element complexed 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) has been undertaken utilizing time-resolved electron pulse radiolysis/transient absorption spectroscopy and high-level quantum mechanical calculations. Lanthanide ion complexed species, [Ln((HEH[EHP]) 2 ) 3 ], exhibited vastly increased reactivity (over 10× faster) in comparison to the non-complexed ligand in n-dodecane solvent, whose rate coefficient was k = (4.66 ± 0.22) × 10 9 M -1 s -1 . Similar reactivity enhancement was also observed for the corresponding americium ion complex, k = (5.58 ± 0.30) × 10 10 M -1 s -1 . The vastly increased reactivity of these f-element complexes was not due to simple increased diffusion-control of these reactions; rather, enhanced hole transfer mechanisms for the complexes were calculated to become energetically more favourable. Interestingly, the observed reactivity trend with lanthanide ion size was not linear; instead, the rate coefficients showed an initial increase (Lu to Yb) followed by a decrease (Tm to Ho), followed by another increase (Dy to La). This behaviour was excellently predicted by the calculated reaction volumes of these complexes. In conclusion, complementary cobalt-60 gamma irradiations for select lanthanide complexes demonstrated that the measured kinetic differences translated to increased ligand degradation at steady-state timescales, affording ~38% increase in ligand loss of a 1:1 [La((HEH[EHP]) 2 ) 3 ]: HEH[EHP] ratio system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The effect of breaking gravity waves on the dynamics and chemistry of the mesosphere and lower thermosphere (invited review)

The influence of breaking gravity waves on the dynamics and chemical composition of the 60 to 110 km region is investigated with a two dimensional model that includes a parameterization of gravity wave momentum deposition and diffusion. The dynamical model is described by Garcia and Solomon (1983) and Solomon and Garcia (1983) and includes a complete chemical scheme for the mesosphere and lower thermosphere. The parameterization of Lindzen (1981) is used to calculate the momentum deposited and the turbulent diffusion produced by the gravity waves. It is found that wave momentum deposition drives a very vigorous mean meridional circulation, produces a very cold summer mesopause and reverse the zonal wind jets above about 85 km. The seasonal variation of the turbulent diffusion coefficient is consistent with the behavior of mesospheric turbulences inferred from MST radar echoes. The large degree of consistency between model results and various types of dynamical and chemical data supports very strongly the hypothesis that breaking gravity waves play a major role in determining the zonally-averaged dynamical and chemical structure of the 60 to 110 km region of the atmosphere.

Garcia, R. R.↗

Reducing dielectric loss and enhancing electrical insulation for multilayer polymer films by nanoconfined ion transport under high poling electric fields

High temperature polar polymers have demonstrated potential for good thermal stability and high dielectric constant at the same time. However, polarization of contaminated impurity ions in polar polymers, even at the ppm level, can significantly increase the dielectric loss at high temperature and low frequencies. One effective strategy to mitigate this problem is to multilayer them with a high temperature nonpolar dielectric polymer to confine impurity ion transport at the nanometer scale. In this study, confined ion transport in high temperature polycarbonate (HTPC)/poly(vinylidene fluoride) (PVDF) multilayer films under high AC electric fields was studied using a direct analytical simulation method. Different from the ion transport under low fields, the ion diffusion model failed to describe the ion transport under high electric fields. An exponential ion distribution profile, which was observed for the DC poling situation, was employed to implement the direct analytical simulation. Confined impurity ion transport under high AC fields was quantitatively understood. As the AC field increased, the mobile ion concentration decreased whereas the diffusion coefficient increased. The decrease of mobile ion concentration was explained by the blockage of impurity ions by the HTPC layers. This knowledge helped in the determination of optimal conditions to polarize impurity ions from the PVDF layers into the HTPC layers. After cooling below the glass transition temperature of HTPC, polarized impurity ions were locked inside the HTPC layers. As a result, increased discharge efficiency and enhanced electrical insulation (i.e., increased dielectric breakdown strength) were achieved for the polarized multilayer films.

36 MATERIALS SCIENCE↗

Ionic Peltier effect in Li-ion electrolytes

The coupled transport of charge and heat provide fundamental insights into the microscopic thermodynamics and kinetics of materials. We describe a sensitive ac differential resistance bridge that enables measurements of the temperature difference on two sides of a coin cell with a resolution of better than 10 μK. We use this temperature difference metrology to determine the ionic Peltier coefficients of symmetric Li-ion electrochemical cells as a function of Li salt concentration, solvent composition, electrode material, and temperature. The Peltier coefficients Π are negative, i.e., heat flows in the direction opposite to the drift of Li ions in the applied electric field, large, –Π > 30 kJ mol –1 , and increase with increasing temperature at T > 300 K. The Peltier coefficient is approximately constant on time scales that span the characteristic time for mass diffusion across the thickness of the electrolyte, suggesting that heat of transport plays a minor role in comparison to the changes in partial molar entropy of Li at the interface between the electrode and electrolyte. Here, our work demonstrates a new platform for studying the non-equilibrium thermodynamics of electrochemical cells and provides a window into the transport properties of electrochemical materials through measurements of temperature differences and heat currents that complement traditional measurements of voltages and charge currents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Construction of generalized quasilinear diffusion coefficient using neural networks with physical restrictions

The quasilinear diffusion coefficient (D QL ) derived from our machine learning framework shows comparable trends with the ground truth D QL obtained from GENRAY-CQL3D simulations. Additionally, for the strong absorption cases, the radial current drive profiles generated using the D QL from our model exhibit consistent behavior with those obtained from the original simulation. These findings indicate the potential of our surrogate modeling approach with physical restrictions to replicate key wave–plasma interaction characteristics while reducing computational costs. Traditionally, calculating D QL for wave–particle interactions relies on computationally intensive wave simulations coupled with Fokker–Planck solvers. To address this challenge, we developed a machine learning-based surrogate model with physical restrictions derived from cold plasma theory and bounce-averaged damping effects. First, we establish the propagation domain of Lower Hybrid Waves in the (N∥, ρ) space by identifying the accessibility limit and determining the upper and lower bounds of N∥ using the Potential Power Deposition (PPD) method. Subsequently, leveraging a database constructed using Latin hypercube sampling alongside the underlying physical restrictions (e.g. PPD), machine learning methods including U-Net and Recurrent Neural Networks are employed to design a physics-restricted machine learning framework capable of reconstructing D QL .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The altitude distribution of the Venus ultraviolet nightglow and implications on vertical transport

The altitude distribution of the nitric oxide nightglow was measured with an ultraviolet spectrometer on board Pioneer Venus, in order to study the effects of the distribution on the Venus nightside lower thermosphere transport properties. Limb profiles were obtained with an 8 ms integration period on several orbits near periapsis. The observations were made between P minus 2 min and P plus 4 min, where altitude ranges between 150 and 350 km, and latitude varies from 24 degrees N to 9 degrees S. A method independent of the spacecraft attitude data was used to fit the observed limb profiles, and to find the altitude of the maximum of the layer (115 plus or minus 2 km), and the topside scale height (about 3 km). It is shown that downward transport by diffusion alone is not sufficient, and if vertical motion is parameterized by eddy diffusion, an eddy diffusion coefficient is deduced from the altitude of the layer.

Gerard, J. C.↗

Wind-US Results for the AIAA 2nd Propulsion Aerodynamics Workshop

This presentation contains Wind-US results presented at the 2nd Propulsion Aerodynamics Workshop. The workshop was organized by the American Institute of Aeronautics and Astronautics, Air Breathing Propulsion Systems Integration Technical Committee with the purpose of assessing the accuracy of computational fluid dynamics for air breathing propulsion applications. Attendees included representatives from government, industry, academia, and commercial software companies. Participants were encouraged to explore and discuss all aspects of the simulation process including the effects of mesh type and refinement, solver numerical schemes, and turbulence modeling. The first set of challenge cases involved computing the thrust and discharge coefficients for a 25deg conical nozzle for a range of nozzle pressure ratios between 1.4 and 7.0. Participants were also asked to simulate two cases in which the 25deg conical nozzle was bifurcated by a solid plate, resulting in vortex shedding (NPR=1.6) and shifted plume shock (NPR=4.0). A second set of nozzle cases involved computing the discharge and thrust coefficients for a convergent dual stream nozzle for a range of subsonic nozzle pressure ratios. The workshop committee also compared the plume mixing of these cases across various codes and models. The final test case was a serpentine inlet diffuser with an outlet to inlet area ratio of 1.52 and an offset of 1.34 times the inlet diameter. Boundary layer profiles, wall static pressure, and total pressure at downstream rake locations were examined.

Wind-US↗

Electric Field Effects on Water and Ion Structure and Diffusion at the Orthoclase (001)–Water Interface

Understanding the electrochemical properties of mineral–water interfaces tends to rely upon electrical double layer (EDL) models, but these models are based on the assumption that electrostatic equilibrium is constantly maintained. In reality, interfacial reactions, ion diffusion, and their electrochemical signatures are based in nonequilibrium conditions of locally or globally imbalanced electrical fields where current EDL models have limited purview. In this work, we performed molecular dynamics (MD) simulations of the orthoclase (001) surface in contact with a 1 M NaCl aqueous solution under various electric fields, to explore the interplay between EDL structure and dynamics when perturbed by electric fields of different direction and strength, by confinement, and by different distributions of structural surface charge. The simulations showed that confinement between two opposing (001) surfaces led to the development of an induced field when the applied field was perpendicular to the surfaces and, as a result, to ionic diffusion coefficients that were independent of electric field strength. In contrast, when the applied field was parallel to the surfaces, confinement resulted in ionic diffusion coefficients that were more strongly dependent on the magnitude of the electric field than in bulk water. Differences in the density and distribution of aluminol groups on the two surfaces had a significant impact on how the interfacial structure and dynamics varied in the presence of an electric field. Notably, these differences resulted in an electro-osmotic flow with opposite directions at the two surfaces under parallel applied electric field. Overall, the MD simulations highlighted the importance of considering atomic-level structure and heterogeneities when developing models of the electrochemical properties of mineral–water interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗

Transport coefficients for magnetic-field evolution in inviscid magnetohydrodynamics

The magnetized resistivity and electrothermal tensors when substituted into the induction equation lead to electrothermal magnetic field generation, resistive magnetic diffusion, and magnetic field advection due to resistivity gradients, temperature gradients and currents. Here, the advection terms driven by temperature gradient and current have cross-field components (perpendicular to both the magnetic field and the driving term) that depend on significantly modified versions of Braginskii’s transport coefficients. The improved fits to Braginskii’s coefficients given by Epperlein and Haines and Ji and Held give physically incorrect results for cross-field advection at small Hall parameters (product of cyclotron frequency and collision time). The errors in Epperlein and Haines’ fits are particularly severe, giving increasing advection velocities below a Hall parameter of one when they should decrease linearly to zero. Epperlein and Haines’ fits can also give erroneous advection terms due to variations in effective atomic number. The only serious error in Braginskii’s fits is an overestimate in advection due to perpendicular resistivity. New fits for the cross-field advection terms are obtained from a direct numerical solution of the Fokker-Planck equation and Ji and Held’s higher order expansion approach that are continuous functions of the effective atomic number.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Dripping Handrail Model: Transient Chaos in Accretion Systems

We define and study a simple dynamical model for accretion systems, the "dripping handrail" (DHR). The time evolution of this spatially extended system is a mixture of periodic and apparently random (but actually deterministic) behavior. The nature of this mixture depends on the values of its physical parameters - the accretion rate, diffusion coefficient, and density threshold. The aperiodic component is a special kind of deterministic chaos called transient chaos. The model can simultaneously exhibit both the quasiperiodic oscillations (QPO) and very low frequency noise (VLFN) that characterize the power spectra of fluctuations of several classes of accretion systems in astronomy. For this reason, our model may be relevant to many such astrophysical systems, including binary stars with accretion onto a compact object - white dwarf, neutron star, or black hole - as well as active galactic nuclei. We describe the systematics of the DHR's temporal behavior, by exploring its physical parameter space using several diagnostics: power spectra, wavelet "scalegrams," and Lyapunov exponents. In addition, we note that for large accretion rates the DHR has periodic modes; the effective pulse shapes for these modes - evaluated by folding the time series at the known period - bear a resemblance to the similarly- determined shapes for some x-ray pulsars. The pulsing observed in some of these systems may be such periodic-mode accretion, and not due to pure rotation as in the standard pulsar model.

Young, Karl↗

Reaction diffusion in the NiCrAl and CoCrAl systems

The paper assesses the effect of overlay coating and substrate composition on the kinetics of coating depletion by interdiffusion. This is accomplished by examining the constitution, kinetics and activation energies for a series of diffusion couples primarily of the NiCrAl/Ni-10Cr or CoCrAl/Ni-10Cr type annealed at temperatures in the range 1000-1205 C for times up to 500 hr. A general procedure is developed for analyzing diffusion in multicomponent multiphase systems. It is shown that by introducing the concept of beta-source strength, which can be determined from appropriate phase diagrams, the Wagner solution for consumption of a second phase in a semiinfinite couple is successfully applied to the analysis of MCrAl couples. Thus, correlation of beta-recession rate constants with couple composition, total and diffusional activation energies, and interdiffusion coefficients are determined.

Levine, S. R.↗

Effect of Sulfonation Level on the Percolated Morphology and Proton Conductivity of Hydrated Fluorine-Free Copolymers: Experiments and Simulations

Using all-atom molecular dynamics simulations and a variety of experimental methods, we previously reported on a linear polyethylene with pendant phenyl sulfonated groups precisely on every fifth carbon along the backbone. With increasing relative humidity this fluorine-free polymer self-assembled to form nanoscale water channels and exhibited exceptional proton conductivity. Expanding upon those findings, here we explore partially sulfonated random copolymers, referred to as p 5PhSH-Y. Using either acetyl sulfate or sulfuric acid, a wide range of sulfonation levels were prepared ( Y = 34−98%) corresponding to ion-exchange capacities (IEC) of 2.0−4.4 mmol/g. Combining experimental techniques and all-atom molecular dynamics simulations, we study the effect of Y on water uptake, nanoscale morphology, and the proton/water transport properties of p5PhSH- Y . The proton conductivity of p 5PhSH- Y increases with relative humidity and with Y and achieves values in excess of 0.1 S/cm. These high conductivities are attributed to high IEC and welldeveloped nanoscale percolated hydrophilic domains made possible by the flexible backbone. We quantitatively describe the nature of the water channels using the characteristic distance, channel width distribution, the area per sulfonate group at the hydrophilic/ hydrophobic interface, and the fractal dimension. Notably, the channel widths and the areas per sulfonate group are nominally independent of the level of sulfonation, while depending significantly on the level of hydration. The fractal dimension of the water channels correlates strongly with the water diffusion coefficients calculated from the molecular dynamics (MD) simulations. These findings demonstrate that the p 5PhSH- Y hydrocarbon copolymers can be modified to tune properties, particularly proton conductivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Full-scale testing of an Ogee tip rotor

Full scale tests were utilized to investigate the effect of the ogee tip on helicopter rotor acoustics, performance, and loads. Two facilities were used: the Langley whirl tower and a UH-1H helicopter. The text matrix for hover on the whirl tower involved thrust values from 0 to 44 480 N (10,000 lb) at several tip Mach numbers for both standard and Ogee rotors. The full scale testing on the UH-1H encompassed the major portion of the flight envelope for that aircraft. Both near field acoustic measurements and far field flyover data were obtained for both the ogee and standard rotors. Data analysis of the whirl tower test shows that the ogee tip does significantly diffuse the tip vortex while providing some improvement in hover performance at low and moderate thrust coefficients. Flight testing of both rotors indicates that the strong impulsive noise signature of the standard rotor can be reduced with the ogee rotor. Analysis of the spectra indicates a reduction in energy in the 250 Hz and 1000 Hz range for the ogee rotor. Forward flight performance was significantly improved with the ogee configuration for a large number of flight conditions. Further, rotor control loads were reduced through use of this advanced tip rotor.

Mantay, W. R.↗

Effect of area ratio on the performance of a 5.5:1 pressure ratio centrifugal impeller

A centrifugal impeller which was initially designed for a pressure ratio of approximately 5.5 and a mass flow rate of 0.959 kg/sec was tested with a vaneless diffuser for a range of design point impeller area ratios from 2.322 to 2.945. The impeller area ratio was changed by successively cutting back the impeller exit axial width from an initial value of 7.57 mm to a final value of 5.97 mm. In all, four separate area ratios were tested. For each area ratio a series of impeller exit axial clearances was also tested. Test results are based on impeller exit surveys of total pressure, total temperature, and flow angle at a radius 1.115 times the impeller exit radius. Results of the tests at design speed, peak efficiency, and an exit tip clearance of 8 percent of exit blade height show that the impeller equivalent pressure recovery coefficient peaked at a design point area ratio of approximately 2.748 while the impeller aerodynamic efficiency peaked at a lower value of area ratio of approximately 2.55. The variation of impeller efficiency with clearance showed expected trends with a loss of approximately 0.4 points in impeller efficiency for each percent increase in exit axial tip clearance for all impellers tested.

Schumann, L. F.↗

Effect of area ratio on the performance of a 5.5:1 pressure ratio centrifugal impeller

A centrifugal impeller which was initially designed for a pressure ratio of approximately 5.5 and a mass flow rate of 0.959 kg/sec was tested with a vaneless diffuser for a range of design point impeller area ratios from 2.322 to 2.945. The impeller area ratio was changed by successively cutting back the impeller exit axial width from an initial value of 7.57 mm to a final value of 5.97 mm. In all, four separate area ratios were tested. For each area ratio a series of impeller exit axial clearances was also tested. Test results are based on impeller exit surveys of total pressure, total temperature, and flow angle at a radius 1.115 times the impeller exit radius. Results of the tests at design speed, peak efficiency, and an exit tip clearance of 8 percent of exit blade height show that the impeller equivalent pressure recovery coefficient peaked at a design point area ratio of approximately 2.748 while the impeller aerodynamic efficiency peaked at a lower value of area ratio of approximately 2.55. The variation of impeller efficiency with clearance showed expected trends with a loss of approximately 0.4 points in impeller efficiency for each percent increase in exit axial tip clearance for all impellers tested.

Schumann, L. F.↗