SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Can High Current Density HTS Magnets be Quench Protected Using Methods Used to Protect High Current Density LTS Magnets?

In the 1970s the Lawrence Berkeley Laboratory (LBL) was testing high current density magnets where the EJ 2 limit was up to two orders of magnitude higher than the quench protection methods of the period would permit at the time. LBL was designing a high current density (J) thin solenoid with a large volume and stored energy (E). The accepted EJ 2 limit at the time was ~10 23 J A 2 m -4 . In order to exceed the EJ 2 limit by two orders of magnitude, one had to greatly increase the quench protection voltages (V) and magnet current (I). In 1977, LBL tested four quench protection methods on a 2-meter diameter high current density solenoid with a stored energy of 2 MJ. Three of these quench protection methods worked well, but the fourth didn't appear to work. These methods will be discussed as to whether they can be applied to HTS or MgB 2 high current density magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Transient Stability-Enhancing Method for Grid-Forming Inverters Under Current Limiting

The primary controller of a grid-forming (GFM) inverter governs the internal reference voltage and angle, which enables the inverter to naturally synchronize and share power with the connected grid; however, during disturbances on the grid caused by, e.g., line faults, overloading, or frequency and phase shifts, a GFM inverter can be pushed into a current-limited operation to protect the device against thermal damage. This can lead to a windup of the internal reference frequency generated by the primary controller and cause a total loss of synchronism with the grid. Here, to make the primary controller of GFM inverters more robust against transient instability, this paper proposes a method that manipulates the reference signals governed by the primary controller by introducing the concept of fictitious power. The fictitious power counteracts the integrator windup in the primary controller, which enables the current-limited inverter to conceive a stable operating point. This method is benchmarked against five state-of-the-art GFM inverter stabilization methods to evaluate its stability-enhancing performance. Full-order electromagnetic-transient simulations and experimental results of a hardware inverter test setup validate the benefits of the method to improve the transient stability of GFM inverters in the face of grid disturbances.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Alignment Method for Synchronized Phase Angle Measurement With Presence of Practical Time Shift

Synchronized Phase Angle Measurements (SPAMs) are widely used in power systems in the applications of situational awareness. However, the practical time shifts can lead to unexpected angle differences among Phasor Measurement Units (PMUs) manufactured from various vendors. Moreover, since most of PMUs calculate phase angle via Discrete Fourier Transform based approaches, this issue becomes even worse under the off-nominal frequency condition. To mitigate the impact of practical timing shift, this letter presents a fast method to rectify the SPAM for appropriate alignment. To verify the performance of the proposed method, laboratory and field tests have been conducted by implementing the method in PMUs and phasor data concentrators, respectively. The results demonstrate that the standard deviations of the phase angle differences have significantly decreased to 0.1$^{\circ }$ order with the adoption of the proposed method.

42 ENGINEERING↗

DOE OSTI2024

Estimation and Visualization of Isosurface Uncertainty from Linear and High-Order Interpolation Methods

Isosurface visualization is fundamental for exploring and analyzing 3D volumetric data. Marching cubes (MC) algorithms with linear interpolation are commonly used for isosurface extraction and visualization. Although linear interpolation is easy to implement, it has limitations when the underlying data is complex and high-order, which is the case for most real-world data. Linear interpolation can output vertices at the wrong location. Its inability to deal with sharp features and features smaller than grid cells can lead to an incorrect isosurface with holes and broken pieces. Despite these limitations, isosurface visualizations typically do not include insight into the spatial location and the magnitude of these errors. We utilize high-order interpolation methods with MC algorithms and interactive visualization to highlight these uncertainties. Our visualization tool helps identify the regions of high interpolation errors. It also allows users to query local areas for details and compare the differences between isosurfaces from different interpolation methods. In addition, we employ high-order methods to identify and reconstruct possible features that linear methods cannot detect. We showcase how our visualization tool helps explore and understand the extracted isosurface errors through synthetic and real-world data.

Ouermi, Timbwaoga↗

DOE OSTI2020

A Method to Represent a Well in a Three‐Dimensional Discrete Fracture Network Model

Abstract In discrete fracture network (DFN) modeling, fractures are randomly generated and placed in the model domain. The rock matrix is considered impermeable. Small fractures and isolated fractures are often ignored to reduce computational expense. As a result, the rock matrix between fractures could be large and intersections may not be found between a well introduced in the model and the hydraulically connected fracture networks (fracture backbones). To overcome this issue, this study developed a method to conceptualize a well in a three‐dimensional (3D) DFN using two orthogonal rectangular fractures oriented along the well's axis. Six parameters were introduced to parameterize the well screen and skin zone, and to control the connectivity between the well and the fracture backbones. The two orthogonal fractures were discretized using a high‐resolution mesh to improve the quality of flow and transport simulations around and along the well. The method was successfully implemented within dfnWorks 2.0 (Hyman et al. 2015) to incorporate a well in a 3D DFN and to track particles leaving an injection well and migrating to a pumping well. Verification of the method against MODFLOW/MODPATH found a perfect match in simulated hydraulic head and particle tracking. Using three examples, the study showed that the method ensured the connectivity between wells and fracture backbones, and honored the physical processes of flow and transport along and around wells in DFNs. Recommendations are given for estimating the values of the six introduced well parameters in a real‐world case study.

Pham, Hai (ORCID:0000000310164394)↗

DOE OSTI2020

An Assessment of Homogeneous Mixture Method Cavitation Models in Predicting Cavitation in Nozzle Flow

The homogeneous mixture method (HMM) is a popular class of models used in the computational prediction of cavitation. Several cavitation models have been developed for use with HMM to govern the development and destruction of vapor in a fluid system. Two models credited to Kunz and Schnerr–Sauer are studied in this paper. The goal of this work is to provide an assessment of the two cavitation submodels in their ability to predict cavitation in nozzle flow. Validation data were obtained via experiments which employ both passive cavitation detection, (PCD) via acoustic sensing and optical cavitation detection (OCD) via camera imaging. The experiments provide quantitative information on cavitation inception and qualitative information on the vapor in the nozzle. The results show that initial vapor formation is not predicted precisely but within reason. A sensitivity analysis of the models to input parameters shows that the Schnerr–Sauer method does not depend upon the estimation of nuclei size and number density. Small changes in the vapor formation rate but not the total vapor volume can be seen when weighting parameters are modified. In contrast, changes to the input parameters for the Kunz model greatly change the final total vapor volume prediction. The assessment also highlights the influence of vapor convection within the method. Finally, the analysis shows that if the fluid and nozzle walls do not support nuclei larger than 40 μm, the methods would still predict cavitation when indeed there would be none in practice.

Engineering↗

DOE OSTI2024

Development of a Method for Shape Optimization for a Gas Turbine Fuel Injector Design Using Metal-Additive Manufacturing

Adjoint shape optimization has enabled physics-based optimal designs for aerodynamic surfaces. Additive manufacturing (AM) makes it possible to manufacture complex shapes. However, there has been a gap between optimal and manufacturable surfaces due to the inherent limitations of commercial computational fluid dynamics (CFD) codes to implement geometric constraints during adjoint computation. In such cases, the design sensitivities are exported and used to perform constrained shape modifications using parametric information stored in computer aided design (CAD) files to satisfy manufacturability constraints. However, modifying the design using adjoint methods in CFD solvers and performing constrained shape modification in CAD can lead to inconsistencies due to different shape parameterization schemes. This paper describes a method to enable the simultaneous optimization of the fluid domain and impose AM manufacturability constraints, resolving one of the key issues of geometry definition for isogeometric analysis. Similar to a grid convergence study, the proposed method verifies the consistencies between shape parameterization techniques present within commercial CAD and CFD software during mesh movement as a part of the adjoint shape optimization routine. By identifying the appropriate parameters essential to a shape optimization study, the error metric between the different parameterization techniques converges to demonstrate sufficient consistencies for justifiable exchange of data between CAD and CFD. For the identified shape optimization parameters, the error metric to measure the deviation between the two parameterization schemes lies within the AM laser-powder bed fusion (L-PBF) process tolerance. Additionally, comparison for subsequent objective function calculations between iterations of the optimization loop showed acceptable differences within 1% variation between the modified geometries obtained using the two parameterization schemes. This method provides justification for the use of multiphysics guided adjoint design sensitivities computed in CFD software to perform shape modifications in CAD to incorporate AM manufacturability constraints during the shape optimization loop such that optimal designs are also additively manufacturable.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

High Order Strong Stability Preserving MultiDerivative Implicit and IMEX Runge--Kutta Methods with Asymptotic Preserving Properties

In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A New Semistructured Algebraic Multigrid Method

Multigrid methods are well suited to large massively parallel computer architectures because they are mathematically optimal and display good parallelization properties. Since current architecture trends are favoring regular compute patterns to achieve high performance, the ability to express structure has become much more important. The hypre software library provides high-performance multigrid preconditioners and solvers through conceptual interfaces, including a semistructured interface that describes matrices primarily in terms of stencils and logically structured grids. This paper presents a new semistructured algebraic multigrid (SSAMG) method built on this interface. The numerical convergence and performance of a CPU implementation of this method are evaluated for a set of semistructured problems. In conclusion, SSAMG achieves significantly better setup times than hypre’s unstructured AMG solvers and comparable convergence. In addition, the new method is capable of solving more complex problems than hypre’s structured solvers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

A Fourth-Order Embedded Boundary Finite Volume Method for the Unsteady Stokes Equations with Complex Geometries

A fourth-order finite volume embedded boundary (EB) method is presented for the unsteady Stokes equations. The algorithm represents complex geometries on a Cartesian grid using EB, employing a technique to mitigate the ``small cut-cell"" problem without mesh modifications, cell merging, or state redistribution. Spatial discretizations are based on a weighted least-squares technique that has been extended to fourth-order operators and boundary conditions, including an approximate projection to enforce the divergence-free constraint. Solutions are advanced in time using a fourth-order additive implicit-explicit Runge-Kutta method, with the viscous and source terms treated implicitly and explicitly, respectively. Formal accuracy of the method is demonstrated with several grid convergence studies, and results are shown for an application with a complex bio-inspired material. In conclusion, the developed method achieves fourth-order accuracy and is stable despite the pervasive small cells arising from complex geometries.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Explicit Runge–Kutta Methods that Alleviate Order Reduction

Explicit Runge–Kutta (RK) methods are susceptible to a reduction in the observed order of convergence when applied to an initial boundary value problem with time-dependent boundary conditions. We study conditions on explicit RK methods that guarantee high order convergence for linear problems; we refer to these conditions as weak stage order conditions. We prove a general relationship between the method’s order, weak stage order, and number of stages. Furthermore, we derive explicit RK methods with high weak stage order and demonstrate, through numerical tests, that they avoid the order reduction phenomenon up to any order for linear problems and up to order three for nonlinear problems.

explicit Runge–Kutta↗

DOE OSTI2026

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗

DOE OSTI2023

Performance Analysis and Optimal Node-aware Communication for Enlarged Conjugate Gradient Methods

Krylov methods are a key way of solving large sparse linear systems of equations but suffer from poor strong scalability on distributed memory machines. Furthermore, this is due to high synchronization costs from large numbers of collective communication calls alongside a low computational workload. Enlarged Krylov methods address this issue by decreasing the total iterations to convergence, an artifact of splitting the initial residual and resulting in operations on block vectors. In this article, we present a performance study of an enlarged Krylov method, Enlarged Conjugate Gradients (ECG), noting the impact of block vectors on parallel performance at scale. Most notably, we observe the increased overhead of point-to-point communication as a result of denser messages in the sparse matrix-block vector multiplication kernel. Additionally, we present models to analyze expected performance of ECG, as well as motivate design decisions. Most importantly, we introduce a new point-to-point communication approach based on node-aware communication techniques that increases efficiency of the method at scale.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Methods—Kintsugi Imaging of Battery Electrodes: Distinguishing Pores from the Carbon Binder Domain using Pt Deposition

The mesostructure of porous electrodes used in lithium-ion batteries strongly influences cell performance. Accurate imaging of the distribution of phases in these electrodes would allow this relationship to be better understood through simulation. However, imaging the nanoscale features in these components is challenging. While scanning electron microscopy is able to achieve the required resolution, it has well established difficulties imaging porous media. This is because the flat imaging planes prepared using focused ion beam milling will intersect with the pores, which makes the images hard to interpret as the inside walls of the pores are observed. It is common to infiltrate porous media with resin prior to imaging to help resolve this issue, but both the nanoscale porosity and the chemical similarity of the resins to the battery materials undermine the utility of this approach for most electrodes. In this study, a technique is demonstrated which uses in situ infiltration of platinum to fill the pores and thus enhance their contrast during imaging. Reminiscent of the Japanese art of repairing cracked ceramics with precious metals, this technique is referred to as the kintsugi method. The images resulting from applying this technique to a conventional porous cathode are presented and then segmented using a multi-channel convolutional method. We show that while some cracks in active material particles were empty, others appear to be filled (perhaps with the carbon binder phase), which will have implications for the rate performance of the cell. Energy dispersive X-ray spectroscopy was used to validate the distribution of phases resulting from image analysis, which also suggested a graded distribution of the binder relative to the carbon additive. The equipment required to use the kintsugi method is commonly available in major research facilities and so we hope that this method will be rapidly adopted to improve the imaging of electrode materials and porous media in general.

25 ENERGY STORAGE↗

DOE OSTI2021

A Spatiotemporal Prediction Model for Regional Scheduling of Shared Bicycles Based on the INLA Method

Dock-less bicycle-sharing programs have been widely accepted as an efficient mode to benefit health and reduce congestions. And modeling and prediction has always been a core proposition in the field of transportation. Most of the existing demand prediction models for shared bikes take regions as research objects; therefore, a POI-based method can be a beneficial complement to existing research, including zone-level, OD-level, and station-level techniques. Point of interest (POI) is the location description of spatial entities, which can reflect the cycling route characteristics for both commuting and noncommuting trips to a certain extent, and is also the main generating point and attraction point of shared-bike travel flow. In this study, we make an effort to model a POI-level cycling demand with a Bayesian hierarchical method. The proposed model combines the integrated nested Laplace approximation (INLA) and random partial differential equation (SPDE) to cope with the huge computation in the modeling process. In particular, we have adopted the dock-less bicycle-sharing rental records of Mobike as a case study to validate our method; the study area was one of the fastest growing urban districts in Shanghai in August 2016. The operation results show that the method can help better understand, measure, and characterize spatiotemporal patterns of bike-share ridership at the POI level and quantify the impact of the spatiotemporal effect on bicycle-sharing use.

Yu, Zhuoran↗

DOE OSTI2021

Recognition Method of Vehicle Cluster Situation Based on Set Pair Logic considering Driver’s Cognition

The recognition of vehicle cluster situations is one of the critical technologies of advanced driving, such as intelligent driving and automated driving. The accurate recognition of vehicle cluster situations is helpful for behavior decision‐making safe and efficient. In order to accurately and objectively identify the vehicle cluster situation, a vehicle cluster situation model is proposed based on the interval number of set pair logic. The proposed model can express the traffic environment’s knowledge considering each vehicle’s characteristics, grouping relationships, and traffic flow characteristics in the target vehicle’s interest region. A recognition method of vehicle cluster situation is designed to infer the traffic environment and driving conditions based on the connection number of set pair logic. In the proposed model, the uncertainty of the driver’s cognition is fully considered. In the recognition method, the relative uncertainty and relative certainty of driver’s cognition, traffic information, and vehicle cluster situation are fully considered. The verification results show that the proposed recognition method of vehicle cluster situations can realize accurate and objective recognition. The proposed anthropomorphic recognition method could provide a basis for vehicle autonomous behavior decision‐making.

Liu, Shijie↗

DOE OSTI2022

Verification and Validation of Two Hydrodynamic Methods for Simulations of High Energy Density Physics Problems

A 3D verification and validation suite of test problems is presented and used to evaluate hydrodynamic methods within a radiation hydrodynamics code, xRAGE. These test problems exercise different levels of complexity, building towards ICF problems which in addition to hydrodynamics also include three temperature plasma physics, thermal conduction, and radiation diffusion. Among the problems in the test suite are the Kidder ball problem, the Verney shell problem, and a 5-material compression problem, which exercise different purely hydrodynamic methods implemented within xRAGE. There is excellent agreement between 2D and 3D XRAGE simulation results and between the xRAGE results and the benchmark solutions. Two 3D ICF test problems are also presented, based on an OMEGA direct drive capsule experiment and on a NIF indirect drive capsule experiment. It is demonstrated that the newer unsplit hydrodynamic method in xRAGE produces more vorticity relative to the older default method. For the indirect drive capsule, the 3D simulations are in reasonable agreement with the experimental values of ion temperature and neutron production.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Method of Moving Asymptotes Optimization Algorithm in Python

Python implementation of the Method of Moving Asymptotes optimization algorithm described in[Svanberg, K., The method of moving asymptotes- a new method for structural optimization. International journal for numerical methods in engineering, 1987. 24(2): p.359 (https://onlinelibrary.wiley.com/doi/abs/10.1002/nme.1620240207). Originally implemented in[GetDP](https://ieeexplore.ieee.org/stamp/stamp.jsp?arnumber=717799) [project](https://gitlab.onelab.info/getdp/getdp), although it has been deleted in their repository.

Salazar De Troya, Miguel↗