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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Microstructural Modeling of the α + β Phase in Ti-6Al-4V: A Diffusion-Based Approach

Complex heat treatment operations and advanced manufacturing processes such as laser or electron-beam welding will see the metallic workpiece experience a considerable range of temperatures and heating/cooling rates. These intrinsic conditions will have a significant bearing upon the microstructure of the material, and in turn upon the thermo-mechanical properties. In this work, a diffusion-based approach to model the growth and shrinkage of precipitates in the alpha + beta field of the common titanium alloy Ti-6Al-4V is established. Further, the numerical model is extended using a JMA-type approach to explore the dependency of the beta-transus temperature on extremely high heating rates, whereby dissolution alone is too slow to accurately describe the alpha to beta-phase transformation. Experimental heat treatments at rates of 5, 50, and 500 °C/s were performed, and metallographic analysis of the samples was used to validate the numerical modeling framework predictions for lamellar shrinkage, while data from the literature has been used to evaluate the numerical modeling framework for the growth of equiaxed microstructures. The agreement between measurements and numerical predictions was found to be good.

36 MATERIALS SCIENCE↗
NASA NTRS1999

Numerical Simulations of Acoustically Driven, Burning Droplets

This computational study focuses on understanding and quantifying the effects of external acoustical perturbations on droplet combustion. A one-dimensional, axisymmetric representation of the essential diffusion and reaction processes occurring in the vicinity of the droplet stagnation point is used here in order to isolate the effects of the imposed acoustic disturbance. The simulation is performed using a third order accurate, essentially non-oscillatory (ENO) numerical scheme with a full methanol-air reaction mechanism. Consistent with recent microgravity and normal gravity combustion experiments, focus is placed on conditions where the droplet is situated at a velocity antinode in order for the droplet to experience the greatest effects of fluid mechanical straining of flame structures. The effects of imposed sound pressure level and frequency are explored here, and conditions leading to maximum burning rates are identified.

Kim, H.-C.↗
NASA NTRS1959

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗
DOE OSTI2018

Gas Flow and Ion Transfer in Heated ESI Capillary Interfaces

Transfer capillaries are the preferred means to transport ions, generated by electrospray ionization, from ambient conditions to vacuum. During the transfer of ions through the narrow, long tubes into vacuum, substantial losses are typical. However, recently it was demonstrated that these losses can be avoided altogether. To understand the experimental observation and provide a general model for the ion transport, here, we investigate the ion transport through capillaries by numerical simulation of interacting ions. The simulation encompasses all relevant factors, such as space charge, diffusion, gas flow, and heating. Special attention is paid to the influence of the gas flow on the transmission and especially the change imposed by heating. The gas flow is modeled by a one-dimensional gas dynamics description. A large number of ions are treated as point particles in this gas flow. This allows to investigate the influence of the capillary heating on the gas flow and by this on the ion transport. The results are compared with experimental findings. .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

HomPINNs: Homotopy physics-informed neural networks for learning multiple solutions of nonlinear elliptic differential equations

Physics-informed neural networks (PINNs) based machine learning is an emerging framework for solving nonlinear differential equations. However, due to the implicit regularity of neural network structure, PINNs can only find the flattest solution in most cases by minimizing the loss functions. In this paper, we combine PINNs with the homotopy continuation method, a classical numerical method to compute isolated roots of polynomial systems, and propose a new deep learning framework, named homotopy physics-informed neural networks (HomPINNs), for solving multiple solutions of nonlinear elliptic differential equations. The implementation of an HomPINN is a homotopy process that is composed of the training of a fully connected neural network, named the starting neural network, and training processes of several PINNs with different tracking parameters. The starting neural network is to approximate a starting function constructed by the trivial solutions, while other PINNs are to minimize the loss functions defined by boundary condition and homotopy functions, varying with different tracking parameters. These training processes are regraded as different steps of a homotopy process, and a PINN is initialized by the well-trained neural network of the previous step, while the first starting neural network is initialized using the default initialization method. Finally, several numerical examples are presented to show the efficiency of our proposed HomPINNs, including reaction-diffusion equations with a heart-shaped domain.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pressure falloff testing to characterize CO 2 plume and dry-out zone during CO 2 injection in saline aquifers

The fluid mobilities and the lateral extent of the CO 2 plume and the dry-out zone are among important unknowns required to effectively manage CO 2 storage in saline aquifers projects. These unknowns can be derived from pressure falloff testing which has been widely used to characterize the subsurface. This paper provides a three-region analytical solution for falloff pressure during the infinite-acting period of CO 2 injection in saline aquifers. The solution is derived using Laplace transformation technique of the governing diffusivity equations and the accompanying initial and boundary conditions. The exact solution is validated against numerical results obtained from a black-oil reservoir simulator for two cases. The solution is approximated into simple forms in real-time domain considering different time periods. Then, a graphical semi-log analysis technique based on the approximate solutions is provided as an interpretation tool of falloff tests where the fluid mobilities and the extents of dry-out zone and CO 2 plume can be inferred. The interpretation technique is applied to synthetic falloff test data to investigate its potential and reliability.

02 PETROLEUM↗
NASA NTRS1972

Experimental determination of the turbulence in a liquid rocket combustion chamber

The intensity of turbulence and the Lagrangian correlation coefficient for a liquid rocket combustion chamber were determined experimentally using the tracer gas diffusion method. The results indicate that the turbulent diffusion process can be adequately modeled by the one-dimensional Taylor theory; however, the numerical values show significant disagreement with previously accepted values. The intensity of turbulence is higher by a factor of about two, while the Lagrangian correlation coefficient which was assumed to be unity in the past is much less than unity.

Hara, J.↗
NASA NTRS1975

Thermodynamics and kinetics of reactions in protective coating systems

A study of the aluminization of Ni from packs containing various percentages of unalloyed Al confirmed that the surface aluminum content of specimens aluminized tends to decrease with time and consequently a simple parabolic law for the weight-gain vs. time relationship is not obeyed. The diffusivity-composition relationship in NiAl was examined, and a set of curves is presented. A numerical method for the calculation of coating dissolution rates was developed and applied to NiAl-Ni3Al type of coatings.

Gupta, B.↗
NASA NTRS1979

On the acceleration of the auroral plasma of the Earth due to conic instability

The role of electrostatic conic instability in the diffusion of spilled particles and the acceleration of the auroral plasma of the Earth are numerically investigated by means of quasi-linear equations. Linear increments to the conic instability under magnetospheric conditions are introduced. The quasi-linear equation is averaged along tubes of force taking into account the drift across the magnetic field lines. It is shown that the existence of a conic instability leads to small spills and isotropic distribution of particles, but is also accompanied by significant acceleration of electrons.

Churayev, R. S.↗
NASA NTRS1979

A limit to the X-ray luminosity of nearby normal galaxies

Emission is studied at luminosities lower than those for which individual discrete sources can be studied. It is shown that normal galaxies do not appear to provide the numerous low luminosity X-ray sources which could make up the 2-60 keV diffuse background. Indeed, upper limits suggest luminosities comparable with, or a little less than, that of the galaxy. This is consistent with the fact that the average optical luminosity of the sample galaxies within approximately 20 Mpc is slightly lower than that of the galaxy. An upper limit of approximately 1% of the diffuse background from such sources is derived.

Worrall, D. M.↗
NASA NTRS1979

Effects of turbulence model selection on the prediction of complex aerodynamic flows

Numerical simulations of viscous transonic flow over a circular-arc airfoil and in a diffuser are described. The simulations are made with a new computer program designed to serve as a tool in the development of improved turbulence models for complex flows. The program incorporates zero-, one-, and two-equation eddy viscosity models and includes a variety of subsonic and supersonic boundary conditions. The airfoil flow contains a shock-separated boundary-layer interaction that has resisted previous attempts at simulation. The diffuser flow also contains a shock-boundary-layer interaction, which has not been simulated previously. Calculations using standard turbulence models, developed originally for incompressible unseparated flows, are described. Results indicate that although there are interesting differences in predictions between the various models, none of them predict the flows accurately. Suggestions for improved turbulence models are discussed.

Coakley, T. J.↗
NASA NTRS1983

An EBIC equation for solar cells

When an electron beam of a scanning electron microscope (SEM) impinges on an N-P junction, the generation of electron-hole pairs by impact ionization causes a characteristic short circuit current I(sc) to flow. The I(sc), i.e., EBIC (electron beam induced current) depends strongly on the configuration used to investigate the cell's response. In this paper the case where the plane of the junction is perpendicular to the surface is considered. An EBIC equation amenable to numerical computations is derived as a function of cell thickness, source depth, surface recombination velocity, diffusion length, and distance of the junction to the beam-cell interaction point for a cell with an ohmic contact at its back surface. It is shown that the EBIC equation presented here is more general and easier to use than those previously reported. The effects of source depth, ohmic contact, and diffusion length on the normalized EBIC characteristic are discussed.

Luke, K. L.↗
NASA NTRS1995

Dynamic response of a pulsed Burke-Schumann diffusion flame

Turbulent flames are often envisioned as an ensemble of random vortices interacting with the combustion process. A better understanding of the vortex-flame interactions therefore would be useful in improving the modeling of turbulent diffusion flames. Substantial simplification may be made by investigating controlled interactions in a laminar flame, as opposed to random interactions in a turbulent flame. The general goals of the research project are to improve our understanding of (1) the influence of buoyancy on co-flow diffusion flames and (2) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. As a first step toward objective (2), we conducted a joint experimental and numerical investigation of the vortex-flame interaction. Vortices were produced by mechanically pulsing the fuel flow at a low frequency, e.g., 10 Hz. Experiments were conducted using a nonflickering Burke-Schumann flame in both microgravity (mu-g) and normal gravity (1g) as a means of varying the buoyant force without modification of the pressure (i.e., density). The effects of buoyant convection may then be determined by a comparison of the mu-g and 1g results. The mu-g results may also reveal the important mechanisms which are masked or overwhelmed by buoyant convection in 1g. A numerical investigation was conducted using a validated, time-accurate numerical code to study the underlying physics during the flame interaction and to assist the interpretation of the experimental results.

Sheu, Jyh-Cherng↗
NASA NTRS1998

MAMA Detector

Work carried out under this grant led to fundamental discoveries and over one hundred publications in the scientific literature. Fundamental developments in instrumentation were made including all the instrumentation on the EUVE satellite, the invention of a whole new type of grazing instrument spectrometer and the development of fundamentally new photon counting detectors including the Wedge and Strip used on EUVE and many other missions and the Time Delay detector used on OREFUS and FUSE. The Wedge and Strip and Time Delay detectors were developed under this grant for less than two million dollars and have been used in numerous missions most recently for the FUSE mission. In addition, a fundamentally new type of diffuse spectrometer has been developed under this grant which has been used in instrumentation on the MMSAT spacecraft and the Lewis spacecraft. Plans are underway to use this instrumentation on several other missions as well.

Bowyer, Stuart↗
NASA NTRS1999

Pulse Phase Thermography for Defect Detection and Visualization

Pulse Phase Thermography (PPT) has been reported as a novel powerful technique of the thermal NDE. It employs application of the Discrete Fourier Transform (DFT) to thermal images obtained following flash heating of the front surface of a specimen. The computed phasegrams are excellent for defect visualization in a wide range of materials. This is in part due to their low sensitivity to uneven heating. This work is an attempt to analyze advantages and limitations of PPT. Results of application of the DFT to simulated temperature decays are presented. The temperature evolution on a surface has been simulated based on an analytical solution of the one-dimensional heat diffusion problem. A more sophisticated study has been done for different sizes of defects using numerical solution of the three-dimensional mathematical model. Capabilities of PPT for in-depth scanning and for monitoring of the material loss are discussed. The recommendations for the practical application of the PPT are presented. Experimental results obtained following these recommendations are reported.

Marinetti, Sergio↗
NASA NTRS2008

Vibroacoustic Response Data of Stiffened Panels and Cylinders

NASA has collected vibroacoustic response data on a variety of complex, aerospace structures to support research into numerical modeling of such structures. This data is being made available to the modeling community to promote the development and validation of analysis methods for these types of structures. Existing data from two structures is described, as well as plans for a data set from a third structure. The first structure is a 1.22 m by 1.22 m stiffened aluminum panel, typical of a commercial aircraft sidewall section. The second is an enclosed, stiffened aluminum cylinder, approximately 3.66 m long and 1.22 m in diameter, constructed to resemble a small aircraft fuselage with no windows and a periodic structure. The third structure is a filament-wound composite cylinder with composite stiffeners. Numerous combinations of excitation and response variables were measured on the structures, including: shaker excitation; diffuse acoustic field; velocity response from a laser vibrometer; intensity scans; and point acceleration.

Cabell, Randolph↗
NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗
DOE OSTI2023

Learning effective stochastic differential equations from microscopic simulations: Linking stochastic numerics to deep learning

We identify effective stochastic differential equations (SDEs) for coarse observables of fine-grained particle- or agent-based simulations; these SDEs then provide useful coarse surrogate models of the fine scale dynamics. We approximate the drift and diffusivity functions in these effective SDEs through neural networks, which can be thought of as effective stochastic ResNets. The loss function is inspired by, and embodies, the structure of established stochastic numerical integrators (here, Euler–Maruyama and Milstein); our approximations can thus benefit from backward error analysis of these underlying numerical schemes. They also lend themselves naturally to “physics-informed” gray-box identification when approximate coarse models, such as mean field equations, are available. Existing numerical integration schemes for Langevin-type equations and for stochastic partial differential equations can also be used for training; we demonstrate this on a stochastically forced oscillator and the stochastic wave equation. Our approach does not require long trajectories, works on scattered snapshot data, and is designed to naturally handle different time steps per snapshot. We consider both the case where the coarse collective observables are known in advance, as well as the case where they must be found in a data-driven manner.

97 MATHEMATICS AND COMPUTING↗