Results for “Numerical robustness” — Page 3

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

NASA NTRS2022

Accelerance Decoupling: An Approach for Removing the Influence of the Test Stand from the Integrated Modal Test

The main objective for launch vehicle (LV) modal testing is to quantify the LV’s modal properties in the free-free state (post pad separation). However, given the size of most LV systems, free-free testing is a challenge and often not feasible. With this, a test stand, typically the launch pad itself, is introduced as the means of support. This shifts the challenge to developing robust numerical methods for removing the influence of the launch pad from the integrated system modal test. The Space Launch System (SLS) is no exception where the mobile launcher (ML) is used to support the vehicle for the integrated modal test (IMT). For the IMT, it is well understood from pre-test analysis with finite element models (FEMs) of the SLS and SLS coupled to ML that the ML has a significant influence on the SLS modal properties especially in the lower frequency range where the primary SLS bending modes exist. An accelerance decoupling (AD) method has been formulated for the purpose of “subtracting out” the influence of the ML from the IMT results. With AD, the SLS decoupled frequency response functions (FRFs) are directly extracted from the IMT FRFs. The subject approach is aimed to utilize measured data only and achieve a robust FRF decoupling scheme. AD is derived from a widely used coupling technique called “Receptance Coupling” (RC). The AD core equation reverses the RC process and utilizes a pair of auxiliary equations that enable the core equation to be resolved based on measured data only. In AD, the decoupled component FRFs are extracted from the coupled system FRFs with a transformation to remove the contribution of the “subtractive component”. This paper addresses the AD’s operational flexibility to resolve SLS free-free modal properties from coupled system measured data but also the possibility to include data from FEM if there is enough confidence in the FEM or if it is asserted that the effect to the final outcome is reasonable.

Accelerance decoupling↗

DOE OSTI2026

Fast and Accurate Intersections on a Sphere

We introduce a fast, high-precision algorithm for calculating intersections between great circle arcs and lines of constant latitude on the unit sphere. We first propose a simplified intersection point formula with improved speed and numerical robustness over the ones traditionally implemented in geoscience software. We then show how algorithms based on the concept of error-free transformations (EFT) can be applied to evaluate this formula within a relative error bound that is on the order of machine precision. Here, we demonstrate that, with a vectorized and parallelized implementation, this enhanced accuracy is achieved with no compute time overhead compared to a direct calculation in hardware floating point, making our algorithm suitable for performance-sensitive applications like regridding of high-resolution climate data. In contrast, evaluating our formula using high-precision data types like quadruple precision and arbitrary precision, or using the robust intersection computation routines from the Computational Geometry Algorithms Library, leads to significant computational overhead, especially since these alternatives inhibit vectorization. More generally, our work demonstrates how EFT techniques can be combined and extended to implement nontrivial geometric calculations with high accuracy and speed.

Environmental sciences↗

DOE OSTI2015

Brownian dynamics of confined rigid bodies

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗

NASA NTRS1996

Artificial Boundary Conditions for Computation of Oscillating External Flows

In this paper, we propose a new technique for the numerical treatment of external flow problems with oscillatory behavior of the solution in time. Specifically, we consider the case of unbounded compressible viscous plane flow past a finite body (airfoil). Oscillations of the flow in time may be caused by the time-periodic injection of fluid into the boundary layer, which in accordance with experimental data, may essentially increase the performance of the airfoil. To conduct the actual computations, we have to somehow restrict the original unbounded domain, that is, to introduce an artificial (external) boundary and to further consider only a finite computational domain. Consequently, we will need to formulate some artificial boundary conditions (ABC's) at the introduced external boundary. The ABC's we are aiming to obtain must meet a fundamental requirement. One should be able to uniquely complement the solution calculated inside the finite computational domain to its infinite exterior so that the original problem is solved within the desired accuracy. Our construction of such ABC's for oscillating flows is based on an essential assumption: the Navier-Stokes equations can be linearized in the far field against the free-stream back- ground. To actually compute the ABC's, we represent the far-field solution as a Fourier series in time and then apply the Difference Potentials Method (DPM) of V. S. Ryaben'kii. This paper contains a general theoretical description of the algorithm for setting the DPM-based ABC's for time-periodic external flows. Based on our experience in implementing analogous ABC's for steady-state problems (a simpler case), we expect that these boundary conditions will become an effective tool for constructing robust numerical methods to calculate oscillatory flows.

Tsynkov, S. V.↗

DOE OSTI2018

TopoMS: Comprehensive topological exploration for molecular and condensed-matter systems

We introduce TopoMS, a computational tool enabling detailed topological analysis of molecular and condensed-matter systems, including the computation of atomic volumes and charges through the quantum theory of atoms in molecules, as well as the complete molecular graph. With roots in techniques from computational topology, and using a shared-memory parallel approach, TopoMS provides scalable, numerically robust, and topologically consistent analysis. TopoMS can be used as a command-line tool or with a GUI (graphical user interface), where the latter also enables an interactive exploration of the molecular graph. This paper presents algorithmic details of TopoMS and compares it with state-of-the-art tools: Bader charge analysis v1.0 (Arnaldsson et al., 01/11/17) and molecular graph extraction using Critic2 (Otero-de-la-Roza et al., Comput. Phys. Commun. 2014, 185, 1007). TopoMS not only combines the functionality of these individual codes but also demonstrates up to 4× performance gain on a standard laptop, faster convergence to fine-grid solution, robustness against lattice bias, and topological consistency. TopoMS is released publicly under BSD License.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Preconditioners for multiphase poromechanics with strong capillarity

This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.

42 ENGINEERING↗

DOE OSTI2017

Multi-loop positivity of the planar $ \mathcal{N} $ = 4 SYM six-point amplitude

We study the six-point NMHV ratio function in planar $ \mathcal{N} $ = 4 SYM theory in the context of positive geometry. The Amplituhedron construction of the integrand for the amplitudes provides a kinematical region in which the integrand was observed to be positive. It is natural to conjecture that this property survives integration, i.e. that the final result for the ratio function is also positive in this region. Establishing such a result would imply that preserving positivity is a surprising property of the Minkowski contour of integration and it might indicate some deeper underlying structure. We find that the ratio function is positive everywhere we have tested it, including analytic results for special kinematical regions at one and two loops, as well as robust numerical evidence through five loops. There is also evidence for not just positivity, but monotonicity in a “radial” direction. We also investigate positivity of the MHV six-gluon amplitude. While the remainder function ceases to be positive at four loops, the BDS-like normalized MHV amplitude appears to be positive through five loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗

DOE OSTI2013

A smoothed particle hydrodynamics model for droplet and film flow on smooth and rough fracture surfaces

Flow on fracture surfaces has been identified by many authors as an important flow process in unsaturated fractured rock formations. Given the complexity of flow dynamics on such small scales, robust numerical methods have to be employed in order to capture the highly dynamic interfaces and flow intermittency. In this work we present microscale free-surface flow simulations using a three-dimensional multiphase Smoothed Particle Hydrodynamics (SPH) code. Pairwise solid-fluid and fluid-fluid interaction forces are used to control the wetting behavior and cover a wide range of static and transient contact angles as well as Reynolds numbers encountered in droplet flow on rock surfaces. We validate our model via comparison with existing empirical and semi-analyical solutions for droplet flow. We use the model to investigate the occurence of adsorbed trailing films of droplets under various flow conditions and its importance for the flow dynamics when films and droplets coexist. We show that flow velocities are higher on prewetted surfaces covered by a thin film which is qualitatively attributed to the enhanced dynamic wetting and dewetting at the trailing and advancing contact line.

unsaturated flow, fractured rocks, smoothed partic↗

DOE OSTI2015

Full quantum mechanical analysis of atomic three-grating Mach–Zehnder interferometry

Atomic three-grating Mach–Zehnder interferometry constitutes an important tool to probe fundamental aspects of the quantum theory. There is, however, a remarkable gap in the literature between the oversimplified models and robust numerical simulations considered to describe the corresponding experiments. Consequently, the former usually lead to paradoxical scenarios, such as the wave–particle dual behavior of atoms, while the latter make difficult the data analysis in simple terms. Here these issues are tackled by means of a simple grating working model consisting of evenly-spaced Gaussian slits. As is shown, this model suffices to explore and explain such experiments both analytically and numerically, giving a good account of the full atomic journey inside the interferometer, and hence contributing to make less mystic the physics involved. More specifically, it provides a clear and unambiguous picture of the wavefront splitting that takes place inside the interferometer, illustrating how the momentum along each emerging diffraction order is well defined even though the wave function itself still displays a rather complex shape. To this end, the local transverse momentum is also introduced in this context as a reliable analytical tool. The splitting, apart from being a key issue to understand atomic Mach–Zehnder interferometry, also demonstrates at a fundamental level how wave and particle aspects are always present in the experiment, without incurring in any contradiction or interpretive paradox. On the other hand, at a practical level, the generality and versatility of the model and methodology presented, makes them suitable to attack analogous problems in a simple manner after a convenient tuning. - Highlights: • A simple model is proposed to analyze experiments based on atomic Mach–Zehnder interferometry. • The model can be easily handled both analytically and computationally. • A theoretical analysis based on the combination of the position and momentum representations is considered. • Wave and particle aspects are shown to coexist within the same experiment, thus removing the old wave-corpuscle dichotomy. • A good agreement between numerical simulations and experimental data is found without appealing to best-fit procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2013

Solution of the nonlinear Poisson-Boltzmann equation: Application to ionic diffusion in cementitious materials

A robust numerical solution of the nonlinear Poisson-Boltzmann equation for asymmetric polyelectrolyte solutions in discrete pore geometries is presented. Comparisons to the linearized approximation of the Poisson-Boltzmann equation reveal that the assumptions leading to linearization may not be appropriate for the electrochemical regime in many cementitious materials. Implications of the electric double layer on both partitioning of species and on diffusive release are discussed. The influence of the electric double layer on anion diffusion relative to cation diffusion is examined.

36 MATERIALS SCIENCE↗

DOE OSTI2019

A comprehensive modeling procedure for estimating statistical properties of forced ignition

In this work, a comprehensive modeling procedure for estimating the probability of ignition with application to high altitude relights of aircraft combustors is developed. In these configurations, an ignitor is used to introduce high-enthalpy discharge into a fuel-laden but stratified flow. Because of the inherent variabilities inflow conditions, kernel discharge process, and the chaotic turbulent flow, ignition cannot be described deterministically, but only as a probabilistic measure. The proposed modeling framework consists of three components. The turbulent flow is represented using the large eddy simulation (LES) framework. The ignition process is modeled using a manifold approach, where the initial kernel evolution is represented using a homogeneous reactor while the latter part of the evolution is represented as a competition between diffusion and chemical reactions using a flamelet-type mapping. A combined lookup table that can track the evolution of the kernel through these two distinct reaction stages is developed. The table lookup variables are solved within the LES framework. The resulting simulation tool is then embedded within an uncertainty quantification approach, where variations in the turbulent flow, as well as operating and kernel properties, are simulated using a Monte-Carlo-based sampling approach. Techniques to reduce computational cost are used to obtain a robust, numerically accurate, and physically representative model for engine relight. The method is validated using experimental data for ignition of methane/air mixtures. Due to the comprehensive nature of the modeling procedure, it is found that the simulation tool reproduces experimentally observed ignition probabilities over a wide range of operating conditions.

42 ENGINEERING↗

DOE OSTI2025

Jacobian-based model diagnostics and application to equation oriented modeling of a carbon capture system

It can be difficult to identify the specific variables or equations responsible for convergence issues in large mathematical programming models. The Institute for the Design of Advanced Energy Systems Integrated Platform (IDAES-IP) contains a tool to identify poorly scaled constraints and variables by searching for rows and columns of the Jacobian matrix with small L2-norms. A singular value decomposition is then performed to identify degenerate sets of equations and remaining scaling issues. Here, this work presents a flowsheet developed for post-combustion carbon capture using a monoethanolamine (MEA) solvent system as a case study. This work takes the reader through the entire process of model diagnostics and reformulation, from a basic introduction to the mathematics behind these model diagnostics to the reformulations necessary to make the model numerically robust, including a significantly modified enhancement factor model.

IDAES↗

DOE OSTI2026

Phenomena-based graph representations and applications to chemical process simulation

Rapid and robust simulation of chemical processes is critical to conduct process design, optimization, techno-economic analysis, and sustainability analysis. Yet, efficiently solving simulation models remains a challenge due to the highly coupled and nonlinear nature of the underlying algebraic equations that capture the physical phenomena taking place in the process (e.g., material and energy conservation, phase equilibrium, reactions). In this work, we show that graph-theoretic representations of the physical phenomena within unit operations can help navigate and decompose equations to systematically identify alternative approaches for fast and robust numerical solutions. Specifically, we present a graph-theoretic abstraction that captures the connectivity between the model variables/equations and use this abstraction to group variables/equations into fundamental phenomena. We show that phenomena-based decomposition of the underlying equations can help decouple nonlinearities and enforce material/energy conservation at the process level to accelerate convergence. The proposed decomposition approach differs from the more traditional sequential modular simulation approach, in which equations are grouped and decomposed by unit operations. We implemented the phenomena-based decomposition in BioSTEAM—an open-source process simulation platform in Python—and demonstrated that this approach can converge a variety of separation process models. Compared to sequential modular simulation, the phenomena-based approach can converge idealized systems faster, but it can be slower for (or even fail to converge) highly coupled and nonideal process systems.

Convergence↗

DOE OSTI2020

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗