Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,387 records · Page 77

The Implications of Experimentally Controlled Gravitational Accelerations for Combustion Science

An overview of basic combustion problems which have been investigated under the condition of reduced gravity is presented to identify promising research directions. Attention is given to the broad categories of gas-jet diffusion flames, droplet combustions, particle clouds, flame spreading over liquid pools, smoldering, and flame spreading over solid fuels. Fire safety in spacecraft is the primary application that is addressed by the studies of combustion under microgravity. The need for more complete testing of the issues discussed in orbiting spacecraft is identified in the light of limited earth-based testing. Attention is also directed toward the need for advanced diagnostic methods for in-flight and other combustion investigations.

Sacksteder, Kurt R.↗

Influence of test configuration on the combustion characteristics of polymers as ignition sources

The experimental evaluation of polymers as ignition sources for metals was accomplished at the NASA White Sands Test Facility (WSTF) using a standard promoted combustion test. These tests involve the transient burning of materials in high-pressure oxygen environments. They have provided data from which design decisions can be made; data include video recordings of ignition and non-ignition for specific combinations of metals and polymers. Other tests provide the measured compositions of combustion products for polymers at select burn times and an empirical basis for estimating burn rates. With the current test configuration, the detailed analysis of test results requires modeling a three-dimensional, transient convection process involving fluid motion, thermal conduction and convection, the diffusion of chemical species, and the erosion of sample surface. At the high pressure extremes, it even requires the analysis of turbulent, transient convection where the physics of the problem are not well known and the computation requirements are not practical at this time. An alternative test configuration that can be analyzed with a relatively-simple convection model was developed during the summer period. The principal change constitutes replacing a large-diameter polymer disk at the end of the metal test rod with coaxial polymer cylinders that have a diameter nearer to that of the metal rod. The experimental objective is to assess the importance of test geometries on the promotion of metal ignition by testing with different lengths of the polymer and, with an extended effort, to analyze the surface combustion in the redesigned promoted combustion tests through analytical modeling of the process. The analysis shall use the results of cone-calorimeter tests of the polymer material to model primary chemical reactions and, with proper design of the promoted combustion test, modeling of the convection process could be conveniently limited to a quasi-steady boundary layer analysis where the economical solution of parabolic equations is involved. The products for the summer period are: (1) a conceptual-level redesign of the test apparatus, and (2) the development and use of an approximate integral boundary layer analysis to demonstrate the influence of geometry changes prior to testing. A computer code STAN5, an accurate numerical boundary layer model whose earlier versions were developed for the NASA Lewis Research Center by the Fellow, also was installed and validated on the WSTF and New Mexico State University computer systems as a starting point in the development of a more detailed fluid mechanics and combustion model.

Julien, Howard L.↗

Advanced Electromagnetic Methods for Aerospace Vehicles

The imminent destructive threats of Lightning on helicopters and other airborne systems has always been a topic of great interest to this research grant. Previously, the lightning induced currents on the surface of the fuselage and its interior were predicted using the finite-difference time-domain (FDTD) method as well as the NEC code. The limitations of both methods, as applied to lightning, were identified and extensively discussed in the last meeting. After a thorough investigation of the capabilities of the FDTD, it was decided to incorporate into the numerical method a subcell model to accurately represent current diffusion through conducting materials of high conductivity and finite thickness. Because of the complexity of the model, its validity will be first tested for a one-dimensional FDTD problem. Although results are not available yet, the theory and formulation of the subcell model are presented and discussed here to a certain degree. Besides lightning induced currents in the interior of an aircraft, penetration of electromagnetic fields through apertures (e.g., windows and cracks) could also be devastating for the navigation equipment, electronics, and communications systems in general. The main focus of this study is understanding and quantifying field penetration through apertures. The simulation is done using the FDTD method and the predictions are compared with measurements and moment method solutions obtained from the NASA Langley Research Center. Cavity-backed slot (CBS) antennas or slot antennas in general have many applications in aircraft-satellite type of communications. These can be flushmounted on the surface of the fuselage and, therefore, they retain the aerodynamic shape of the aircraft. In the past, input impedance and radiation patterns of CBS antennas were computed using a hybrid FEM/MoM code. The analysis is now extended to coupling between two identical slot antennas mounted on the same structure. The predictions are performed using both the hybrid FEM/MoM and the FDTD NEWS code. The results are compared with each other as well as with measurements performed in the ElectroMagnetic Anechoic Chamber (EMAC) of ASU. In addition to self and mutual impedances versus frequency, the comparisons include mutual coupling S(sub 12) as a function of distance for various slot orientations.

Balanis, Constantine A.↗

Towards robust autonomous impedance spectroscopy analysis: A calibrated hierarchical Bayesian approach for electrochemical impedance spectroscopy (EIS) inversion

Distribution-based analyses, such as the distribution of relaxation times (DRT) and the distribution of diffusion times (DDT), present model-free alternatives to equivalent circuit modeling for analysis of electrochemical impedance spectroscopy (EIS) data. However, reconstructing such distributions from noisy impedance data is an ill-posed problem that must be solved with specialized inversion algorithms, requiring careful control and tuning. Furthermore, most inversion algorithms developed to date can only solve problems of limited complexity. Herein, we present a new hierarchical Bayesian method for EIS inversion, leveraging efficient algorithms for optimization and Hamiltonian Monte Carlo (HMC) sampling to solve models of arbitrary complexity. We overcome the challenge of ad-hoc parameter tuning by encoding intrinsic characteristics of the DRT and DDT into flexible prior distributions and “pre-calibrating” the model to simulated data. This approach is versatile, highly robust to noise, and provides quantitative estimates of both the error structure of the data and the uncertainty in the recovered distributions. The model is validated with simulated data to demonstrate accurate recovery of the DRT and the DDT. The method also shows promise for simultaneous recovery of multiple distributions, raising the intriguing possibility of semi-autonomous EIS analysis and ad-hoc model construction. Finally, the practical utility of the method is illustrated with experimental data. Throughout, we draw comparisons to several recently published EIS inversion methodologies.

36 MATERIALS SCIENCE↗

Physics-Informed Generative Adversarial Networks for Stochastic Differential Equations

We developed a new class of physics-informed generative adversarial networks (PI-GANs) to solve forward, inverse, and mixed stochastic problems in a unified manner based on a limited number of scattered measurements. Unlike standard GANs relying solely on data for training, here we encoded into the architecture of GANs the governing physical laws in the form of stochastic differential equations (SDEs) using automatic differentiation. In particular, we applied Wasserstein GANs with gradient penalty (WGAN-GP) for its enhanced stability compared to vanilla GANs. We first tested WGAN-GP in approximating Gaussian processes of different correlation lengths based on data realizations collected from simultaneous reads at sparsely placed sensors. We obtained good approximation of the generated stochastic processes to the target ones even if there is a mismatch between the input noise dimensionality and the effective dimensionality of the target stochastic processes. We also studied the overfitting issue for both the discriminator and the generator, and we found that overfitting occurs also in the generator in addition to the discriminator as previously reported. Subsequently, we considered the solution of elliptic SDEs requiring approximations of three stochastic processes, namely the solution, the forcing, and the diffusion coefficient. Here again, we assumed data collected from simultaneous reads at a limited number of sensors for the multiple stochastic processes. Furthermore, three generators were used for the PI-GANs: two of them were feed forward deep neural networks (DNNs), while the other one was the neural network induced by the SDE. For the case where we have one group of data, we employed one feed forward DNN as the discriminator, while for the case of multiple groups of data we employed multiple discriminators in PI-GANs. We solved forward, inverse, and mixed problems without changing the framework of PI-GANs, obtaining both the means and the standard deviations of the stochastic solution and the diffusion coefficient in good agreement with benchmarks. In this work, we have demonstrated the effectiveness of PI-GANs in solving SDEs for about 120 dimensions. In principle, PI-GANs could tackle very high dimensional problems given more sensor data with low-polynomial growth in computational cost.

97 MATHEMATICS AND COMPUTING↗

Analysis of Integrated Fuselage Composite Seat for Small Aircraft in Crashworthiness Applications

The research paper is addressed to develop and validate an aircraft seat cushion to be used for the passenger seat directly integrated in the fuselage. For very light aircraft, these systems result to be a good choice in terms of weight reduction and space optimization, but can represent a challenging problem for crashworthiness design and a key energy-absorbing function for human protection from injuries so to improve the level of the safety. The first issue was dedicated to investigate the mechanical properties of four commercial foams, diffused in different industrial applications, which present nonlinear and time-dependent viscoelastic behavior, and finally to determine whether typical seat cushion materials are rate sensitive or not. Then, numerical simulations were validated with measurements to simulate the mechanical properties of foams and to evaluate the contact pressures, stress distribution, and the dependence of stress on the strain rate in order to provide adequate levels of both comfort and energy absorption.

36 MATERIALS SCIENCE↗

Asymptotic solution of the diffusion equation in slender impermeable tubes of revolution. I. The leading-term approximation

The anisotropic 3D equation describing the pointlike particles diffusion in slender impermeable tubes of revolution with cross section smoothly depending on the longitudinal coordinate is the object of our study. We use singular perturbations approach to find the rigorous asymptotic expression for the local particles concentration as an expansion in the ratio of the characteristic transversal and longitudinal diffusion relaxation times. The corresponding leading-term approximation is a generalization of well-known Fick-Jacobs approximation. This result allowed us to delineate the conditions on temporal and spatial scales under which the Fick-Jacobs approximation is valid. A striking analogy between solution of our problem and the method of inner-outer expansions for low Knudsen numbers gas kinetic theory is established. With the aid of this analogy we clarify the physical and mathematical meaning of the obtained results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Understanding Detrimental and Beneficial Grain Boundary Effects in Halide Perovskites

Abstract Grain boundaries play a key role in the performance of thin‐film optoelectronic devices and yet their effect in halide perovskite materials is still not understood. The biggest factor limiting progress is the inability to identify grain boundaries. Noncrystallographic techniques can misidentify grain boundaries, leading to conflicting literature reports about their influence; however, the gold standard – electron backscatter diffraction (EBSD) – destroys halide perovskite thin films. Here, this problem is solved by using a solid‐state EBSD detector with 6000 times higher sensitivity than the traditional phosphor screen and camera. Correlating true grain size with photoluminescence lifetime, carrier diffusion length, and mobility shows that grain boundaries are not benign but have a recombination velocity of 1670 cm s −1 , comparable to that of crystalline silicon. Amorphous grain boundaries are also observed that give rise to locally brighter photoluminescence intensity and longer lifetimes. This anomalous grain boundary character offers a possible explanation for the mysteriously long lifetime and record efficiency achieved in small grain halide perovskite thin films. It also suggests a new approach for passivating grain boundaries, independent of surface passivation, to lead to even better performance in optoelectronic devices.

Adhyaksa, Gede W. P.↗

Nonlinear Polarization Effects in Dielectrics with Hydrogen Bonds

Using quasiclassical kinetic theory, nonlinear phenomena under the formation of space-charge polarization in crystals with hydrogen bonds (HBC) are investigated. From the solution of the nonlinear system of Fokker-Planck and Poisson equations with blocking electrodes, it is established that for the mathematical description of the relaxation polarization in HBCs in weak fields (100-1000 kV/m) and at high temperatures (T> 350 K), it is sufficient to use a linear approximation of perturbation theory. Theoretically, it was found that at the first odd frequency of the alternating field, nonlinear effects caused by the interaction of relaxation modes of various frequency orders start to appear. The diffusion and mobility coefficients are calculated taking into account both mechanisms of the protons transitions (thermally activated and tunneling transitions) through a parabolic potential barrier. A recurrence expression is constructed for calculating complex amplitudes of relaxation modes generated at an arbitrary frequency (multiple to the fundamental frequency) of the alternating field in an arbitrary approximation of perturbation theory. The proposed scheme for solving the kinetic equation can be applied to other crystals with ionic conductivity similar to HBCs in the type and properties of the crystal lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The case for electron re-acceleration at galaxy cluster shocks

On the largest scales, the Universe consists of voids and filaments making up the cosmic web. Galaxy clusters are located at the knots in this web, at the intersection of filaments. Clusters grow through accretion from these large-scale filaments and by mergers with other clusters and groups. In a growing number of galaxy clusters, elongated Mpc-sized radio sources have been found. Also known as radio relics, these regions of diffuse radio emission are thought to trace relativistic electrons in the intracluster plasma accelerated by low-Mach-number shocks generated by cluster–cluster merger events. A long-standing problem is how low-Mach-number shocks can accelerate electrons so efficiently to explain the observed radio relics. Here, we report the discovery of a direct connection between a radio relic and a radio galaxy in the merging galaxy cluster Abell 3411–3412 by combining radio, X-ray and optical observations. This discovery indicates that fossil relativistic electrons from active galactic nuclei are re-accelerated at cluster shocks. Lastly, it also implies that radio galaxies play an important role in governing the non-thermal component of the intracluster medium in merging clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

A collective elastic fluctuation mechanism for decoupling and stretched relaxation in glassy colloidal and molecular liquids

We propose a microscopic theory for the decoupling of self-diffusion and structural relaxation in glass-forming liquids within the Elastically Collective Nonlinear Langevin Equation (ECNLE) activated dynamics framework. Our central hypothesis is that the heterogeneity relevant to this problem is static fluctuations of local density on the scale of 3–4 particle diameters and how this changes local packing correlations. These fluctuations modify the degree of dynamical cage expansion that mechanistically couples intracage large amplitude hopping and longer range collective elasticity in ECNLE theory. Decoupling only emerges in the deeply supercooled regime where the strongly temperature dependent elastic barrier becomes non-negligible relative to its noncooperative local analog. The theory makes predictions for various aspects of the decoupling phenomenon, including apparent fractional power law Stokes-Einstein behavior, that appear to be consistent with experiments and simulations on hard sphere fluids and molecular liquids. Of central importance is a microscopic connection between the barrier fluctuation variance and most probable barrier height. Sensible results are also obtained for the nonexponential stretching of a generic relaxation time correlation function and its temperature evolution. Nonuniversality can arise from the relative importance of the local and collective barriers (related to fragility) and the precise magnitude of the length scale that defines the transition from local cage to elastic physics. Comparison is made with a traplike model based on a Gaussian distribution of barriers.

Xie, Shi-Jie (ORCID:0000000333174720)↗

A Numerical Tool for the Design of Enhanced Geothermal Reservoirs based on Modeling of Thermo-Hydro-Mechanical Processes and Explicit Representation of 3D Discrete Fracture Networks

This project has developed a coupled thermo-hydro-mechanical numerical modeling tool for prediction of energy production from Enhanced Geothermal Reservoirs (EGS). The tool is unique because it can explicitly represent thousands of pre-existing natural fractures. The capabilities developed in this project were built on the already existing expertise of Itasca arising from thirty years of research and development in discrete element modeling of fractured rocks. A thermal module suitable for analysis of different heat transfer mechanisms in fractured rocks percolated by fluids was developed and coupled with the already existing hydro-mechanical calculation in 3DEC, a commercial software for analysis of jointed rock based on the discrete element method (DEM). The heat transfer module simulates heat conduction, heat advection by moving fluid in rock fractures, and convective heat exchange at the fracture walls between the rock and the moving fluid. Techniques for spatial discretization of heat transfer in EGS-type problems were evaluated, and a boundary layer meshing approach was developed. Also, given the significant differences in time scale of physical processes involved, temporal integration schemes for coupled advection-diffusion processes and coupled thermal-hydro-mechanical processes were developed by taking advantage of both implicit and explicit time integration and also optimizing the coupling intervals. These improvements were necessary for modeling EGS reservoirs. Also, the DFN capacities in 3DEC were extended to include a larger number of fractures. Research was carried out on approaches to simplify the DFN so that the model can meet computational demands, while preserving thermo-hydro-mechanical characteristics of the DFN. The developed modules were coupled with the already existing capabilities of 3DEC. The results were presented as a submitted conference paper (Riahi et al., 2018a). A series of benchmark and test examples showed the developed capabilities. Finally, the developed tool was tested for a large-scale reservoir model for longer periods of injection and production. The thermal draw-down curves were generated and presented. The project successfully met all the milestones laid out in the proposal. Results demonstrated the feasibility of developing a fast and reliable software for the EGS community to use for assessment and prediction of reservoir performance. Some of its applications include assessment of stimulated volume and produced temperature and power. It can also be used for improved engineering, design and optimization of EGS reservoirs as well as other unconventional resources such as shale gas. Finally, the 3D analysis presented in this study showed that simulating EGS reservoirs with tens of thousands of fractures is feasible, and confirmed that the developed tool can serve toward meeting the goals and objectives of Geothermal Technologies Office in overcoming technical barriers associated with commercial EGS development.

15 GEOTHERMAL ENERGY↗

Thermal radiation heat transfer.

A comprehensive discussion of heat transfer by thermal radiation is presented, including the radiative behavior of materials, radiation between surfaces, and gas radiation. Among the topics considered are property prediction by electromagnetic theory, the observed properties of solid materials, radiation in the presence of other modes of energy transfer, the equations of transfer for an absorbing-emitting gas, and radiative transfer in scattering and absorbing media. Also considered are radiation exchange between black isothermal surfaces, radiation exchange in enclosures composed of diffuse gray surfaces and in enclosures having some specularly reflecting surfaces, and radiation exchange between nondiffuse nongray surfaces. The use of the Monte Carlo technique in solving radiant-exchange problems and problems of radiative transfer through absorbing-emitting media is explained.

Siegel, R.↗

The chemistry and diffusion of aircraft exhausts in the lower stratosphere during the first few hours after fly-by

An analysis of the hydrogen-nitrogen-oxygen reaction systems in the lower stratosphere as they are initially perturbed by individual aircraft engine exhaust plumes was conducted in order to determine whether any significant chemical reactions occur, either among exhaust chemical species, or between these species and the environmental ozone, while the exhaust products are confined to intact plume segments at relatively high concentrations. The joint effects of diffusive mixing and chemical kinetics on the reactions were also studied, using the techniques of second-order closure diffusion/chemistry models. The focus of the study was on the larger problem of the potential depletion of ozone by supersonic transport aircraft exhaust materials emitted into the lower stratosphere.

Hilst, G. R.↗

Additions and corrections to the absorption coefficients of CO2 ice - Applications to the Martian south polar cap

Reflectance spectra of carbon dioxide frosts were calculated using the optical constants provided by Warren (1986) for the wavelength region 2-6 microns. In comparing these calculated spectra to spectra of frosts observed in the laboratory and on the surface of Mars, problems in the optical constants presented by Warren (1986) became apparent. Absorption coefficients for CO2 ice have been derived using laboratory reflectance measurements and the Hapke (1981) model for calculating diffuse reflectance. This provides approximate values in regions where no data were previously available and indicates where corrections to the compilation by Warren (1986) are required. Using these coefficients to calculate the reflectance of CO2 ice at varying grain sizes indicates that a typical Mariner polar cap spectrum is dominated by absorptions due to CO2 frost or ice at grain sizes that are quite large, probably of the order of millimeters to centimeters. There are indications of contamination of water frost or dust, but confirmation will require more precise absorption coefficients for solid CO2 than can be obtained from the method used here.

Calvin, Wendy M.↗

Electron heating in quasi-perpendicular shocks - A Monte Carlo simulation

To study the problem of electron heating in quasi-perpendicular shocks, under the combined effects of 'reversible' motion, in the shock electric potential and magnetic field, and wave-particle interactions a diffusion equation was derived, in the drift (adiabatic) approximation and it was solved by using a Monte Carlo method. The results show that most of the observations can be explained within this framework. The simulation has also definitively shown that the electron parallel temperature is determined by the dc electromagnetic field and not by any wave particle induced heating. Wave-particle interactions are effective in smoothing out the large gradients in phase space produced by the 'reversible' motion of the electrons, thus producing a 'cooling' of the electrons. Some constraints on the wave-particle interaction process may be obtained from a detailed comparison between the simulation and observations. In particular, it appears that the adiabatic approximation must be violated in order to explain the observed evolution of the perpendicular temperature.

Veltri, Pierluigi↗

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results.

Glicksman, M. E.↗

Renormalization Group Tutorial

Complex physical systems sometimes have statistical behavior characterized by power- law dependence on the parameters of the system and spatial variability with no particular characteristic scale as the parameters approach critical values. The renormalization group (RG) approach was developed in the fields of statistical mechanics and quantum field theory to derive quantitative predictions of such behavior in cases where conventional methods of analysis fail. Techniques based on these ideas have since been extended to treat problems in many different fields, and in particular, the behavior of turbulent fluids. This lecture will describe a relatively simple but nontrivial example of the RG approach applied to the diffusion of photons out of a stellar medium when the photons have wavelengths near that of an emission line of atoms in the medium.

Bell, Thomas L.↗