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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Solutions and reductions for radiative energy transport in laser-heated plasma

A full symmetry classification is given for models of energy transport in radiant plasma when the mass density is spatially variable and the diffusivity is nonlinear. A systematic search for conservation laws also leads to some potential symmetries and to an integrable nonlinear model. Classical point symmetries, potential symmetries, and nonclassical symmetries are used to effect variable reductions and exact solutions. The simplest time-dependent solution is shown to be stable and relevant to a closed system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Light scattering by a rough surface of human skin. 2. Diffuse reflectance

Based on the previously calculated luminance factors, we have investigated the integral characteristics of light reflection from a rough surface of the skin with large-scale inhomogeneities under various conditions of the skin illumination. Shadowing of incident and scattered beams by relief elements is taken into account. Diffuse reflectances by the Gaussian and the quasi-periodic surfaces are compared and, in general, both these roughness models are shown to give similar results. We have studied the effect of the angular structure of radiation multiply scattered deep in the tissue and the refraction of rays as they propagate from the dermis to the surface of the stratum corneum on the reflection characteristics of the skin surface. The importance of these factors is demonstrated. The algorithms constructed can be included in the schemes of calculation of the light fields inside and outside the medium in solving various direct and inverse problems of optics of biological tissues. (biophotonics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Electrochemical Amyloid β Immunosensor Based on Ti 3 C 2 T x MXene Nanosheets

Alzheimer’s disease (AD) is the most common neurodegenerative disease and a major health-care problem. More than 55 million people in the world, including 6.7 million Americans, are living with AD. As AD is characterized by production and deposition of β-amyloid peptide (Aβ) in the brain, Aβ is considered as a potential biomarker for diagnosing and monitoring the progression of AD. Nanomaterials-based electrochemical immunosensor is one of the most promising methods for early diagnosis and therapeutic of AD, due to its highly selective, sensitive, and label-free advantages. The graphene-like two-dimensional (2D) MXene nanosheets have shown large surface area, good conductivity, biocompatibility, high ion transport, and low diffusion barrier, attracting attention in recent years to use it as the nanomaterial substrate for biosensor fabrication. In this paper, we presented a Ti 3 C 2 T x MXene nanosheets-based electrochemical immunosensor to detect trace Aβ for early diagnosis of AD.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A Hybrid Finite-Volume, Discontinuous Galerkin Discretization for the Radiative Transport Equation

In this report we propose a hybrid spatial discretization for the radiative transport equation that combines a second-order discontinuous Galerkin (DG) method and a second-order finite-volume (FV) method. The strategy relies on a simple operator splitting that has been used previously to combine different angular discretizations. Unlike standard FV methods with upwind fluxes, the hybrid approach is able to accurately simulate problems in scattering dominated regimes. However, it requires less memory and yields a faster computational time than a uniform DG discretization. In addition, the underlying splitting allows naturally for hybridization in both space and angle. Numerical results are given to demonstrate the efficiency of the hybrid approach in the context of discrete ordinate angular discretizations and Cartesian spatial grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

The Extinction Properties of and Distance to the Highly Reddened Type IA Supernova 2012cu

Correcting Type Ia Supernova brightnesses for extinction by dust has proven to be a vexing problem. Here we study the dust foreground to the highly reddened SN 2012cu, which is projected onto a dust lane in the galaxy NGC 4772. The analysis is based on multi-epoch, spectrophotometric observations spanning from 3300–9200 Å, obtained by the Nearby Supernova Factory. Phase-matched comparison of the spectroscopically twinned SN 2012cu and SN 2011fe across 10 epochs results in the best-fit color excess of (E(B−V), RMS) = (1.00, 0.03) and total-to-selective extinction ratio of (R {sub V}, RMS) = (2.95, 0.08) toward SN 2012cu within its host galaxy. We further identify several diffuse interstellar bands and compare the 5780 Å band with the dust-to-band ratio for the Milky Way (MW). Overall, we find the foreground dust-extinction properties for SN 2012cu to be consistent with those of the MW. Furthermore, we find no evidence for significant time variation in any of these extinction tracers. We also compare the dust extinction curve models of Cardelli et al., O’Donnell, and Fitzpatrick, and find the predictions of Fitzpatrick fit SN 2012cu the best. Finally, the distance to NGC4772, the host of SN 2012cu, at a redshift of z = 0.0035, often assigned to the Virgo Southern Extension, is determined to be 16.6 ± 1.1 Mpc. We compare this result with distance measurements in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A regularized deconvolution method for turbulent closure modeling in implicitly filtered large-eddy simulation

Turbulence closure models that are based on a regularized deconvolution method (RDM) are proposed for application to implicitly filtered Large-Eddy Simulation (LES). This method reconstructs sub-grid scale contributions by approximately recovering filtered quantities through an optimization approach. Physical principles on the boundedness and conservation of reconstructed quantities are preserved by constraining the optimization problem. Further, to account for effects of scales that are not resolved in implicitly filtered LES, a reconstruction method is proposed that combines sub-grid projection and energy-similarity to enrich unresolved scales. With this method, closure models are developed for turbulent scalar fluxes and turbulence–chemistry interaction that are intrinsically consistent with the LES formulation. These closure models are evaluated in application to a piloted turbulent jet flame with inhomogeneous inlet composition, in which combustion is represented by a two-scalar manifold model. The RDM-formulation is benchmarked against simulations obtained using a presumed PDF method and eddy diffusivity closure model. In partially-premixed flame regions, it is observed that RDM improves the prediction of turbulent mixing, and faster grid convergence is observed for results with RDM compared to simulations employing a presumed PDF method.

42 ENGINEERING↗
DOE OSTI2025

Sustainable upcycling of polyethylene waste to compatibilizers and valuable chemicals

Controllable functionalization of polyethylene (PE) waste could generate new polymeric materials that are generally difficult to manufacture sustainably while also addressing the growing plastics waste problem. However, these modifications remain challenging due to the inherent stability of the PE backbone. Non-thermal atmospheric plasma enables molecular activation under mild conditions while utilizing renewable energy but is primarily employed for surface modification, as plasmas do not penetrate the bulk of materials. Herein, controllable bulk oxidative functionalization of PE wax (PEW) and low-density PE (LDPE) of varying molecular weights was achieved, with up to 6 mol% oxygen incorporation, by manipulating melt viscosity. This functionalization was accomplished either through temperature adjustment or by introducing a melt viscosity modifier, removable via simple extraction methods, to reduce LDPE viscosity, enhance diffusion and chain mobility, and enable bulk oxidation. The oxidized LDPE induces compatibilization in blends of poly(lactic acid) (PLA) and LDPE with improved interfacial adhesion and mechanical properties, such as a 70% increase in elongation-at-break values vs. the control. These findings pave the way for catalyst-free upcycling of direct plastics waste and plastics waste-derived products, enabling the creation of high-value products across various markets.

Plastic upcycling↗
DOE OSTI2015

IMPULSIVE SPOT HEATING AND THERMAL EXPLOSION OF INTERSTELLAR GRAINS REVISITED

The problem of the impulsive heating of dust grains in cold, dense interstellar clouds is revisited theoretically with the aim of better understanding the leading mechanisms of the explosive desorption of icy mantles. We rigorously show that if the heating of a reactive medium occurs within a sufficiently localized spot (e.g., the heating of mantles by cosmic rays (CRs)), then the subsequent thermal evolution is characterized by a single dimensionless number λ. This number identifies a bifurcation between two distinct regimes: when λ exceeds a critical value (threshold), the heat equation exhibits the explosive solution, i.e., the thermal (chemical) explosion is triggered. Otherwise, thermal diffusion causes the deposited heat to spread over the entire grain—this regime is commonly known as whole-grain heating. The theory allows us to find a critical combination of physical parameters that govern the explosion of icy mantles due to impulsive spot heating. In particular, our calculations suggest that heavy CR species (e.g., iron ions) colliding with dust are able to trigger the explosion. Based on recently calculated local CR spectra, we estimate the expected rate of explosive desorption. The efficiency of the desorption, which in principle affects all solid species independent of their binding energy, is shown to be comparable to other CR desorption mechanisms typically considered in the literature. Also, the theory allows us to estimate the maximum abundances of reactive species that may be stored in the mantles, which provides important constraints on the available astrochemical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Scalable hierarchical PDE sampler for generating spatially correlated random fields using nonmatching meshes: Scalable hierarchical PDE sampler using nonmatching meshes

This work describes a domain embedding technique between two nonmatching meshes used for generating realizations of spatially correlated random fields with applications to large-scale sampling-based uncertainty quantification. The goal is to apply the multilevel Monte Carlo (MLMC) method for the quantification of output uncertainties of PDEs with random input coefficients on general and unstructured computational domains. We propose a highly scalable, hierarchical sampling method to generate realizations of a Gaussian random field on a given unstructured mesh by solving a reaction–diffusion PDE with a stochastic right-hand side. The stochastic PDE is discretized using the mixed finite element method on an embedded domain with a structured mesh, and then, the solution is projected onto the unstructured mesh. This work describes implementation details on how to efficiently transfer data from the structured and unstructured meshes at coarse levels, assuming that this can be done efficiently on the finest level. We investigate the efficiency and parallel scalability of the technique for the scalable generation of Gaussian random fields in three dimensions. An application of the MLMC method is presented for quantifying uncertainties of subsurface flow problems. Here, we demonstrate the scalability of the sampling method with nonmatching mesh embedding, coupled with a parallel forward model problem solver, for large-scale 3D MLMC simulations with up to 1.9·109 unknowns.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗
DOE OSTI2019

High-speed 3D shape measurement using the optimized composite fringe patterns and stereo-assisted structured light system

In this study, we propose a high-speed 3D shape measurement technique based on the optimized composite fringe patterns and stereo-assisted structured light system. Stereo phase unwrapping, as a new-fashioned method for absolute phase retrieval based on the multi-view geometric constraints, can eliminate the phase ambiguities and obtain a continuous phase map without projecting any additional patterns. However, in order to ensure the stability of phase unwrapping, the period of fringe is generally around 20, which limits the accuracy of 3D measurement. To solve this problem, we develop an optimized method for designing the composite pattern, in which the speckle pattern is embedded into the conventional 4-step phase-shifting fringe patterns without compromising the fringe modulation, and thus the phase measurement accuracy. We also present a simple and effective evaluation criterion for the correlation quality of the designed speckle pattern in order to improve the matching accuracy significantly. When the embedded speckle pattern is demodulated, the periodic ambiguities in the wrapped phase can be eliminated by combining the adaptive window image correlation with geometry constraint. Finally, some mismatched regions are further corrected based on the proposed regional diffusion compensation technique (RDC). These proposed techniques constitute a complete computational framework that allows to effectively recover an accurate, unambiguous, and distortion-free 3D point cloud with only 4 projected patterns. Experimental results verify that our method can achieve high-speed, high-accuracy, robust 3D shape measurement with dense (64-period) fringe patterns at 5000 frames per second.

42 ENGINEERING↗
DOE OSTI2014

The isotropy problem of sub-ankle ultra high energy cosmic rays

We study the time dependent propagation of sub-ankle ultra high energy cosmic rays (UHECRs) originating from point-like Galactic sources. We show that drift in the Galactic magnetic field (GMF) may play an important role in the propagation of UHECRs and their measured anisotropy, particularly when the transport is anisotropic. To fully account for the discreteness of UHECR sources in space and time, a Monte Carlo method is used to randomly place sources in the Galaxy. The low anisotropy measured by Auger is not generally characteristic of the theoretical models, given that the sources are distributed in proportion to the star formation rate, but it can possibly be understood as (1) intermittency effects due to the discrete nature of the sources or, with extreme parameters, (2) a cancellation of drift current along a current sheet with outward radial diffusive flux. We conclude that it is possible to interpret the Galactic sub-ankle CR flux as being due entirely to intermittent discrete Galactic sources distributed in proportion to star formation, but only with a probability of roughly 35%, of which the spectrum is in accord with observations about 30% of the time. An alternative explanation for the low anisotropy may be that they are mostly extragalactic and/or heavy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Implementation and verification of a user-defined element (UEL) for coupled thermal-hydraulic-mechanical-chemical (THMC) processes in saturated geological media

Efficient and accurate modeling of the coupled thermal-hydraulic-mechanical-chemical (THMC) processes in various rock formations is indispensable for designing energy geo-structures such as underground repositories for high-level nuclear wastes. Here, this work focuses on developing and verifying an implicit finite element solver for generic coupled THMC problems in geological settings. Starting from the mass, momentum, and energy balance laws, a specialized set of governing equations and a thermoporoelastic constitutive model is derived. This system is then solved by an implicit finite element (FE) scheme. Specifically, the residuals and the Jacobians are scripted in a user-defined element (UEL) subroutine which is then combined with the general-purpose FE software Abaqus Standard to solve initial-boundary value problems. Considering the complexity of the system, the UEL development follows a stepwise manner by first solving the coupled hydraulic-mechanical (HM) and thermal-hydraulic-mechanical (THM) equations before moving on to the full THMC problem. Each implementation step consists of at least one verification test by comparing computed results with closed-form analytical solutions to ensure that the various coupling effects are correctly realized. To demonstrate the robustness of the algorithm and to validate the UEL, a three-dimensional case study is performed with reference to the in-situ heating test of ATLAS at Belgium in 1980s. A hypothetical radionuclide leakage event is then simulated by activating the chemical-concentration degree of freedom and prescribing a constant high concentration at the heater's surface. The model predicts a limited contaminated regime after six years considering both diffusion and advection effects on species transport.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Study of continuous-wave domain fluorescence diffuse optical tomography for quality control on agricultural produce

The importance of monitoring the quality of vegetables and fruits is prosperity by giving a competitive advantage for producer and providing a more healthy food for consumer. Diffuse Optical Tomography (DOT) is offering the possibility to detect the internal defects of the agricultural produce quality. Fluorescence diffuse optical tomography (FDOT) is the development of DOT, offering the possibilities to improve spatial resolution and to contrast image. The purpose of this research is to compare FDOT and DOT in forward analysis with continuous wave approach. The scattering and absorbing parameters of potatoes are used to represent the real condition. The object was illuminated by the NIR source from some positions on the boundary of object. A set of NIR detector are placed on the peripheral position of the object to measure the intensity of propagated or emitted light. In the simulation, we varied a condition of object then we analyzed the sensitivity of forward problem. The result of this study shows that FDOT has a better sensitivity than DOT and a better potential to monitor internal defects of agricultural produce because of the contrast value between optical and fluorescence properties of agricultural produce normal tissue and defects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A study of deformation and strain induced in bulk by the oxide layers formation on a Fe-Cr-Al alloy in high-temperature liquid Pb-Bi eutectic

At elevated temperatures, heavy liquid metals and their alloys are known to create a highly corrosive environment that causes irreversible degradation of most iron-based materials. In this paper, it has been found that an appropriate concentration of oxygen in the liquid alloy can significantly reduce this issue by creating a passivating oxide scale that controls diffusion, especially if Al is present in Fe-based materials (by Al-oxide formation). However, the increase of the temperature and of oxygen content in liquid phase leads to the increase of oxygen diffusion into bulk, and to promotion of the internal Al oxidation. This can cause a strain in bulk near the oxide layer, due either to mismatch between the thermal expansion coefficients of the oxides and bulk material, or to misfit of the crystal lattices (bulk vs. oxides). This work investigates the strain induced into proximal bulk of a Fe-Cr-Al alloy by oxide layers formation in liquid lead-bismuth eutectic utilizing synchrotron X-ray Laue microdiffraction. Finally, it is found that internal oxidation is the most likely cause for the strain in the metal rather than thermal expansion mismatch as a two-layer problem.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Entropy–viscosity method for the single material Euler equations in Lagrangian frame

In this work, a new finite element method for solving the Euler equations in Lagrangian coordinates is proposed. The method is stabilized by adding artificial diffusion terms compatible with positivity of mass and internal energy, a minimum principle on the specific entropy, and all generalized entropy inequalities. Two options of first-order artificial diffusion are considered. One is in the spirit of the Eulerian based method (Guermond et al., 2011 [23, 22]; Zingan et al., 2013) and the other is similar to existing viscosity stabilizations in Lagrangian frame, e.g., Dobrev et al. (2012). The method is verified to be high-order for smooth solutions even with active viscosity terms. This is achieved by using high-order finite element spaces and an entropy-based viscosity stabilization that degenerates the first-order viscous terms. This stabilization automatically distinguishes smooth and singular regions. The formal accuracy and convergence properties of the proposed methods are tested on a series of benchmark problems. This is the first result extending the entropy–viscosity methodology to the Lagrangian hydrodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Nonlinear phase field model for electrodeposition in electrochemical systems

A nonlinear phase-field model has been developed for describing the electrodeposition process in electrochemical systems that are highly out of equilibrium. Main thermodynamic driving forces for the electrode-electrolyte interface (EEI) evolution are limited to local variations of overpotential and ion concentration. Application of the model to Li-ion batteries describes the electrode interface motion and morphology change caused by charge mass transfer in the electrolyte, an electrochemical reaction at the EEI and cation deposition on the electrode surface during the charging operation. The Li electrodeposition rate follows the classical Butler-Volmer kinetics with exponentially and linearly depending on local overpotential and cation concentration at the electrode surface, respectively. Simulation results show that the Li deposit forms a fiber-like shape and grows parallel to the electric field direction. The longer and thicker deposits are observed both for higher current density and larger rate constant where the surface reaction rate is expected to be high. The proposed diffuse interface model well captures the metal electrodeposition phenomena in plenty of non-equilibrium electrochemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

HIGH-TEMPERATURE PROCESSING OF SOLIDS THROUGH SOLAR NEBULAR BOW SHOCKS: 3D RADIATION HYDRODYNAMICS SIMULATIONS WITH PARTICLES

A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H{sub 2} is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ∼few× 10{sup –8} L{sub ☉}. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.

79 ASTRONOMY AND ASTROPHYSICS↗