SEARCH · Search NASA
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.
Local emissivity of vertex point of diffuse conical or V-groove cavity, noting dependence on wall material, angle factor and agreement with theoretical solution
We study the dynamical aspects of dark energy in the context of a non-minimally coupled scalar field with curvature and torsion. Whereas the scalar field acts as the source of the trace mode of torsion, a suitable constraint on the torsion pseudo-trace provides a mass term for the scalar field in the effective action. In the equivalent scalar-tensor framework, we find explicit cosmological solutions representing dark energy in both Einstein and Jordan frames. We demand the dynamical evolution of the dark energy to be weak enough, so that the present-day values of the cosmological parameters could be estimated keeping them within the confidence limits set for the standard LCDM model from recent observations. For such estimates, we examine the variations of the effective matter density and the dark energy equation of state parameters over different redshift ranges. In spite of being weakly dynamic, the dark energy component differs significantly from the cosmological constant, both in characteristics and features, for e.g. it interacts with the cosmological (dust) fluid in the Einstein frame, and crosses the phantom barrier in the Jordan frame. We also obtain the upper bounds on the torsion mode parameters and the lower bound on the effective Brans-Dicke parameter. The latter turns out to be fairly large, and in agreement with the local gravity constraints, which therefore come in support of our analysis.
A data-driven ab initio generalized Langevin equation (AIGLE) approach is developed to learn and simulate high-dimensional, heterogeneous, coarse-grained (CG) conformational dynamics. Constrained by the fluctuation–dissipation theorem, the approach can build CG models in dynamical consistency (DC) with all-atom molecular dynamics. Here, we also propose practical criteria for AIGLE to enforce long-term DC. Case studies of a toy polymer, with 20 CG sites, and the alanine dipeptide, with two dihedral angles, elucidate why one should adopt AIGLE or its Markovian limit for modeling CG conformational dynamics in practice.
Energetic electron intensities radial flow diffusion in outer radiation belt near geomagnetic equator
We present a family of multi-stage preconditioners for coupled thermal-compositional-reactive reservoir simulation problems. The most common preconditioner used in industrial practice, the Constrained Pressure Residual (CPR) method, was designed for isothermal models and does not offer a specific strategy for the energy equation. For thermal simulations, inadequate treatment of the temperature unknown can cause severe conver-gence degradation. When strong thermal diffusion is present, the energy equation exhibits significant elliptic behavior that cannot be accurately corrected by CPR's second stage. In this work, we use Schur-complement decompositions to extract a temperature subsystem and apply an Algebraic MultiGrid (AMG) approximation as an additional preconditioning stage to improve the treatment of the energy equation. We present results for several two-dimensional hot air injection problems using an extra heavy oil, including challenging reactive In-Situ Combustion (ISC) cases. We show improved performance and robustness across different thermal regimes, from advection dominated (high Péclet number) to diffusion dominated (low Péclet number). The number of linear iterations is reduced by 40–85% compared to standard CPR for both homogeneous and heterogeneous media, and the new methods exhibit almost no sensitivity to the thermal regime.
An existing method which solves the multi-dimensional ensemble-averaged compressible time-dependent Navier-Stokes equations in conjunction with mixing length turbulence model and shock capturing technique has been extended to include the shock-tracking adaptive grid systems. The numerical scheme for solving the governing equations is based on a linearized block implicit approach. The effects of grid-motion and grid-distribution on the calculated flow solutions have been studied in relative detail and this is carried out in the context of physically steady, shocked flows computed with non-stationary grids. Subsequently, the unsteady dynamics of the flows occurring in a supercritically operated transonic diffuser and a mixed compression supersonic inlet have been investigated with the adaptive grid systems by solving the Navier-Stokes equations.
The validity of imposing a zero chemical potential for neutrinos in hydrodynamic calculations of collapsing supernovae is investigated in the diffusion approximation of neutrino transport. A coupled system of equations is solved for neutrino and energy diffusion fluxes as well as lepton diffusion in a collapsing supernovae ambient medium, and the results indicate a substantial growth in the neutrino chemical potential for densities greater than 10 to the 12th power gm/cu cm. The rate of energy transport is shown to be significantly affected by increases in Fermi integrals and chemical-potential gradients accompanied by decreases in temperature, and the extent of neutrino particle/antiparticle reactions is found also to affect energy diffusion rates. It is concluded that the photon-like behavior usually assumed for neutrinos may be incorrect and that an extension of the Sn transport approximation to include lepton characteristics is necessary for a definitive answer to the question of neutrino transport in supernovae.
We propose to re-interpret the cosmic microwave background spectral distortions as solutions to the Boltzmann equation. This approach makes it possible to solve the second order Boltzmann equation explicitly, with the spectral y distortion and the momentum independent second order temperature perturbation, while generation of μ distortion cannot be explained even at second order in this framework. We also extend our method to higher order Boltzmann equations systematically and find new type spectral distortions, assuming that the collision term is linear in the photon distribution functions, namely, in the Thomson scattering limit. As an example, we concretely construct solutions to the cubic order Boltzmann equation and show that the equations are closed with additional three parameters composed of a cubic order temperature perturbation and two cubic order spectral distortions. The linear Sunyaev-Zel'dovich effect whose momentum dependence is different from the usual y distortion is also discussed in the presence of the next leading order Kompaneets terms, and we show that higher order spectral distortions are also generated as a result of the diffusion process in a framework of higher order Boltzmann equations. The method may be applicable to a wider class of problems and has potential to give a general prescription to non-equilibrium physics.
Heat capacity expression due to chemical reaction, relating composition variations and binary diffusion coefficient with Lewis number, presenting integral equation for system enthalpy
We calculate the full probability density function (PDF) of inflationary curvature perturbations, even in the presence of large quantum backreaction. Making use of the stochastic-δ N formalism, two complementary methods are developed, one based on solving an ordinary differential equation for the characteristic function of the PDF, and the other based on solving a heat equation for the PDF directly. In the classical limit where quantum diffusion is small, we develop an expansion scheme that not only recovers the standard Gaussian PDF at leading order, but also allows us to calculate the first non-Gaussian corrections to the usual result. In the opposite limit where quantum diffusion is large, we find that the PDF is given by an elliptic theta function, which is fully characterised by the ratio between the squared width and height (in Planck mass units) of the region where stochastic effects dominate. We then apply these results to the calculation of the mass fraction of primordial black holes from inflation, and show that no more than ∼ 1 e -fold can be spent in regions of the potential dominated by quantum diffusion. We explain how this requirement constrains inflationary potentials with two examples.
Cardiovascular implantable devices alter the biofluid dynamics and biochemistry of the blood in which they are placed. These perturbations can lead to thrombus formation which may or may not be desired, depending on the application. In this work, a computational model is developed that couples biofluid dynamics and biochemistry to predict the clotting response of blood to such devices. The model consists of 28 advection–diffusion–reaction partial differential equations to track proteins in the blood involved in clotting and utilizes boundary flux terms to model the initiation of the intrinsic clotting pathway at thrombogenic device surfaces. We use this model to simulate the transient clot growth within a 2D idealized bifurcation aneurysm filled with various distributions of bare metal coils with similar packing densities. The clot model predicts initial clot formation to occur in areas along coil surfaces where flow is minimal and where time-averaged shear rates are the smallest. Among the six coil-filled aneurysm cases simulated, maximum thrombus occlusion ranged between 80.8 and 92.2% of the post-treatment aneurysm volume and was achieved 325–450 s after treatment. Finally, with further refinement and validation, the computational clotting model will be a valuable engineering tool for evaluating and comparing the relative performance of cardiovascular implantable devices.
Protein rotation in viscous environments can be measured by fluorescence depletion anisotropy (FDA) which combines long lifetimes of chromophore triplet states with the sensitivity of fluorescence excitation and detection. FDA achieves sensitivity well beyond that attainable by the more common technique of time-resolved phosphorescence anisotropy (TPA). We have now combined benefits of both time-domain and frequency-domain FDA into a single continuous technique (CFDA). Intensity and polarization of a single laser beam are modulated continuously according to a complex, repeating waveform. Fluorescence signals excited from triplet-forming fluorescent probes are digitized over recurring waveform periods by a high-speed signal averager. CFDA experiments typically involve substantial ground state depletion. Thus signals, unlike those of TPA, are not linear in the exciting light intensity and simple data analysis based on such linearity is not appropriate. An exact solution of the coupled diffusion and triplet production/decay equation describing CFDA within individual data points has been combined with simulated annealing optimization to extract triplet and anisotropy decay kinetics from experimental data. Related calculations compare possible excitation waveforms with respect to rotational information provided per fluorescence photon. We present CFDA results for the model system of eosin conjugates of carbonic anhydrase, BSA and immunoglobulin G in 90% glycerol at various temperatures and initial cellular results on eosin-IgE bound to 2H3 cell Type I Fcε receptors. We explore how CFDA reflects rotational parameters of heterogeneous systems and discuss challenges of extending this method to single cell microscopic measurements.
This research article describes, an eco-friendly activated carbon prepared from the Gracilaria corticata seaweeds which was employed for the preparation of biodegradable polymeric beads for the efficient removal of crystal violet dye in an aqueous solution. The presence of chemical functional groups in the adsorbent material was detected using FTIR spectroscopy. The morphology and physical phases of the adsorbent materials were analyzed using SEM and XRD studies respectively. Batch mode dye adsorption behavior of the activated carbon/Zn/alginate polymeric beads was investigated as a function of dosage, solution pH, contact time, initial dye concentration and temperature. Maximum dye removal was observed at a pH of 5.0, 1 g of adsorbent dosage with 60 mg/L dye concentration, 50 min of contact time and at 30 °C. The equilibrium modeling studies were analyzed with Freundlich and Langmuir adsorption isotherms and the adsorption dynamics was predicted with Lagergren’s pseudo-first order, pseudo-second order equations and intra particle diffusion models. The process of dye removal followed a pseudo second-order kinetics rather than pseudo first order. The thermodynamic parameters like standard Gibbs free energy change (ΔG°), enthalpy change (ΔH°) and entropy change (ΔS°) were determined and the results imply that the adsorption process was spontaneous, endothermic and increases the randomness between the adsorbent and adsorbate. The results from the experimental and correlation data reveal that the Gracilaria corticata activated carbon/Zn/alginate polymeric beads have proved to be an excellent adsorbent material for the removal of CV dye.
Diffusion-driven charge transport (DDCT) in III–V light-emitting diodes (LEDs) can enable unconventional optoelectronic devices and functionality by fundamentally changing device design and the current injection principle. In our recent study, an AlGaAs/GaAs DDCT–LED consisting of an array of lateral heterojunctions was studied for large-area applications at high powers. Here, we investigate the current spreading and recombination uniformity of a modulation doped GaInP/GaAs DDCT–LED. In particular, we analyze how the background doping of the lower GaInP cladding layer (CL) and the GaAs substrate changes the carrier distribution within the active region of the device. Our charge transport simulations based on the drift-diffusion current and continuity equations predict that modulation doping by a p-doped CL provides much higher recombination uniformity at high powers compared to an n-doped CL. Most importantly, improved current spreading is achieved while maintaining excellent device performance.
In deep subsurface environments, where porous rocks are saturated with fluids, hydromechanical coupling plays a critical role in the formation and propagation of localized deformation zones, such as compaction bands. In this study, we numerically investigate the influence of fluid-rock interactions on the development of compaction localization in fluid-saturated rocks. For this purpose, it uses a numerical framework that adapts the heat equation to model the diffusion of pore pressure caused by deformation. This is achieved with specialized subroutines for material behavior and heat transfer formulated within the ABAQUS platform and by allowing changes in rock permeability during strain localization. The results highlight that hydromechanical coupling significantly alters the compaction behavior compared to dry conditions. Under high confining stress, excess pore pressure inhibits the progression of compaction fronts, while in time-dependent processes like rock creep, pore pressure transients accelerate deformation by creating feedback loops between mechanical forces and fluid flow. These findings underscore the importance of saturation in controlling the dynamics of compaction localization and provide insights into the interplay between pore pressure, deformation, and permeability. The implications of this study extend to applications in geological storage, resource extraction, and subsurface geomechanics, offering a deeper understanding of coupled hydromechanical processes in porous rocks.
Ammonia is carbon-free and is regarded as a potential fuel to address global warming issues. Here, in this work, three-dimensional direct numerical simulations (DNS) of ammonia/air turbulent premixed flames were performed to explore the influence of turbulence and equivalence ratio on the flame structure and NO formation characteristics. Two equivalence ratios were considered, i.e. Φ=0.9 and Φ=1.1. The general flame structures were presented and species distributions were examined. The NO mass fraction was found to be the highest in the product for the lean case and in the reaction zone for the rich case. The conditional mean values of species mass fractions and reaction rates were compared with those of the unstrained and strained laminar premixed flames to explore how well the laminar flame structures can approach those of the turbulent flame. The budget analysis of the species transport equations showed that turbulent diffusion plays an important role in species transport. The turbulent diffusivity D T was estimated using the gradient hypothesis based on the DNS data. Various laminar flame simulations with different diffusivities were carried out. It was shown that the conditional means of the DNS agree well with those of the laminar flames including D T in the transport property calculation. The global and local NO formation characteristics were investigated. It was found that the mean NO production conditioned on the progress variable is lower compared with the corresponding laminar flame in the rich case. However, the relative contributions from various NO pathways are rarely affected by turbulence. The NO mass fraction is higher in negative curvature regions compared with positive curvature regions of the flame surface for the rich case, which is due to the preferential diffusion of H 2 and other radicals and the enhanced NO pathways in negatively curved regions.
Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.
Tabulated chemistry reduces the cost of modeling reacting flows by transporting fewer scalars. It also has the benefit of introducing fewer terms that must be closed for Reynolds-averaged Navier-Stokes and large eddy simulation modeling. In this study of premixed lean hydrogen–air flames, an extension to an existing tabulated chemistry model was developed by including differential diffusion and thermal diffusion of multiple species, using a mixture-averaged diffusion model. The transport equations for the tabulation variables (a progress variable and a mixture fraction-like variable) have been formulated, and curvature effects are intrinsic to the transport equations derived from first principles. The mathematical derivation is general, such that other fuel/air mixtures could be examined. The new tabulated chemistry model was then validated using a range of configurations, including flames with curvature, extinction, and hot spots. The model validation indicated excellent agreement with detailed chemistry simulations for lean premixed hydrogen–air mixtures. Conditional means of the progress variable source term were shown to agree well with detailed chemistry simulations of two-dimensional freely propagating flames and three-dimensional turbulent flames. Finally, progress variable source terms three times higher than one-dimensional flat flame results were identified and validated with detailed chemistry, and super-adiabatic hot spots were also predicted in the turbulent flame configuration.