Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Systems, methods and computer program products for collaborative agent control

Systems, methods and unmanned agents for collaboratively controlling agents in a collaborative network by one or more agents continuously simulating numeric models of one or more other agents in the network to dramatically reduce the computational bandwidth required between agents, and improve the quality of shared estimates of the agent locations as well as the locations and characteristics of other objects of interest, e.g. targets. Bandwidth is reduced by using the models to intelligently filter data before communicating.

Buerger, Stephen↗
DOE OSTI2022

Contribution to the IRPHEP handbook since Physor 2018

The Organisation for Economic Co-operation and Development (OECD), the Nuclear Energy Agency (NEA), the International Reactor Physics Experiment Evaluation Project (IRPhEP) continue to preserve and evaluate integral reactor physics experiment data to support nuclear energy and technology needs. International contributions are collated within the International Handbook of Evaluated Reactor Physics Benchmark Experiments (IRPhEP Handbook). The 2021 edition of the IRPhEP Handbook now contains data for 57 unique nuclear facilities with benchmark specifications for 169 experimental series, of which four are draft evaluations. Distributed with the IRPhEP Handbook and available online is the IRPhEP Database and Analysis Tool (IDAT), allowing users to search and interrogate the data. Since Physor 2018, there have been 11 minor revisions to existing benchmark evaluations, significant revisions to two evaluations, one draft evaluation, and nine new evaluations. The corona virus pandemic has delayed international activities supporting the production of additional content for the IRPhEP Handbook. However, efforts continue for numerous endeavors to complete new reactor physics benchmark evaluations to support current and future validation needs. Those interested in serving as evaluators and/or reviewers are strongly encouraged to participate. This paper provides a summary of IRPhEP content and activities since Physor 2018. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1987

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗
NASA NTRS1982

Numerical comparisons of nonlinear convergence accelerators

As part of a continuing program of numerical tests of convergence accelerators, the iterated Aitken's Delta-squared method, Wynn's epsilon algorithm, Brezinski's theta algorithm, and Levin's u transform are compared on a broad range of test problems: linearly convergence alternating, monotone, and irregular-sign series, logarithmically convergent series, power method and Bernoulli method sequences, alternating and monotone asymptotic series, and some perturbation series arising in applications. In each category either the epsilon algorithm or the u transform gives the best results of the four methods tested. In some cases differences among methods are slight, and in others they are quite striking.

Smith, D. A.↗
NASA NTRS2019

Modern Aerocapture Guidance to Enable Reduced-Lift Vehicles at Neptune

Aerocapture is covered extensively in the literature as means of achieving orbital insertion with dramatic mass-saving results compared to fully-propulsive systems. One of the primary obstacles facing aerocapture is the inherent uncertainty associated with passing through a planet’s upper atmosphere. In-flight dispersions due to delivery errors, environment variables, and aerodynamic performance impose a large flight envelope. System studies for aerocapture often select high lift-to-drag ratios to compensate for these uncertainties. However, modern predictor-corrector guidance strategies have shown promise in recent years to provide robust control schemes in-situ. These algorithms do not rely on a pre-calculated reference trajectory and instead employ a numerical optimizer to continuously solve nonlinear equations of motion each guidance cycle. Numerical predictor-corrector strategies may provide considerable accuracy over heritage guidance schemes. The goal of this study is reproduce a landmark study of Neptune aerocapture and apply modern guidance to illustrate relative performance improvements and cost-saving potential. Capture constraints based on the theoretical corridor width are considered. Results indicate that heritage vehicles with moderate lift-to-drag ratios, lower than previous studies have indicated, may prove viable for aerocapture at Neptune.

Heidrich, C. R.↗
DOE OSTI2023

Comparative verification of hydro-mechanical fracture behavior: Task G of international research project DECOVALEX–2023

Numerical simulations become a necessity when experimental approaches cannot cover the required physical and time scale of interest. One of such area is a simulation of long-term host rock behaviors for nuclear waste disposal and simulation tools involved in the assessment must go through rigorous validation tests. The DECOVALEX project (Development of COupled models and their VAlidation against EXperiments) is dedicated to this purpose by international participants. This work is part of the ongoing phase DECOVALEX–2023 (D–2023, Task G) particularly aiming to simulate fracture behaviors under various conditions. Here, in this study, we cross-verified a variety of numerical methods including continuous and discontinuous approaches against four benchmark exercises with emphasis on numerical accuracy and parameterization of the various numerical approaches. The systematic inter-comparisons of test cases highlight advantages and disadvantages of the different numerical models. Numerical details on discretization effects (e.g. mesh density and orientation) and domain size were investigated in detail for practical applications. It became evident that meticulous attention to mesh resolution and domain size is imperative for achieving accurate numerical simulations, even for static cracks. Moreover, when comparing numerical methods to closed-form solutions for static cracks, all models successfully reproduced the maximum crack opening but encountered challenges near the crack tips. Finally, the paper discusses how to convert between and therefore compare parameters of various numerical approaches. Our benchmark studies reveal that each model necessitates a distinct number of parameters, even in simple scenarios like static crack aperture benchmarks. It is generally more practical to employ fewer parameters to mitigate model over-parameterization and enhance experimental feasibility.

58 GEOSCIENCES↗
DOE OSTI2021

Reference document for LANL stack sampling and ANSI N13.1 (Article) Gielow RL and McNamee MR 1993. Numerical Flue Gas Flow Modeling for Continuous Emissions Monitoring Applications. EPRI CEM Users Group Meeting. Baltimore. RP1961-13

American National Standard N13.1 “sets forth guidelines and performance criteria for sampling the emissions of airborne radioactive substances in the air discharge ducts and stacks of nuclear facilities. Emphasis is on extractive sampling from a location in a stack or duct where the contaminant is well mixed. At such a location, sampling may be conducted at a single point. This standard provides performance-based criteria for the use of air sampling probes, transport lines, sample collectors, sample monitoring instruments, and gas flow measuring methods. This standard also covers sampling program objectives, quality assurance issues, developing air sampling action levels, system optimization, and system performance verification. Workplace, containment, and environmental air monitoring are not addressed. Specific sample analysis methods and the reporting or interpreting of results are also not addressed.” (HPS 2011).

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2021

Analytic continuation of noisy data using Adams Bashforth residual neural network

We propose a data-driven learning framework for the analytic continuation problem in numerical quantum many-body physics. Designing an accurate and efficient framework for the analytic continuation of imaginary time using computational data is a grand challenge that has hindered meaningful links with experimental data. The standard Maximum Entropy (MaxEnt)-based method is limited by the quality of the computational data and the availability of prior information. Also, the MaxEnt is not able to solve the inversion problem under high level of noise in the data. Here we introduce a novel learning model for the analytic continuation problem using a Adams-Bashforth residual neural network (AB-ResNet). Additionally, the advantage of this deep learning network is that it is model independent and, therefore, does not require prior information concerning the quantity of interest given by the spectral function. More importantly, the ResNet-based model achieves higher accuracy than MaxEnt for data with higher level of noise. Finally, numerical examples show that the developed AB-ResNet is able to recover the spectral function with accuracy comparable to MaxEnt where the noise level is relatively small.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Noise Robustness and Experimental Demonstration of a Quantum Generative Adversarial Network for Continuous Distributions

Abstract The potential advantage of machine learning in quantum computers is a topic of intense discussion in the literature. Theoretical, numerical, and experimental explorations will most likely be required to understand its power. There have been different algorithms proposed to exploit the probabilistic nature of variational quantum circuits for generative modeling. In this paper, a hybrid architecture for quantum generative adversarial networks (QGANs) is employed and their robustness in the presence of noise is studied. A simple way of adding different types of noise to the quantum generator circuit is devised, and the noisy hybrid QGANs (HQGANs) are simulated numerically to learn continuous probability distributions, and to show that the performance of HQGANs remains unaffected. The effect of different parameters on the training time is also investigated to reduce the computational scaling of the algorithm and simplify its deployment on a quantum computer. The training on Rigetti's Aspen‐4‐2Q‐A quantum processing unit is then performed, and the results from the training are presented. The authors' results pave the way for experimental exploration of different quantum machine learning algorithms on noisy intermediate‐scale quantum devices.

Anand, Abhinav↗
NASA NTRS2024

Transfer-AE: A novel autoencoder-based impact detection model for structural digital twin

Accurately detecting the location and intensity of impacts is crucial for ensuring structural safety. Currently, AI-based structural impact detection methods are widely used for their excellent detection accuracy. However, their generalization capability is limited by the scenarios present in the training data. Many complex and dangerous impact scenarios are difficult to conduct real-world experiments on to collect sufficient samples. To capture all impact scenarios and fully leverage the advantages of AI-based detection technologies, advanced methods involve combining real-world structural monitoring data with corresponding numerical models to construct digital twins. These methods continuously refine the created numerical models with limited real-world data and provide diverse impact scenarios through numerical model simulations. However, there are inevitable differences between digital models and physical models that are challenging to correct through mechanical means. This discrepancy in data distribution between the two models significantly hinders the application of digital twin technology in impact/event identification tasks. To address this challenge, this study proposes a novel model based on autoencoders, named Transfer-AE. Transfer-AE encodes the common features of digital twins in the latent space to bridge the uncertainty gap at a macro scale between numerical models and physical models and synchronously fits the magnitude and location of the impact load in the decoder. This enables consistent detection results for the same impact event, whether the sample comes from the numerical model or the physical model. Transfer-AE includes two operating modes: Mode 1 has a fixed computational complexity with stable inference speed, but the training cost and difficulty increase with data distribution. Mode 2's computational complexity increases with data distribution, but it has a fixed training cost and speed. In both cases involving the geodesic dome structure simulating a deep space habitat and the IASC-ASCE benchmark structure, Transfer-AE demonstrated the best performance in impact localization and quantification tasks compared to mainstream domain-adaptive transfer models.

Chengjia Han↗
NASA NTRS2011

Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods

A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.

Atkins, Harold L.↗
NASA NTRS1994

A Continuing Search for a Near-Perfect Numerical Flux Scheme

While enjoying demonstrated improvement in accuracy, efficiency, and robustness over existing schemes, the Advection Upstream Splitting Scheme (AUSM) was found to have some deficiencies in extreme cases. This recent progress towards improving the AUSM while retaining its advantageous features is described. The new scheme, termed AUSM+, features: unification of velocity and Mach number splitting; exact capture of a single stationary shock; and improvement in accuracy. A general construction of the AUSM+ scheme is layed out and then focus is on the analysis of the a scheme and its mathematical properties, heretofore unreported. Monotonicity and positivity are proved, and a CFL-like condition is given for first and second order schemes and for generalized curvilinear co-ordinates. Finally, results of numerical tests on many problems are given to confirm the capability and improvements on a variety of problems including those failed by prominent schemes.

Liou, Meng-Sing↗
NASA NTRS2005

Simulations of Evaporating Multicomponent Fuel Drops

A paper presents additional information on the subject matter of Model of Mixing Layer With Multicomponent Evaporating Drops (NPO-30505), NASA Tech Briefs, Vol. 28, No. 3 (March 2004), page 55. To recapitulate: A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The model is used to perform direct numerical simulations in continuing studies directed toward understanding the behaviors of sprays of liquid petroleum fuels in furnaces, industrial combustors, and engines. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and drops, respectively, and incorporates a concept of continuous thermodynamics, according to which the chemical composition of a fuel is described by use of a distribution function. In this investigation, the distribution function depends solely on the species molar weight. The present paper reiterates the description of the model and discusses further in-depth analysis of the previous results as well as results of additional numerical simulations assessing the effect of the mass loading. The paper reiterates the conclusions reported in the cited previous article, and states some new conclusions. Some new conclusions are: 1. The slower evaporation and the evaporation/ condensation process for multicomponent-fuel drops resulted in a reduced drop-size polydispersity compared to their single-component counterpart. 2. The inhomogeneity in the spatial distribution of the species in the layer increases with the initial mass loading. 3. As evaporation becomes faster, the assumed invariant form of the molecular- weight distribution during evaporation becomes inaccurate.

Bellan, Josette↗
NASA NTRS1990

Continuous assimilation of Geosat altimetric sea level observations into a numerical synoptic ocean model of the California Current

The Geosat altimetric sea level observations for the period from January to December 1987 were continuously assimilated into a realistic wind-driven numerical synoptic ocean model of the California Current in order to evaluate the effectiveness of using a realistic synoptic ocean model for interpolating (dynamically) real altimetric sea level observations onto a regular grid. The accuracy of dynamical interpolation was tested by comparing the gridded sea level residuals to ones estimated from in situ observations (by expendable bathythermographs) collected in the California Current region during the same period. The comparison yielded nearly exact agreement at low frequency (i.e., semiannual cycle), but less agreement on month-to-month time scales of variability, possibly due to the unfiltered nature of the in situ estimates.

White, Warren B.↗
DOE OSTI2021

Recent Thermal and Hydraulic Stimulation Results at Raft River, ID EGS Site

The Raft River geothermal system, located 160 km northwest of Salt Lake City, Utah is the site of a Department of Energy Enhanced Geothermal System (EGS) stimulation demonstration project. The target well RRG-9 ST-1, is roughly 1.6 km (1 mile) south of the power plant. This well has undergone a series of thermal and hydraulic stimulation treatments. Furthermore, the well passes through approximately 1,524 m MD (5,000 ft. MD) of laterally discontinuous Quaternary and Tertiary volcanoclastic and volcanic rocks that unconformably overlie the Precambrian basement. The Elba Quartzite formation located in the Precambrian basement is the target zone for the stimulation program. Fluid injected into the well exits through a fracture zone in the Elba Quartzite at 1,567 m TVD (5,140 ft. TVD). In June 2013 low-rate stimulation began at RRG-9 ST-1 using injection water from the geothermal plant. The average injection rate, temperature, and well head pressure of the low-rate stimulation was 163 Lpm (43 gpm), 34°C, and 1,931kPa (280 psig), respectively. This injection continued through late August of 2013 when wellhead pressures were intentionally increased to approximately 5,516 kPa (800 psig) at approximately 984 Lpm (260 gpm). From September 12 to September 24, 2013 cold well water was injected. The average surface injection temperature was approximately 12°C. On September 25, 2013 injection of the warmer plant water was resumed. This was continued through the rest of 2013 and into 2014. Between April 1, 2014 and April 4, 2014, a higher rate hydraulic stimulation was conducted. Colder water (~10°C) was injected at rates of 20 bpm to 30 bpm (3,178 Lpm to 4,770 Lpm) at wellhead pressures of 5,861 kPa (850 psig) over approximately 24 hours. Microseismic activity was recorded near the Tertiary-Precambrian contact. Subsequently, injection of water from the plant was resumed. As of mid-January, 2015, over 666 million liters (176 million gallons) have been injected into the well. The fate of this water is the subject of ongoing numerical analyses and continuous site monitoring. Various analyses have been implemented to monitor and assess the ongoing stimulation program. Modified Hall plot analysis shows increasing conductivity and/or a decreasing skin factor around the wellbore. Since April, 2014, the injectivity index has increased significantly from 0.28 Lpm/kPa (0.5 gpm/psig) to 0.93 Lpm/kPa (1.7 gpm/psig).

enhanced geothermal system↗
DOE OSTI2024

Developing a predictive science of the biosphere requires the integration of scientific cultures

Increasing the speed of scientific progress is urgently needed to address the many challenges associated with the biosphere in the Anthropocene. Consequently, the critical question becomes: How can science most rapidly progress to address large, complex global problems? We suggest that the lag in the development of a more predictive science of the biosphere is not only because the biosphere is so much more complex, or because we do not have enough data, or are not doing enough experiments, but, in large part, because of unresolved tension between the three dominant scientific cultures that pervade the research community. We introduce and explain the concept of the three scientific cultures and present a novel analysis of their characteristics, supported by examples and a formal mathematical definition/representation of what this means and implies. The three cultures operate, to varying degrees, across all of science. However, within the biosciences, and in contrast to some of the other sciences, they remain relatively more separated, and their lack of integration has hindered their potential power and insight. Our solution to accelerating a broader, predictive science of the biosphere is to enhance integration of scientific cultures. The process of integration—Scientific Transculturalism—recognizes that the push for interdisciplinary research, in general, is just not enough. Unless these cultures of science are formally appreciated and their thinking iteratively integrated into scientific discovery and advancement, there will continue to be numerous significant challenges that will increasingly limit forecasting and prediction efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Engineering Privacy at the Edge: A Practical Guide to Differential Privacy in System Architectures

The rapid expansion of distributed and edge computing platforms—spanning autonomous vehicles, IoT sensors, and healthcare monitors—has heightened concerns about data privacy. Differential Privacy (DP) offers a rigorous mathematical framework to protect sensitive information while retaining analytical utility. This tutorial introduces the foundations of DP for both numerical and categorical datasets and extends the discussion to correlation-aware techniques tailored for structured and high-dimensional data. Hands-on demonstrations will begin with the PETINA (Privacy prEservaTIoN Algorithms) package for numerical data and continue with MIC-DP (Maximum Information Correlated Differential Privacy) for tabular data. Designed for researchers and practitioners in secure systems, embedded architectures, and AI accelerators, the tutorial emphasizes practical and scalable methods for integrating DP into real-world system designs.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗
DOE OSTI2022

Structure-Informed Graph Learning of Networked Dependencies for Online Prediction of Power System Transient Dynamics

Online transient analysis plays an increasingly important role in dynamic power grids as the renewable generation continues growing. Traditional numerical methods for transient analysis not only are computationally intensive but also require precise contingency information as input, and therefore, are not suitable for online applications. Existing online transient assessment studies focus on the determination of post-contingency system stability or stability margin. Here, this paper develops a novel graph-learning framework, Deep-learning Neural Representation or DNR, for online prediction, of the time-series trajectories of the system states using initial system responses that can be measured by phasor measurement units (PMUs). The proposed DNR framework consists of two sequential modules: a Network Constructor that captures network dependencies among generators, and a Dynamics Predictor that predicts the system trajectories. The key to improved prediction performance is the introduction of the spatio-temporal message-passing operations into graph neural networks with structural knowledge. Its effectiveness and scalability are validated through comparative studies, demonstrating the prediction performance under different contingency scenarios for systems of different sizes. This framework provides a solution to online predicting post-fault system dynamics based on real-time PMU measurements. Additionally, it can also be applied to facilitate the offline transient simulation without simulating the entire trajectories.

24 POWER TRANSMISSION AND DISTRIBUTION↗