SEARCH · Engineering Papers
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.
We present and discuss the general case of a fractional nonlinear partial differential equation using similarity reductions and recover results associated with Harry Dym-type equations. Particular cases are discussed and some plots illustrate the results. The fractional derivative is considered in the Caputo sense.
Incompressible two dimensional calculations are reported for the impulsively started lid driven cavity with aspect ratio two. The algorithm is based on the time dependent streamfunction equation, with a Crank-Nicolson differencing scheme for the diffusion terms, and with an Adams-Bashforth scheme for the convection terms. A multigrid method is used to solve the linear implicit equations at each time step. Periodic asymptotic solutions have been found for Re = 10000 and for Re = 5000. The Re = 5000 results are validated by grid refinement calculations. The solutions are shown to be precisely periodic, and care is taken to demonstrate that asymptotic states were reached. A discussion is included about the indicators that are used to show that an asymptotic state was reached, and to show that the asymptotic state is indeed periodic.
Next, a robust and efficient framework for simulating NSPD in multiple dimensions is developed. The reactive Navier-Stokes equations are extended to include a drift-diffusion plasma-fluid model with a local field approximation (LFA) in a finite-volume solver, which uses an adaptive mesh refinement (AMR) strategy to address the wide separation of length scales in the problem. A two-way coupling strategy is used whereby the plasma-fluid model and reactive Navier-Stokes equations are integrated simultaneously. The oxidation of ethylene/air mixtures mediated by NSPD is simulated in a pin-to-pin configuration. All phases of the plasma discharge are simulated explicitly (including streamer ignition, propagation, and connection, as well as the subsequent spark phase), along with the evolution of the plasma during the inter-pulse period. Temporally and spatially-resolved results are presented, with an emphasis on the analysis of heating and energy deposition, as well as of the evolution of the concentration of active particles generated during the NSPD and their influence on ignition.
In biomembranes lipid mobility exhibits deviations from the classical diffusive behavior of Brownian particles, i.e. ‘anomalous’ diffusion. The question arises as to how this anomalous diffusive behavior varies in gel, ripple, and fluid biomembrane phases. Here, in this study, we perform all-atom molecular dynamics (MD) simulations of dimyristoylphosphatidylcholine (DMPC) bilayers in the three different phases and analyze the results using the framework of the Generalized Langevin Equation (GLE). This analysis emphasizes sub-diffusive behavior on the relatively short, picosecond-nanosecond timescales, capturing local molecular constraints and transient caging effects during the crossover of atomic dynamics from vibrational to incipient anharmonic motion. The ripple and gel phases are found to exhibit strong transient caging and prolonged memory effects resulting in distinct sub-diffusive behavior. The role of hydrogen bonding in lipid confinement is also examined, demonstrating its influence on phase-dependent molecular ordering and on short-time diffusional constraints. These findings demonstrate the GLE framework’s utility in characterizing molecular transport and lipid dynamics, with implications for longer timescale membrane dynamics.
Canopy layers control momentum and solute transport to and from the overlying water surface layer. These transfer mechanisms strongly dependent on canopy geometry, affect the amount of solute in the river, the hydrological retention and availability of dissolved solutes to organisms located in the vegetated layers, and are critical to improve water quality. In this work, we consider steady state transport in a vegetated channel under fully developed flow conditions. Under the hypothesis that the canopy layer can be described as an effective porous medium with prescribed properties, i.e., porosity and permeability, we model solute transport above and within the vegetated layer with an advection-dispersion equation with a spatially variable dispersion coefficient (diffusivity). By means of the Generalized Integral Transform Technique , we derive a semianalytical solution for the concentration field in submerged vegetated aquatic systems. We show that canopy layer’s permeability affects the asymmetry of the concentration profile, the effective vertical spreading behavior, and the magnitude of the peak concentration. Due to its analytical features, the model has a low computational cost. The proposed solution successfully reproduces previously published experimental data.
This paper presents an approach for characterizing and estimating statistical dependence between a large number of observables in a composite material system. Conditional regression is carried out using the estimated joint density function, permitting a systematic exploration of interdependence between fine scale and coarse observables that can be used for both prognosis and design of complex material systems. An example demonstrates the integration of experimental data with a computational database. The statistical approach is based on the probabilistic learning on manifolds recently developed by the authors. This approach leverages intrinsic structure detected through diffusion on graphs with projected stochastic differential equations to generate samples constrained to that structure.
We study how the mass transfer of a chemical species in shear flow is suppressed by the production of a second species that partially reacts away the first species. The second species is produced at the surface with the first species as a reactant. Our work is directly motivated by electrochemical CO 2 reduction in aqueous flow reactors, where OH– molecules generated by the CO evolution reaction react away CO 2 molecules, ultimately inhibiting the mass transfer of CO 2 to the cathode surface. We derive a simple approximation of the Sherwood number, a dimensionless measure of the mass flux into the surface, as a function of the Péclet number, surface Damköhler number, and bulk Damköhler number, and we compare our approximation to numerical solutions of the governing equations. We find that in the diffusion-limited regime, the Sherwood number is well described by the classical Graetz-Lévêque result reduced by a reaction factor due to the competing bulk reaction; with the stoichiometry relevant to electrochemical CO 2 reduction, this reaction factor is 1/2. While the model problem we solve provides insight into how OH– production affects CO 2 mass transfer, experimental systems often involve more complex chemistry. We thus also show how a common buffered electrolyte, KHCO 3 , affects our results.
In this paper, we systematically derive a model for turbulence spreading from the basic kinetic equation. The model contains explicit nonlocal nonlinear diffusion and nonlocal growth. When the nonlocality scale parameter δ b (banana width) vanishes, this model reduces to the usual turbulence spreading model. We elucidate the mechanisms of nonlinear saturation and nonlocal growth. Results show that nonlocal effects, especially the nonlocal growth, thicken the turbulence spreading front and increase the speed of front propagation. More turbulence intensity penetrates the stable region when δ b increases. The penetration depth Δ p is proportional to $ \delta_b/L_T $, therefore the fraction of turbulence in the unstable region scales as 1 – δ b* . Lastly, the transport scales the same way.
Quantum fluid droplets made of helium-3 ( 3 He) or helium-4 ( 4 He) isotopes have long been considered as ideal cryogenic nanolabs, enabling unique ultracold chemistry and spectroscopy applications. The droplets were believed to provide a homogeneous environment in which dopant atoms and molecules could move and react almost as in free space but at temperatures close to absolute zero. Here, we report ultrafast x-ray diffraction experiments on xenon-doped 3 He and 4 He nanodroplets, demonstrating that the unavoidable rotational excitation of isolated droplets leads to highly anisotropic and inhomogeneous interactions between the host matrix and enclosed dopants. Superfluid 4 He droplets are laced with quantum vortices that trap the embedded particles, leading to the formation of filament-shaped clusters. In comparison, dopants in 3 He droplets gather in diffuse, ring-shaped structures along the equator. The shapes of droplets carrying filaments or rings are direct evidence that rotational excitation is the root cause for the inhomogeneous dopant distributions.
The nano/micro-scale geometry of polymer electrolyte fuel cell (PEFC) catalyst layers critically affects cell performance. The small length scales and complex structure of these composite layers make it challenging to analyze cell performance and physics at the particle scale by experiment. We present a computational method to simulate transport and chemical reaction phenomena at the pore/particle-scale and apply it to a PEFC cathode with platinum group metal free (PGM-free) catalyst. Here, we numerically solve the governing equations for the physics with heterogeneous oxygen diffusion coefficient and proton conductivity evaluated using the actual electrode structure and ionomer distribution obtained using nano-scale resolution X-ray computed tomography (nano-CT). Using this approach, the oxygen concentration and electrolyte potential distributions imposed by the oxygen reduction reaction are solved and the impact of the catalyst layer structure on performance is evaluated.
Transient gas-flow process in multilayer insulation systems during evacuation predicted by equations in conjunction with measured permeabilities and diffusion coefficients
Laminar flame propagation equations in stagnation flow using flat flame diffusion burner
Pitch angle diffusion of trapped electrons in terrestrial radiation zones, discussing diffusion theory based on Fokker-Planck equation, atmospheric Coulomb scattering mechanism, etc
Near IR radiation diurnal variations in corn canopy accounted for by Duntley equations for unidirectionally incident light propagation through diffusing medium
The influence of molecular transport is included in the computation by treating viscous and thermal diffusion terms in the governing partial differential equations as correction terms in the method of characteristics scheme. The development of a production type computer program is reported which is capable of calculating the flow field in a variety of axisymmetric mixed-compression aircraft inlets. The results agreed well with those produced by the two-dimensional method characteristics when axisymmetric flow fields are computed. For three-dimensional flow fields, the results agree well with experimental data except in regions of high viscous interaction and boundary layer removal.
Using the finite difference method with overrelaxation, numerical solutions of the steady-state vorticity transport equation were obtained for a continuous flow diffusion chamber of the Hudson-Squires type. The calculation neglected the effects due to temperature, gravity, and saturation. The size and shape of the manifold used to inject the aerosol laden flow were varied to obtain a design which would improve the performance of the chamber from strictly low Reynolds number (less than 20) fluid dynamical considerations.
Experimental data have been obtained in an incompressible turbulent flow over a rearward-facing step in a diverging channel flow. Mean velocities, Reynolds stresses, and triple products that were measured by a laser Doppler velocimeter are presented for two cases of tunnel wall divergence. Eddy viscosities, production, convection, turbulent diffusion, and dissipation (balance of kinetic energy equation) terms are extracted from the data. These data are compared with various eddy-viscosity turbulence models. Numerical calculations incorporating the k-epsilon and algebraic-stress turbulence models are compared with the data. When determining quantities of engineering interest, the modified algebraic-stress model (ASM) is a significant improvement over the unmodified ASM and the unmodified k-epsilon model; however, like the others, it dramatically overpredicts the experimentally determined dissipation rate.
Simulations of the electron acoustic instability for a finite-size electron beam system are performed with a particle-in-cell code to investigate the heating phenomena associated with the instability and the width of the heating region. The results show that the beam radiates electrostatic electron acoustic waves whose decay time outside the beam agrees with the spatial decay length derived from the linear dispersion equation. The ambient cold electrons in a diffusion region surrounding the beam are heated to a higher temperature by absorbing the radiated electron acoustic waves, with the heating occurring mainly in the parallel direction. In the heat diffusion region, the temperature of the cold electrons decreases with distance from the beam with a temperature gradient length smaller than the decay length of the wave energy. The results are discussed with respect to the DE 1 plasma and wave observations in the polar cusp region.