Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Permeability Variations Associated With Fault Reactivation in a Claystone Formation Investigated by Field Experiments and Numerical Simulations

We studied the relation between rupture and changes in permeability within a fault zone intersecting the Opalinus Clay formation at 300 m depth in the Mont Terri Underground Research Laboratory (Switzerland). A series of water injection experiments were performed in a borehole straddle interval set within the damage zone of the main fault. A three-component displacement sensor allowed an estimation of the displacement of a minor fault plane reactivated during a succession of step rate pressure tests. The experiment reveals that the fault hydromechanical (HM) behavior is different from one test to the other with varying pressure levels needed to trigger rupture and different slip behavior under similar pressure conditions. Numerical simulations were performed to better understand the reason for such different behavior and to investigate the relation between rupture nucleation, permeability change, pressure diffusion, and rupture propagation. Here, our main findings are as follows: (i) a rate frictional law and a rate-and-state permeability law can reproduce the first test, but it appears that the rate constitutive parameters must be pressure dependent to reproduce the complex HM behavior observed during the successive injection tests; (ii) almost similar ruptures can create or destroy the fluid diffusion pathways; (iii) a too high or too low diffusivity created by the main rupture prevents secondary rupture events from occurring whereas “intermediate” diffusivity favors the nucleation of a secondary rupture associated with the fluid diffusion. However, because rupture may in certain cases destroy permeability, this succession of ruptures may not necessarily create a continuous hydraulic pathway.

58 GEOSCIENCES↗
NASA NTRS1993

A numerical analysis of high-temperature heat pipe startup from the frozen state

Continuum and rarefied vapor flows co-exist along the heat pipe length for most of the startup period. A two-region model is proposed in which the vapor flow in the continuum region is modeled by the compressible Navier-Stokes equations, and the vapor flow in the rarefied region is simulated by a self-diffusion model. The two vapor regions are linked with appropriate boundary conditions, and heat pipe wail, wick, and vapor flow are solved as a conjugate problem. The numerical solutions for the entire heat pipe startup process from the frozen state are compared with the corresponding experimental data with good agreement.

Cao, Y.↗
DOE OSTI2026

Carbon fiber design improvements based on economic models

Circular fiber geometries are predominant in commercial carbon fiber material systems, but the use of this fiber shape has numerous limitations. Circular geometries have continuous symmetry, which is helpful for various processing considerations, but also have the largest possible maximum diffusion thickness for a given fiber area. This characteristic means that circular carbon fibers always have the highest material processing cost and lowest production throughput compared to any other fiber shape with the same area and tow count. To quantify material cost and other benefits for non-circular carbon fiber geometries, process models for polyacrylonitrile based carbon fiber production are developed in relationship to the carbon fiber shape, size, and tow count. For a given fiber shape, precursor production costs are shown to favor maximizing fiber size while conversion costs are minimized by the smallest fiber size. These competing cost trends result in a numerically optimal fiber size for a given shape and tow count, while both cost components are decreased by increasing tow count. Cost–performance tradeoffs for three lobe fiber geometries are studied by supplementing the cost trends with a numerical failure model to predict compressive strength for discrete shape variants. The shape selection is shown to be more sensitive to variations in cost than compressive strength while suboptimal shape designs can improve manufacturing robustness and achievable fiber volume fractions. Finally, an optimal three lobe carbon fiber is identified that balances the set of considerations while reducing costs and embodied energy and increasing production throughput compared to a commercial carbon fiber.

Carbon fiber↗
DOE OSTI2018

Quantitative soil zymography: Mechanisms, processes of substrate and enzyme diffusion in porous media

Soil membrane zymography enables 2D mapping of enzyme activities on the surface of soil samples. The method is according to diffusion of components of enzymatically-mediated reactions to/from membrane, and, hence, reflects the distribution of enzyme activities at the intact soil surface. Zymography has been already successfully implemented in numerous soil ecology applications. Here we identify two methodological aspects for further improvement and expansion of the method at micro and macro scales: first, accounting for the area of contact between the soil surface and the zymography membranes and, second, accounting for diffusion effects during the zymography procedure. We tested three methods, namely, laser-scanning, staining with a fluorescent product (e.g. MUF: 4-methylumbelliferone), and X-ray computed micro-tomography, for assessing the area of the soil surface in contact with the membranes. We quantified diffusion of MUF, enzymes and substrate between the substrate-saturated membrane and soil as well as diffusion processes during membrane zymography via HP2 software. Diffusion of the substrate from the membrane and of the MUF-product to the membrane was detected, while there were no clear evidence of enzyme diffusion to/in the membrane. According to the model simulations, the enzyme activities detected via 2D zymography probably represent only a small portion, about 20%, of the actual reactions within the soil volume that is in both direct contact and in hydrological contact with zymography membranes. This is a result of omnidirectional diffusion of reaction products. The membrane contact with the soil surface estimated by three methods ranged from 3.4 to 36.5% further signifying that only a fraction of enzymes activity is detectable in a course of 2D soil zymography.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1994

On the cross-field diffusion of ions in one- and two-dimensional hybrid simulations of collisionless shocks

It can be demonstrated analytically that under certain geometries used in numerical simulations of collisionless shocks in which there is at least one ignorable spatial coordinate, the transport of particles across the magnetic field is essentially zero. This notion is tested using one- and two-dimensional hybrid simulations (kinetic ions/fluid electrons). We find, as the theorem predicts, the particles treated kinetically are tied to the same field line on which they start.

Giacalone, Joe↗
DOE OSTI2025

Generalizable Porous Aromatic Framework‐Included Polymer Membranes for Diffusion‐Enhanced Gas Separations

Industrial separation processes account for 10-15% of global energy consumption. Membrane-based processes are less energy-intensive than traditional gas separation technologies; however, enhanced material separation performance and stability for numerous gas mixtures are needed for widespread industrial adoption. This work presents a generalizable strategy for preparing mixed-matrix gas separation membranes exceeding the performance upper bounds of existing polymer membranes for a wide variety of industrial gases. By incorporating robust porous aromatic framework (PAF) particles into various dense commercial polymer matrices, gas diffusivity and solubility can be enhanced. For diverse gas mixtures (e.g., CO2/N2, O2/N2, He/CH4, H2/N2, and C2H4/C2H6), the resulting composite membranes exhibit enhanced gas permeabilities-by as much as 520%-and largely unchanged selectivities even after 6 years of aging under simulated flue gas conditions. These improvements arise from the ultrahigh porosity, excellent chemical compatibility, and unique physicochemical properties of the embedded PAF particles. Functionalizing the PAFs with polyamines also enables composite membranes that achieve among the highest reported performances against plasticization, a common obstacle in commercializing gas separation membranes. Significantly, the PAF-1 particles are readily dispersible in various common membrane casting solvents, suggesting their broader utility as a filler for designing high-performance membranes for many industrial gas separations.

Uliana, Adam A↗
NASA NTRS1986

Astrophysical data on 5 eV to 1 keV radiation from the radiative decay of fundamental particles - Current limits and prospects for improvement

Line emission from the decay of fundamental particles, integrated over cosmological distances, can give rise to detectable spectral features in the diffuse astronomical background between 5 eV and 1 keV. Spectroscopic observations may allow these features to be separated from line emission from the numerous local sources of radiation. The current observational status and existing evidence for such features are reviewed. No definitive detections of nongalactic line features have been made. Several local sources of background mask the features at many wavelengths and confuse the interpretation of the data. No systematic spectral observations have been carried out to date. Upcoming experiments which can be expected to provide significantly better constraints on the presence of spectral features in the diffuse background from 5 eV to 1 keV are reviewed.

Bowyer, Stuart↗
DOE OSTI2017

The equilibrium-diffusion limit for radiation hydrodynamics

The equilibrium-diffusion approximation (EDA) is used to describe certain radiation-hydrodynamic (RH) environments. When this is done the RH equations reduce to a simplified set of equations. The EDA can be derived by asymptotically analyzing the full set of RH equations in the equilibrium-diffusion limit. Here, we derive the EDA this way and show that it and the associated set of simplified equations are both first-order accurate with transport corrections occurring at second order. Having established the EDA’s first-order accuracy we then analyze the grey nonequilibrium-diffusion approximation and the grey Eddington approximation and show that they both preserve this first-order accuracy. Further, these approximations preserve the EDA’s first-order accuracy when made in either the comoving-frame (CMF) or the lab-frame (LF). And while analyzing the Eddington approximation, we found that the CMF and LF radiation-source equations are equivalent when neglecting O(β 2 ) terms and compared in the LF. Of course, the radiation pressures are not equivalent. It is expected that simplified physical models and numerical discretizations of the RH equations that do not preserve this first-order accuracy will not retain the correct equilibrium-diffusion solutions. As a practical example, we show that nonequilibrium-diffusion radiative-shock solutions devolve to equilibrium-diffusion solutions when the asymptotic parameter is small.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Formation and evolution of protostellar accretion discs – I. Angular-momentum budget, gravitational self-regulation, and numerical convergence

Here, we investigate the formation and early evolution of a protostellar disc from a magnetized prestellar core using non-ideal magnetohydrodynamic (MHD) simulations including ambipolar diffusion and Ohmic dissipation. The dynamical contraction of the prestellar core ultimately leads to the formation of a first hydrostatic core, after ambipolar diffusion decouples the magnetic field from the predominantly neutral gas. The hydrostatic core accumulates angular momentum from the infalling material, evolving into a rotationally supported torus; this ‘first hydrostatic torus’ then forms an accreting protostar and a rotationally supported disc. The disc spreads out by gravitational instability, reaching ~30 au in diameter at ~3 kyr after protostar formation. The total mass and angular momentum of the protostar–disc system are determined mainly by accretion of gas from an infalling pseudo-disc, which has low specific angular momentum because of magnetic braking; their removal from the protostar–disc system by outflow and disc magnetic braking are negligible, in part because the magnetic field is poorly coupled there. The redistribution of angular momentum within the protostar–disc system is facilitated mainly by gravitational instability; this allows formation of relatively large discs even when the specific angular momentum of infalling material is low. We argue that such discs should remain marginally unstable as they grow (with Toomre Q ~ 1–2), an idea that is broadly consistent with recent observational estimates for Class 0/I discs. We discuss the numerical convergence of our results, and show that properly treating the inner boundary condition is crucial for achieving convergence at an acceptable computational cost.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

Positional Dependence of SNPP VIIRS Solar Diffuser BRDF Change Factor: An Empirical Approach

The Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership satellite regularly performs on-orbit radiometric calibration of its reflective solar bands (RSBs), primarily through observations of an onboard sunlit solar diffuser (SD). The on-orbit change of the SD bidirectional reflectance distribution function (BRDF) value, quantified by a numerical factor called the H-factor, is determined by the onboard SD stability monitor. Our previous study showed that the H-factor is solar angle and view direction dependent. In this study, we determine the dependence of the H-factor on the detector SD view footprint location. We fit an empirical model to the NASA Collection 1 SNPP VIIRS Level 1B (L1B) spectral reflectance difference across the detectors in an RSB over uniform Earth scenes of the Libya 4 desert and deep convective clouds. We apply the model predicted SD positional dependent H-factor to calibrate the RSBs. Under this new calibration scheme, the original unreal striping is removed from the homogeneous Libya 4 desert and the deep convective cloud images, as well as the original unreal striping from the Dunhuang desert image. The SD positional dependent H-factor has been used to calculate the SNPP VIIRS RSB radiometric correction factor for the NASA Collection 2.0 SNPP VIIRS L1B products.

SNPP VIIRS↗
NASA NTRS2004

Characteristics of Electromagnetic Pulse Propagation in Metal

It is well known that the solution of the diffusion equation for an electromagnetic field with a time harmonic term, e(sup iwt), is in the form of a traveling wave whose amplitude attenuates over distance into a conducting medium. As the attenuation is an increasing function of frequency, the high frequency components attenuate more rapidly than those of low ones upon entering a well conducting object. At the same time, the phase velocity of an individual component is also an increasing function of frequency causing a broadening of the pulse traveling inside a conductor. In the results of our previous study of numerical simulations, the problem of using a gaussian input pulse was immediately clear. First, having the dominant frequency components distributed around zero, the movement of the peak was not well defined. Second, with the amplitude of fourier components varying slowly over a wide range, the dispersion-induced blurring of the peak position was seen to be severe. For the present study, we have used a gaussian modulated single frequency sinusoidal wave, i. e., the carrier, as an input pulse in an effort to improve the issues related to the unclear movement of peak and dispersion as described above. This was based on the following two anticipated advantages: First, the packet moves in a conductor at the group velocity calculated at the carrier frequency, which means it is well controllable. Second, the amplitude of frequency components other than that of the carrier can be almost negligible, such that the effect of dispersion can be significantly reduced. A series of experiments of transmitting electromagnetic pulses through aluminum plates of various thickness was performed to test the validity of the above points. The results of numerical simulation based on wave propagation are discussed with respect to the experimental results. Finally, a simple simulation was performed based on diffusion of a continuous sine wave input and the results are compared with those of a single frequency sinusoidal wave observed over time at difference locations inside a conductor.

Namkung, M.↗
NASA NTRS2004

Damped Windows for Aircraft Interior Noise Control

Windows are a significant path for structure-borne and air-borne noise transmission into aircraft. To improve the acoustical performance, damped windows were fabricated using two or three layers of plexiglas with transparent viscoelastic damping material sandwiched between the layers. In this paper, numerical and experimental results are used to evaluate the acoustic benefits of damped windows. Tests were performed in the Structural Acoustic Loads and Transmission Facility at NASA Langley Research Center to measure the transmission loss for diffuse acoustic excitation and radiated sound power for point force excitation. Comparisons between uniform and damped plexiglas windows showed increased transmission loss of 6 dB at the first natural frequency, 6 dB at coincidence, and 4.5 dB over a 50 to 4k Hz range. Radiated sound power was reduced up to 7 dB at the lower natural frequencies and 3.7 dB over a 1000 Hz bandwidth. Numerical models are presented for the prediction of radiated sound power for point force excitation and transmission loss for diffuse acoustic excitation. Radiated sound power and transmission loss predictions are in good agreement with experimental data. A parametric study is presented that evaluates the optimum configuration of the damped plexiglas windows for reducing the radiated sound power.

Buehrle, Ralph D.↗
DOE OSTI2021

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗
DOE OSTI2025

Nonintrusive projection-based reduced order modeling using stable learned differential operators

Nonintrusive projection-based reduced order models (ROMs) are essential for dynamics prediction in multi-query applications where underlying governing equations are known but the access to the source of the underlying full order model (FOM) is unavailable; that is, FOM is a glass-box. This article proposes a learn-then-project approach for nonintrusive model reduction. In the first step of this approach, high-dimensional stable sparse learned differential operators (S-LDOs) are determined using the generated data. In the second step, the ordinary differential equations, comprising these S-LDOs, are used with suitable dimensionality reduction and low-dimensional subspace projection methods to provide equations for the evolution of reduced states. This approach allows easy integration into the existing intrusive ROM framework to enable nonintrusive model reduction while allowing the use of Petrov–Galerkin projections. The applicability of the proposed approach is demonstrated for Galerkin and LSPG projection-based ROMs through four numerical experiments: 1-D scalar advection, 1-D Burgers, 2-D scalar advection and 1-D scalar advection–diffusion–reaction equations. In conclusion, the results indicate that the proposed nonintrusive ROM strategy provides accurate and stable dynamics prediction.

42 ENGINEERING↗
DOE OSTI2020

SubTSBR to tackle high noise and outliers for data-driven discovery of differential equations

Data-driven discovery of differential equations has been an emerging research topic. We propose a novel algorithm subsampling-based threshold sparse Bayesian regression (SubTSBR) to tackle high noise and outliers. The subsampling technique is used for improving the accuracy of the Bayesian learning algorithm. It has two parameters: subsampling size and the number of subsamples. When the subsampling size increases with fixed total sample size, the accuracy of our algorithm goes up and then down. When the number of subsamples increases, the accuracy of our algorithm keeps going up. We demonstrate how to use our algorithm step by step and compare our algorithm with threshold sparse Bayesian regression (TSBR) for the discovery of differential equations. We show that our algorithm produces better results. We also discuss the merits of discovering differential equations from data and demonstrate how to discover models with random initial and boundary condition as well as models with bifurcations. The numerical examples are: (1) predator-prey model with noise, (2) shallow water equations with outliers, (3) heat diffusion with random initial and boundary condition, and (4) fish-harvesting problem with bifurcations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Microstructure evolution and mechanical performance of nickel based superalloy C1023 at elevated temperatures

Highlights: • Growth and strengthening mechanisms of carbide in superalloy C1023 were investigated. • Carbide dendrite primary arm grows selectively based on the vertex of blocky carbonitride. • Monolayer or multilayer carbide dendrite arm has certain crystal relationship with γ′ and γ. • Plastic deformation of γ hinders interior crack propagating from MC to matrix at 650 °C. • Oxidation of GBs significantly weakens the tensile strength of the superalloy at 850 °C. - Abstract: Microstructure evolution of polycrystalline nickel based superalloy C1023 casting was investigated by comparing as cast superalloy and heat treated superalloy in the paper. Influences of strengtheners and carbides on the tensile behaviors of heat treated superalloy at different temperatures (20 °C, 650 °C, and 850 °C) were studied. Resulting data from deep corrosion suggested that growth mode of scripted carbide included nucleation on the vertex of blocky carbonitride, expansion of dendrite arm in C-rich area, and enlargement of dendrite cap. Decomposition of MC carbide generated numerous octahedral M{sub 23}C{sub 6} carbide as well as nanoscale Cr-Mo-rich nuclei. Diffusion of Mo from γ to precipitates such as γ', MC, and M{sub 23}C{sub 6}, occurred during the heat treatment while irregular γ' (4 nm–500 nm) transformed into cuboidal particles. At 650 °C, no significant deterioration of tensile strength was found while plastic deformation ability of γ base was improved. Oxidation of grain boundaries (GBs) dominated the crack behavior of the superalloy at 850 °C.

77 NANOSCIENCE AND NANOTECHNOLOGY↗