Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “anisotropic diffusion”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quasi-one-dimensional Ising-like antiferromagnetism in the rare-earth perovskite oxide TbScO 3

The rare-earth perovskite TbScO 3 has been widely used as a substrate for the growth of epitaxial ferroelectric and multiferroic thin films, while its detailed low-temperature magnetic properties were rarely reported. In this paper, we performed detailed magnetization, specific heat, and single crystal neutron scattering measurements, along with the crystalline electric field calculations to study the low-temperature magnetic properties of TbScO 3. All our results suggest the magnetic Tb 3+ has an Ising-like pseudo-doublet ground state at low temperatures. Due to the constrain of local point symmetry, these Tb 3+ Ising moments are confined in the ab plane with a tilt angle of φ=±48° to the a axis. In zero field, the system undergoes an antiferromagnetic phase transition at T N =2.53 K, and forms a G x A y noncollinear magnetic structure below T N . Here, we find the dipole-dipole interactions play an important role to determine the magnetic ground state, which are also responsible for the quasi-one-dimensional magnetism in TbScO 3 . The significant anisotropic diffuse scatterings further confirm the quasi-one-dimensional magnetism along the c axis. The magnetic phase diagram with the field along the easy b axis is well established. In addition to the G x A y antiferromagnetic state, there is an exotic field-induced phase emerged near the critical field B c ≃0.7 T, where three-dimensional magnetic order is suppressed but strong one-dimensional correlations may still exist.

36 MATERIALS SCIENCE↗
DOE OSTI2014

MO-G-17A-03: MR-Based Cortical Bone Segmentation for PET Attenuation Correction with a Non-UTE 3D Fast GRE Sequence

Purpose: To determine the feasibility of identifying cortical bone on MR images with a short-TE 3D fast-GRE sequence for attenuation correction of PET data in PET/MR. Methods: A water-fat-bone phantom was constructed with two pieces of beef shank. MR scans were performed on a 3T MR scanner (GE Discovery™ MR750). A 3D GRE sequence was first employed to measure the level of residual signal in cortical bone (TE{sub 1}/TE{sub 2}/TE{sub 3}=2.2/4.4/6.6ms, TR=20ms, flip angle=25°). For cortical bone segmentation, a 3D fast-GRE sequence (TE/TR=0.7/1.9ms, acquisition voxel size=2.5×2.5×3mm{sup 3}) was implemented along with a 3D Dixon sequence (TE{sub 1}/TE{sub 2}/TR=1.2/2.3/4.0ms, acquisition voxel size=1.25×1.25×3mm{sup 3}) for water/fat imaging. Flip angle (10°), acquisition bandwidth (250kHz), FOV (480×480×144mm{sup 3}) and reconstructed voxel size (0.94×0.94×1.5mm{sup 3}) were kept the same for both sequences. Soft tissue and fat tissue were first segmented on the reconstructed water/fat image. A tissue mask was created by combining the segmented water/fat masks, which was then applied on the fast-GRE image (MRFGRE). A second mask was created to remove the Gibbs artifacts present in regions in close vicinity to the phantom. MRFGRE data was smoothed with a 3D anisotropic diffusion filter for noise reduction, after which cortical bone and air was separated using a threshold determined from the histogram. Results: There is signal in the cortical bone region in the 3D GRE images, indicating the possibility of separating cortical bone and air based on signal intensity from short-TE MR image. The acquisition time for the 3D fast-GRE sequence was 17s, which can be reduced to less than 10s with parallel imaging. The attenuation image created from water-fat-bone segmentation is visually similar compared to reference CT. Conclusion: Cortical bone and air can be separated based on intensity in MR image with a short-TE 3D fast-GRE sequence. Further research is required to optimize the strategy to reduce Gibbs artifacts.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

WE-G-BRD-07: Automated MR Image Standardization and Auto-Contouring Strategy for MRI-Based Adaptive Brachytherapy for Cervix Cancer

Purpose: MRI-based adaptive brachytherapy (ABT) is an emerging treatment modality for patients with gynecological tumors. However, MR image intensity non-uniformities (IINU) can vary from fraction to fraction, complicating image interpretation and auto-contouring accuracy. We demonstrate here an automated MR image standardization and auto-contouring strategy for MRI-based ABT of cervix cancer. Methods: MR image standardization consisted of: 1) IINU correction using the MNI N3 algorithm, 2) noise filtering using anisotropic diffusion, and 3) signal intensity normalization using the volumetric median. This post-processing chain was implemented as a series of custom Matlab and Java extensions in MIM (v6.4.5, MIM Software) and was applied to 3D T2 SPACE images of six patients undergoing MRI-based ABT at 3T. Coefficients of variation (CV=σ/µ) were calculated for both original and standardized images and compared using Mann-Whitney tests. Patient-specific cumulative MR atlases of bladder, rectum, and sigmoid contours were constructed throughout ABT, using original and standardized MR images from all previous ABT fractions. Auto-contouring was performed in MIM two ways: 1) best-match of one atlas image to the daily MR image, 2) multi-match of all previous fraction atlas images to the daily MR image. Dice’s Similarity Coefficients (DSCs) were calculated for auto-generated contours relative to reference contours for both original and standardized MR images and compared using Mann-Whitney tests. Results: Significant improvements in CV were detected following MR image standardization (p=0.0043), demonstrating an improvement in MR image uniformity. DSCs consistently increased for auto-contoured bladder, rectum, and sigmoid following MR image standardization, with the highest DSCs detected when the combination of MR image standardization and multi-match cumulative atlas-based auto-contouring was utilized. Conclusion: MR image standardization significantly improves MR image uniformity. The combination of MR image standardization and multi-match cumulative atlas-based auto-contouring produced the highest DSCs and is a promising strategy for MRI-based ABT for cervix cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-C-207B-05: Tissue Segmentation of Computed Tomography Images Using a Random Forest Algorithm: A Feasibility Study

Purpose: Develop an automated Random Forest algorithm for tissue segmentation of CT examinations. Methods: Seven materials were classified for segmentation: background, lung/internal gas, fat, muscle, solid organ parenchyma, blood/contrast, and bone using Matlab and the Trainable Weka Segmentation (TWS) plugin of FIJI. The following classifier feature filters of TWS were investigated: minimum, maximum, mean, and variance each evaluated over a pixel radius of 2n, (n = 0–4). Also noise reduction and edge preserving filters, Gaussian, bilateral, Kuwahara, and anisotropic diffusion, were evaluated. The algorithm used 200 trees with 2 features per node. A training data set was established using an anonymized patient’s (male, 20 yr, 72 kg) chest-abdomen-pelvis CT examination. To establish segmentation ground truth, the training data were manually segmented using Eclipse planning software, and an intra-observer reproducibility test was conducted. Six additional patient data sets were segmented based on classifier data generated from the training data. Accuracy of segmentation was determined by calculating the Dice similarity coefficient (DSC) between manual and auto segmented images. Results: The optimized autosegmentation algorithm resulted in 16 features calculated using maximum, mean, variance, and Gaussian blur filters with kernel radii of 1, 2, and 4 pixels, in addition to the original CT number, and Kuwahara filter (linear kernel of 19 pixels). Ground truth had a DSC of 0.94 (range: 0.90–0.99) for adult and 0.92 (range: 0.85–0.99) for pediatric data sets across all seven segmentation classes. The automated algorithm produced segmentation with an average DSC of 0.85 ± 0.04 (range: 0.81–1.00) for the adult patients, and 0.86 ± 0.03 (range: 0.80–0.99) for the pediatric patients. Conclusion: The TWS Random Forest auto-segmentation algorithm was optimized for CT environment, and able to segment seven material classes over a range of body habitus and CT protocol parameters with an average DSC of 0.86 ± 0.04 (range: 0.80–0.99).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Multilevel Convergence Analysis of Multigrid-Reduction-in-Time

This study presents a multilevel convergence framework for multigrid-reduction-in-time (MGRIT) as a generalization of previous two-grid estimates. The framework provides a priori upper bounds on the convergence of MGRIT V- and F-cycles, with different relaxation schemes, by deriving the respective residual and error propagation operators. The residual and error operators are functions of the time-stepping operator, analyzed directly and bounded in the norm, both numerically and analytically. We present various upper bounds of different computational cost and varying sharpness. These upper bounds are complemented by proposing analytic formulae for the approximate convergence factor of V-cycle algorithms that take the number of fine grid time points, the temporal coarsening factors, and the eigenvalues of the time-stepping operator as parameters. The paper concludes with supporting numerical investigations of parabolic (anisotropic diffusion) and hyperbolic (wave equation) model problems. We assess the sharpness of the bounds and the quality of the approximate convergence factors. Observations from these numerical investigations demonstrate the value of the proposed multilevel convergence framework for estimating MGRIT convergence a priori and for the design of a convergent algorithm. We further highlight that observations in the literature are captured by the theory, including that two-level Parareal and multilevel MGRIT with F-relaxation do not yield scalable algorithms and the benefit of a stronger relaxation scheme. An important observation is that with increasing numbers of levels MGRIT convergence deteriorates for the hyperbolic model problem, while constant convergence factors can be achieved for the diffusion equation. The theory also indicates that L-stable Runge--Kutta schemes are more amendable to multilevel parallel-in-time integration with MGRIT than A-stable Runge--Kutta schemes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A two-level GPU-accelerated incomplete LU preconditioner for general sparse linear systems

This paper presents a parallel preconditioning approach based on incomplete LU (ILU) factorizations in the framework of Domain Decomposition (DD) for general sparse linear systems. We focus on distributed memory parallel architectures, specifically, those that are equipped with graphic processing units (GPUs). In addition to block-Jacobi, we present general purpose two-level ILU Schur complement-based approaches, where different strategies are presented to solve the coarse-level reduced system. These strategies are combined with modified ILU methods in the construction of the coarse-level operator, in order to effectively remove smooth errors by targeting an algebraically smooth vector. We leverage available GPU-based sparse matrix kernels to accelerate the setup and the solve phases of the proposed ILU preconditioner. We evaluate the efficiency of the proposed methods as a smoother for algebraic multigrid (AMG) and as a preconditioner for Krylov subspace methods on challenging anisotropic diffusion problems and a collection of general sparse matrices.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Updates to the Generation of Physics Data Inputs for MAMMOTH Simulations of the Transient Reactor Test Facility - FY2016

The INL is currently evolving the modeling and simulation (M&S) capability that will enable improved core operation as well as design and analysis of TREAT experiments. This M&S capability primarily uses MAMMOTH, a reactor physics application being developed under Multi-physics Object Oriented Simulation Environment (MOOSE) framework. MAMMOTH allows the coupling of a number of other MOOSE-based applications. This second year of work has been devoted to the generation of a deterministic reference solution for the full core, the preparation of anisotropic diffusion coefficients, the testing of the SPH equivalence method, and the improvement of the control rod modeling. In addition, this report includes the progress made in the modeling of the M8 core configuration and experiment vehicle since January of this year.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

GLOBAL SIMULATIONS OF GALACTIC WINDS INCLUDING COSMIC-RAY STREAMING

Galactic outflows play an important role in galactic evolution. Despite their importance, a detailed understanding of the physical mechanisms responsible for the driving of these winds is lacking. In an effort to gain more insight into the nature of these flows, we perform global three-dimensional magnetohydrodynamical simulations of an isolated Milky Way-size starburst galaxy. We focus on the dynamical role of cosmic rays (CRs) injected by supernovae, and specifically on the impact of the streaming and anisotropic diffusion of CRs along the magnetic fields. We find that these microphysical effects can have a significant effect on the wind launching and mass loading factors, depending on the details of the plasma physics. Due to the CR streaming instability, CRs propagating in the interstellar medium scatter on self-excited Alfvén waves and couple to the gas. When the wave growth due to the streaming instability is inhibited by some damping process, such as turbulent damping, the coupling of CRs to the gas is weaker and their effective propagation speed faster than the Alfvén speed. Alternatively, CRs could scatter from “extrinsic turbulence” that is driven by another mechanism. We demonstrate that the presence of moderately super-Alfvénic CR streaming enhances the efficiency of galactic wind driving. Cosmic rays stream away from denser regions near the galactic disk along partially ordered magnetic fields and in the process accelerate more tenuous gas away from the galaxy. For CR acceleration efficiencies broadly consistent with the observational constraints, CRs reduce the galactic star formation rates and significantly aid in launching galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1971

Transport of solar flare protons: Comparison of a new analytic model with spacecraft measurements

An analytic solution was obtained to the complete Fokker-Planck equation including the effects of convection, interplanetary deceleration and acceleration, corotation, and anisotropic diffusion. The solution yields a time to maximum for the particle flux of about 10 h and an exponential decay time of about 5 h. Several solar flare particle events were observed with the solar and galactic cosmic ray experiment on OGO-6. Comparisons of the calculated time dependence of the fluxes with these observations of 1 to 70 MeV protons show that the model adequately describes both the rise and decay times.

Lupton, J. E.↗
NASA NTRS1971

The propagation of solar cosmic-ray bursts.

Description of a model showing analytically the three phases of anisotropy which occur during solar cosmic-ray events in the 7.5 to 21 MeV kinetic-energy interval and reported by McCracken et al. (1971): (1) a highly anisotropic, near field-aligned, initial phase, (2) a convective phase, and (3) a late-time phase in which the anisotropy is approximately perpendicular to the mean interplanetary magnetic field. The model is based on the cosmic-ray particles being convectively transported out from the sun, undergoing anisotropic diffusion along the interplanetary magnetic-field lines, and losing energy by adiabatic deceleration or by collision processes. The event is seen simply as a pulse moving outward from the sun after a cosmic-ray burst with a negative density-gradient in front of it and a positive gradient behind.

Ng, C. K.↗
NASA NTRS1981

Temporal Aperture Modulation

The two types of modulation techniques useful to X-ray imaging are reviewed. The use of optimum coded temporal aperature modulation is shown, in certain cases, to offer an advantage over a spatial aperture modulator. Example applications of a diffuse anisotropic X-ray background experiment and a wide field of view hard X-ray imager are discussed.

Proctor, R. J.↗
NASA NTRS1985

The expected cosmic ray density and stream distributions at the heliolatitudinal asymmetry of solar wind

The results of the spatial distribution of cosmic ray density, gradients, and anisotropy obtained on the basis of the numerical solution of the anisotropic diffusion equation with an account of solar wind velocity change depending on the latitudinal angle theta of the form U=u sub oe sup alpha theta and the diffusion coefficient depending on the spatial coordinates and the particle rigidity are presented. It is shown that the increase of the solar wind velocity and the diffusion coefficient with heliolatitude leads to gradient distributions that are in accord with experimental data observed in space. The results of the energetic spectrum of 11 and 22-year cosmic ray variations obtained with an account of direction of the general magnetic field of the Sun are presented are given.

Alania, M. V.↗
NASA NTRS1985

Spatial distribution of cosmic rays in magnetic cycles of the Sun

The space-time distribution of cosmic ray (CR) density is determined from the solution of the anisotropic diffusion equation. The boundary problem is solved in a quasistationary axisymmetric case, taking into account the delay of cosmic rays relative to the cyclic variations of solar activity. The specificty of the modulation model is the use of the empirical dependence of the free path of particles for scattering on the observed latitudinal distribution of solar activity. Another specific feature of the model is an account of an additional effect due to the drift of charged particles towards the helio-equator plane or backward, depending on the direction of the total magnetic field of the Sun.

Alania, M. V.↗
NASA NTRS1985

The modulation features of the long-period cosmic ray variations in connection with the sign change of the general magnetic field of the Sun

On the basis of the model and experimental investigations, the spatial distribution of cosmic ray anisotropy for the different epochs of solar activity is studied. A solution is offered to the anisotropic diffusion equation with regard to the electromagnetic conditions for the periods of minimum and maximum solar activity and cosmic ray particle drift in a regular interplanetary magnetic field. It is shown that the long period changes amplitude and phase of the diurnal variations of cosmic rays is limited not only by convection and diffusion of particles but also by the drift effect before and after the sing change of the general magnetic field of the Sun. The calculated model is compared with the results obtained on the basis of an analysis of the experimental data from the neutron super monitor station, and it is shown that for the periods when the lines of magnetic force of the Sun come from the Northern Hemisphere the phase of the first harmonic diurnal variation is shifted forwards to an earlier time.

Iskra, K.↗
NASA NTRS1985

The theoretical and experimental investigation of cosmic ray Forbush-effects

The theoretical results of analysis of the expected spatial distributions of density, gradients and anisotropy of cosmic rays, obtained on the basis of the numerical solution of anisotropic diffusion equation in the presence of the disturbances of shock wave type in the interplanetary space are presented. The theoretical calculations on the definition of the energy spectrum and anisotropy of galactic cosmic rays during Forbush effect are compared with the experimental data.

Alania, M. V.↗
DOE OSTI2018

Formation of superlattice structure in irradiated materials: a phase field study

The formation and self-organization of fission gas bubbles under irradiation is a unique phenomenon, which affects the swelling behavior and mechanical stability of nuclear materials. We develop a phase field model to investigate the formation mechanism of fission gas bubbles, and associated material properties under different irradiation condition. The thermodynamic and kinetic origin of superlattice formation is systematically analyzed. In particular, the symmetry of the superlattice is dictated by the anisotropic diffusion of interstitial atoms, while the lattice parameter is determined by the instability of vacancy and gas concentration waves. The spatial distribution of gas bubbles strongly depends on the nucleation, growth and coarsening process, which leads to the formation of self-organized superlattice. This study provides a new insight into designing target experiments to tailor desired microstructure under irradiation.

36 MATERIALS SCIENCE↗