Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Implicit-explicit Runge-Kutta for radiation hydrodynamics I: Gray diffusion

Radiation hydrodynamics are a challenging multiscale and multiphysics set of equations. To capture the relevant physics of interest, one typically must time step on the hydrodynamics timescale, making explicit integration the obvious choice. On the other hand, the coupled radiation equations have a scaling such that implicit integration is effectively necessary in non-relativistic regimes. A first-order Lie-Trotter-like operator split is the most common time integration scheme used in practice, alternating between an explicit hydrodynamics step and an implicit radiation solve and energy deposition step. However, such a scheme is limited to first-order accuracy, and nonlinear coupling between the radiation and hydrodynamics equations makes a more general additive partitioning of the equations non-trivial. Here, we develop a new formulation and partitioning of radiation hydrodynamics with gray diffusion that allows us to apply (linearly) implicit-explicit Runge-Kutta time integration schemes. In conclusion, we prove conservation of total energy in the new framework, and demonstrate 2nd-order convergence in time on multiple radiative shock problems, achieving error 3–5 orders of magnitude smaller than the first-order Lie-Trotter operator split at the hydrodynamic CFL, even when Lie-Trotter applies a 3rd-order TVD Runge-Kutta scheme to the hydrodynamics equations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Optimizing long-term monitoring of radiation air-dose rates after the Fukushima Daiichi Nuclear Power Plant

Radiation air dose rates near the Fukushima Daiichi Nuclear Power Plant (FDNPP) have been steadily decreasing over the past eight years since the release of radioactive elements in March 2011. Currently, the radiation monitoring program is expected to transition to long-term monitoring after most of the remediation activities are completed. The main long-term monitoring objectives are to (1) confirm the continuing reduction of contaminant and hazard levels, (2) provide assurance for the public, (3) accumulate the basic datasets for scientific knowledge and future preparation, and (4) detect changes or anomalies in contaminant mobility (if they occur), or any unexpected processes or events. In this work, we have developed a methodology for optimizing the monitoring locations of radiation air dose-rate monitoring. Our approach consists of three steps in order to determine monitoring locations in a systematic manner: (1) prioritizing the critical locations, such as schools or regulatory requirement locations, (2) diversifying locations that cover the key environmental controls that are known to influence contaminant mobility and distributions, and (3) capturing the heterogeneity of radiation air-dose rates across the domain. Therefore, for the second step, we use a Gaussian mixture model to identify the representative locations among multiple environmental variables, such as elevation and land-cover types. For the third step, we use a Gaussian process model to capture and estimate the heterogeneity of air-dose rates across the domain. Employing an integrated dose-rate map derived from Bayesian geostatistical methods as a reference map, we distribute the monitoring locations in such a way as to capture the heterogeneity of the reference map. Our results have shown that this approach allows us to select monitoring locations in a systematic manner such that the heterogeneity of air dose rates is captured by the minimal number of monitoring locations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Enhancing isotope mixing in U-10Mo downblend castings with electromagnetic stirring

Uranium alloyed with 10 wt percent molybdenum (U-10Mo) is the proposed fuel for use in the United States’ high-performance research reactors. The U-10Mo fuel is fabricated by a two-step casting process that downblends highly enriched uranium with depleted uranium (DU) and/or natural uranium (NU) and subsequently alloys the resulting high-assay low-enriched uranium (HALEU) with Mo. Currently, a two-step casting process is needed to meet the ingots’ U-235 enrichment homogeneity specifications. This work demonstrates the ability to provide more homogeneous U-235 distributions during downblend casting by using zirconia for crucibles rather than graphite. When graphite is used as a crucible, the electromagnetic field produced by the induction heater couples directly with the graphite. Zirconia is nonconducting in this study, and therefore does not couple with the field. The induction field couples directly with the metal, and causes electromagnetic stirring (EMS) in the molten metal pool. Eight downblend casting experiments were carried out in this work. Four were performed with zirconia crucibles and four with graphite crucibles. The U-235 enrichment was measured at nine discrete points in each casting using laser ablation multi-collector inductively coupled mass spectrometry (LA-MC-ICP-MS). The enrichment homogeneity was approximately eight times better in zirconia crucible castings than in graphite crucible castings as measured by the plate enrichment range and enrichment coefficient of variation. The results show a single casting step could be used to meet USHPRR enrichment specifications if electromagnetic stirring is present during casting. A related investigation was carried out to determine whether it is appropriate to measure enrichment on as-cast ingots or whether heat treated specimens should be used to accurately characterize enrichment. The results show that the as-cast microstructure plays a significant role in enrichment homogeneity. It is recommended that all enrichment measurements be done after the homogenization heat treatment typical of U-10Mo processing.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Emulation of neutron damage with proton irradiation and its effects on microstructure and microchemistry of Zircaloy-4

This work assesses the potential of proton irradiation to simulate the neutron damage to the matrix and laves phase Zr(Fe,Cr) 2 precipitates in Zircaloy-4. Isothermal proton irradiation has been performed on Zircaloy-4 samples at irradiation temperatures ranging from 250 to 350 °C. Two-step proton irradiation was also performed to enhance the amorphization of and iron loss from the laves phase Zr(Fe,Cr) 2 precipitates. The irradiated microstructures, including dislocation loops and rafts near SPPs, were observed in proton irradiated Zircaloy-4, which are consistent with neutron irradiated material at a similar damage level. The amount of irradiation-induced hardening after proton irradiation was similar to post neutron irradiated data. The significant amorphization of the SPPs and concurrent Fe redistribution observed on neutron irradiated materials can be effectively emulated using a two-step proton irradiation on Zircaloy-4. Hence, the neutron irradiation effect on Zircaloy-4 can be mostly captured using the two-step proton irradiation described in this study.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using dynamic mode decomposition to predict the dynamics of a two-time non-equilibrium Green’s function

Computing the numerical solution of the Kadanoff–Baym equations, a set of nonlinear integral differential equations satisfied by the two-time Green's functions derived from many-body perturbation theory for a quantum many-body system away from equilibrium, is a challenging task. Recently, we have successfully applied dynamic mode decomposition (DMD) to construct a data driven reduced order model that can be used to extrapolate the time-diagonal of a two-time Green's function from numerical solutions of the KBE within a small time window. In this paper, we extend the previous work and use DMD to predict off-diagonal elements of the two-time Green's function. We partition the two-time Green's function into a number of one-time functions along the diagonal and subdiagonals of the two-time window as well as in horizontal and vertical directions. We use DMD to construct separate reduced order models to predict the dynamics of these one-time functions in a two-step procedure. We extrapolate along diagonal and several subdiagonals within a subdiagonal band of a two-time window in the first step. In the second step, we use DMD to extrapolate the Green's function outside of the sub-diagonal band. In conclusion, we demonstrate the efficiency and accuracy of this approach by applying it to a two-band Hubbard model problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Proton transfer reactions in donor site mutants of ESR, a retinal protein from Exiguobacterium sibiricum

Light-driven proton transport by microbial retinal proteins such as archaeal bacteriorhodopsin involves carboxylic residues as internal proton donors to the catalytic center which is a retinal Schiff base (SB). The proton donor, Asp96 in bacteriorhodopsin, supplies a proton to the transiently deprotonated Schiff base during the photochemical cycle. Subsequent proton uptake resets the protonated state of the donor. Here this two step process became a distinctive signature of retinal based proton pumps. Similar steps are observed also in many natural variants of bacterial proteorhodopsins and xanthorhodopsins where glutamic acid residues serve as a proton donor. Recently, however, an exception to this rule was found. A retinal protein from Exiguobacterium sibiricum, ESR, contains a Lys residue in place of Asp or Glu, which facilitates proton transfer from the bulk to the SB. Lys96 can be functionally replaced with the more common donor residues, Asp or Glu. Proton transfer to the SB in the mutants containing these replacements (K96E and K96D/A47T) is much faster than in the proteins lacking the proton donor (K96A and similar mutants), and in the case of K96D/A47T, comparable with that in the wild type, indicating that carboxylic residues can replace Lys96 as proton donors in ESR. We show here that there are important differences in the functioning of these residues in ESR from the way Asp96 functions in bacteriorhodopsin. Reprotonation of the SB and proton uptake from the bulk occur almost simultaneously during the M to N transition (as in the wild type ESR at neutral pH), whereas in bacteriorhodopsin these two steps are well separated in time and occur during the M to N and N to O transitions, respectively.

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

Precipitation of gadolinium from magnetic resonance imaging contrast agents may be the Brass tacks of toxicity

The formation of gadolinium-rich nanoparticles in multiple tissues from intravenous magnetic resonance imaging contrast agents may be the initial step in rare earth metallosis. The mechanism of gadolinium-induced diseases is poorly understood, as is how these characteristic nanoparticles are formed. Gadolinium deposition has been observed with all magnetic resonance imaging contrast agent brands. Aside from endogenous metals and acidic conditions, little attention has been paid to the role of the biological milieu in the degradation of magnetic resonance imaging contrast agents into nanoparticles. Herein, we describe the decomposition of the commercial magnetic resonance imaging contrast agents Omniscan and Dotarem in the presence of oxalic acid, a well-known endogenous compound. Omniscan dechelated rapidly and preluded measurement by the means available, while Dotarem underwent a two-step decomposition process. The decomposition of both magnetic resonance imaging contrast agents by oxalic acid formed gadolinium oxalate (Gd 2 [C 2 O 4 ] 3 , Gd 2 Ox 3 ). Furthermore, both observed steps of the Dotarem reaction involved the associative addition of oxalic acid. Adding protein (bovine serum albumin) increased the rate of dechelation. Displacement reactions could occur at lysosomal pH. Through these studies, we have demonstrated that magnetic resonance imaging contrast agents can be dissociated by endogenous molecules, thus illustrating a metric by which gadolinium-based contrast agents (GBCAs) might be destabilized in vivo.

Gadodiamide↗
DOE OSTI2022

Preparation and characterization of isotopically pure Mo targets for nuclear science measurements

Isotopically pure molybdenum targets are essential in various nuclear science measurements. For studies of cross sections of capture reactions on molybdenum, uniform targets with areal densities of hundreds of g/cm 2 are desired. Mechanical rolling of enriched metal pieces or vacuum sputtering methods used in previous studies to prepare molybdenum targets are challenging and involve labor- and time-intensive processes. Therefore, reliable and straightforward methods to produce molybdenum targets are crucial. We have developed a new double-step method to produce targets meeting these needs. The first step consists of evaporating isotopically enriched MoO 3 onto gold backings. The second step involves the reduction of the oxide layer in a 5% H 2 -Ar environment to obtain molybdenum targets. X-ray fluorescence (XRF) spectroscopy, X-ray Diffraction (XRD), a scanning electron microscope/focused ion beam milling, and energy-dispersive X-ray spectroscopy (EDS) were used to characterize the targets. These methods show that the targets exhibit a uniform coverage of Mo over the gold backings. Furthermore, significant volume shrinkage occurring during the reduction of MoO 3 creates a porous Mo layer where the gold diffuses due to the high processing temperatures, resulting in mixed Mo–Au layers.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Persistent global greening over the last four decades using novel long-term vegetation index data with enhanced temporal consistency

Advanced Very High-Resolution Radiometer (AVHRR) satellite observations have provided the longest global daily records from 1980s, but the remaining temporal inconsistency in vegetation index datasets has hindered reliable assessment of vegetation greenness trends. To tackle this, we generated novel global long-term Normalized Difference Vegetation Index (NDVI) and Near-Infrared Reflectance of vegetation (NIRv) datasets derived from AVHRR and Moderate Resolution Imaging Spectroradiometer (MODIS). We addressed residual temporal inconsistency through three-step post processing including cross-sensor calibration among AVHRR sensors, orbital drifting correction for AVHRR sensors, and machine learning-based harmonization between AVHRR and MODIS. After applying each processing step, we confirmed the enhanced temporal consistency in terms of detrended anomaly, trend and interannual variability of NDVI and NIRv at calibration sites. Our refined NDVI and NIRv datasets showed a persistent global greening trend over the last four decades (NDVI: 0.0008 yr -1 ; NIRv: 0.0003 yr -1 ), contrasting with those without the three processing steps that showed rapid greening trends before 2000 (NDVI: 0.0017 yr -1 ; NIRv: 0.0008 yr -1 ) and weakened greening trends after 2000 (NDVI: 0.0004 yr -1 ; NIRv: 0.0001 yr -1 ). These findings highlight the importance of minimizing temporal inconsistency in long-term vegetation index datasets, which can support more reliable trend analysis in global vegetation response to climate changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Modeling evacuation demand during no-notice emergency events: Tour formation behavior

Disastrous events have been drastically increasing – both in frequency and destructive capacity – over the past few years. While advance-notice events have received a great deal of attention in the literature of disaster management, not much attention so far has been given to the no-notice events mainly because of the scarcity of data. As an attempt to address this critical gap, the current study proposes a disaggregate evacuation demand framework to understand evacuees’ travel behavior in case of no-notice emergency events. Here, the proposed framework comprises four main steps of evacuation decision, evacuation planning, tour formation, and activity schedule update. This article is dedicated to the introduction of the framework structure and elaboration on the tour formation step. In this step, we first estimate the total number of intermediate stops, travel time, and distance of the evacuation tours for those who decide to evacuate through a joint modeling structure and then, determine the type of each intermediate stop (if any). It is found that a broad range of factors including evacuees’ demographic profiles, built-environment attributes, and characteristics of the disastrous event plays a significant role in people’s evacuation behavior during no-notice emergency events. The findings of this study can assist responsible agencies in understanding evacuees’ complex behavior, and consequently, in devising effective strategies to alleviate economic damages and casualties resulted by such events.

99 GENERAL AND MISCELLANEOUS↗