Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Computational and experimental study of laminar flames

During the past three years, our research has centered on an investigation of the effects of complex chemistry and detailed transport on the structure and extinction of hydrocarbon flames in coflowing axisymmetric configurations. We have pursued both computational and experimental aspects of the research in parallel on both steady-state and time-dependent systems. The computational work has focused on the application of accurate and efficient numerical methods for the solution of the steady-state and time-dependent boundary value problems describing the various reacting systems. Detailed experimental measurements were performed on axisymmetric coflow flames using two-dimensional imaging techniques. Previously, spontaneous Raman scattering, chemiluminescence, and laser-induced fluorescence were used to measure the temperature, major and minor species profiles. Particle image velocimetry (PIV) has been used to investigate velocity distributions and for calibration of time-varying flames. Laser-induced incandescence (LII) with an extinction calibration was used to determine soot volume fractions, while soot surface temperatures were measured with three-color optical pyrometry using a color digital camera. A blackbody calibration of the camera allows for determination of soot volume fraction as well, which can be compared with the LII measurements. More recently, we have concentrated on a detailed characterization of soot using a variety of techniques including time-resolved LII (TiRe-LII) for soot primary particles sizes, multi-angle light scattering (MALS) for soot radius of gyration, and spectrally-resolved line of sight attenuation (spec-LOSA). Combining the information from all of these soot measurements can be used to determine the soot optical properties, which are observed to vary significantly depending on spatial location and fuel dilution. Our goal has been to obtain a more fundamental understanding of the important fluid dynamic and chemical interactions in these flames so that this information can be used effectively in combustion modeling.

02 PETROLEUM↗
DOE OSTI2024

Open‐Source Anaerobic Digestion Modeling Platform, Anaerobic Digestion Model No. 1 Fast (ADM1F)

An open‐source modeling platform, called Anaerobic Digestion Model No. 1 Fast (ADM1F), is introduced to achieve fast and numerically stable simulations of anaerobic digestion processes. ADM1F is compatible with an iPython interface to facilitate model configuration, simulation, data analysis, and visualization. Faster simulations and more stable results are accomplished by implementing an advanced open‐source library of numerical methods called Portable Extensive Toolkit for Scientific Computation (PETSc) to solve the ADM1 system of equations. Leveraging PETSc, ADM1F can consistently complete a steady‐state simulation under 0.2 s, over 99% faster than a benchmark ADM1 model implemented with MATLAB while achieving agreement of model outputs within 1% of those obtained with the benchmark model. For dynamic simulations, however, ADM1F has a computational speed advantage only when the influent characteristics update more frequently than every 4 h. The ability of ADM1F to be useful as a tool to study anaerobic digestion systems is demonstrated through two example implementations of ADM1F: (1) a two‐phase co‐digestion scenario evaluating the impact of the organic loading rate and the substrate composition on reactor performance and stability, and (2) a conventional digester scenario assessing the effectiveness of recovery strategies after disruptions that led to instability. These examples demonstrate how the high simulation speed and the convenience of the iPython interface allow ADM1F to complete complex analyses within minutes, much faster than computational strategies currently reported in the literature.

anaerobic co-digestion↗
DOE OSTI2017

Comparison of Spatial Discretization Methods for Solving the S N Equations Using a Three-Dimensional Method of Manufactured Solutions Benchmark Suite with Escalating Order of Nonsmoothness

For this study, a comparison of the accuracy and computational efficiency of spatial discretization methods of the three-dimensional S N equations is conducted, including discontinuous Galerkin finite element methods, the arbitrarily high-order transport method of nodal type (AHOTN), the linear-linear method, the linear-nodal (LN) method, and the higher-order diamond difference method. For this purpose, we have developed a suite of method of manufactured solutions benchmarks that provides an exact solution of the S N equations even in the presence of scattering. Most importantly, our benchmark suite permits the user to set an arbitrary level of smoothness of the exact solution across the singular characteristics. Our study focuses on the computational efficiency of the considered spatial discretization methods. Numerical results indicate that the best-performing method depends on the norm used to measure the discretization error. We employ discrete Lp norms and integral error norms in this work. For configurations with continuous exact angular flux, high-order AHOTNs perform best under L p error norms, while the LN method performs best when measured by integral error norms. Lastly, when the angular flux is discontinuous, a new singular-characteristic tracking method for three-dimensional geometries performs best among the considered methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Review of recent activities with MOOSE, an open-source finite element & finite volume multi-fidelity simulation framework

Modeling and simulation are an increasing part of engineering. This is undoubtedly driven by the high costs of constructing experimental facilities, but also enabled by the exponential increase in computing powers over the last decades, which allows computational models to be closer than ever to reality. One of the main drivers for the development of MOOSE is supporting advanced nuclear reactor simulations. A challenging aspect of modeling advanced nuclear reactors is the plurality of physics involved, including neutronics, thermal hydraulics and fuel performance. These physics are all coupled to some extent and are generally solved in a sequential but iterative fashion. The United States (U.S.) national laboratories have been developing MOOSE, an open source multiphysics framework since its inception at the Idaho National Laboratory (INL) in 2008. This framework enables seamless coupling of multiphysics simulations and facilitates the implementation of new physics and material governing laws. It is continuously expanded with novel numerical methods and new pre-implemented physics module. Numerous applications, developed within the Department of Energy (DOE) laboratories, academia, and industry, including outside of nuclear engineering, have been developed to study specialized physics problems. International collaborations are welcome on this open-source modeling and simulation project.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Design and convergence analysis of the conforming virtual element method for polyharmonic problems

In the recent years, there has been a tremendous interest to numerical methods that approximate partial differential equations (PDEs) on computational meshes with arbitrarily-shaped polygonal/polyhedral (polytopal, for short) elements. A nonexhaustive list of such methods include the Mimetic Finite Difference method [4, 9, 21, 24, 28, 33, 35, 62, 64–67], the Polygonal Finite Element Method [68, 73, 74], the polygonal Discontinuous Galerkin Finite Element Methods [5, 7, 8, 10, 13–15, 18, 36, 38, 39, 41, 50], the Hybridizable Discontinuous Galerkin and Hybrid High-Order Methods [46–48, 51], the Gradient Discretization method [49, 53, 54], the Finite Volume Method [52], and the BEM-based FEM [71]. An alternative approach that also proved to be very successful is provided by the Virtual Element method (VEM), which was originally proposed in [19] for the numerical treatment of second-order elliptic problems [43, 44], and readily extended to linear and nonlinear elasticity [22, 25, 56], plate bending problems [37], Cahn-Hilliard equation [3], Stokes equations [2], Laplace-Beltrami equation [55], Darcy-Brinkam equation [76], discrete topology optimization problems [6], fracture networks problems [29], eigenvalue problems [59, 75]. The mixed virtual element formulation was proposed in [20, 34]; the nonconforming Virtual element formulations was proposed for second-order elliptic problems in [16], and later extended to general advection-reaction-diffusion problems, Stokes equation, the biharmonic problems, the iegenvalue problem, and the Schrodinger equation in [11, 30, 42, 45, 58, 80]. The p- and hp-version of the VEM were developed in [26, 27] and efficient multigrid methods for the resulting linear system of equations in [12]. A posteriori error estimates can be found in [40].

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

On the regularization method for Fredholm integral equations with odd weakly singular kernel

In this paper, we propose a numerical method to approach the solution of a Fredholm integral equation with a weakly singular kernel by applying the convolution product as a regularization operator and the Fourier series as a projection. Preliminary numerical results show that the order of convergence of the method is better than the one of conventional projection methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A Posteriori Error Analysis of Two-Stage Computation Methods with Application to Efficient Discretization and the Parareal Algorithm

We consider numerical methods for initial value problems that employ a two stage approach consisting of solution on a relatively coarse discretization followed by solution on a relatively fine discretization. Examples include adaptive error control, parallel-in-time solution schemes, and efficient solution of adjoint problems for computing a posteriori error estimates. We describe a general formulation of two stage computations then perform a general a posteriori error analysis based on computable residuals and solution of an adjoint problem. The analysis accommodates various variations in the two stage computation and in formulation of the adjoint problems. We apply the analysis to compute “dual-weighted” a posteriori error estimates, to develop novel algorithms for efficient solution that take into account cancellation of error, and to the Parareal Algorithm. We test the various results using several numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Methods to measure contact angles in scCO 2 -brine-sandstone systems

Numerous methods are used to measure contact angles ( θ ) in multiphase systems. The wettability and θ are primary controls on CO 2 residual trapping during Geologic Carbon Storage (GCS) and determining these values within rock pores is paramount to increasing storage efficiency. One traditional experimental approach is the sessile drop method which involves θ measurements on a single image of droplets. More recent developments utilize X-ray micro-computed tomography (CT) scans which provide the resolutions necessary to image in situ θ of fluids at representative conditions; however, experimental micro-CT data is limited and varied. To further examine θ distributions in supercritical-CO 2 -brine-sandstone systems, a combination of manual and automated θ measurement methods were utilized to measure θ using both sessile drop and micro-CT images of two sandstone cores. Here, the purpose of this work was threefold: (1) compare two current and two new θ measuring methods using micro-CT images of scCO 2 -brine-sandstone systems; (2) determine how traditional experimental method (sessile drop) θ results compare to in situ θ results (micro-CT); and (3) determine if the Matlab Contact Angle Algorithm (MCAA) from Klise et al. (2016) can be used to measure θ scCO 2 -brine-sandstone systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A high-order multiscale finite-element method for time-domain elastic wave modeling in strongly heterogeneous media

Efficient and accurate numerical methods for elastic wave modeling in complex media have many important applications. However, it is fairly challenging to model elastic wave propagation in strongly heterogeneous media with high computational efficiency and high-order accuracy simultaneously. We develop a novel high-order multiscale finite-element method to model elastic wave propagation in strongly heterogeneous media in the time domain. The most important feature of our method is a generalization of standard multiscale finite element method by using high-order multiscale finite-element basis functions to capture the fine-scale heterogeneities on the coarse mesh, in contrast to conventional finite-element basis functions that are merely determined by the order of polynomials. These multiscale basis functions leads to a system matrix with significantly reduced dimension, thus enable us to solve the elastic wave equation on the coarse mesh with high-order accuracy and very low computational time cost. In conclusion, we use 2D and 3D numerical examples to demonstrate the superior efficiency and accuracy of our new modeling method compared with the conventional spectral-element method.

58 GEOSCIENCES↗
DOE OSTI2014

A harmonic polynomial cell (HPC) method for 3D Laplace equation with application in marine hydrodynamics

We propose a new efficient and accurate numerical method based on harmonic polynomials to solve boundary value problems governed by 3D Laplace equation. The computational domain is discretized by overlapping cells. Within each cell, the velocity potential is represented by the linear superposition of a complete set of harmonic polynomials, which are the elementary solutions of Laplace equation. By its definition, the method is named as Harmonic Polynomial Cell (HPC) method. The characteristics of the accuracy and efficiency of the HPC method are demonstrated by studying analytical cases. Comparisons will be made with some other existing boundary element based methods, e.g. Quadratic Boundary Element Method (QBEM) and the Fast Multipole Accelerated QBEM (FMA-QBEM) and a fourth order Finite Difference Method (FDM). To demonstrate the applications of the method, it is applied to some studies relevant for marine hydrodynamics. Sloshing in 3D rectangular tanks, a fully-nonlinear numerical wave tank, fully-nonlinear wave focusing on a semi-circular shoal, and the nonlinear wave diffraction of a bottom-mounted cylinder in regular waves are studied. The comparisons with the experimental results and other numerical results are all in satisfactory agreement, indicating that the present HPC method is a promising method in solving potential-flow problems. The underlying procedure of the HPC method could also be useful in other fields than marine hydrodynamics involved with solving Laplace equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Shock waves simulated using the dual domain material point method combined with molecular dynamics

Here in this work we combine the dual domain material point method with molecular dynamics in an attempt to create a multiscale numerical method to simulate materials undergoing large deformations with high strain rates. In these types of problems, the material is often in a thermodynamically nonequilibrium state, and conventional constitutive relations or equations of state are often not available. In this method, the closure quantities, such as stress, at each material point are calculated from a molecular dynamics simulation of a group of atoms surrounding the material point. Rather than restricting the multiscale simulation in a small spatial region, such as phase interfaces, or crack tips, this multiscale method can be used to consider nonequilibrium thermodynamic effects in a macroscopic domain. This method takes the advantage that the material points only communicate with mesh nodes, not among themselves; therefore molecular dynamics simulations for material points can be performed independently in parallel. The dual domain material point method is chosen for this multiscale method because it can be used in history dependent problems with large deformation without generating numerical noise as material points move across cells, and also because of its convergence and conservation properties. In conclusion, to demonstrate the feasibility and accuracy of this method, we compare the results of a shock wave propagation in a cerium crystal calculated using the direct molecular dynamics simulation with the results from this combined multiscale calculation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Code-verification techniques for the method-of-moments implementation of the combined-field integral equation

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. Here, we demonstrate the effectiveness of these approaches for cases with and without coding errors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Integrated simulation of vertical fracture propagation induced by water injection and its borehole electromagnetic responses in shale gas systems

We analyze fracture propagation induced by hydraulic fracturing with water injection and examine their detectability with crosswell electromagnetic (EM) geophysical methods. For rigorous 3D coupled flow-geomechanical modeling, we employ a numerical method that can model failure by tensile and shear stresses, dynamic nonlinear permeability, dual continuum approach, and thermo-poro-mechanical effects. From numerical simulation, we find that the fracture propagation is not the same as propagation of the water front, because fracturing is governed by geomechanics whereas water saturation is determined by multiphase flow. At early times, the water front is almost identical to the fracture tip, suggesting that the fracture is mostly filled with the injected water. However, at late times, movement of the water front is retarded compared to fracture propagation, yielding a significant gap between the water front and the fracture top, filled with reservoir gas. During fracture propagation, the coupled flow-geomechanical models are transformed via a rock-physics model into electrical conductivity models. We employ a full 3D finite-element EM geophysical simulator to evaluate the sensitivity of the crosswell EM method to fracture propagation. It is shown that anomalous distribution of electrical conductivity is closely related to the injected water saturation, but not closely related to newly created unsaturated fractures. Our numerical modeling experiments demonstrate that the crosswell EM method can be highly sensitive to electrical conductivity changes that directly indicate the migration pathways of the injected fluid. Accordingly, the EM method can serve as an effective monitoring tool for monitoring the injected fluids during hydraulic fracturing operations.

58 GEOSCIENCES↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2015

Relativistic opacities for astrophysical applications

Here, we report on the use of the Los Alamos suite of relativistic atomic physics codes to generate radiative opacities for the modeling of astrophysically relevant plasmas under local thermodynamic equilibrium (LTE) conditions. The atomic structure calculations are carried out in fine-structure detail, including full configuration interaction. Three example applications are considered: iron opacities at conditions relevant to the base of the solar convection zone, nickel opacities for the modeling of stellar envelopes, and samarium opacities for the modeling of light curves produced by neutron star mergers. In the first two examples, comparisons are made between opacities that are generated with the fully and semi-relativistic capabilities in the Los Alamos suite of codes. As expected for these highly charged, iron-peak ions, the two methods produce reasonably similar results, providing confidence that the numerical methods have been correctly implemented. However, discrepancies greater than 10% are observed for nickel and investigated in detail. In the final application, the relativistic capability is used in a preliminary investigation of the complicated absorption spectrum associated with cold lanthanide elements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Quantitative study of particle size distribution in an in-situ grown Al–TiB{sub 2} composite by synchrotron X-ray diffraction and electron microscopy

Synchrotron X-ray diffraction and transmission electron microscopy (TEM) were applied to quantitatively characterize the average particle size and size distribution of free-standing TiB{sub 2} particles and TiB{sub 2} particles in an in-situ grown Al–TiB{sub 2} composite. The detailed evaluations were carried out by X-ray line profile analysis using the restricted moment method and multiple whole profile fitting procedure (MWP). Both numerical methods indicate that the formed TiB{sub 2} particles are well crystallized and free of crystal defects. The average particle size determined from different Bragg reflections by the restricted moment method ranges between 25 and 55 nm, where the smallest particle size is determined using the 110 reflection suggesting the highest lateral-growth velocity of (110) facets. The MWP method has shown that the in-situ grown TiB{sub 2} particles have a very low dislocation density (~ 10{sup 11} m{sup −} {sup 2}) and their size distribution can be described by a log-normal distribution. Good agreement was found between the results obtained from the restricted moment and MWP methods, which was further confirmed by TEM. - Highlights: • Accurate quantitative characterization of in-situ grown T{sub i}B{sub 2} particles has been achieved. • Particle size anisotropy was revealed indicating 110 facets being largest during T{sub i}B{sub 2} growth. • A wide size distribution was observed for T{sub i}B{sub 2} particles with a dominant size smaller than 100 nm.

36 MATERIALS SCIENCE↗