Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytical methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

RAPID

Parallel computer code for the simulator for dynamics of power systems which has the capability to initiate the system and create different faults for the dynamic analysis. The code is based on time-parallel method (Parareal) with Adaptive Method Reduction (AMR). The coarse solvers for the Parareal algorithm include several Semi Analytical Solution methods. Also, Integrated simulation of coupled transmission and distribution systems can be studied.

Simunovic, Srdjan [Oak Ridge National Lab. (ORNL),↗
DOE OSTI2021

Advancing lightweight mirror design: a paradigm shift in mirror preforms by utilizing design for additive manufacturing

Additive manufacturing is a disruptive technology that can be leveraged by the redesign of components in most engineering fields. Fundamental engineering resources for lightweight mirrors were developed more than 30 years ago with a main design limitation, state of the art manufacturing. In this study, we present two design methodologies for the design of lightweight mirrors. The first method utilizes analytical expressions to design a traditional isogrid mirror, which provided the foundation for most lightweight mirrors to date. The second method employs a combination of topology optimization, lattice infill, and analytical estimation to develop an advanced lightweight mirror designed for additive manufacturing. The advanced mirror design outperforms the traditional design for each functional requirement, including a 94% reduction in predicted surface quilting and a higher specific stiffness. The manufacturing of the advanced mirror is only possible with an additive manufacturing process.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Determination of Phenolic Groups in Bio-Oils Using Revised Folin-Ciocalteu Methods: Single Cuvette and Plate Reader. Laboratory Analytical Procedure (LAP), Issue Date: May 12, 2022

Phenol components are ubiquitous in wood-derived bio-oils and biocrudes. Their reaction with other functional groups (e.g., aldehydes) may contribute to the formation of carbonaceous species and the expected thermal instability of pyrolysis oils. During hydrotreating, phenols can be recalcitrant species, requiring higher reaction temperatures than other oxygen-containing functional groups. Phenols are also present in upgraded products and have been shown to lead to catalyst deactivation during hydrotreating. They also are the first oxygenated functional group to re-appear in the upgraded product, signaling catalyst deactivation. This procedure covers the determination of phenolic compounds in fast pyrolysis oils. This Laboratory Analytical Procedure (LAP) includes two methods, the first allowing for shorter analysis time at increased reaction temperature, while the second employs a longer analysis time but at room temperature. Additionally, the use of both a single cuvette and of a plate reader are also presented.

09 BIOMASS FUELS↗
DOE OSTI2022

Determination of asymptotic normalization coefficients for the channel 16 O→ α+ 12 C: excited state 16 O(0 + ; 6.05 MeV)

Asymptotic normalization coefficients (ANC) determine the overall normalization of cross sections of peripheral radiative capture reactions. In the present paper, we treat the ANC C for the virtual decay 16 O(0 + ; 6.05 MeV)→ α+ 12 C(g.s.), the known values of which are characterized by a large spread (0.29 - 1.65) × 10 3 fm -1/2 . The ANC C is found by analytic continuation in energy of the α 12 C s-wave scattering amplitude, known from the phase-shift analysis of experimental data, to the pole corresponding to the 16O bound state and lying in the unphysical region of negative energies. To determine C, two different methods of analytic continuation are used. In the first method, the scattering data are approximated by the sum of polynomials in energy in the physical region and then extrapolated to the pole. The best way of extrapolation is chosen on the basis of the exactly solvable model. Within the second approach, the ANC C is found by solving the Schrödinger equation for the two-body α 12 C potential, the parameters of which are selected from the requirement of the best description of the phase-shift analysis data at a fixed experimental binding energy of 16 O(0 + ; 6.05 MeV) in the α+ 12 C channel. The values of the ANC C obtained within these two methods lie in the interval (886–1139) fm -1/2 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Systematic and statistical uncertainties of the hilbert-transform based high-precision FID frequency extraction method

Pulsed nuclear magnetic resonance (NMR) is widely used in high-precision magnetic field measurements. The absolute value of the magnetic field is determined from the precession frequency of nuclear magnetic moments. The Hilbert transform is widely used to extract the phase function from the observed free induction decay (FID) signal and then its frequency. In this paper, a detailed implementation of a Hilbert-transform based FID frequency extraction method is described. How artifacts and noise level in the FID signal affect the extracted phase function are derived analytically. A method of mitigating the artifacts in the extracted phase function of an FID is discussed. Correlations between noises of the phase function samples are studied for different noise spectra. We discovered that the error covariance matrix for the extracted phase function is nearly singular and improper for constructing the $\chi^2$ used in the fitting routine. A down-sampling method for fixing the singular covariance matrix has been developed, so that the minimum $\chi^2$-fit yields properly the statistical uncertainty of the extracted frequency. Other practical methods of obtaining the statistical uncertainty are also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Combining shallow-water and analytical wake models for tidal array micro-siting

For tidal-stream energy to become a competitive renewable energy source, clustering multiple turbines into arrays is paramount. Array optimisation is thus critical for achieving maximum power performance and reducing cost of energy. However, ascertaining an optimal array layout is a complex problem, subject to specific site hydrodynamics and multiple inter-disciplinary constraints. In this work, we present a novel optimisation approach that combines an analytical-based wake model, FLORIS, with an ocean model, Thetis. The approach is demonstrated through applications of increasing complexity. By utilising the method of analytical wake superposition, the addition or alteration of turbine position does not require re-calculation of the entire flow field, thus allowing the use of simple heuristic techniques to perform optimisation at a fraction of the computational cost of more sophisticated methods. Using a custom condition-based placement algorithm, this methodology is applied to the Pentland Firth for arrays with turbines of 3.05 m/s rated speed, demonstrating practical implications whilst considering the temporal variability of the tide. For a 24-turbine array case, micro-siting using this technique delivered an array 15.8% more productive on average than a staggered layout, despite flow speeds regularly exceeding the rated value. Performance was evaluated through assessment of the optimised layout within the ocean model that treats turbines through a discrete turbine representation. Used iteratively, this methodology could deliver improved array configurations in a manner that accounts for local hydrodynamic effects.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Online Parameter Estimation Methods for Adaptive Cruise Control Systems

Modeling Adaptive Cruise Control (ACC) vehicles enables the understanding of the impact of these vehicles on traffic flow. In this work, two online methods are used to provide real time system identification of ACC enabled vehicles. The first technique is a recursive least squares (RLS) approach, while the second method solves a nonlinear joint state and parameter estimation problem via particle filtering (PF). We provide a parameter identifiability analysis for both methods to analytically show that the model parameters are not identifiable using equilibrium driving. The accuracy and computational runtime of the online methods are compared to a commonly used offline simulation-based optimization (i.e., batch optimization) approach. The methods are tested on synthetic data as well as on empirical data collected directly from a 2019 model year ACC vehicle using data from sensors that are part of the stock ACC system. The online methods are scalable and provide comparable accuracy to the batch method. RLS runs in real time and is two orders of magnitude faster than the batch method for modest sized (e.g., 15 min) datasets. The particle filter also runs in real- time, and is also suitable in streaming applications in which the datasets can grow arbitrarily large.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

DPC Disposal Thermal Scoping Analysis

This is a progress report on thermal modeling for dual-purpose canister (DPCs) direct disposal that covers several available calculation methods and addresses creep and temperature-dependent properties in a salt repository. Three modeling approaches are demonstrated: A semi-analytical calculation method that uses linear solutions with superposition and imaging, to represent a central waste package in a larger array; A finite difference model of coupled thermal creep, implemented in FLAC2D; and An integrated finite difference thermal-hydrologic modeling approach for repositories in different generic host media, implemented in PFLOTRAN. These approaches are at different levels of maturity, and future work is expected to add refinements and establish the best applications for each.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Comparative Analysis of Standard and Advanced USL Methodologies [Slides]

This presentation is on the comparative analysis of standard and advanced USL methodologies. Each method has a different way to address/account for experimental biases and uncertainties. This study showed why bias and uncertainty predictions by these popular methods differ significantly in some cases. In order to establish confidence to our validation efforts, especially with new reactor concepts, in designing new critical experiments and claiming experimental coverage for our existing applications, we need to understand capabilities of our methods further.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Analysis of Seismic Signals Generated by Vehicle Traffic with Application to Derivation of Subsurface Q-Values

Correct identification and modeling of anthropogenic sources of ground motion are of considerable importance for many studies, including detection of small earthquakes and imaging seismic properties below the surface. Here, to understand signals generated by common vehicle traffic, we use seismic data recorded by closely spaced geophones normal to roads at two sites on San Jacinto fault zone. To quantify the spatiotemporal and frequency variations of the recorded ground motions, we develop a simple analytical solution accounting for propagation and attenuation of surface waves. The model reproduces well-observed bell-shaped spectrograms of car signals recorded by geophones close to roads, and it can be used to estimate frequency-dependent Q-values of the subsurface materials. The data analysis indicates Q-values of 3–40, for frequencies up to 150 Hz for road-receiver paths at the two examined sites. The derived Q-values are consistent with attenuation factors of surface waves previously obtained with other methods. The analytical results and analysis procedure provide a highly efficient method for deriving Q-values of shallow subsurface materials.

58 GEOSCIENCES↗
DOE OSTI2020

Simultaneous material, shape and topology optimization

Using three design fields we develop an optimization environment that can simultaneously optimize material, shape and topology. We use the implicit representation of the boundaries with level-set functions that define the shape and topology. Differentiable R-functions allow us to combine these shapes and topology descriptions with Boolean operations. Additionally, we incorporate design dependent-stiffness materials with another design field. Notably, this framework accommodates design dependent loads, has the ability to introduce holes, and ensures the satisfaction of optimality criteria. It builds upon the fictitious domain, ersatz material, material interpolation and level-set methods. Additionally, it also borrows from parameterized density-based topology optimization methods. Since analytical sensitivities can be computed, we use efficient nonlinear programming algorithms to update the design instead of the Hamilton–Jacobi’s scheme of level-set methods. We illustrate the features of our framework by designing a cantilever beam with octet truss microlattice, a dam with design-dependent loads, and a composite clevis plate.

42 ENGINEERING↗
DOE OSTI2023

Deterministic-Monte Carlo Hybrid Methods for Eigenvalue Sensitivity Coefficient Calculations [Slides]

Hybrid method was developed based on a need to generate accurate sensitivities for specified systems with CLUTCH. This new method provides improved sensitivities with HMF-028-001, specifically with 238 U in the large reflector region. With the importance of each voxel predetermined with the adjoint flux, the hybrid method is able to generate more accurate sensitivities. More testing is needed for other types of systems and materials (i.e., thermal and intermediate energy ranges and different moderators and reflectors). Initial results are very promising and continual development of the new hybrid method is currently in progress.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Counterdiabatic Driving with Performance Guarantees

Counterdiabatic (CD) driving has the potential to speed up adiabatic quantum state preparation by suppressing unwanted excitations. However, existing approaches either require intractable classical computations or are based on approximations that do not have performance guarantees. We propose and analyze a nonvariational, system-agnostic CD expansion method and analytically show that it converges exponentially quickly in the expansion order. In finite systems, the required resources scale inversely with the spectral gap, which we argue is asymptotically optimal. To extend our method to the thermodynamic limit and suppress errors stemming from high-frequency transitions, we leverage finite-time adiabatic protocols. In particular, we show that a time determined by the quantum speed limit is sufficient to prepare the desired ground state, without the need to optimize the adiabatic trajectory. Numerical tests of our method on the quantum Ising chain show that our method can outperform state-of-the-art variational CD approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗
DOE OSTI2022

Development and validation of a quantitative proton NMR method for the analysis of pregnenolone

Pregnenolone (PREG) is an endogenous steroid frequently sold as an over-the-counter dietary supplement touted to promote neurological and immunological health. While the PREG dietary supplement is added to the diet for health benefits, there are no FDA approved PREG drugs. However, compounded PREG drug products are available to U.S. patients. The FDA works with state regulatory authorities on the oversight of compounding activities, including developing 503A and 503B lists of bulk substances that compounders are permitted to use. PREG is one of the substances publicly nominated to be included on the 503B list. Compounded hormone therapies such as those using PREG are of interest given the lack of standardization in compounded drug products which may increase the possibility of underdosing, overdosing, or contamination. However, no USP monograph currently exists to evaluate the quality of PREG drug substance or product. To address knowledge gaps and assist in quality control, a simple and rapid quantitative proton nuclear magnetic resonance spectroscopy (qNMR) method for the identification and assay of PREG in different types of PREG products was developed and validated. PREG samples were characterized using 1D 1 H and 2D 1 H– 13 C HSQC NMR spectra. The qNMR assay method (taking approximately 10 min per NMR spectrum) was validated for precision, accuracy, specificity, robustness and linearity per ICH Q2(R1) guidance. The method was validated in a range from 0.032 to 3.2 mg/mL. As a proof of concept, seven PREG bulk substance samples, three tablet and two capsule PREG dietary supplements were assessed by the qNMR analytical procedure. NMR data from all tested samples met the expected criteria for identification and assay. The results demonstrate the potential of qNMR for the quality assessment of different types of PREG samples.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Sum-of-Fractions Methodology for Actinides in Water- and Polyethylene-Moderated and -Reflected Systems

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232 U, 233 U, 234 U, 235 U, 237 Np, 236 Pu, 238 Pu, 239 Pu, 240 Pu, 241 Pu, 242 Pu, 241 Am, 242 m Am, 243 Am, 242 Cm, 243 Cm, 244 Cm, 245 Cm, 246 Cm, 247 Cm, 249 Cf, and 251 Cf. The methodology uses available benchmark data for mixtures of 233 U, 235 U, and 239 Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

07 ISOTOPE AND RADIATION SOURCES↗