Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical 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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

WE-AB-204-07: Spatiotemporal Distribution of the FDG PET Tracer in Solid Tumors: Contributions of Diffusion and Convection Mechanisms

Purpose: In this study, a mathematical model is utilized to simulate FDG distribution in tumor tissue. In contrast to conventional compartmental modeling, tracer distributions across space and time are directly linked together (i.e. moving beyond ordinary differential equations (ODEs) to utilizing partial differential equations (PDEs) coupling space and time). The diffusion and convection transport mechanisms are both incorporated to model tracer distribution. We aimed to investigate the contributions of these two mechanisms on FDG distribution for various tumor geometries obtained from PET/CT images. Methods: FDG transport was simulated via a spatiotemporal distribution model (SDM). The model is based on a 5K compartmental model. We model the fact that tracer concentration in the second compartment (extracellular space) is modulated via convection and diffusion. Data from n=45 patients with pancreatic tumors as imaged using clinical FDG PET/CT imaging were analyzed, and geometrical information from the tumors including size, shape, and aspect ratios were classified. Tumors with varying shapes and sizes were assessed in order to investigate the effects of convection and diffusion mechanisms on FDG transport. Numerical methods simulating interstitial flow and solute transport in tissue were utilized. Results: We have shown the convection mechanism to depend on the shape and size of tumors whereas diffusion mechanism is seen to exhibit low dependency on shape and size. Results show that concentration distribution of FDG is relatively similar for the considered tumors; and that the diffusion mechanism of FDG transport significantly dominates the convection mechanism. The Peclet number which shows the ratio of convection to diffusion rates was shown to be of the order of 10−{sup 3} for all considered tumors. Conclusion: We have demonstrated that even though convection leads to varying tracer distribution profiles depending on tumor shape and size, the domination of the diffusion phenomenon prevents these factors from modulating FDG distribution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Quantum simulations of hydrodynamics via the Madelung transformation

Developing numerical methods to simulate efficiently nonlinear fluid dynamics on universal quantum computers is a challenging problem. In this paper, a generalization of the Madelung transform is defined to solve quantum relativistic charged fluid equations interacting with external electromagnetic forces via the Dirac equation. The Dirac equation is discretized into discrete-time quantum walks which can be efficiently implemented on universal quantum computers. A variant of this algorithm is proposed to implement simulations using current noisy intermediate scale quantum (NISQ) devices in the case of homogeneous external forces. High resolution (up to N=2 17 grid points) numerical simulations of relativistic and nonrelativistic hydrodynamical shocks on current IBM NISQs are performed with this algorithm. Here, this paper demonstrates that fluid dynamics can be simulated on NISQs, and opens the door to simulating other fluids, including plasmas, with more general quantum walks and quantum automata.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

A TWO-DIMENSIONAL METHOD OF MANUFACTURED SOLUTIONS BENCHMARK SUITE BASED ON VARIATIONS OF LARSEN'S BENCHMARK WITH ESCALATING ORDER OF SMOOTHNESS OF THE EXACT SOLUTION

The quantification of the discretization error associated with the spatial discretization of the Discrete Ordinate(DO) equations in multidimensional Cartesian geometries is the central problem in error estimation of spatial discretization schemes for transport theory as well as computer code verification. Traditionally fine mesh solutions are employed as reference, because analytical solutions only exist in the absence of scattering. This approach, however, is inadequate when the discretization error associated with the reference solution is not small compared to the discretization error associated with the mesh under scrutiny. Typically this situation occurs if the mesh of interest is only a couple of refinement levels away from the reference solution or if the order of accuracy of the numerical method (and hence the reference as well) is lower than expected. In this work we present a Method of Manufactured Solutions (MMS) benchmark suite with variable order of smoothness of the underlying exact solution for two-dimensional Cartesian geometries which provides analytical solutions aver- aged over arbitrary orthogonal meshes for scattering and non-scattering media. It should be emphasized that the developed MMS benchmark suite first eliminates the aforementioned limitation of fine mesh reference solutions since it secures knowledge of the underlying true solution and second that it allows for an arbitrary order of smoothness of the underlying ex- act solution. The latter is of importance because even for smooth parameters and boundary conditions the DO equations can feature exact solution with limited smoothness. Moreover, the degree of smoothness is crucial for both the order of accuracy and the magnitude of the discretization error for any spatial discretization scheme.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

ThinCurr: An open-source 3D thin-wall eddy current modeling code for the analysis of large-scale systems of conducting structures

In this paper we present a new thin-wall eddy current modeling code, ThinCurr, for studying inductively-coupled currents in 3D conducting structures -- with primary application focused on the interaction between currents flowing in coils, plasma, and conducting structures of magnetically-confined plasma devices. The code utilizes a boundary finite element method on an unstructured, triangular grid to accurately capture device structures. The new code, part of the broader Open FUSION Toolkit, is open-source and designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. Scalability to large models is enabled through use of hierarchical off-diagonal low-rank compression of the inductance matrix, which is otherwise dense. Ease of handling large models of complicated geometry is further supported by automatic determination of supplemental elements through a greedy homology approach. Here, a detailed description of the numerical methods of the code and verification of the implementation of those methods using cross-code comparisons against the VALEN code and Ansys commercial analysis software is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Wave propagation speeds and source term influences in single and integral porosity shallow water equations

In urban flood modeling, so-called porosity shallow water equations (PSWEs), which conceptually account for unresolved structures, e.g., buildings, are a promising approach to addressing high CPU times associated with state-of-the-art explicit numerical methods. The PSWE can be formulated with a single porosity term, referred to as the single porosity shallow water model (SP model), which accounts for both the reduced storage in the cell and the reduced conveyance, or with two porosity terms: one accounting for the reduced storage in the cell and another accounting for the reduced conveyance. The latter form is referred to as an integral or anisotropic porosity shallow water model (AP model). The aim of this study was to analyze the differences in wave propagation speeds of the SP model and the AP model and the implications of numerical model results. First, augmented Roe-type solutions were used to assess the influence of the source terms appearing in both models. It is shown that different source terms have different influences on the stability of the models. Second, four computational test cases were presented and the numerical models were compared. It is observed in the eigenvalue-based analysis as well as in the computational test cases that the models converge if the conveyance porosity in the AP model is close to the storage porosity. If the porosity values differ significantly, the AP model yields different wave propagation speeds and numerical fluxes from those of the BP model. In this study, the ratio between the conveyance and storage porosities was determined to be the most significant parameter.

58 GEOSCIENCES↗
DOE OSTI2016

Time-optimal excitation of maximum quantum coherence: Physical limits and pulse sequences

Here we study the optimum efficiency of the excitation of maximum quantum (MaxQ) coherence using analytical and numerical methods based on optimal control theory. The theoretical limit of the achievable MaxQ amplitude and the minimum time to achieve this limit are explored for a set of model systems consisting of up to five coupled spins. In addition to arbitrary pulse shapes, two simple pulse sequence families of practical interest are considered in the optimizations. Compared to conventional approaches, substantial gains were found both in terms of the achieved MaxQ amplitude and in pulse sequence durations. For a model system, theoretically predicted gains of a factor of three compared to the conventional pulse sequence were experimentally demonstrated. Motivated by the numerical results, also two novel analytical transfer schemes were found: Compared to conventional approaches based on non-selective pulses and delays, double-quantum coherence in two-spin systems can be created twice as fast using isotropic mixing and hard spin-selective pulses. Also it is proved that in a chain of three weakly coupled spins with the same coupling constants, triple-quantum coherence can be created in a time-optimal fashion using so-called geodesic pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Modeling of the oxygen reduction reaction for dense LSM thin films

In this study, the oxygen reduction reaction mechanism is investigated using numerical methods on a dense thin (La 1-x Sr x ) y MnO 3±δ film deposited on a YSZ substrate. This 1-D continuum model consists of defect chemistry and elementary oxygen reduction reaction steps coupled via reaction rates. The defect chemistry model contains eight species including cation vacancies on the A- and B-sites. The oxygen vacancy is calculated by solving species transportation equations in multiphysics simulations. Due to the simple geometry of a dense thin film, the oxygen reduction reaction was reduced to three elementary steps: surface adsorption and dissociation, incorporation on the surface, and charge transfer across the LSM/YSZ interface. The numerical simulations allow for calculation of the temperature- and oxygen partial pressure-dependent properties of LSM. The parameters of the model are calibrated with experimental impedance data for various oxygen partial pressures at different temperatures. The results indicate that surface adsorption and dissociation is the rate-determining step in the ORR of LSM thin films. With the fine-tuned parameters, further quantitative analysis is performed. The activation energy of the oxygen exchange reaction and the dependence of oxygen non-stoichiometry on oxygen partial pressure are also calculated and verified using the literature results.

36 MATERIALS SCIENCE↗
DOE OSTI2010

Scalable direct Vlasov solver with discontinuous Galerkin method on unstructured mesh.

This paper presents the development of parallel direct Vlasov solvers with discontinuous Galerkin (DG) method for beam and plasma simulations in four dimensions. Both physical and velocity spaces are in two dimesions (2P2V) with unstructured mesh. Contrary to the standard particle-in-cell (PIC) approach for kinetic space plasma simulations, i.e., solving Vlasov-Maxwell equations, direct method has been used in this paper. There are several benefits to solving a Vlasov equation directly, such as avoiding noise associated with a finite number of particles and the capability to capture fine structure in the plasma. The most challanging part of a direct Vlasov solver comes from higher dimensions, as the computational cost increases as N{sup 2d}, where d is the dimension of the physical space. Recently, due to the fast development of supercomputers, the possibility has become more realistic. Many efforts have been made to solve Vlasov equations in low dimensions before; now more interest has focused on higher dimensions. Different numerical methods have been tried so far, such as the finite difference method, Fourier Spectral method, finite volume method, and spectral element method. This paper is based on our previous efforts to use the DG method. The DG method has been proven to be very successful in solving Maxwell equations, and this paper is our first effort in applying the DG method to Vlasov equations. DG has shown several advantages, such as local mass matrix, strong stability, and easy parallelization. These are particularly suitable for Vlasov equations. Domain decomposition in high dimensions has been used for parallelization; these include a highly scalable parallel two-dimensional Poisson solver. Benchmark results have been shown and simulation results will be reported.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2014

Thermally-activated non-Schmid glide of screw dislocations in W using atomistically-informed kinetic Monte Carlo simulations

Thermally-activated 1/2<111> screw dislocation motion is the controlling plastic mechanism at low temperatures in body-centered cubic (bcc) crystals. Dislocation motion proceeds by nucleation and propagation of atomic-sized kink pairs in close-packed planes. The atomistic character of kink pairs can be studied using techniques such as molecular dynamics (MD). However, MD’s natural inability to properly sample thermally-activated processes as well as to capture {110} screw dislocation glide calls for the development of other methods capable of overcoming these limitations. Here we develop a kinetic Monte Carlo (kMC) approach to study single screw dislocation dynamics from room temperature to 0.5T m and at stresses 0 < σ < 0.9σ P , where Tm and σ P are the melting point and the Peierls stress. The method is entirely parameterized with atomistic simulations using an embedded atom potential for tungsten. To increase the physical fidelity of our simulations, we calculate the deviations from Schmid’s law prescribed by the interatomic potential used and we study single dislocation kinetics using both projections. We calculate dislocation velocities as a function of stress, temperature, and dislocation line length. We find that considering non-Schmid effects has a strong influence on both the magnitude of the velocities and the trajectories followed by the dislocation. We finish by condensing all the calculated data into effective stress and temperature dependent mobilities to be used in more homogenized numerical methods.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Circumventing mesh bias by r- and h-adaptive techniques for variational eigenfracture

This work introduces and compares mesh r- and h-adaptivity for the eigenfracture model originally proposed in Schmidt et al. (Multiscale Model Simul 7:1237–1266, 2009), Pandolfi and Ortiz (J Numer Methods Eng 92:694–714, 2012), with the goal of suppressing potential mesh bias due to the element deletion. In the r-adaptive approach, we compute the configurational force at each incremental step and move nodes near the crack tip parallel to the normalized configurational forces field such that the crack propagation direction can be captured more accurately within each incremental step. In the h-adapative approach, we introduce mesh refinement via a quad-tree algorithm to introduce more degrees of freedom within the nonlocal ϵ-neighborhood such that a more refined crack path can be reproduced with a higher mesh resolution. Our numerical examples indicate that the r-adaptive approach is able to replicate curved cracks and complex geometrical features, whereas the h-adaptive approach is advantageous in simulating sub-scale fracture when the nonlocal regions are smaller than the un-refined coarse mesh.

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

Simulation of sharp interface multi-material flows involving an arbitrary number of components through an extended five-equation model

In this paper, we present an anti-diffusive method dedicated to the simulation of interface flows on Cartesian grids involving an arbitrary number m of compressible components. Our work is two-fold: first, we introduce a m-component flow model that generalizes a classic two material five-equation model. In that way, interfaces are localized using color function discontinuities and a pressure equilibrium closure law is used to complete this new model. The resulting model is demonstrated to be hyperbolic under simple assumptions and consistent. Second, we present a discretization strategy for this model relying on a Lagrange–Remap scheme. Here, the projection step involves an anti-dissipative mechanism allowing to prevent numerical diffusion of the material interfaces. The proposed solver is built ensuring consistency and stability properties but also that the sum of the color functions remains equal to one. The resulting scheme is first order accurate and conservative for the mass, momentum, energy and partial masses. Furthermore, the obtained discretization preserves Riemann invariants like pressure and velocity at the interfaces. Finally, validation computations of this numerical method are performed on several tests in one and two dimensions. The accuracy of the method is also compared to results obtained with the upwind Lagrange–Remap scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

The global version of the gyrokinetic turbulence code GENE

The understanding and prediction of transport due to plasma microturbulence is a key open problem in modern plasma physics, and a grand challenge for fusion energy research. Ab initio simulations of such small-scale, low-frequency turbulence are to be based on the gyrokinetic equations, a set of nonlinear integro-differential equations in reduced (five-dimensional) phase space. In the present paper, the extension of the well-established and widely used gyrokinetic code GENE [F. Jenko, W. Dorland, M. Kotschenreuther, B.N. Rogers, Electron temperature gradient driven turbulence, Phys. Plasmas 7 (2000) 1904-1910] from a radially local to a radially global (nonlocal) version is described. The necessary modifications of both the basic equations and the employed numerical methods are detailed, including, e.g., the change from spectral methods to finite difference and interpolation techniques in the radial direction and the implementation of sources and sinks. In addition, code verification studies and benchmarks are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Bayesian inference for high confidence signal validation and sensor calibration assessment - 328

Online Monitoring (OLM) systems has been investigated for applications in the nuclear industry for the past two decades. These systems offer economic and technological advantages over the currently implemented time based testing methods for sensors and associated instrumentations. The advantages offered not only augment economic viability of current nuclear reactors but will be an inevitable change called for to suit the technologies of the advanced reactors (Gen III+). Irrespective, the acceptance of OLM into the regulatory structure in the U.S. is still hindered by technical and structural clarity on uncertainty quantification used by these systems. A Bayesian inference technique is implemented in the form of a residual model to quantify uncertainty in OLM predictions. The developed system used Gaussian processes to represent the OLM predictions and model inadequacy, the stationary representation of model based uncertainty. In addition, model selection techniques are used in conjunction with numerical methods such as Genetic Algorithm (GA) and Gibb's sampler to add robustness to the estimation model inadequacy, particularly in the presence of drifts and faulted conditions. The prediction model developed herein find applicability in fault detection, obtaining real values of processes after the onset of faults etc. Further development looks into the extension of current method for a multiple output scale implementation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2013

A diffuse interface model with immiscibility preservation

A new, simple, and computationally efficient interface capturing scheme based on a diffuse interface approach is presented for simulation of compressible multiphase flows. Multi-fluid interfaces are represented using field variables (interface functions) with associated transport equations that are augmented, with respect to an established formulation, to enforce a selected interface thickness. The resulting interface region can be set just thick enough to be resolved by the underlying mesh and numerical method, yet thin enough to provide an efficient model for dynamics of well-resolved scales. A key advance in the present method is that the interface regularization is asymptotically compatible with the thermodynamic mixture laws of the mixture model upon which it is constructed. It incorporates first-order pressure and velocity non-equilibrium effects while preserving interface conditions for equilibrium flows, even within the thin diffused mixture region. We first quantify the improved convergence of this formulation in some widely used one-dimensional configurations, then show that it enables fundamentally better simulations of bubble dynamics. Demonstrations include both a spherical-bubble collapse, which is shown to maintain excellent symmetry despite the Cartesian mesh, and a jetting bubble collapse adjacent a wall. Comparisons show that without the new formulation the jet is suppressed by numerical diffusion leading to qualitatively incorrect results.

42 ENGINEERING↗
DOE OSTI2013

Simulation of the conditions of wave excitation for the optimization of the lower hybrid current drive in the Globus-M spherical tokamak

At present, the method of current drive by means of lower hybrid waves is not applied to low-aspect-ratio tokamaks, because, in the traditional approach, it would be necessary to use waves with a very high slowing-down factor. However, studies of new transparency regions for waves in a nonuniform magnetized plasma, performed earlier at the Ioffe Physical Technical Institute, Russian Academy of Sciences, made it possible to develop an approach in which slow waves are excited in the poloidal (rather than toroidal) direction. In this approach, moderately slowed-down waves first propagate in the poloidal direction, but then turn in the toroidalal direction and get into the dense plasma. In this work, this approach is further developed using numerical methods. In particular, the influence of the density profile in the edge plasma on the efficiency of wave excitation for given antenna parameters is studied in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Numerical Polynomial Homotopy Continuation Method and String Vacua

Finding vacua for the four-dimensional effective theories for supergravity which descend from flux compactifications and analyzing them according to their stability is one of the central problems in string phenomenology. Except for some simple toy models, it is, however, difficult to find all the vacua analytically. Recently developed algorithmic methods based on symbolic computer algebra can be of great help in the more realistic models. However, they suffer from serious algorithmic complexities and are limited to small system sizes. In this paper, we review a numerical method called the numerical polynomial homotopy continuation (NPHC) method, first used in the areas of lattice field theories, which by construction finds all of the vacua of a given potential that is known to have only isolated solutions. The NPHC method is known to suffer from no major algorithmic complexities and is embarrassingly parallelizable , and hence its applicability goes way beyond the existing symbolic methods. We first solve a simple toy model as a warm-up example to demonstrate the NPHC method at work. We then show that all the vacua of a more complicated model of a compactified M theory model, which has an S U ( 3 ) structure, can be obtained by using a desktop machine in just about an hour, a feat which was reported to be prohibitively difficult by the existing symbolic methods. Finally, we compare the various technicalities between the two methods.

Mehta, Dhagash↗