Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse domain approach”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Optimizing Sub-Sampled STEM Imaging for Beam Sensitive Materials and Dynamic Processes

The goal of this research is to quantify how sub-sampling and inpainting (broadly classified as methods in Artificial Intelligence) can be used to obtain the highest resolution images under any experimental conditions, and to use this to reproducibly control dynamic in-situ processes in gases and liquids on the nanoscale. As atomic diffusion under non-extreme conditions is typically in the sub-ms time domain, the ability to use this sub-sampling approach to reduce the effect of the electron beam and to speed-up the imaging process represents the optimum combination of spatial and temporal resolution achievable in a conventional microscope. This LDRD funding demonstrated the feasibility of imaging dynamics in complex energy systems using this methodology by establishing the optimum experimental sampling conditions (optimizing the number of pixels, speed, contrast and resolution) for each experiment. For in-situ liquid cell experiments this analysis shows that simply by distributing the same dose and dose rate maximally in space and time, the kinetics for radiolysis products can be modified significantly. The final stage of this work, i.e. the experimental demonstration of its use for nucleation and growth, was severely curtailed by issues with COVID19. However, all the experiments are designed and currently underway, with the goal to complete the final demonstration and publication as soon as safety concerns for laboratory access are met.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

The 3-D numerical study of airflow in the compressor/combustor prediffuser and dump diffuser of an industrial gas turbine

This paper describes the 3D computations of a flow field in the compressor/combustor diffusers of an industrial gas turbine. The geometry considered includes components such as the combustor support strut, the transition piece and the impingement sleeve with discrete cooling air holes on its surface. Because the geometry was complex and 3D, the airflow path was divided into two computational domains sharing an interface region. The body-fitted grid was generated independently in each of the two domains. The governing equations for incompressible Navier-Stokes equations were solved using the finite volume approach. The results show that the flow in the prediffuser is strongly coupled with the flow in the dump diffuser and vice versa. The computations also revealed that the flow in the dump diffuser is highly nonuniform.

Agrawal, Ajay K.↗
DOE OSTI2020

Dirichlet absorbing boundary conditions for classical and peridynamic diffusion-type models

Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of the PD diffusion models proposed so far are applied to bounded domains only. In this study, we propose an effective way to handle unbounded domains both with PD and classical diffusion models. For the former, we employ a meshfree discretization, whereas for the latter the finite element method (FEM) is employed. The proposed ABCs are time-dependent and Dirichlet-type, making the approach easy to implement in the available models. The performance of the approach, in terms of accuracy and stability, is illustrated by numerical examples in 1D, 2D, and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Effects of local order parameter dependent transport coefficient in diblock copolymers under applied electric fields

In this work, we present an approach for constructing thermodynamically consistent time-dependent models relevant to thin films of diblock copolymers in applied electric fields. The approach is based on the principles of linear irreversible thermodynamics, and, in this work, it is applied to study the effects of electric fields on thin films of incompressible diblock copolymers. Enforcement of local incompressibility constraint at all times leads to a local order parameter dependent transport coefficient in the model for the diblock copolymers. The dependence of the transport coefficient on the local order parameter is used to relate it with the diffusion constant of Rouse chains and leads to sensitivity of the model to initial conditions. In addition, transient behavior is found to be affected when compared with an ad hoc model assuming a constant transport coefficient. Numerical results such as electric field induced alignment of lamellae domains due to the field are found to be in qualitative agreement with experiments. This approach opens up a systematic way of developing kinetic models for simulating effects of electrolytes added to thin films containing diblock copolymers in the presence of applied electric fields.

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

Sub‐Nanosecond Reconfiguration of Ferroelectric Domains in Bismuth Ferrite

Abstract Domain switching is crucial for achieving desired functions in ferroic materials that are used in various applications. Fast control of domains at sub‐nanosecond timescales remains a challenge despite its potential for high‐speed operation in random‐access memories, photonic, and nanoelectronic devices. Here, ultrafast laser excitation is shown to transiently melt and reconfigure ferroelectric stripe domains in multiferroic bismuth ferrite on a timescale faster than 100 picoseconds. This dynamic behavior is visualized by picosecond‐ and nanometer‐resolved X‐ray diffraction and time‐resolved X‐ray diffuse scattering. The disordering of stripe domains is attributed to the screening of depolarization fields by photogenerated carriers resulting in the formation of charged domain walls, as supported by phase‐field simulations. Furthermore, the recovery of disordered domains exhibits subdiffusive growth on nanosecond timescales, with a non‐equilibrium domain velocity reaching up to 10 m s −1 . These findings present a new approach to image and manipulate ferroelectric domains on sub‐nanosecond timescales, which can be further extended into other complex photoferroic systems to modulate their electronic, optical, and magnetic properties beyond gigahertz frequencies. This approach could pave the way for high‐speed ferroelectric data storage and computing, and, more broadly, defines new approaches for visualizing the non‐equilibrium dynamics of heterogeneous and disordered materials.

36 MATERIALS SCIENCE↗
NASA NTRS1985

Solution method for a hovering helicopter rotor using the Euler equations

A method for the calculation of the flow field of a helicopter rotor blade in hover is presented. The approach uses a finite volume solution of the three dimensional Euler equations for the blade near field. In the Euler solver the velocity field is decomposed into two parts. One is the induced velocity of the vortex wake extending below the blade, found from a free wake calculation procedure. The other part is the unknown additional velocity field of the rotor blade. This approach eliminates numerical diffusion of the rolled up wake vorticity due to truncation error and artificial viscosity. Also, the effects of the far wake are included in the limited computational domain. Solutions are presented for an isolated wing and a model helicopter rotor and compared to experiment.

Roberts, T. W.↗
DOE OSTI2023

Role of Intra-Domain Heterogeneity on Ion and Polymer Dynamics in Block Polymer Electrolytes: An Approach for Spatially Resolving Dynamics and Ion Transport

The design of safe and high-performance, nanostructured, block polymer (BP) electrolytes for lithium-ion batteries requires a thorough understanding of the key parameters that govern local structure and dynamics. Yet, the interfaces between microphase separated domains can introduce complexities in this local behavior that can be challenging to quantify. Herein, the local polymer, cation (Li + ), and anion dynamics were described in salt-doped polystyrene block poly(oligo-oxyethylene methyl ether methacrylate) (PS-b-POEM) through a quantitative framework that considered the effects of polymer architecture, segmental mixing, chain stretching, and confinement on polymer mobility and ion transport. This framework was validated through nuclear magnetic resonance (NMR) spectroscopy measurements on solid (dry) polymer electrolyte samples. Notably, a mobility transition temperature (T mobility ) was identified through NMR spectroscopy that captured the local dynamics more accurately than the thermal glass transition temperature. Additionally, the approach quantitatively described the mobility gradient across a domain when segmental mixing effects were combined with chain stretching and confinement information, especially at higher segregation strengths – facilitating the assessment of local ion diffusion and conductivity. As a result, spatially averaged local ion diffusion predictions quantitatively matched NMR-measured ion diffusivities in the BP samples, while spatially summed ionic conductivity predictions across a domain qualitatively captured trends in the measured ionic conductivities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Critical fluid thermal equilibration experiment (19-IML-1)

Gravity sometimes blocks all experimental techniques of making a desired measurement. Any pure fluid possesses a liquid-vapor critical point. It is defined by a temperature, pressure, and density state in thermodynamics. The critical issue that this experiment attempts to understand is the time it takes for a sample to reach temperature and density equilibrium as the critical point is approached; is it infinity due to mass and thermal diffusion, or do pressure waves speed up energy transport while mass is still under diffusion control. The objectives are to observe: (1) large phase domain homogenization without and with stirring; (2) time evolution of heat and mass after temperature step is applied to a one phase equilibrium sample; (3) phase evolution and configuration upon going two phase from a one phase equilibrium state; (4) effects of stirring on a low g two phase configuration; (5) two phase to one phase healing dynamics starting from a two phase low g configuration; and (6) effects of shuttle acceleration events on spatially and temporally varying compressible critical fluid dynamics.

Wilkinson, R. Allen↗
NASA NTRS1997

A Model Study of Zonal Forcing in the Equatorial Stratosphere by Convectively Induced Gravity Waves

A two-dimensional cloud-resolving model is used to examine the possible role of gravity waves generated by a simulated tropical squall line in forcing the quasi-biennial oscillation (QBO) of the zonal winds in the equatorial stratosphere. A simulation with constant background stratospheric winds is compared to simulations with background winds characteristic of the westerly and easterly QBO phases, respectively. In all three cases a broad spectrum of both eastward and westward propagating gravity waves is excited. In the constant background wind case the vertical momentum flux is nearly constant with height in the stratosphere, after correction for waves leaving the model domain. In the easterly and westerly shear cases, however, westward and eastward propagating waves, respectively, are strongly damped as they approach their critical levels, owing to the strongly scale-dependent vertical diffusion in the model. The profiles of zonal forcing induced by this wave damping are similar to profiles given by critical level absorption, but displaced slightly downward. The magnitude of the zonal forcing is of order 5 m/s/day. It is estimated that if 2% of the area of the Tropics were occupied by storms of similar magnitude, mesoscale gravity waves could provide nearly 1/4 of the zonal forcing required for the QBO.

Alexander, M. J.↗
DOE OSTI2024

Visualizing the Internal Nanocrystallinity of Calcite Due to Nonclassical Crystallization by 3D Coherent X–Ray Diffraction Imaging

The internal crystallinity of calcite is investigated for samples synthesized using two approaches: precipitation from solution and the ammonium carbonate diffusion method. Scanning electron microscopy (SEM) analyses reveal that the calcite products precipitated using both approaches have a well–defined rhombohedron shape, consistent with the euhedral crystal habit of the mineral. The internal structure of these calcite crystals is characterized using Bragg coherent diffraction imaging (BCDI) to determine the 3D electron density and the atomic displacement field. BCDI reconstructions for crystals synthesized using the ammonium carbonate diffusion approach have the expected euhedral shape, with internal strain fields and few internal defects. In contrast, the crystals synthesized by precipitation from solution have very complex external shapes and defective internal structures, presenting null electron density regions and pronounced displacement field distributions. These heterogeneities are interpreted as multiple crystalline domains, created by a nonclassical crystallization mechanism, where smaller nanoparticles coalescence into the final euhedral particles. The combined use of SEM, X–ray diffraction (XRD), and BCDI allows for structurally differentiating calcite crystals grown with different approaches, opening new opportunities to understand how grain boundaries and internal defects alter calcite reactivity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Long-time integration of parametric evolution equations with physics-informed DeepONets

Ordinary and partial differential equations (ODEs/PDEs) play a paramount role in analyzing and simulating complex dynamic processes across all corners of science and engineering. In recent years machine learning tools are aspiring to introduce new effective ways of simulating such equations, however existing approaches are not able to reliably return stable and accurate predictions across long temporal horizons. We aim to address this challenge by introducing an effective framework for learning evolution operators that map random initial conditions to associated ODE/PDE solutions within a short time interval. Such operators can be parametrized by deep neural networks that are trained in an entirely self-supervised manner without requiring one to generate any paired input-output observations. Global long-time predictions across a range of initial conditions can be then obtained by iteratively evaluating the trained model using each prediction as the initial condition for the next evaluation step. Here, this introduces a new approach to temporal domain decomposition that is shown to be effective in performing accurate long-time simulations for a wide range of parametric ODE and PDE systems, from wave propagation, to reaction-diffusion dynamics and stiff chemical kinetics, introducing a new way of rapidly emulating non-equilibrium processes in science and engineering.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Optimal Experimental Design for Inverse Problems in the Presence of Observation Correlations

Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.

Correlated observations↗
DOE OSTI2021

An efficient reconstruction algorithm for diffusion on triangular grids using the nodal discontinuous Galerkin method

High-energy-density (HED) hydrodynamics studies such as those relevant to inertial confinement fusion and astrophysics require highly disparate densities, temperatures, viscosities, and other diffusion parameters over relatively short spatial scales. This presents a challenge for high-order accurate methods to effectively resolve the hydrodynamics at these scales, particularly in the presence of highly disparate diffusion. A significant volume of engineering and physics applications use an unstructured discontinuous Galerkin (DG) method developed based on the finite element mesh generation and algorithmic framework. This work discusses the application of an affine reconstructed nodal DG method for unstructured grids of triangles. Solving the diffusion terms in the DG method is non-trivial due to the solution representations being piecewise continuous. Hence, the diffusive flux is not defined on the interface of elements. The proposed numerical approach reconstructs a smooth solution in a parallelogram that is enclosed by the quadrilateral formed by two adjacent triangle elements. The interface between these two triangles is the diagonal of the enclosed parallelogram. Similar to triangles, the mapping of parallelograms from a physical domain to a reference domain is an affine mapping, which is necessary for an accurate and efficient implementation of the numerical algorithm. Thus, all computations can still be performed on the reference domain, which promotes efficiency in computation and storage. This reconstruction does not make assumptions on choice of polynomial basis. Reconstructed DG algorithms have previously been developed for modal implementations of the convection–diffusion equations. However, to the best of the authors’ knowledge, this is the first practical guideline that has been proposed for applying the reconstructed algorithm on a nodal discontinuous Galerkin method with a focus on accuracy and efficiency. As a result, the algorithm is demonstrated on a number of benchmark cases as well as a challenging substantive problem in HED hydrodynamics with highly disparate diffusion parameters.

Computational efficiency↗
NASA NTRS1988

Boundary identification for 2-D parabolic problems arising in thermal testing of materials

Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in an output least square approach. Parameter estimation techniques based on the method of mappings are discussed, and approximation schemes are developed based on a finite-element Galerkin approach. Theoretical convergence results for computational techniques are given, and the results are applied to the identification of two kinds of boundary shapes.

Banks, H. T.↗
NASA NTRS1989

Boundary estimation problems arising in thermal tomography

Problems on the identification of two-dimensional spatial domains arising in the detection and characterization of structural flaws in materials are considered. For a thermal diffusion system with external boundary input, observations of the temperature on the surface are used in a output least squares approach. Parameter estimation techniques based on the method of mappings are discussed and approximation schemes are developed based on a finite element Galerkin approach. Theoretical convergence results for computational techniques are given and the results are applied to experimental data for the identification of flaws in the thermal testing of materials.

Banks, H. T.↗
DOE OSTI2022

Monte-Carlo modeling of phonon thermal transport using DFT-based anisotropic dispersion relations over the full Brillouin zone

In this work, we present a Monte Carlo (MC) approach to solve the phonon Boltzmann transport equation (BTE) in which the anisotropic phonon dispersion relations over the full Brillouin zone (BZ) are used. In this approach, the discretization of the BZ used to compute the phonon relaxation time places constraints on the direction of scattered phonons in the real-space simulation domain. The phonon dispersion and phonon relaxation times are calculated using the density functional theory (DFT) approach. The modified MC approach is validated by a close examination of its ability to simulate phonon transport in both the ballistic and diffusive regimes for multiple materials including GaAs, InAs, ThO 2 , and α-U. In doing so, the phonon thermal conductivities from 100 K to 1000 K are calculated and compared with traditional non-transport solution of the phonon BTE. It is found that the phonon thermal conductivities of α-U and ThO 2 obtained from MC simulations using isotropic dispersion are larger than the values obtained using anisotropic phonon dispersion relations over the full BZ. The effect of phonon-defect scattering on the thermal conductivity of ThO 2 is also studied as an application of the current MC approach and found to agree with previously computed values in the literature. The MC solver developed here has been parallelized as a step to demonstrate its potential to solving computationally intensive phonon thermal transport problems at the mesoscale.

36 MATERIALS SCIENCE↗