Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference method”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Experimental Characterization and Numerical Modeling of the Interaction Between Carbon Fiber Composite Prepregs During a Preforming Process

We report carbon fiber reinforced composites have received growing attention because of their superior performance and high potential for lightweight systems. An economic method to manufacture the parts made of these composites is a sequence of forming followed by a compression molding. The first step in this sequence is called preforming that forms the prepreg, which is the fabric impregnated with the uncured resin, to the product geometry, while the molding process cures the resin. Slip between different prepreg layers is observed in the preforming step, and it is believed to have a non-negligible impact on the resulting geometry. This paper reports a method to characterize the interaction between different prepreg layers, which should be valuable for future predictive modeling and design optimization. An experimental device was built to evaluate the interactions with respect to various industrial production conditions. The experimental results were analyzed for an in-depth understanding about how temperature, relative sliding speed, and fiber orientation affect the tangential interaction between two prepreg layers. Moreover, a hydro-lubricant model was introduced to study the relative motion mechanism of this fabric-resin-fabric system, and the results agreed well with the experiment data. Lastly, the interaction factors obtained from this research will be implemented in a preforming process finite element simulation model.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Stabilization of numerical interchange in spectral-element magnetohydrodynamics

In this study, auxiliary numerical projections of the divergence of flow velocity and vorticity parallel to magnetic field are developed and tested for the purpose of suppressing unphysical interchange instability in magnetohydrodynamic simulations. The numerical instability arises with equal-order C 0 finite- and spectral-element expansions of the flow velocity, magnetic field, and pressure and is sensitive to behavior at the limit of resolution. The auxiliary projections are motivated by physical field-line bending, and coercive responses to the projections are added to the flow-velocity equation. Their incomplete expansions are limited to the highest-order orthogonal polynomial in at least one coordinate of the spectral elements. Cylindrical eigenmode computations show that the projections induce convergence from the stable side with first-order ideal-MHD equations during h-refinement and p-refinement. Hyperbolic and parabolic projections and responses are compared, together with different methods for avoiding magnetic divergence error. Lastly, the projections are also shown to be effective in linear and nonlinear time-dependent computations with the NIMROD code [C. R. Sovinec, et al., J. Comput. Phys. 195 (2004) 355-386], provided that the projections introduce numerical dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The phase diagram of quantum chromodynamics in one dimension on a quantum computer

The quantum chromodynamics (QCD) phase diagram, which reveals the state of strongly interacting matter at different temperatures and densities, is key to answering open questions in physics, ranging from the behaviour of particles in neutron stars to the conditions of the early universe. However, classical simulations of QCD face significant computational barriers, such as the sign problem at finite matter densities. Quantum computing offers a promising solution to overcome these challenges. Here, we take an important step toward exploring the QCD phase diagram with quantum devices by preparing thermal states in one-dimensional non-Abelian gauge theories. We experimentally simulate the thermal states of SU(2) and SU(3) gauge theories at finite densities on a trapped-ion quantum computer using a variational method. This is achieved by introducing two features: Firstly, we add motional ancillae to the existing qubit register to efficiently prepare thermal probability distributions. Secondly, we introduce charge-singlet measurements to enforce colour-neutrality constraints. This work pioneers the quantum simulation of QCD at finite density and temperature for two and three colours, laying the foundation to explore QCD phenomena on quantum platforms.

Quantum information↗
DOE OSTI2016

WE-DE-201-12: Thermal and Dosimetric Properties of a Ferrite-Based Thermo-Brachytherapy Seed

Purpose: The novel thermo-brachytherapy (TB) seed provides a simple means of adding hyperthermia to LDR prostate permanent implant brachytherapy. The high blood perfusion rate (BPR) within the prostate motivates the use of the ferrite and conductive outer layer design for the seed cores. We describe the results of computational analyses of the thermal properties of this ferrite-based TB seed in modelled patient-specific anatomy, as well as studies of the interseed and scatter (ISA) effect. Methods: The anatomies (including the thermophysical properties of the main tissue types) and seed distributions of 6 prostate patients who had been treated with LDR brachytherapy seeds were modelled in the finite element analysis software COMSOL, using ferrite-based TB and additional hyperthermia-only (HT-only) seeds. The resulting temperature distributions were compared to those computed for patient-specific seed distributions, but in uniform anatomy with a constant blood perfusion rate. The ISA effect was quantified in the Monte Carlo software package MCNP5. Results: Compared with temperature distributions calculated in modelled uniform tissue, temperature distributions in the patient-specific anatomy were higher and more heterogeneous. Moreover, the maximum temperature to the rectal wall was typically ∼1 °C greater for patient-specific anatomy than for uniform anatomy. The ISA effect of the TB and HT-only seeds caused a reduction in D90 similar to that found for previously-investigated NiCu-based seeds, but of a slightly smaller magnitude. Conclusion: The differences between temperature distributions computed for uniform and patient-specific anatomy for ferrite-based seeds are significant enough that heterogeneous anatomy should be considered. Both types of modelling indicate that ferrite-based seeds provide sufficiently high and uniform hyperthermia to the prostate, without excessively heating surrounding tissues. The ISA effect of these seeds is slightly less than that for the previously-presented NiCu-based seeds.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

SPIDER beam dump as diagnostic of the particle beam

The beam power produced by the negative ion source for the production of ion of deuterium extracted from RF plasma is mainly absorbed by the beam dump component which has been designed also for measuring the temperatures on the dumping panels for beam diagnostics. A finite element code has been developed to characterize, by thermo-hydraulic analysis, the sensitivity of the beam dump to the different beam parameters. The results prove the capability of diagnosing the beam divergence and the horizontal misalignment, while the entity of the halo fraction appears hardly detectable without considering the other foreseen diagnostics like tomography and beam emission spectroscopy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

Hidden Mn magnetic-moment disorder and its influence on the physical properties of medium-entropy NiCoMn solid solution alloys

In this work, the ab initio Korringa-Kohn-Rostoker method combined with the coherent potential approximation (CPA) was employed to investigate the electronic, magnetic, and transport properties of medium-entropy face-centered-cubic (fcc) NiCoMn solid solution alloys. By comparing the CPA electronic structure with that from supercell calculations, we uncovered an unconventional CPA ground state, which correctly distinguishes two equally populated Mn CPA components—with large spin moments but opposite orientations. Using the spin spiral calculations, we further demonstrated that this ground state is most energetically favorable in the presence of spin noncollinearity, and no significant longitudinal spin fluctuation is observed, justifying the applicability of the Heisenberg model. The finite-temperature magnetism was further studied using different approximations based on the Heisenberg model, and we found the Mn moments to be fully disordered at low temperature due to a small net effective Weiss field on Mn. In addition, the magnetic effect on the electron scattering at finite temperatures was evaluated and compared with other scattering mechanisms. Since the magnetization-induced electron scattering is almost saturated in the ground state, (full) spin disorder only yields a small addition to the resistivity, whereas the thermal displacements increase it modestly. Finally, we elucidate the role of hydrostatic pressure on the magnetic and transport properties. In conclusion, these findings reflect the importance of the magnetic signatures on the physical properties of alloys, and they provide a window into magnetism-controlled electronic structure and energy dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Investigation of assumptions and approximations in the virtual fields method for a viscoplastic material model

The Virtual Fields Method (VFM) is an inverse technique used for parameter estimation and calibration of constitutive models. Many assumptions and approximations—such as plane stress, incompressible plasticity, and spatial and temporal derivative calculations—are required to use VFM with full-field deformation data, for example, from Digital Image Correlation (DIC). Here, this work presents a comprehensive discussion of the effects of these assumptions and approximations on parameters identified by VFM for a viscoplastic material model for 304L stainless steel. We generated synthetic data from a Finite-Element Analysis (FEA) in order to have a reference solution with a known material model and known model parameters, and we investigated four cases in which successively more assumptions and approximations were included in the data. We found that VFM is tolerant to small deviations from the plane stress condition in a small region of the sample, and that the incompressible plasticity assumption can be used to estimate thickness changes with little error. A local polynomial fit to the displacement data was successfully employed to compute the spatial displacement gradients. The choice of temporal derivative approximation (i.e., backwards difference versus central difference) was found to have a significant influence on the computed rate of deformation and on the VFM results for the rate-dependent model used in this work. Finally, the noise introduced into the displacement data from a stereo-DIC simulator was found to have negligible influence on the VFM results. Evaluating the effects of assumptions and approximations using synthetic data is a critical first step for verifying and validating VFM for specific applications. Finally, the results of this work provide the foundation for confidently using VFM for experimental data.

304L stainless steel↗
DOE OSTI2021

Momentum distribution of the uniform electron gas at finite temperature: Effects of spin polarization

We carry out extensive direct path integral Monte Carlo (PIMC) simulations of the uniform electron gas (UEG) at finite temperature for different values of the spin-polarization ξ. This allows us to unambiguously quantify the impact of spin effects on the momentum distribution function n(k) and related properties. We find that interesting physical effects like the interaction-induced increase in the occupation of the zero-momentum state n(0) substantially depend on ξ. Our results further advance the current understanding of the UEG as a fundamental model system, and are of practical relevance for the description of transport properties of warm dense matter in an external magnetic field. All PIMC results are freely available online and can be used as a benchmark for the development of methods and applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

(U) Estimating the Photonics Budget, Resolution, and Signal Requirements for a Multi-Monochromatic X-ray Imager

This document examines the performance of a generic flat-mirror multimonochromatic imager (MMI), with special emphasis on existing instruments at NIF and Omega. We begin by deriving the standard equation for the mean number of photons detected per resolution element. The pinhole energy bandwidth is a contributing factor; this is dominated by the finite size of the source and may be considerable. The most common method for estimating the spatial resolution of such a system (quadrature addition) is, technically, mathematically invalid for this case. However, under the proper circumstances it may produce good estimates compared to a rigorous calculation based on the convolution of point-spread functions. Diffraction is an important contribution to the spatial resolution. Common approximations based on Fraunhofer (farfield) diffraction may be inappropriate and misleading, as the instrument may reside in multiple regimes depending upon its configuration or the energy of interest. It is crucial to identify the correct diffraction regime; Fraunhofer and Fresnel (near-field) diffraction profiles are substantially different, the latter being considerably wider. Finally, we combine the photonics and resolution analyses to derive an expression for the minimum signal level such that the resulting images are not dominated by photon statistics. This analysis is consistent with observed performance of the NIF MMI.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

A micro-scale cutting model for UD CFRP composites with thermo-mechanical coupling

Cutting a unidirectional carbon fiber-reinforced polymer (UD CFRP) structure is the basic unit for CFRP machining, which is a complex thermal-mechanically coupled process. To reveal the deformation mechanism and predict cutting force in UD CFRP micro cutting, a micro-scale fracture model for UD CFRP cutting with thermal-mechanical coupling is demonstrated in this paper, which captures the failure modes for fibers, matrix and the interface based on a micro-level RVE using a relatively simple damage based fracture method. The thermal-mechanical coupling model at the micro scale is developed on the basis of the plastic energy dissipation and frictional heating during cutting. Failure models for the fiber, matrix and interface region are applied depending on the material properties of each of these three phases. Numerical simulations based on the above model with different fiber orientations were performed to predict the deformation and forces of different components in UD CFRP. Cutting experiments with the same fiber orientations as considered in the simulations were carried out to validate the force and deformation results. The predicted force and deformation patterns match well with evidence from our experiments. In general, the cutting force is larger than the thrust force regardless of fiber orientation. The cutting force reaches a maximum as the fiber orientation approaches 90 , but thrust forces do not vary substantially across cases. When the fiber orientation is acute, the deformation of fibers is much smaller than when the cutting angle is obtuse. Surface roughness follows the same trend with cutting angle as fiber deformation.

42 ENGINEERING↗
DOE OSTI2020

Fast GPU 3D diffeomorphic image registration

3D image registration is one of the most fundamental and computationally expensive operations in medical image analysis. Here, we present a mixed-precision, Gauss–Newton–Krylov solver for diffeomorphic registration of two images. Our work extends the publicly available CLAIRE library to GPU architectures. Despite the importance of image registration, only a few implementations of large deformation diffeomorphic registration packages support GPUs. Our contributions are new algorithms to significantly reduce the run time of the two main computational kernels in CLAIRE: calculation of derivatives and scattered-data interpolation. Additionally, we deploy (i) highly-optimized, mixed-precision GPU-kernels for the evaluation of scattered-data interpolation, (ii) replace Fast-Fourier-Transform (FFT)-based first-order derivatives with optimized 8th-order finite differences, and (iii) compare with state-of-the-art CPU and GPU implementations. As a highlight, we demonstrate that we can register clinical images in less than 6 s on a single NVIDIA Tesla V100. This amounts to over 20 speed-up over the current version of CLAIRE and over 30 speed-up over existing GPU implementations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

EMPIRE-PIC: A Performance Portable Unstructured Particle-in-Cell Code

In this study we introduce EMPIRE-PIC, a finite element method particle-in-cell (FEM-PIC) application developed at Sandia National Laboratories. The code has been developed in C++ using the Trilinos library and the Kokkos Performance Portability Framework to enable running on multiple modern compute architectures while only requiring maintenance of a single codebase. EMPIRE-PIC is capable of solving both electrostatic and electromagnetic problems in two- and three-dimensions to second-order accuracy in space and time. In this paper we validate the code against three benchmark problems — a simple electron orbit, an electrostatic Langmuir wave, and a transverse electromagnetic wave propagating through a plasma. We demonstrate the performance of EMPIRE-PIC on four different architectures: Intel Haswell CPUs, Intel's Xeon Phi Knights Landing, ARM Thunder-X2 CPUs, and NVIDIA Tesla V100 GPUs attached to IBM POWER9 processors. This analysis demonstrates scalability of the code up to more than two thousand GPUs, and greater than one hundred thousand CPUs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Automated property optimization via ab initio O(N) elongation method: Application to (hyper-)polarizability in DNA

An automated property optimization method was developed based on the ab initio O(N) elongation (ELG) method and applied to the optimization of nonlinear optical (NLO) properties in DNA as a first test. The ELG method mimics a polymerization reaction on a computer, and the reaction terminal of a starting cluster is attacked by monomers sequentially to elongate the electronic structure of the system by solving in each step a limited space including the terminal (localized molecular orbitals at the terminal) and monomer. The ELG-finite field (ELG-FF) method for calculating (hyper-)polarizabilities was used as the engine program of the optimization method, and it was found to show linear scaling efficiency while maintaining high computational accuracy for a random sequenced DNA model. Furthermore, the self-consistent field convergence was significantly improved by using the ELG-FF method compared with a conventional method, and it can lead to more feasible NLO property values in the FF treatment. The automated optimization method successfully chose an appropriate base pair from four base pairs (A, T, G, and C) for each elongation step according to an evaluation function. From test optimizations for the first order hyper-polarizability (β) in DNA, a substantial difference was observed depending on optimization conditions between “choose-maximum” (choose a base pair giving the maximum β for each step) and “choose-minimum” (choose a base pair giving the minimum β). In contrast, there was an ambiguous difference between these conditions for optimizing the second order hyper-polarizability (γ) because of the small absolute value of γ and the limitation of numerical differential calculations in the FF method. It can be concluded that the ab initio level property optimization method introduced here can be an effective step towards an advanced computer aided material design method as long as the numerical limitation of the FF method is taken into account.

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

Experimentally Validated Computational Modeling of Creep and Creep-Cracking for Nuclear Concrete Structures

In a Nuclear Power Plant, one of the most important components is the concrete nuclear reactor cavity, which serves both a structural and protective function as the biological radiation shield. Given that creep has been identified as a major knowledge gap in the assessment of nuclear structures (NUREG/CR-7153), this work helps to further the understanding of creep behavior of massive concrete containment structures for decades to enable safe and long-term operation of these facilities. This project has developed a robust, experimentally validated model to predict creep in nuclear concrete structures for up 60 years using short-term creep data thereby enabling a longer service life of critical facilities and early detection of structural failure. The work presented in this report is a pairing of computational and experimental methods. For the first time, the time temperature superposition (TTS) principle was successfully used to generate a uniaxial creep compliance master curve to predict mortar creep response for up to 22,500 days (nearly 60 years) at a reference temperature of 20°C. These data were used as input into finite element analysis (FEA) codes that use highly realistic random, 3D concrete microstructures from reconstructed coarse limestone aggregates. Finite element analysis performed provides the ability to quickly upscale mortar viscoelastic behavior to long-term concrete creep/relaxation data. A master creep compliance curve, constructed from the TTS principle, spanning 27 years, was used to validate two and a half decades of simulated concrete creep. Concurrently, three different simulated wall specimens were designed to mimic the behavior of post-tensioned concrete nuclear containment facility vessel walls over time as a result of concrete creep. The specimens were designed with different thicknesses, transverse and longitudinal reinforcement ratios, and level of post-tensioning stress. Each specimen contained various instrumentation to measure internal concrete temperature, concrete strain, and post-tensioning strain hourly for over 3 years. The concrete creep model developed in this project, based on the FEA concrete simulations, was applied to simulate the structural-scale experiments of prestressed concrete walls conducted in this project using the Grizzly code. These models can represent the effects of reinforcing and prestressing. Although there are some discrepancies with the experimental data, the model can predict the overall trends of the creep response in these experiments. One of these experimental models was also applied to an extended time to demonstrate how the findings from this study can be used to predict the behavior of actual structures of interest that have been in service for extended periods of time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Accurate and precise determination of critical properties from Gibbs ensemble Monte Carlo simulations

Since the seminal paper by Panagiotopoulos [Mol. Phys. 61, 813 (1997)], the Gibbs ensemble Monte Carlo (GEMC) method has been the most popular particle-based simulation approach for the computation of vapor–liquid phase equilibria. However, the validity of GEMC simulations in the near-critical region has been questioned because rigorous finite-size scaling approaches cannot be applied to simulations with fluctuating volume. Valleau [Mol. Simul. 29, 627 (2003)] has argued that GEMC simulations would lead to a spurious overestimation of the critical temperature. More recently, Patel et al. [J. Chem. Phys. 134, 024101 (2011)] opined that the use of analytical tail corrections would be problematic in the near-critical region. To address these issues, we perform extensive GEMC simulations for Lennard-Jones particles in the near-critical region varying the system size, the overall system density, and the cutoff distance. For a system with N = 5500 particles, potential truncation at 8σ and analytical tail corrections, an extrapolation of GEMC simulation data at temperatures in the range from 1.27 to 1.305 yields T{sub c} = 1.3128 ± 0.0016, ρ{sub c} = 0.316 ± 0.004, and p{sub c} = 0.1274 ± 0.0013 in excellent agreement with the thermodynamic limit determined by Potoff and Panagiotopoulos [J. Chem. Phys. 109, 10914 (1998)] using grand canonical Monte Carlo simulations and finite-size scaling. Critical properties estimated using GEMC simulations with different overall system densities (0.296 ≤ ρ{sub t} ≤ 0.336) agree to within the statistical uncertainties. For simulations with tail corrections, data obtained using r{sub cut} = 3.5σ yield T{sub c} and p{sub c} that are higher by 0.2% and 1.4% than simulations with r{sub cut} = 5 and 8σ but still with overlapping 95% confidence intervals. In contrast, GEMC simulations with a truncated and shifted potential show that r{sub cut} = 8σ is insufficient to obtain accurate results. Additional GEMC simulations for hard-core square-well particles with various ranges of the attractive well and for n-decane molecules represented by the TraPPE force field yield data that support the trends observed for Lennard-Jones particles. The finite-size dependence of the critical properties obtained from GEMC simulations is significantly smaller than those from grand-canonical ensemble simulations. Thus, when resources are not available for a rigorous finite-size scaling study, GEMC simulations provide a straightforward route to determine fairly accurate critical properties using relatively small system sizes.

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

On convergent finite difference schemes for variational–PDE-based image processing

We study an adaptive anisotropic Huber functional-based image restoration scheme. Using a combination of L2–L1 regularization functions, an adaptive Huber functional-based energy minimization model provides denoising with edge preservation in noisy digital images. We study a convergent finite difference scheme based on continuous piecewise linear functions and use a variable splitting scheme, namely the Split Bregman (In: Goldstein and Osher, SIAM J Imaging Sci 2(2):323–343, 2009) algorithm, to obtain the discrete minimizer. Experimental results are given in image denoising and comparison with additive operator splitting, dual fixed point, and projected gradient schemes illustrates that the best convergence rates are obtained for our algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Symmetric spectrum current drive due to finite radial drift effects

We report the drift surfaces of minority heated ions differ from flux surfaces due to finite poloidal gyroradius effects. As the minority poloidal gyroradius approaches radial scale lengths in the plasma, the difference between drift and flux surfaces can modify the heating and lead to a symmetric spectrum minority counter-current being driven. In response, a corresponding overall net co-current of comparable size is driven. This beneficial symmetric spectrum current drive in a tokamak is due to the parallel velocity asymmetry in the drift departure from a flux surface. As this new source of driven current is a side effect of minority heating it comes without any additional economic cost to reactor power balance. The symmetric spectrum current driven for near Maxwellian minorities is evaluated by an adjoint method and found to be modest. However, minority heating typically results in strong non-Maxwellian features on minority distributions so it may be possible to drive a significantly larger co-current. A related evaluation is performed for alpha particles in a deuterium minority heated plasma with a tritium majority. The low density of the alphas tends to keep this driven symmetric spectrum current small, but at very high heating levels a significant co-current might be driven. Other mechanisms to drive co-current with a symmetric spectrum are discussed and estimated, including asymmetric electron drag and focusing of the applied minority heating radio frequency fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗