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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Seismic modeling with radial basis function-generated finite differences (RBF-FD) – a simplified treatment of interfaces

In a previous study of seismic modeling with radial basis function-generated finite differences (RBF-FD), we outlined a numerical method for solving 2-D wave equations in domains with material interfaces between different regions. The method was applicable on a mesh-free set of data nodes. It included all information about interfaces within the weights of the stencils (allowing the use of traditional time integrators), and was shown to solve problems of the 2-D elastic wave equation to 3rd-order accuracy. In the present paper, we discuss a refinement of that method that makes it simpler to implement. It can also improve accuracy for the case of smoothly-variable model parameter values near interfaces. We give several test cases that demonstrate the method solving 2-D elastic wave equation problems to 4th-order accuracy, even in the presence of smoothly-curved interfaces with jump discontinuities in the model parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

On the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers: Selective mode-type excitation and specific mode enhancement

Here, we investigate the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers. Non-uniform polarization is represented as a non-uniform electromechanical coupling coefficient using a polarization function. Governing equations are derived for the forced in-plane vibrations of a thin wafer under the assumption of generalized plane stress. The effect of non-uniform polarization is explicitly obtained in the forcing terms of the governing equations. These equations are then recast into a variational weak form that is then solved using the finite element method to obtain the displacement fields for different modes. The electromechanical response of a non-uniformly polarized piezoelectric wafer is derived in terms of the out-of-plane displacement profile on the surface of the wafer. Using the derived analytical expression, a necessary and sufficient condition for the presence/absence of a vibrational mode in the electromechanical impedance spectrum is obtained. Based on this condition, criteria for selective mode-type excitation and specific mode enhancement of vibrational modes in the electromechanical impedance spectrum are postulated. Selective mode-type excitation of in-plane extensional and shear modes is demonstrated for a square wafer and that of in-plane bending modes is demonstrated for a rectangular wafer. Specific mode enhancement is demonstrated for both square and rectangular wafers. In addition, it is also demonstrated how the criteria can be used to suppress specific vibrational modes in the electromechanical impedance spectrum. The proposed methodology of using non-uniformly polarized piezoelectric wafers finds application in the design of single element transducers with multi-frequency operation, frequency-tuned receivers/sensors, acoustic holograms, designing acoustic beams of prescribed shape/lobes, and other non-traditional applications such as information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Model for describing plasmonic nanolasers using Maxwell-Liouville equations with finite-difference time-domain calculations

In this work, we present a theoretical study of lasing action when plasmonic metallic structures that show lattice plasmon resonances are embedded in a gain medium. Our model combines classical electrodynamics for arrays of gold nanoparticles with a four-level quantum Liouville model of the laser dye photophysics. A numerical solution was implemented using finite-difference time-domain calculations coupled with a finite-difference solution to the Liouville equation. A particular focus of this work is the influence of dephasing in the quantum dynamics on the emission intensity at the threshold for lasing. We find that dephasing in the quantum system leads to reduced lasing emission, but with little effect on the long-term population inversion. Both electronic and vibrational dephasing is considered, but only electronic dephasing is significant, with the fully dephased result appearing for dephasing times comparable to plasmon dephasing (~10 fs) while fully coherent results involve >100 ps dephasing times as determined by the rate of stimulated emission. There are factor-of-2 differences between the Maxwell-Liouville results (greater emission intensities and narrower widths) compared to the corresponding results of rate-equation models of the dye states, which indicates the importance of using the Maxwell-Liouville approach in modeling these systems. We also examine rate-equation models with and without constraints arising from the Pauli exclusion principle, and we find relatively small effects.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Enhancing two-photon spontaneous emission in rare earths using graphene and graphene nanoribbons

The enhancement of two-photon spontaneous emission (2PSE) from trivalent and divalent rare earth ions in proximity to graphene and graphene nanoribbons is calculated for achievable experimental conditions using a combination of finite difference time domain simulations and direct computation of transition rates between energy levels in rare earths. For Er 3+ , we find that the 2PSE rate is initially 8 orders lower than the single-photon spontaneous emission rate but that, with enhancement, 2PSE can reach 2.5% of the overall decay. When graphene nanoribbons are used, we also show that the emission of free-space photon pairs from Er 3+ at 3–3.2 μ⁢m via 2PSE can be increased by ~400. Our calculations show significantly less relative graphene-enhanced 2PSE than previous works, and we attribute this variation to differences in emitter size and assumed graphene mobility. Here, we also show that the internal energy structure of the ion can have an impact on the degree of 2PSE enhancement achievable and find that divalent rare earths are more favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Instability of flux flow and production of vortex-antivortex pairs by current-driven Josephson vortices in layered superconductors

We report numerical simulations of the nonlinear dynamics of Josephson vortices driven by strong dc currents in layered superconductors. Dynamic equations for interlayer phase differences in a stack of coupled superconducting layers were solved to calculate a drag coefficient η(J) of the vortex as a function of the perpendicular dc current density J. It is shown that Cherenkov radiation produced by a moving vortex causes significant radiation drag increasing η(v) at high vortex velocities v and striking instabilities of driven Josephson vortices moving faster than a terminal velocity v c . The steady-state flux flow breaks down at v > v c as the vortex starts producing a cascade of expanding vortex-antivortex pairs evolving into either planar macrovortex structures or branching flux patterns propagating both along and across the layers. This vortex-antivortex pair production triggered by a rapidly moving vortex is most pronounced in a stack of underdamped planar junctions where it can occur at J > J s well below the interlayer Josephson critical current density. Both v c and J s were calculated as functions of the quasiparticle damping parameter, and the dc magnetic field was applied parallel to the layers. The effects of vortex interaction on the Cherenkov instability of moving vortex chains and lattices in annular stacks of Josephson junctions were considered. It is shown that a vortex driven by a current density J > J s in a multilayer of finite length excites self-sustained large-amplitude standing waves of magnetic flux, resulting in temporal oscillations of the total magnetic moment. We evaluated a contribution of this effect to the power W radiated by the sample and showed that W increases strongly as the number of layers increases. Furthermore, these mechanisms can result in nonlinearity of the c-axis electromagnetic response and contribute to THz radiation from the layered cuprates at high dc current densities flowing perpendicular to the ab planes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Mechanism of structural colors in binary mixtures of nanoparticle-based supraballs

Inspired by structural colors in avian species, various synthetic strategies have been developed to produce noniridescent, saturated colors using nanoparticle assemblies. Nanoparticle mixtures varying in particle chemistry and size have additional emergent properties that affect the color produced. For complex multicomponent systems, understanding the assembled structure and a robust optical modeling tool can empower scientists to identify structure-color relationships and fabricate designer materials with tailored color. Here, we demonstrate how we can reconstruct the assembled structure from small-angle scattering measurements using the computational reverse-engineering analysis for scattering experiments method and use the reconstructed structure in finite-difference time-domain calculations to predict color. We successfully, quantitatively predict experimentally observed color in mixtures containing strongly absorbing nanoparticles and demonstrate the influence of a single layer of segregated nanoparticles on color produced. The versatile computational approach that we present is useful for engineering synthetic materials with desired colors without laborious trial-and-error experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Parametric study of optical transmission through plasmonic hole arrays modulated by the phase transition of vanadium dioxide

We have performed comprehensive electromagnetic simulations and preliminary experiments to explore the effects of geometrical and material parameters on the extraordinary optical transmission (EOT) through periodic arrays of subwavelength holes in a bilayer stack consisting of a gold or silver film atop a vanadium dioxide film (Au/Ag + VO 2 ), where the latter undergoes a semiconductor-to-metal phase transition. Using the finite-difference time-domain (FDTD) and finite-element methods (FEM), we vary iteratively the array periodicity, VO 2 film thickness and hole diameters, as well as the refractive index inside the VO 2 -layer holes and the VO 2 optical constants. For each variation, we compare the metallic-to-semiconducting ratios of the zero-order transmission ( T 00 ) peaks and find sharp maxima in these ratios within narrow parameter ranges. The maxima arise from Fabry-Perot and Fano-type resonances that minimize T 00 in the semiconducting phase of the perforated bilayers. At a fixed array period, the primary factors controlling the VO 2 -enabled EOT modulation are the VO 2 thickness, diameter of the VO 2 -layer holes, and absorption in the two VO 2 phases. Besides uncovering the origins of the higher metallic-phase T 00 , this study provides a protocol for optimizing the performance of the bilayer hole arrays for potential uses as dynamically tunable nano-optical devices.

Donev, Eugenii U. (ORCID:000000030555219X)↗
DOE OSTI2020

ATDM AD projects to ECP Report (Q3 FY2020)

The SPARC team has completed much of the workflow development needed for the coupled and uncoupled analyses needed to support our flight test validation work. The SPARC team has continued to make progress on development activities to support unsteady, full reentry vehicle aero analysis, which have focused on turbulence modeling, uniform mesh refinement, in-situ visualization (with the Paraview/Catalyst team), and performance analysis. The SPARC team is making solid progress towards our Q4 goal of demonstrating an improved Wall Modeled Large Eddy Simulation (WMLES) capability for standard geometries using high-order finite difference, discontinuous Galerkin, and low-dissipation finite volume methods. The SPARC team has met the minimum completion criteria for parts of the L1, including documenting performance and scaling of SPARC on Trinity, Sierra and Astra, and performing runs for the flight test validation case.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Ultrafast Laser Material Damage Simulation—A New Look at an Old Problem

The chirped pulse amplification technique has enabled the generation of pulses of a few femtosecond duration with peak powers multi-Tera and Peta–Watt in the near infrared. Its implementation to realize even shorter pulse duration, higher energy, and higher repetition rate laser systems relies on overcoming the limitations imposed by laser damage of critical components. In particular, the laser damage of coatings in the amplifiers and in post-compression optics have become a bottleneck. The robustness of optical coatings is typically evaluated numerically through steady-state simulations of electric field enhancement in multilayer stacks. However, this approach cannot capture crucial characteristics of femtosecond laser induced damage (LID), as it only considers the geometry of the multilayer stack and the optical properties of the materials composing the stack. This approach neglects that in the interaction of an ultrashort pulse and the materials there is plasma generation and associated material modifications. Here, we present a numerical approach to estimate the LID threshold of dielectric multilayer coatings based on strong field electronic dynamics. In this dynamic scheme, the electric field propagation, photoionization, impact ionization, and electron heating are incorporated through a finite-difference time-domain algorithm. We applied our method to simulate the LID threshold of bulk fused silica, and of multilayer dielectric mirrors and gratings. The results are then compared with experimental measurements. The salient aspects of our model, such as the implementation of the Keldysh photoionization model, the impact ionization model, the electron collision model for ‘low’-temperature, dense plasma, and the LID threshold criterion for few-cycle pulses are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of representative ground motion level on seismic PSA with the boundary between overestimation and underestimation

One commonly used approach in seismic probabilistic safety assessment (PSA) is the discrete method. This method follows the standard PSA framework and can be applied to various models, such as multi-unit models, while reducing computational costs using standard software. However, due to the inability to subdivide intervals infinitely, the discrete method approximates with a finite number of subintervals. In practice, different numbers of subintervals are applied, and the representative ground motion level is selected based on expert judgment. When employing a smaller number of subintervals, it is important to take caution to prevent underestimation. This study analyzes the impact of the representative ground motion level on seismic risk. It confirms that underestimation can occur with a small number of subintervals depending on the representative ground motion level. This study also proposes a method for determining the boundary of underestimation and overestimation. The method is demonstrated through examples, providing a mathematical foundation for selecting appropriate representative ground motion levels. By avoiding underestimation, this research helps prevent the oversight of significant risk contributors and enhances the understanding of seismic risk.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2020

Nonoverlapping Localized Exponential Time Differencing Methods for Diffusion Problems

In this paper, we propose nonoverlapping localized exponential time differencing (ETD) methods for diffusion problems. The model time-dependent diffusion equation is first reformulated on subdomains based on the nonoverlapping domain decomposition, in which Neumann boundary conditions are imposed on the interfaces for the subdomain problems and Dirichlet type conditions are enforced to form a space-time interface problem. After spatial discretization by standard central finite differences and temporal integration with the first or second order ETD methods, the fully discrete interface problem is obtained. Such an interface problem is then solved iteratively either at each time step or over the whole time interval: the former involves the solution of stationary problems in each subdomain at each iteration while the latter involves the solution of time-dependent subdomain problems at each iteration. For both approaches, we prove that localized ETD solutions conserve mass exactly and converge in time to the exact space semidiscrete solution. Numerical experiments in two dimensions are also presented to illustrate the performance of the proposed methods.

Mass conservation↗
DOE OSTI2025

Boundary Corrections for Kernel Approximation to Differential Operators

The kernel-based approach to operator approximation for partial differential equations has been shown to be unconditionally stable for linear PDEs and numerically exhibit unconditional stability for non-linear PDEs. These methods have the same computational cost as an explicit finite difference scheme but can exhibit order reduction at boundaries. In previous work on periodic domains, order reduction was addressed, yielding high-order accuracy. The issue addressed in this work is the elimination of order reduction of the kernel-based approach for a more general set of boundary conditions. Further, we consider the case of both first and second order operators. To demonstrate the theory, we provide not only the mathematical proofs but also experimental results by applying various boundary conditions to different types of equations. The results agree with the theory, demonstrating a systematic path to high order for kernel-based methods on bounded domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Efficient three-dimensional Laplace-Fourier domain acoustic-wave simulations using discontinuous finite-difference meshes with embedded boundaries

© 2018 SEG. We develop embedded boundary and discontinuous mesh methods to handle arbirarily shaped topography and accurately simulate acoustic seismic wave propagation in Laplace-Fourier domain. The purpose of the embedded boundary method is to enhance accurate wave simulation near the surface and the discontinuous mesh method is used to achieve considerable savings in both computation time and memory savings relative to fixed mesh schemes.

AlSalem, Hussain↗
DOE OSTI2018

Effect of laser intensity and exposure time on photothermal therapy with nanoparticles heated by a 793-nm diode laser and tissue optical clearing

Laser-induced thermotherapy is a promising method for cancer treatment, the outcome of which is affected by the exposure time. An inappropriate exposure time and laser intensity cause incomplete tumour destruction, tumour regrowth, and metastasis. Also possible is irreversible damage, i.e. death of healthy cells, and so numerical models are necessary to provide an optimised laser intensity and exposure time for different cancerous tumours. In this study, a model based on finite element method (FEM) is used for solving the bio-heat transfer equation and the Arrhenius equation describing tissue damage. The cancerous tumour is considered as a perfect cylinder with a diameter of 20 mm and a thickness of 2, 3, 4, and 5 mm, filled up by highly absorbing nanoparticles and surrounded by healthy cylindrical tissue with a diameter of 40 mm and a length of 10 mm, which ahs a low scattering coefficient due to optical clearing. The results show that 243 s is a safe and appropriate exposure time when a diode laser with a wavelength of 793 nm and intensity of 0.75 W cm{sup −2} together with gold nanorods of concentration 0.0001 % is used for the treatment of a 3-mm-thick tumour. Then, the results are developed and extensive numerical simulations are used to reveal mathematical relationships between two critical parameters, input power and optimised exposure time, for a series of different tumour thicknesses. Treatment protocols are presented.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

THERMAL/STRUCTURAL ANALYSIS OF THE AXION QUANTUM METROLOGY CAVITY AND ITS COMPONENTS

This research was centered around maximizing the capability to cool dielectric material within a containment unit, or Photonic Band Gap (PBG) cavity, designed for detecting axion dark matter and identifying the unit’s thermal properties. There are multiple types of PBG cavities, but the latest version that axion researchers wish to use has been theorized to contain possible issues related to its thermal properties. Thermal conductivity is an issue with the dielectric material because it is made from alumina which is highly insulative. This is important since the research is being done in a cryogenic environment and the thermal noise affects the quantum bit used for detecting the axion to photon conversion process. Therefore, any improvements to this unit should be justified and implemented but are not entirely limited to thermal contact related aspects of PBG cavities. A prospect of using sapphire in place of alumina also exists, but this is a more expensive, less tested, and more elusive material to justify the creation of a full dielectric structure out of sapphire. Thermal aspects of the cavity were analyzed using finite element method (FEM) and an experiment designed to test different thermal joint materials. FEM was used to check the contraction of the cavity during cooling, the contact quality between the dielectric material and cavity wall, theheat flow rate through the assembly, and helps visualize the cavity’s reaction to different design changes. A simple comparison between thermal conductivity curves justified the usage of sapphire over alumina from a thermal property point of view. Thus, the method of experimentation had an additional dielectric material, sapphire, to test with and compare to alumina. The thermal test identified which material is best to use for a thermal joint but also simultaneously found the conductance of the joint and dielectric material as well as an estimate for what the temperature is inside the larger scale cavity since a thermometer cannot be placed inside the PBG when it cools. Initially, the cavity simulation was tested for structural deformation properties and stress distribution since it is made of copper which contracts heavily in comparison to most other materials. The reduction in volume of the cavity was less than anticipated which gave more room for a possible usage of a modified thermal joint for additional contact area and higher conductivity.

43 PARTICLE ACCELERATORS↗