Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 883 records · Page 49

Plasticity of irradiated materials at the nano and micro-scales

Here, we review here our recent work on plastic deformation in irradiated materials at the nano- and micro-scales, as revealed by Discrete Dislocation Dynamics (DDD) simulations. Two methods of including irradiation effects in the DDD framework are presented. The first directly captures the atomistic interaction mechanisms, while the second can effectively study high-dose irradiation. Computer simulations lead to new understanding of the dynamics of collective dislocation-irradiation defect interactions, as well as the quantitative analysis of the temporal and spatial characteristics associated with plastic instabilities. Based on these insights, theoretical models are developed to predict the critical conditions for dislocation channel formation. A simple probability model is proposed and demonstrated to predict the width of dislocation channels in bulk irradiated materials with good agreement with experimental data. The fundamental understanding of the origins of plastic flow localization in irradiated materials sheds light on the design of future generations of radiation-resistant materials.

36 MATERIALS SCIENCE↗

Discrete breakup arcs and kinetic Alfven waves

Attention is confined to the extremely narrow auroral structures whose characteristic dimensions perpendicular to the magnetic fields are of the order of several hundred meters or even less. The starting point of the model used here is the assumption that discrete breakup arcs are caused by a change in the magnetosphere convection pattern that is localized in space. Particular attention is given to the shear Alfven mode. It is proposed that discrete breakup arcs are a direct manifestation of kinetic interaction between charged particles and an 'mhd' perturbation. It is stressed that the model applies only to narrow structures (thickness of the order of the plasma skin depth) and to short-lived phenomena (time scales shorter than a dispersion time). It is shown that a change of convection pattern in a spatially limited region leads to the emission of a shear-Alfven wave packet which propagates along the magnetic field line with a characteristic velocity.

Goertz, C. K.↗

Particle Size Distributions in Atmospheric Clouds

In this note, we derive a transport equation for a spatially integrated distribution function of particles size that is suitable for sparse particle systems, such as in atmospheric clouds. This is done by integrating a Boltzmann equation for a (local) distribution function over an arbitrary but finite volume. A methodology for evolving the moments of the integrated distribution is presented. These moments can be either tracked for a finite number of discrete populations ('clusters') or treated as continuum variables.

Paoli, Roberto↗

Development and Verification of the Charring, Ablating Thermal Protection Implicit System Simulator

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver (CATPISS) is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method (FEM) with first and second order fully implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton s method, while the linear system is solved via the Generalized Minimum Residual method (GMRES). Verification results from exact solutions and Method of Manufactured Solutions (MMS) are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗

DEM Modeling and Validation of Pebble Bed Packing Using Chrono::GPU

Accurate prediction of pebble packing structure is important for pebble bed reactors because the spatial distribution of void fraction directly affects coolant flow, pressure drop, heat transfer, and neutronic behavior. However, experimentally validated DEM studies that directly evaluate local void-fraction structure in reactor-relevant pebble beds remain limited. In this work, the pebble bed experiment conducted at Missouri University of Science and Technology is simulated using the graphics processing unit (GPU)-based discrete element method (DEM) code Chrono::GPU. The study focuses on evaluating the ability of Chrono::GPU to reproduce the packing arrangement and void-fraction distribution of a randomly packed spherical pebble bed. The DEM results are first verified against established radial void-fraction correlations, including the Mueller and Vortmeyer-Schuster models, to assess the predicted bulk porosity, near-wall behavior, and oscillatory packing structure. The simulation is then verified against reference DEM data and validated against gamma-ray computed tomography (CT) experimental data at three axial locations. The Chrono::GPU results reproduce the main features of the experimental packing, including the high void fraction near the wall, the first near-wall trough, and the damped oscillatory radial profile caused by wall-induced ordering. Quantitative comparison with DEM data and the CT-based radial profiles shows good agreement, with mean absolute errors on the order of 0.07 and root-mean-square errors below 0.09 for the averaged profiles. These results demonstrate that Chrono::GPU can accurately capture the void-fraction structure of spherical pebble beds and provides a reliable DEM framework for future pebble bed reactor packing, recycling, and thermal-hydraulic studies.

97 - MATHEMATICS AND COMPUTING↗

DCT quantization matrices visually optimized for individual images

This presentation describes how a vision model incorporating contrast sensitivity, contrast masking, and light adaptation is used to design visually optimal quantization matrices for Discrete Cosine Transform image compression. The Discrete Cosine Transform (DCT) underlies several image compression standards (JPEG, MPEG, H.261). The DCT is applied to 8x8 pixel blocks, and the resulting coefficients are quantized by division and rounding. The 8x8 'quantization matrix' of divisors determines the visual quality of the reconstructed image; the design of this matrix is left to the user. Since each DCT coefficient corresponds to a particular spatial frequency in a particular image region, each quantization error consists of a local increment or decrement in a particular frequency. After adjustments for contrast sensitivity, local light adaptation, and local contrast masking, this coefficient error can be converted to a just-noticeable-difference (jnd). The jnd's for different frequencies and image blocks can be pooled to yield a global perceptual error metric. With this metric, we can compute for each image the quantization matrix that minimizes the bit-rate for a given perceptual error, or perceptual error for a given bit-rate. Implementation of this system demonstrates its advantages over existing techniques. A unique feature of this scheme is that the quantization matrix is optimized for each individual image. This is compatible with the JPEG standard, which requires transmission of the quantization matrix.

Watson, Andrew B.↗

Origins of the 12-200 micron flux in NGC 6946 - Starlight and continuum dust emission from AN Sc galaxy

Infrared maps and optical CCD images of NGC 6946, a nearly face-on spiral with evidence of active star formation, are presented. Consideration is given, at 1-3 kpc spatial scales, to the relationship of the infrared continuum to spiral and radial structure, tracers of the interstellar medium, and ionizing and nonionizing sources of stellar radiation. Two infrared components were revealed - one with surface brightnesses of 20-60 Jy per 45 arcsec aperture at 160 microns, and the other with a surface brightness of less than 20 Jy. Various continuum distributions show discrete sources coinciding with giant H II regions superposed on a smooth unresolved disk background which decreases exponentially with the radius. The optical and spectral properties of the nuclear source and the disk are described. Dust properties and heating conditions, and far-infrared emitting phases are discussed.

Engargiola, G.↗

Effect of wall temperature spatial variations on aerodynamic heating to an actively cooled aircraft

The design of cooled skin panels, which requires consideration of many parameters such as heat load, skin thickness, panel size, tube spacing, tube size, temperature limits, coolant flow, etc., must minimize weight and thermal stresses. A design procedure for such panels has been developed by Bell Aerospace Company (Helenbrook and Anthony, 1971). This design procedure assumes a uniform heating rate to the panel, i.e., the variation of aerodynamic heating rate due to the cyclic wall temperature distribution, which results from the discrete coolant tube locations in such a panel, has been neglected. This study was made to determine if these variations were sufficiently large to cause variations in aerodynamic heating rate which would require coupling of boundary layer and conductive thermal analysis solutions.

Jones, R. A.↗

Field emission microscopy of carbon nanotube fibers: Evaluating and interpreting spatial emission

In this work, the authors quantify field emission properties of cathodes made from carbon nanotube (CNT) fibers. The cathodes were arranged in different configurations to determine the effect of cathode geometry on the emission properties. Various geometries were investigated including (1) flat cut fiber tip, (2) folded fiber, (3) looped fiber, and (4) fibers wound around a cylinder. The authors employ a custom field emission microscope to quantify I-V characteristics in combination with laterally resolved field-dependent electron emission area. Additionally, they look at the very early emission stages, first when a CNT fiber is turned on for the first time, which is then followed by multiple ramp-up/down runs. Upon the first turn on, all fibers demonstrated limited and discrete emission area. During ramping runs, all CNT fibers underwent multiple (minor and/or major) breakdowns, which improved emission properties in that turn-on field decreased and field enhancement factor and emission area both increased. It is proposed that breakdowns are responsible for removing initially undesirable emission sites caused by stray fibers higher than average. This initial breakdown process gives way to a larger emission area that is created when the CNT fiber subcomponents unfold and align with the electric field. Overall, the authors' results form the basis for careful evaluation of CNT fiber cathodes for dc or low frequency pulsed power systems in which large uniform area emission is required or for narrow beam high frequency applications in which high brightness is a must.

Image processing↗

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗

Bilevel Parameter Optimization for Nonlocal Image Denoising Models

We propose a bilevel optimization approach for the determination of parameters in nonlocal image denoising. We consider both spatial weights in front of the fidelity term, as well as weights within the kernel of the nonlocal operator. In both cases we investigate the differentiability of the solution operator in function spaces and derive a first order optimality system that characterizes local minima. For the numerical solution of the problems, we propose a second-order optimization algorithm in combination with a finite element discretization of the nonlocal denoising models and a computational strategy for the solution of the resulting dense linear systems. Several experiments are run in order to show the suitability of our approach.

97 MATHEMATICS AND COMPUTING↗

Scattering and Diffraction of Electromagnetic Radiation: An Effective Probe to Material Structure

Scattered electromagnetic waves from material bodies of different forms contain, in an intricate way, precise information on the intrinsic, geometrical and physical properties of the objects. Scattering theories, ever deepening, aim to provide dependable interpretation and prediction to the complicated interaction of electromagnetic radiation with matter. There are well-established multiple-scattering formulations based on classical electromagnetic theories. An example is the Generalized Multi-particle Mie-solution (GMM), which has recently been extended to a special version ̶ the GMM-PA approach, applicable to finite periodic arrays consisting of a huge number (e.g., >>106) of identical scattering centers [1]. The framework of the GMM-PA is nearly complete. When the size of the constituent unit scatterers becomes considerably small in comparison with incident wavelength, an appropriate array of such small element volumes may well be a satisfactory representation of a material entity having an arbitrary structure. X-ray diffraction is a powerful characterization tool used in a variety of scientific and technical fields, including material science. A diffraction pattern is nothing more than the spatial distribution of scattered intensity, determined by the distribution of scattering matter by way of its Fourier transform [1]. Since all linear dimensions entered into Maxwell's equations are normalized by wavelength, an analogy exists between optical and X-ray diffraction patterns. A large set of optical diffraction patterns experimentally obtained can be found in the literature [e.g., 2,3]. Theoretical results from the GMM-PA have been scrutinized using a large collection of publically accessible, experimentally obtained Fraunhofer diffraction patterns. As far as characteristic structures of the patterns are concerned, theoretical and experimental results are in uniform agreement; no exception has been found so far. Closely connected with the spatial distribution of scattered intensities are cross sections, such as for extinction, scattering, absorption, and radiation pressure, as a critical type of key quantity addressed in most theoretical and experimental studies of radiative scattering. Cross sections predicted from different scattering theories are supposed to be in general agreement. For objects of irregular shape, the GMM-PA solutions can be compared with the highly flexible Discrete Dipole Approximation (DDA) [4,5] when dividing a target to no more than ~106 unit cells. Also, there are different ways to calculate the cross sections in the GMM-PA, providing an additional means to examine the accuracy of the numerical solutions and to unveil potential issues concerning the theoretical formulations and numerical aspects. To solve multiple scattering by an assembly of material volumes through classical theories such as the GMM-PA, the radiative properties of the component scatterers, the complex refractive index in particular, must be provided as input parameters. When using a PA to characterize a material body, this involves the use of an adequate theoretical tool, an effective medium theory, to connect Maxwell's phenomenogical theory with the atomistic theory of matter. In the atomic theory, one regards matter as composed of interacting particles (atoms and molecules) embedded in the vacuum [6]. However, the radiative properties of atomic-scaled particles are known to be substantially different from bulk materials. Intensive research efforts in the fields of cluster science and nanoscience attempt to bridge the gap between bulk and atom and to understand the transition from classical to quantum physics. The GMM-PA calculations, which place virtually no restriction on the component-particle size, might help to gain certain insight into the transition.

Xu, Yu-Lin↗

Discrete scale invariance of the quasi-bound states at atomic vacancies in a topological material

Recently, log-periodic quantum oscillations have been detected in the topological materials zirconium pentatelluride (ZrTe 5 ) and hafnium pentatelluride (HfTe 5 ), displaying an intriguing discrete scale invariance (DSI) characteristic. In condensed materials, the DSI is considered to be related to the quasi-bound states formed by massless Dirac fermions with strong Coulomb attraction, offering a feasible platform to study the long-pursued atomic-collapse phenomenon. Here, we demonstrate that a variety of atomic vacancies in the topological material HfTe 5 can host the geometric quasi-bound states with a DSI feature, resembling an artificial supercritical atom collapse. The density of states of these quasi-bound states is enhanced, and the quasi-bound states are spatially distributed in the “orbitals” surrounding the vacancy sites, which are detected and visualized by low-temperature scanning tunneling microscope/spectroscopy. By applying the perpendicular magnetic fields, the quasi-bound states at lower energies become wider and eventually invisible; meanwhile, the energies of quasi-bound states move gradually toward the Fermi energy (E F ). These features are consistent with the theoretical prediction of a magnetic field–induced transition from supercritical to subcritical states. The direct observation of geometric quasi-bound states sheds light on the deep understanding of the DSI in quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An immersed interface method for the 2D vorticity-velocity Navier-Stokes equations with multiple bodies

We present an immersed interface method for the vorticity-velocity form of the 2D Navier Stokes equations that directly addresses challenges posed by nonconvex immersed bodies, multiply connected domains, and the calculation of force distributions on immersed surfaces. The immersed interface method is re-interpreted as a polynomial extrapolation of flow quantities and boundary conditions into the immersed solid bodies, reducing computational cost and enabling simulations with nonconvex bodies that could not be discretized with previous immersed interface methods. In the flow, the vorticity transport equation is discretized using a conservative finite difference scheme and explicit Runge-Kutta time integration. The velocity reconstruction problem is transformed to a scalar Poisson equation that is discretized with conservative finite differences, and solved using an FFT-accelerated iterative algorithm. The use of conservative differencing throughout leads to exact enforcement of a discrete Kelvin's theorem, allowing for simulations with multiply connected domains and outflow boundaries that have challenged other immersed interface vortex methods. We also explore novel methods for recovering time-dependent pressure distributions on immersed bodies within a vorticity-based method and present a novel control volume formulation for recovering aerodynamic moments from only the vorticity and velocity fields. The method achieves second order spatial accuracy and third order temporal accuracy, and is validated on a variety of 2D flows in internal and free-space domains.

97 MATHEMATICS AND COMPUTING↗

Matching Curved Lattices to Anisotropic Tangent Planes

Radial quantization would be the ideal formalism for studying strongly-coupled near-conformal quantum field theories but it requires the ability to perform lattice calculations on static, curved manifolds, specifically a very long cylinder whose cross section is a sphere. Smoothly discretizing the surface of a sphere requires a graph with unequal edge lengths. The geometry of such graphs is well understood since 1961 using Regge Calculus. But, lattice quantum field theories are defined in terms of couplings which appear in the action rather than edge lengths and so the relationship between couplings and lengths must be determined dynamically. A simple example is computing the ratio of spatial to temporal lattice spacings in anisotropic lattice QCD. I will discuss our conjecture that computing anisotropic lattice spacing ratios on affine transformations of regular flat lattices is sufficient to determine coupling assignments on curved lattices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computer measurement and representation of the heart in two and three dimensions

Methods for the measurement and display by minicomputer of cardiac images obtained from fluoroscopy to permit an accurate assessment of functional changes are discussed. Heart contours and discrete points can be digitized automatically or manually, with the recorded image in a video, cine, or print format. As each frame is digitized it is assigned a code name identifying the data source, experiment, run, view, and frame, and the images are filed for future reference in any sequence. Two views taken at the same point in the heart cycle are used to compute the spatial position of the ventricle apex and the midpoint of the aortic valve. The remainder of the points on the chamber border are corrected for the linear distortion of the X-rays by projection to a plane containing the chord between the apex and the aortic valve center and oriented so that lines perpendicular to the chord are parallel to the image intensifier face. The image of the chamber surface is obtained by generating circular cross sections with diameters perpendicular to the major chord. The transformed two- and three-dimensional imagery can be displayed in either static or animated form using a graphics terminal.

Rasmussen, D.↗

Steady-State and Transient Boundary Element Methods for Coupled Heat Conduction

Boundary element algorithms for the solution of steady-state and transient heat conduction are presented. The algorithms are designed for efficient coupling with computational fluid dynamic discretizations and feature piecewise linear elements with offset nodal points. The steady-state algorithm employs the fundamental solution approach; the integration kernels are computed analytically based on linear shape functions, linear elements, and variably offset nodal points. The analytic expressions for both singular and nonsingular integrands are presented. The transient algorithm employs the transient fundamental solution; the temporal integration is performed analytically and the nonsingular spatial integration is performed numerically using Gaussian quadrature. A series solution to the integration is derived for the instance of a singular integrand. The boundary-only character of the algorithm is maintained by integrating the influence coefficients from initial time. Numerical results are compared to analytical solutions to verify the current boundary element algorithms. The steady-state and transient algorithms are numerically shown to be second-order accurate in space and time, respectively.

Kontinos, Dean A.↗

Electron Energization and Structure of the Diffusion Region During Asymmetric Reconnection

Results from particle-in-cell simulations of reconnection with asymmetric upstream conditions are reported to elucidate electron energization and structure of the electron diffusion region (EDR). Acceleration of unmagnetized electrons results in discrete structures in the distribution functions and supports the intense current and perpendicular heating in the EDR. The accelerated electrons are cyclotron turned by the reconnected magnetic field to produce the outflow jets, and as such, the acceleration by the reconnection electric field is limited, leading to resistivity without particle-particle or particle-wave collisions. A map of electron distributions is constructed, and its spatial evolution is compared with quantities previously proposed to be EDR identifiers to enable effective identifications of the EDR in terrestrial magnetopause reconnection.

magnetic reconnection↗