Using the Restoring Force Surface Method to Model Modal Coupling.
Abstract not provided.
Roettgen, Daniel↗
SEARCH · Search NASA
Results for “force 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.
Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form
We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.
97 MATHEMATICS AND COMPUTING↗
Reaction analogy based forcing for incompressible scalar turbulence
Here, we present a novel reaction analogy (RA) based forcing method for generating statistically stationary scalar fields in incompressible turbulence. The new method can produce more general scalar probability density functions (PDFs), for example, quasi-double- δ PDF, than current methods, while ensuring that scalar fields remain bounded, unlike existent forcing methodologies that can potentially violate naturally existing bounds. Such features are useful for generating initial fields in nonpremixed combustion, inlet conditions for spatially developing flows, or for studying non-Gaussian scalar turbulence. The RA method mathematically models hypothetical chemical reactions that convert reactants in a mixed state back into its pure unmixed components. Various types of chemical reactions are formulated and the corresponding mathematical expressions derived such that the reaction term is smooth in scalar space and is consistent with mass conservation. For large values of the scalar forcing rate, the method produces statistically stationary quasi-double- δ scalar PDFs. Quasiuniform, Gaussian, and stretched exponential scalar statistics are recovered for smaller values of the scalar forcing rate. The shape of the scalar PDF can be further controlled by changing the stoichiometric coefficients of the reaction. Finally, the ability of the new method to produce fully developed passive scalar fields with quasi-Gaussian PDFs is also investigated, by exploring the convergence of the scalar variance spectrum to the Obukhov-Corrsin scaling and of the third-order mixed structure function to the “four-thirds” Yaglom's law.
42 ENGINEERING↗
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↗
Techniques for improved statistical convergence in quantification of eddy diffusivity moments
While recent approaches, such as the macroscopic forcing method (MFM) or Green's function-based approaches, can be used to compute Reynolds-averaged Navier-Stokes closure operators using forced direct numerical simulations, MFM can also be used to directly compute moments of the effective nonlocal and anisotropic eddy diffusivities. The low-order spatial and temporal moments contain limited information about the eddy diffusivity but are often sufficient for quantification and modeling of nonlocal and anisotropic effects. However, when using MFM to compute eddy diffusivity moments, the statistical convergence can be slow for higher-order moments. In this work, we demonstrate that using the same direct numerical simulation (DNS) for all forced MFM simulations improves statistical convergence of the eddy diffusivity moments. We present its implementation in conjunction with a decomposition method that handles the MFM forcing semianalytically and allows for consistent boundary condition treatment, which we develop for both scalar and momentum transport. We demonstrate that for a two-dimensional Rayleigh-Taylor instability case study, using the same DNS for all forced MFM simulations results in convergence with 𝒪(100) simulations rather than 𝒪(1000) simulations. In conclusion, we then demonstrate the impacts of improved convergence on the quantification of the eddy diffusivity.
general physics↗
A machine learning approach for increased throughput of density functional theory substitutional alloy studies
In this study, a machine learning-based technique is developed to reduce the computational cost required to explore large design spaces of substitutional alloys. The first advancement is based on a neural network (NN) approach to predict the initial position of minority and majority ions prior to DFT relaxation. The second advancement is to allow the NN to predict the total energy for every possibility minority ion position and select the most stable configuration in the absence of relaxing each trial minority configuration. A bismuth oxide materials system, (Bi x La y Yb z ) 2 MoO 6 , was used as the model system to demonstrate the developed methods and quantify the resulting computational speedup. Compared to a brute force method that requires the calculation of every permutation of minority configuration and subsequent DFT relaxation, a 1.3× speedup was realized if the NN predicted the initial configuration of ions prior to relaxation. Implementation of the second advancement allowed the NN to predict the total energy for all possible trial configurations and downselect the most stable configurations prior to relaxation, resulting in a speedup of approximately 37×. Validation was done by comparing position and energy between the NN and DFT predictions. A maximum position vector mean squared error (MSE) of 1.6 × 10 -2 and a maximum energy MSE of 2.3 × 10 -7 was predicted for the worst case configuration. This method demonstrates a significant computational speedup, which has the potential for even greater computational savings for larger compositional design spaces.
36 MATERIALS SCIENCE↗
A propagation-based fault detection and discrimination method and the optimization of sensor deployment
Industrial processes can be affected by faults having a serious impact on operation when not promptly detected and diagnosed. Here in this paper, a propagation-based fault detection and discrimination(PFDD) method is proposed to develop a strategy for fault diagnosis while in the design phase of a system. The PFDD method constructs the system model using the Integrated System Fault Analysis(ISFA) technique. Based on the system model, the propagation of hardware and software faults are simulated qualitatively. Given the results of the simulation, the process by which a fault propagates can be characterized using the qualitative features of system variables including the deviation of the system variables from their expected values, the variation of the system variables over time, and the order in which each variable is influenced during the propagation of the fault. The strategy by which a fault can be detected and discriminated is defined using those features. The PFDD method supports the detection and discrimination of faults in both steady states and transient states. Based on the PFDD method, the optimization of sensor deployment in a system is discussed. A brute force algorithm is developed to examine the system’s capability at diagnosing faults and the cost of sensor deployment for all possible configurations of sensors. The optimal sensor deployment strategy can be derived accordingly. However, the brute force method is only applicable to small-scale systems due to its high computational cost. A genetic algorithm is used to optimize sensor deployment in large-scale systems. The PFDD and sensor deployment optimization methods are applied to the Experimental Breeder Reactor II (EBR-II) for verification.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Discovering the Most Severe K-Point Failure Based on Reinforcement Learning: Preprint
Smart devices are essential to ensure the stability of the power grid and resilience to intermittent energy production. However, smart devices can also be the target of cyber adversaries that may exploit false data injection attacks (FDIAs) to induce unstable grid conditions. A practical consideration of FDIA mitigation approaches is addressed here: given a finite available budget, for which smart device should cyber-threat mitigation be deployed first? In this work, this question is answered by identifying the so-called most-sensitive devices, i.e., the devices that, if compromised, can let an adversary induce the most serious grid instabilities. The method proposed utilizes an adversarial reinforcement learning (RL) framework to identify the k-mostsensitive smart devices (here, smart inverters). The adversarial agent can tamper with the compromised inverters' active and reactive operating power setup points, with the goal of maximizing voltage deviations. Numerical results show that the proposed RL method finds the optimal attack scenarios for 1-point failure and the near-optimal solution for the 2-point case. Additionally, the proposed RL method achieves an 8.8 speed-up ratio in running time compared to the brute force method for the 2-point case.
97 MATHEMATICS AND COMPUTING↗
A demonstration of force estimation and regularization methods for multi-shaker testing.
Abstract not provided.
Schultz, Ryan↗
Reconstructing phase-resolved hysteresis loops from first-order reversal curves
The first order reversal curve (FORC) method is a magnetometry based technique used to capture nanoscale magnetic phase separation and interactions with macroscopic measurements using minor hysteresis loop analysis. This makes the FORC technique a powerful tool in the analysis of complex systems which cannot be effectively probed using localized techniques. However, recovering quantitative details about the identified phases which can be compared to traditionally measured metrics remains an enigmatic challenge. We demonstrate a technique to reconstruct phase-resolved magnetic hysteresis loops by selectively integrating the measured FORC distribution. From these minor loops, the traditional metrics—including the coercivity and saturation field, and the remanent and saturation magnetization—can be determined. In order to perform this analysis, special consideration must be paid to the accurate quantitative management of the so-called reversible features. This technique is demonstrated on three representative materials systems, high anisotropy FeCuPt thin-films, Fe nanodots, and SmCo/Fe exchange spring magnet films, and shows excellent agreement with the direct measured major loop, as well as the phase separated loops.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations
The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.
54 ENVIRONMENTAL SCIENCES↗
Element-specific first order reversal curves measured by magnetic transmission x-ray microscopy
The first-order reversal curve (FORC) method is a macroscopic measurement technique that can be used to extract quantitative and microscopic properties of hysteretic systems. Using magnetic transmission x-ray microscopy (MTXM), local element-specific FORC measurements are performed on a 20 nm thick film of CoTb. The FORCs measured with microscopy reveal a step-by-step domain evolution under the magnetic field cycling protocol and provide a direct visualization of the mechanistic interpretation of FORC diagrams. They are compared with magnetometry FORCs and show good quantitative agreement. Furthermore, the high spatial resolution and element-specific sensitivity of MTXM provide new capabilities to measure FORCs in small regions or specific phases within multicomponent systems, including buried layers in heterostructures. The ability to perform FORCs on very small features is demonstrated with the MTXM-FORC measurement of a rectangular microstructure with vortex-like Landau structures. This work demonstrates the confluence of two uniquely powerful techniques to achieve quantitative insight into nanoscale magnetic behavior.
36 MATERIALS SCIENCE↗
Constraining Galaxy-Halo connection using machine learning
We investigate the potential of machine learning (ML) methods to model small-scale galaxy clustering for constraining Halo Occupation Distribution (HOD) parameters. Our analysis reveals that while many ML algorithms report good statistical fits, they often yield likelihood contours that are significantly biased in both mean values and variances relative to the true model parameters. This highlights the importance of careful data processing and algorithm selection in ML applications for galaxy clustering, as even seemingly robust methods can lead to biased results if not applied correctly. ML tools offer a promising approach to exploring the HOD parameter space with significantly reduced computational costs compared to traditional brute-force methods if their robustness is established. Using our ANN-based pipeline, we successfully recreate some standard results from recent literature. Properly restricting the HOD parameter space, transforming the training data, and carefully selecting ML algorithms are essential for achieving unbiased and robust predictions. Among the methods tested, artificial neural networks (ANNs) outperform random forests (RF) and ridge regression in predicting clustering statistics, when the HOD prior space is appropriately restricted. We demonstrate these findings using the projected two-point correlation function (w p (r p )), angular multipoles of the correlation function (ξ ℓ (r)), and the void probability function (VPF) of Luminous Red Galaxies from Dark Energy Spectroscopic Instrument mocks. Our results show that while combining w p (r p ) and VPF improves parameter constraints, adding the multipoles ξ 0 , ξ 2 , and ξ 4 to w p (r p ) does not significantly improve the constraints.
cosmology↗
Fast calculation of diffraction patterns from an ensemble of aligned molecules
We report an algorithm to calculate electron diffraction patterns for molecules with anisotropic angular distribution, which is significantly faster than existing methods. The algorithm uses a transform to convert the molecular orientation distribution, which is a function of three Euler angles, to the atom-pair distribution functions which depend on the polar and azimuthal angles. The diffraction signal can then be calculated from the atom-pair distributions. We demonstrate the computation method numerically by calculating electron diffraction patterns for a symmetric top molecule (trifluoroiodomethane) and an asymmetric top molecule (formaldehyde) and show that it reduces the calculation time by approximately two orders of magnitude compared to the standard brute-force method. Here, the method can also be applied to the calculation of x-ray diffraction patterns.
74 ATOMIC AND MOLECULAR PHYSICS↗
The statistical defensibility of forensic proteomics
The U.S. Federal court system maintains very high standards for admissibility of scientific evidence. In particular, to allow such evidence into court, a judge must be satisfied it is both reliable, in other words, based on a robust, reproducible, and accurate method, and relevant, meaning it adds weight to the prosecution or defense claims in a case. In 1993, the U.S. Supreme Court issued landmark guidelines for judges considering whether or not to admit scientific testimony. The resulting Daubert criteria changed the forensic sciences in a profound way, leading investigators, researchers, and legal experts to question the scientific validity of established forensic methods that had been widely accepted for years. The Daubert criteria also caused a paradigm shift in the development of new forensic methods, forcing researchers to better prepare new techniques for use in an increasingly adversarial environment. In this chapter, we discuss the current state of forensic proteomics in the context of modern forensics. We present the Daubert criteria, discuss seven elements of a defensible method, and provide guidelines for building statistical defensibility of this emerging discipline.
Jarman, Kristin H.↗
A Method for Determining Impact Force for Single and Tri Axis Resonant Plate Shock Simulations.
Abstract not provided.
Ferri, Brian Anthony↗
Comparison of ozone formation attribution techniques in the northeastern United States
The Integrated Source Apportionment Method (ISAM) has been revised in the Community Multiscale Air Quality (CMAQ) model. This work updates ISAM to maximize its flexibility, particularly for ozone (O 3 ) modeling, by providing multiple attribution options, including products inheriting attribution fully from nitrogen oxide reactants, fully from volatile organic compound (VOC) reactants, equally from all reactants, or dynamically from NO x or VOC reactants based on the indicator gross production ratio of hydrogen peroxide (H 2 O 2 ) to nitric acid (HNO 3 ). The updated ISAM has been incorporated into the most recent publicly accessible versions of CMAQ (v5.3.2 and beyond). This study's primary objective is to document these ISAM updates and demonstrate their impacts on source apportionment results for O 3 and its precursors. Additionally, the ISAM results are compared with the Ozone Source Apportionment Technology (OSAT) in the Comprehensive Air-quality Model with Extensions (CAMx) and the brute-force method (BF). All comparisons are performed for a 4 km horizontal grid resolution application over the northeastern US for a selected 2 d summer case study (9 and 10 August 2018). General similarities among ISAM, OSAT, and BF results add credibility to the new ISAM algorithms. However, some discrepancies in magnitude or relative proportions among tracked sources illustrate the distinct features of each approach, while others may be related to differences in model formulation of chemical and physical processes. Despite these differences, OSAT and ISAM still provide useful apportionment data by identifying the geographical and temporal contributions of O 3 and its precursors. Both OSAT and ISAM attribute the majority of O 3 and NO x contributions to boundary, mobile, and biogenic sources, whereas the top three contributors to VOCs are found to be biogenic, boundary, and area sources.
54 ENVIRONMENTAL SCIENCES↗
Review of thermophysical property methods applied to fueled and un-fueled molten salts
Recent investigation into designs of molten salt reactors (MSR) requires accurate data of candidate salts. Review of the literature shows lacking thermophysical (TP) property data for coolants and fuel-bearing salts of Advanced High Temperature Reactors (AHTR). Measurement uncertainties of current data are even more deficient, and this affects ongoing research where MSR modeling and testing are being performed. Melting temperature, viscosity, density, specific heat capacity, and thermal and electrical conductivity measurements are often not arranged together or discussed in great detail, especially for fuel-bearing salts. Raw data from key findings spanning over 70 years for LiF, LiF-BeF 2 , LiF-NaF-KF, KCl-MgCl 2 , KF-ZrF 4 , and NaF-KF-MgF 2 have been organized and presented here. Brief description of methods used to obtain these results, in addition to the TP data of some of these salts as fuel carriers, has also been discussed. Thermal conductivity data are particularly weak due to older, high-error measurement techniques. More recent measurement systems to evaluate salt thermal conductivity show promise for obtaining consistent absolute values. This work concludes that measurement systems such as transient hot-wire method, forced Rayleigh scattering, and laser flash analysis provide the most promise to fill these data gaps and reduce current uncertainty in thermal conductivity measurements of both fueled and un-fueled molten salts.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗