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 631 records · Page 35

Multilevel Techniques for Compression and Reduction of Scientific Data---The Unstructured Case

Previous work on multilevel techniques for compression and reduction of scientific data is extended to the case of data given on unstructured meshes in two and three dimensions. The centerpiece of the work is a decomposition algorithm which is shown to be optimal, in terms of both storage and operational complexity, applicable to unstructured grids in both two and three dimensions, and which implicitly gives a Riesz basis that can be exploited to reduce the data while maintaining rigorous bounds on the loss incurred. The flexibility of the approach is illustrated by applications to potential flow around an airfoil and the effect of compression on quantities of interest relevant to airfoil design; compression of computational simulation of a nonlinear reaction-diffusion system with special attention given to the problem of time series reduction; and, data from a simulation of magnetically confined plasma in a fusion reactor reduced so as to preserve the electric field computed from the data.

97 MATHEMATICS AND COMPUTING↗

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗

Mercuric iodide (HgI2) growth for nuclear detectors

The purpose of this investigation is to grow more-perfect mercuric iodide crystals in a low-gravity environment by taking advantage of diffusion-controlled growth conditions and by avoiding the problem of strain dislocations produced by the crystal's weight. This crystal has considerable practical importance as a sensitive gamma-ray detector and energy spectrometer that can operate at ambient temperature, as compared to presently available detectors that must be cooled to near liquid nitrogen temperatures. However, the performance of mercuric iodide crystals only rarely approaches the expected performance, presumably because some of the free electrical charges produced within the crystal are not collected at the electrodes, but instead remain trapped or immobilized at crystal defects. An efficient high atomic number semiconductor detector capable of operating at room temperature utilizing single HgI2 crystals offers a greater potential than existing detector technology.

Schnepple, W.↗

Application of the method of lines for solutions of the Navier-Stokes equations using a nonuniform grid distribution

The feasibility of the method of lines for solutions of physical problems requiring nonuniform grid distributions is investigated. To attain this, it is also necessary to investigate the stiffness characteristics of the pertinent equations. For specific applications, the governing equations considered are those for viscous, incompressible, two dimensional and axisymmetric flows. These equations are transformed from the physical domain having a variable mesh to a computational domain with a uniform mesh. The two governing partial differential equations are the vorticity and stream function equations. The method of lines is used to solve the vorticity equation and the successive over relaxation technique is used to solve the stream function equation. The method is applied to three laminar flow problems: the flow in ducts, curved-wall diffusers, and a driven cavity. Results obtained for different flow conditions are in good agreement with available analytical and numerical solutions. The viability and validity of the method of lines are demonstrated by its application to Navier-Stokes equations in the physical domain having a variable mesh.

Abolhassani, J. S.↗

A convective model for turbulent mixing in rotating convection zones

The effects of rotation are included in an analytical model for the convective motions in a plane-parallel layer of an ideal fluid. The turbulent stress tensor, formed by taking products and averages of the various velocity components, is calculated for an arbitrary eddy size and shape. Heuristic formulae presented for determining the size and shape of the dominant eddy then give a fully specified stress tensor. Applications for this stress tensor in problems of stellar internal dynamics, heat flow, scalar diffusion, and dynamo theory are suggested. The resultant stresses tend to produce differential rotation profiles with rapidly rotating equators and interiors. The dynamo activity associated with these convective motions tends to occur near the lower boundary of the convection zone.

Hathaway, D. H.↗

Preliminary study on pressure brazing and diffusion welding of Nb-1Zr to Inconel 718

Future space power systems may include Nb-1Zr/Inconel 718 dissimilar metal joints for operation at 1000 K for 60,000 h. The serviceability of pressure-brazed and diffusion-welded joints was investigated. Ni-based metallic glass foil filler metals were used for brazing. Ni and Fe foils were used as diffusion welding inter-layers. Joint soundness was determined by metallographic examination in the as-brazed and as-welded condition, after aging at 1000 K, and after thermal cycling. Brazed joints thermally cycled in the as-brazed condition and diffusion-welded joints were unsatisfactory because of cracking problems. Brazed joints may meet the service requirements if the joints are aged at 1000 K prior to thermal cycling.

Moore, T. J.↗

New results in the theory of dust grain alignment

Two complementary approaches are used in an attempt to propose an appropriate formulation of the solution to the problem of magnetic alignment of grains in the diffuse and/or the more denser clouds, whatever the mechanism of rotational excitation can be. The interest of such a unified formulation is mainly that the same theoretical expression for polarization can be used everywhere, allowing for easier comparisons between regions where the physical conditions are highly different. The first consists in applying a Monte-Carlo method to a limited number of representative cases, for which all the torques acting on the grain are taken into account: impulsive random torques due to direct collisions with gas atoms, to evaporation of atoms from the surface, and to exo-energetic recombinations forming hydrogen molecules, followed by violent ejections from peculiar sites; magnetic torques. Three characteristic times are associated with these torques: the collisional damping time, the time necessary to change completely the actual sites configuration narrowly bound to the correlation time of the suprathermal torque; and the magnetic damping time. The second approach starts from a heuristic point of view. It consists in a generalization of results (Cugnon, 1983; see also Purcell and Spitzer, 1971; Greenberg, 1978) obtained for thermal alignment to the suprathermal case. It appears indeed that in two extreme cases, the thermal formulation may be used after redefinition of involved times and temperatures.

Cugnon, Pierre↗

Comets: Role and importance to exobiology

The transfer of organic compounds from interstellar space to the outskirts of a protoplanetary disk, their accretion into cometary objects, and the transport of the latter into the inner solar system by orbital diffusion throw a new light on the central problem of exobiology. It suggests the existence of a cosmic mechanism, working everywhere, that can supply prebiotic compounds to ubiquitous rocky planets, in search of the proper environment to start life in many places in the Universe. Under the heading of chemistry of the cometary nucleus, the following topics are covered: radial homogeneity of the nucleus; the dust-to-ice ratio; nature of the dust grains; origin of the dust in comets; nature of the volatile fraction; the CO distribution in comet Halley; dust contribution to the volatile fraction; elemental balance sheet of comet Halley; quantitative molecular analysis of the volatile fraction; and isotopic ratios. Under the heading of exogenous origin of carbon on terrestrial planets the following topics are covered: evidence for a high-temperature phase; from planetesimals to planets; a veneer of volatile and organic material; and cometary contribution.

Delsemme, Armand H.↗

The effects of sooting in droplet combustion

The study of the burning of a single droplet is an ideal problem from which to gain fundamental understanding of diffusion flame characteristics. Droplet combustion is a complex physico-chemical process that involves a chemically-reacting two-phase flow with phase changes and yet simple experiments and analysis can be used to attain important insights into the burning rate, flame dynamics, kinetic extinction and disruption processes. It is a subject that has been actively studied for the past 40 years with most of the fundamental experiments being performed under reduced-gravity conditions for direct comparisons with theoretical/computational analyses that invoke spherical symmetry assumptions. In the earlier studies, the effects of sooting on the overall burning characteristics were not considered. However, recent microgravity investigations performed at the NASA-LeRC droptowers (Droplet Combustion Experiment) and others indicate that effects of soot and sootcloud formation may be significant during the lifetime of the droplet and therefore must be included in the analysis.

Lee, Kyeong-Ook↗

Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers

A four-equation model of turbulence is applied to the numerical simulation of flows with massive separation induced by a sudden expansion. The model constants are a function of the flow parameters, and two different formulations for these functions are tested. The results are compared with experimental data for a high Reynolds-number case and with experimental and DNS data for a low Reynolds-number case. The computations prove that the recovery region downstream of the massive separation is properly modeled only for the high Re case. The problems in this case stem from the gradient diffusion hypothesis, which underestimates the turbulent diffusion.

Michelassi, V.↗

Oxidative Attack of Carbon/Carbon Substrates through Coating Pinholes

A critical issue with oxidation protected carbon/carbon composites used for spacecraft thermal protection is the formation of coating pinholes. In laboratory experiments, artificial pinholes were drilled through SiC-coatings on a carbon/carbon material and the material was oxidized at 600, 1000, and 1400 C at reduced pressures of air. The attack of the carbon/carbon was quantified by both weight loss and a novel cross-sectioning technique. A two-zone, one dimensional diffusion control model was adapted to analyze this problem. Agreement of the model with experiment was reasonable at 1000 and 1400 C; however results at lower temperatures show clear deviations from the theory suggesting that surface reaction control plays a role.

Jacobson, Nathan S.↗

Antimatter/HiPAT Support Services

Techniques were developed for trapping normal matter in the High Performance Antiproton Trap (HiPAT). Situations encountered included discharge phenomena, charge exchange and radial diffusion processes. It is important to identify these problems, since they will also limit the performance in trapping antimatter next year.

Lewis, Raymond A.↗

Extension of Low Dissipative High Order Hydrodynamics Schemes for MHD Equations

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD (magnetohydrodynamic) equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls to limit the amount and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvi-linear grids. The three features of the present MHD scheme over existing schemes in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion magnetized flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike existing schemes for the conservative MHD equations which suffer from ill-conditioned eigen-decompositions, the present scheme makes use of a well-conditioned eigen-decomposition to solve the conservative form of the MHD equations. This is due to, partly. the fact that the divergence of the magnetic field condition is a different type of constraint from its incompressible Navier-Stokes cousin. Third, a new approach to minimize the numerical error of the divergence free magnetic condition for high order scheme is introduced.

Yee, H. C.↗

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows, I: Basic Theory

The objective of this paper is to extend our recently developed highly parallelizable nonlinear stable high order schemes for complex multiscale hydrodynamic applications to the viscous MHD equations. These schemes employed multiresolution wavelets as adaptive numerical dissipation controls t o limit the amount of and to aid the selection and/or blending of the appropriate types of dissipation to be used. The new scheme is formulated for both the conservative and non-conservative form of the MHD equations in curvilinear grids. The four advantages of the present approach over existing MHD schemes reported in the open literature are as follows. First, the scheme is constructed for long-time integrations of shock/turbulence/combustion MHD flows. Available schemes are too diffusive for long-time integrations and/or turbulence/combustion problems. Second, unlike exist- ing schemes for the conservative MHD equations which suffer from ill-conditioned eigen- decompositions, the present scheme makes use of a well-conditioned eigen-decomposition obtained from a minor modification of the eigenvectors of the non-conservative MHD equations t o solve the conservative form of the MHD equations. Third, this approach of using the non-conservative eigensystem when solving the conservative equations also works well in the context of standard shock-capturing schemes for the MHD equations. Fourth, a new approach to minimize the numerical error of the divergence-free magnetic condition for high order schemes is introduced. Numerical experiments with typical MHD model problems revealed the applicability of the newly developed schemes for the MHD equations.

Sjoegreen, Bjoern↗

Updated Computational Model of Cosmic Rays Near Earth

An updated computational model of the galactic-cosmic-ray (GCR) environment in the vicinity of the Earth, Earth s Moon, and Mars has been developed, and updated software has been developed to implement the updated model. This model accounts for solar modulation of the cosmic-ray contribution for each element from hydrogen through iron by computationally propagating the local interplanetary spectrum of each element through the heliosphere. The propagation is effected by solving the Fokker-Planck diffusion, convection, energy-loss boundary-value problem. The Advanced Composition Explorer NASA satellite has provided new data on GCR energy spectra. These new data were used to update the original model and greatly improve the accuracy of prediction of interplanetary GCR.

ONeill, Patrick M.↗

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗

Neural Network Enhanced RKPM for Electrochemical-Mechanical Coupled Damage Modeling of Energy Storage Materials

Energy storage materials undergo significant charge cycling, which makes understanding their reliability and durability fundamental in predicting performance and service life. Strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking, largely along material interfaces and grain boundaries. For Li-ion batteries, for example, lithium moving between electrodes during charging and discharging process causes expansion and contraction of grains, and the strongly anisotropic and nonlinearly [Li]-dependent grain material properties can cause grains to expand into and contract away from each other, leading to chemo-mechanical cracking. In the first part of this work, a RKPM based computational framework for solving the coupled solid-phase lithium conservation with Fickian diffusion and the lithium concentration dependent anisotropic mechanical problem subjected to a highly nonlinear Butler-Volmer boundary condition is introduced. The choice of RKPM completeness conditions for lithium concentration and mechanical deformation fields, and the variational consistency condition for the domain integration of the coupled problem is first determined. In the second part of this work, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1] is leveraged to accurately capture damage and crack propagation throughout the material, by learning the location, orientation, and sharpness of discontinuity while allowing for a coarser nodal distribution than that is necessary for capturing sharp solution transitions using traditional mesh-based methods. NN-RKPM is used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure.

damage modeling↗

Distributed Multi-GPU Community Detection on Exascale Computing Platforms

Community detection is a fundamental operation in graph mining, and by uncovering hidden structures and patterns within complex systems it helps solve fundamental problems pertaining to social networks, such as information diffusion, epidemics, and recommender systems. Scaling graph algorithms for massive networks becomes challenging on modern distributed-memory multi-GPU (Graphics Processing Unit) systems due to limitations such as irregular memory access patterns, load imbalances, higher communication-computation ratios, and cross-platform support. We present a novel algorithm HiPDPL-GPU (Distributed Parallel Louvain) to address these challenges. We conduct experiments involving different partitioning techniques to achieve an optimized performance of HiPDPL-GPU on the two largest supercomputers: Frontier and Summit. Remarkably, HiPDPL-GPU processes a graph with 4.2 billion edges in less than 3 minutes using 1024 GPUs. Qualitatively, the performance of HiPDPL-GPU is similar or better compared to other state-of-the-art CPU- and GPU-based implementations. While prior GPU implementations have predominantly employed CUDA, our first-of-its-kind implementation for community detection is cross-platform, accommodating both AMD and NVIDIA GPUs.

Sattar, Naw Safrin↗