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 505 records · Page 28

[Ribosomal RNA Evolution]

It is generally believed that an RNA World existed at an early stage in the history of life. During this early period, RNA molecules are seen to be potentially involved in both catalysis and the storage of genetic information. Translation presents several interrelated themes of inquiry for exobiology. First, it is essential, for understanding the very origin of life, how peptides and eventually proteins might have come to be made on the early Earth in a template directed manner. Second, it is necessary to understand how a machinery of similar complexity to that found in the ribosomes of modern organisms came to exist by the time of the last common ancestor (as detected by 16S rRNA sequence studies). Third, the ribosomal RNAs themselves likely had a very early origin and studies of their history may be very informative about the nature of the RNA World. Moreover, studies of these RNAs will contribute to a better understanding of the potential roles of RNA in early evolution.During the past year we have ave conducted a comparative study of four completely sequenced bacterial genoames. We have focused initially on conservation of gene order. The second component of the project continues to build on the model system for studying the validity of variant 5S rRNA sequences in the vicinity of the modern Vibrio proteolyticus 5S rRNA that we established earlier. This system has made it possible to conduct a detailed and extensive analysis of a local portion of the sequence space. These core methods have been used to construct numerous mutants during the last several years. Although it has been a secondary focus, this work has continued over the last year such that we now have in excess of 125 V. proteolyticus derived constructs which have been made and characterized. We have also continued high resolution NMR work on RNA oligomers originally initiated by G. Kenneth Smith who was funded by a NASA Graduate Student Researcher's Fellowship Award until May of 1996. Mr. Smith developed synthesis & purification protocols in order to obtain large quantities of RNA oligomers for NMR studies. He succeeded in obtaining preliminary one dimensional spectra for an analog of the helix II/loop B/helix III domain of 5S rRNAs well as a 29mer which includes helix III and the highly characteristic 13 nucleotide loop C. Significant progress in laying the required groundwork for this work has been accomplished during the past year. A plasmid, pLeuS-1, carrying the E. coli leucine tRNA synthetase gene (LeuRS) has been obtained and used to overexpressed the leucine synthetase. A T7 runoff transcription has been established in our laboratory and successfully used to make both minihelix RNAs and tRNAs. In addition, encapsulation studies have been conducted and we have successfully encapsulated tRNA in lipid vesicles. Initial efforts to obtain peptide bond synthesis will be underway over the next several moths.

Source record↗

Energy loss by resonance line photons in an absorbing medium

The mean path length of photons undergoing repeated scatterings in media of large optical thickness is calculated from accurate numerical solutions of the transfer equation including the effect of frequency redistribution characteristic of combined Doppler and natural broadening. Energy loss by continuous absorption processes, such as ionization or dust absorption, is discussed, and asymptotic scaling laws for the energy loss, the mean path length, and the mean number of scatterings are inferred from the numerical data.

Hummer, D. G.↗

Analysis of discretization errors in LES

All numerical simulations of turbulence (DNS or LES) involve some discretization errors. The integrity of such simulations therefore depend on our ability to quantify and control such errors. In the classical literature on analysis of errors in partial differential equations, one typically studies simple linear equations (such as the wave equation or Laplace's equation). The qualitative insight gained from studying such simple situations is then used to design numerical methods for more complex problems such as the Navier-Stokes equations. Though such an approach may seem reasonable as a first approximation, it should be recognized that strongly nonlinear problems, such as turbulence, have a feature that is absent in linear problems. This feature is the simultaneous presence of a continuum of space and time scales. Thus, in an analysis of errors in the one dimensional wave equation, one may, without loss of generality, rescale the equations so that the dependent variable is always of order unity. This is not possible in the turbulence problem since the amplitudes of the Fourier modes of the velocity field have a continuous distribution. The objective of the present research is to provide some quantitative measures of numerical errors in such situations. Though the focus of this work is LES, the methods introduced here can be just as easily applied to DNS. Errors due to discretization of the time-variable are neglected for the purpose of this analysis.

Ghosal, Sandip↗

Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds

Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas. However, exact local energy-momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite the continuous effort. Without such local conservation laws, energy and momentum can be instantaneously transported across spacetime, which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory. The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic fields reside on different manifolds. To overcome this difficulty, we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact, local conservation laws, in particular the energy-momentum conservation laws, for the electromagnetic gyrokinetic system. A weak Euler-Lagrange equation is established to replace the standard Euler-Lagrange equation for the particles. It is discovered that an induced weak Euler-Lagrange current enters the local conservation laws. And it is the new physics captured by the high-order field theory on heterogeneous manifolds. A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy-momentum conservation laws electromagnetic gauge-invariant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Passage of Bolides Through the Atmosphere

Different fragmentation models are applied to a number of events, including the entry of TC3 2008 asteroid in order to reproduce existing observational data. Keywords meteoroid entry fragmentation modeling 1 Introduction Fragmentation is a very important phenomenon which occurs during the meteoroid entry into the atmosphere and adds more drastic effects than mere deceleration and ablation. Modeling of bolide fragmentation (100 106 kg in mass) may be divided into several approaches. Detail fitting of observational data (deceleration and/or light curves) allows the determination of some meteoroid parameters (ablation and shape-density coefficients, fragmentation points, amount of mass loss) (Ceplecha et al. 1993; Ceplecha and ReVelle 2005). Observational data with high accuracy are needed for the gross-fragmentation model (Ceplecha et al. 1993), which is used for the analysis of European and Desert bolide networks data. Hydrodynamical models, which describe the entry of the meteoroid including evolution of its material, are applied mainly for large bodies (>106 kg) (Boslough et al. 1994; Svetsov et al. 1995; Shuvalov and Artemieva 2002, and others). Numerous papers were devoted to the application of standard equations for large meteoroid entry in the attempts to reproduce dynamics and/or radiation for different bolides and to predict meteorite falls. These modeling efforts are often supplemented by different fragmentation models (Baldwin and Sheaffer, 1971; Borovi.ka et al. 1998; Artemieva and Shuvalov, 2001; Bland and Artemieva, 2006, and others). The fragmentation may occur in different ways. For example, few large fragments are formed. These pieces initially interact through their shock waves and then continue their flight independently. The progressive fragmentation model suggests that meteoroids are disrupted into fragments, which continue their flight as independent bodies and may be disrupted further. Similar models were suggested in numerous papers, beginning with Levin (1956) and initial interaction of fragments started to be taken into account after the paper by Passey and Melosh (1980). The progressive fragmentation model with lateral spreading of formed fragments is widely used (Artemieva and Shuvalov, 1996; Nemtchinov and Popova, 1997; Borovi.ka et al. 1998; Bland and Artemieva, 2006).

Popova, O.↗

Time-dependent water vapor desorption isotherm model of hardened cement paste

Highlights: • Short-term and long-term desorption isotherms of hardened cement paste are reviewed. • Water content at RH depends on both thermodynamic and microstructural equilibria. • An extended GAB equation is proposed to fit the cavitation drop during desorption. • A dynamic time-dependent sorption isotherm model is shown and used for simulations. • Numerical results capture anomalous drying in the same RH as seen experimentally. This study aims to enhance moisture transport modeling by elucidating the so-called anomalous water diffusion in cementitious materials. Water desorption isotherms are measured from samples at various drying stages using techniques with different durations to quantify the impact of the drying duration on the water sorption ability. A single water sorption isotherm does not solely give a relationship between the water content and relative humidity but also represents the state of microstructure. The continuous evolution of a desorption isotherm due to drying-induced microstructural rearrangement is demonstrated. For numerical modeling, the microstructural alteration can be explicitly considered through a dynamic desorption isotherm model, which governs the local thermodynamic equilibrium at the capillary meniscus. This approach is implemented into a multiphase transport code, whose ability to predict drying is validated using literature data. Finally, the effect of prolonged drying on the colloidal nature of the calcium-silicate-hydrate gel is discussed.

36 MATERIALS SCIENCE↗

Splitting regions and shrinking islands from higher point constraints

We study constraints from higher-point amplitudes on 2 → 2 scattering in the context of effective field theory (EFT) using the perturbative numerical S-matrix bootstrap. Specifically, we investigate the class of weakly coupled EFTs with amplitudes that obey the hidden zero and split conditions that are known to hold both for Tr(Φ 3 ) theory and for certain string tree amplitudes, including at 4-point the beta function. Requiring the splitting condition for the 5-point amplitude not only fixes nearly all its contact terms, but it also imposes non-linear constraints among the 4-point EFT Wilson coefficients. When included in the bootstrap, the resulting allowed region consistent with positivity is no longer convex but is restricted to a smaller non-convex region — which has a sharp corner near the string beta function! Assuming the absence of an infinite spin tower at the mass gap, the allowed region bifurcates into a trivial region (with states only above a chosen cutoff) and an island that continues to shrink around the string as more constraints are included in the bootstrap. The numerics indicate that in the absence of single-mass infinite spin towers the string beta function is the unique 4-point amplitude compatible with hidden zero and the 5-point splitting constraints. The analysis provides a prototype example for how features of higher-point amplitudes constrain the bootstrap of 4-point amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Three-dimensional fully spectral numerical method for mantle convection with depth-dependent properties

A semi-implicit fully spectral collocation method for the simulation of three-dimensional mantle convection with depth-dependent thermo-dynamic and transport properties is presented. The variable property Navier-Stokes equation expressed in terms of the primitive variable velocity and pressure is solved with the mass continuity and temperature equations. The periodic horizontal boundary conditions allow a Fourier expansion for the two horizontal directions. The stress-free, impermeable isothermal boundary conditions along with the depth dependent coefficients are handled with a Chebyshev expansion in the vertical direction. In the limit of an infinite Prandtl number appropriate to mantle convection, the inertial terms in the momentum equation are unimportant. In this case an explicit solution of a Poisson equation for pressure can be avoided; instead a fourth-order equation for vertical velocity can be solved. Simultaneous imposition of both impermeable and continuity boundary conditions during the vertical velocity evaluation is discussed. The pressure distributions on the top and bottom bounding planes were determined by means of an influence matrix technique. The numerical method employed here avoids time-splitting errors and enforces velocity boundary conditions and continuity over the entire domain, including the boundaries, to machine accuracy. Strongly time-dependent three-dimensional solutions up to a surface Rayleigh number of 1 x 10(exp 7) have been obtained. Strong upwellings, pulsating chaotically, are formed by the collective merging of cylindrical plumes.

Balachandar, S.↗

Numerical Modeling & Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (CRADA Final Report)

This project represents a continuation of the collaboration between iProTech and NLR to simulate, optimize and design the iProTech Pitching Inertial Pump (PIP) device. The objectives of this TEAMER project are twofold: 1. Refining the physical characteristics of the existing iProTech PIP WEC-Sim model to enhance the model’s fidelity and include controllable components. Key model enhancements target the inclusion of Coulomb friction, the introduction of a controllable bypass valve, and the replacement of traditional check valves with advanced motorized ones. 2. Exploring traditional and advanced control algorithms. From traditional methods like latching control to cutting-edge reinforcement learning (RL) algorithms, the goal is to ensure the PIP device's adaptability and optimal performance across a range of ocean conditions. NLR is tasked with augmenting the WEC-Sim model and implementing the control algorithms, culminating in performance comparison analyses. iProTech will update their existing 3D models, advise on model improvements, and determine crucial system metrics. WEC-Sim, developed in MATLAB/SIMULINK with Simscape Multibody, is the main piece of software that will be used in this project. Coupled with the MATLAB RL Toolbox, it offers a robust platform for in-depth simulation and optimization of the iProTech PIP device. Building on previous work to explore the PIP design space and optimize its geometry, mass distribution, center of gravity and other key parameters, this project aims to refine iProTech’s existing numerical models and develop effective control algorithms that can seamlessly integrate into their future hardware testing campaigns.

16 TIDAL AND WAVE POWER↗

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

Second-Order Invariant Domain Preserving ALE Approximation of Euler Equations

Abstract An invariant domain preserving arbitrary Lagrangian-Eulerian method for solving non-linear hyperbolic systems is developed. The numerical scheme is explicit in time and the approximation in space is done with continuous finite elements. The method is made invariant domain preserving for the Euler equations using convex limiting and is tested on various benchmarks.

Guermond, Jean-Luc↗

Continuous-time Monte Carlo renormalization group

We implement the Monte Carlo renormalization group approach in the continuous-time Monte Carlo simulation of a quantum system. In this work, we demonstrate numerically the emergent isotropy between space and time at large distances for the systems that exhibit Lorentz invariance at quantum criticality. Here, this allows us to estimate accurately the sound velocity for these quantum systems. $\textit{Q}$-state Potts models in one and two space dimensions are used to illustrate the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real-Time Modal Analysis of Electric Power Grids - The Need for Dynamic State Estimation

We articulate the reason why dynamic state estimation is needed to push the boundaries of the modal analysis of electric power grids in real-time operation. Then, we demonstrate how to unravel linear and nonlinear modes by using the extended dynamic mode decomposition along with estimates of the synchronous generators’ rotor angles and rotor speed deviations from nominal speed. The estimated modes are associated with electromechanical oscillations that take place continuously in electric power grids because of imbalances between power generation and demand. The numerical simulations are performed on a synthetic, albeit realistic, 2,000-bus network that was designed to resemble the electric power grid of Texas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analytical solutions of the Arrhenius-Semenov problem for constant volume burn

Analytical solutions to the Semenov thermal ignition problem for constant volume burn governed by Arrhenius reaction kinetics are derived. Specifically, an approximate analytical solution technique for the Arrhenius-Semenov differential equation is derived for reaction orders n ϵ R> 0 and exact solutions are also constructed for reaction orders n ϵ N : n ≤ 3. The approximation technique relies on expansion of the respective nondominant terms in the differential equation at the lower and upper bounds of the reaction progress variable in order to create a pair of integrable series. The two integrated series are then connected to create a single continuous analytical solution. Excellent agreement is observed between the analytical approximation and solutions obtained numerically. The presented approximation constitutes a simple and robust strategy for solving the Arrhenius-Semenov problem analytically.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

FY25 Electric Grid Security Annual Report

Sandia’s Electric Grid Security program advances a national vision of energy dominance and accessibility, while applying our national security -emphasis on ensuring of a secure, resilient, and affordable electric system for all users. Our achievements reflect a strategic approach combining technology development; modeling, simulation, and data analytics; and partnered demonstrations and outreach to further the adoption of advanced grid and storage technologies. Our FY25 efforts leverage the strengths of our partnerships—spanning Sandia’s core science and technology competencies as well as external technology leaders—to develop the solutions today which enable the grid of tomorrow. Key accomplishments in this report that support our strategy span our technical program areas and include: • New open-source analytical tools for systems -level planning and optimization, including significant advances to the QuESt analytical environment; • Further advancement of artificial intelligence and machine learning to enhanced grid operations and planning as we rise to the challenge of new large loads; • Development of solid-state power conversion technologies and a new medium-voltage research lab; • New technologies to assess wildfire vulnerabilities and mitigate potential impacts; • Advanced applications of new cybersecurity technologies with industry partners; • Contributions to understanding the impacts of electromagnetic pulses and geomagnetic disturbances on grid components; and • Digital twin development for hybrid microgrids with multiple generators, storage, and loads. This report indicates key areas of research and engagement and summarizes the impact of Sandia’s contributions through notable accomplishments, journal publications, patents, and technical conferences and presentations. It is provided with the hope that readers discover ways we can further team to create our modern grid and apply the outcomes of our efforts. The bulk of work described herein is funded by several offices within the U.S. Department of Energy (USDOE), including the Office of Electricity (OE); Cybersecurity, Energy Security, and Emergency Response (CESER); former offices such as the Office of Energy Efficiency and Renewable Energy (EERE), the Grid Deployment Office (GDO), the Office of Clean Energy Demonstrations (OCED), and other key programs at USDOE. As we continue to state in these annual reports, the contributors to our successes are too numerous to name here, though our team wishes to express our deep gratitude to the numerous program and project sponsors at the US Department of Energy, who often function equally as technical collaborators; our many partners in industry, academia, utilities, and other national labs; and fellow researchers and business partners at Sandia whose leadership and creativity have enabled the accomplishments described herein.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Models for infrared atmospheric radiation

Line and band models for infrared spectral absorption are discussed. Radiative transmittance and integrated absorptance of Lorentz, Doppler, and voigt line profiles were compared for a range of parameters. It was found that, for the intermediate path lengths, the combined Lorentz-Doppler (Voigt) profile is essential in calculating the atmospheric transmittance. Narrow band model relations for absorptance were used to develop exact formulations for total absorption by four wide band models. Several continuous correlations for the absorption of a wide band model were compared with the numerical solutions of the wide band models. By employing the line-by-line and quasi-random band model formulations, computational procedures were developed for evaluating transmittance and upwelling atmospheric radiance. Homogeneous path transmittances were calculated for selected bands of CO, CO2, and N2O and compared with experimental measurements. The upwelling radiance and signal change in the wave number interval of the CO fundamental band were also calculated.

Tiwari, S. N.↗

An improved analysis/synthesis capability based on dual methods - ACCESS 3

Approximation concepts and dual method algorithms are combined to create a new method for minimum weight design of structural systems. Approximation concepts convert the basic mathematical programming statement of the structural synthesis problem into a sequence of explicit primal problems of separable form. These problems are solved by constructing explicit dual functions, which are maximized subject to nonnegativity constraints. The dual method is successfully extended to deal with pure discrete and mixed continuous-discrete design variable problems. The power of the method presented is illustrated with numerical results for example problems, including a thin delta wing with fiber composite skins.

Schmit, L. A.↗

Dual methods and approximation concepts in structural synthesis

Approximation concepts and dual method algorithms are combined to create a method for minimum weight design of structural systems. Approximation concepts convert the basic mathematical programming statement of the structural synthesis problem into a sequence of explicit primal problems of separable form. These problems are solved by constructing explicit dual functions, which are maximized subject to nonnegativity constraints on the dual variables. It is shown that the joining together of approximation concepts and dual methods can be viewed as a generalized optimality criteria approach. The dual method is successfully extended to deal with pure discrete and mixed continuous-discrete design variable problems. The power of the method presented is illustrated with numerical results for example problems, including a metallic swept wing and a thin delta wing with fiber composite skins.

Fleury, C.↗