Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher-Order Modifications”

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 19 records

Higher-Order Approximations for Stabilizing Zero-Energy Modes in Peridynamics Crystal Plasticity Models with Large Horizon Interactions

The non-ordinary state-based peridynamics theory combines non-local dynamic techniques with a desirable correspondence material principle, allowing for the use of continuum mechanics constitutive models. Such an approach presents a unique capability for solving problems involving discontinuities (e.g., strain localization, fracture, and fragmentation). However, the correspondence-based peridynamics models often suffer from zero-energy mode instabilities in numerical implementation, primarily due to the approximations of the non-local deformation gradient tensor. This paper focuses on a computational scheme for eliminating the zero-energy mode oscillations using a choice of influence functions that improve the truncation error in a higher-order Taylor series expansion of the deformation gradient. The novelty here is a tensor-based derivation of the linear constraint equations, which can be used to systematically identify the particle interaction weight functions for various user-specified horizon radii. In this paper, the proposed higher-order stabilization scheme is demonstrated for multi-dimensional examples involving polycrystalline and composite microstructures, along with comparisons against conventional finite element methods. The proposed stabilization scheme is shown to be highly effective in suppressing the spurious zero-energy mode oscillations in all numerical examples while enabling efficient simulations of strain localizations across material interfaces.

Non-Ordinary State-Based Peridynamics↗

The application of finite element techniques to acoustic transmission in lined ducts with flow

The finite element method (FEM) is used to analyze the propagation of sound in two-dimensional nonuniform ducts carrying a compressible subsonic mean flow. Galerkin and residual least squares (RLS) methods with natural and forced boundary conditions are considered. The accuracy of FEM results for the eigenvalue and transmission problems is assessed by comparison with alternative numerical schemes for nonuniform ducts. The results presented and those from associated investigations indicate that modal coupling is a significant feature of the acoustic field, especially at high Mach numbers. A multimodal model therefore appears to be essential if any reliable conclusions are to be drawn in the context of turbofan inlet regions. Improvements to the eigenvalue scheme following the implementation of higher-order Hermitian elements indicate a similar modification for the transmission problem.

Astley, R. J.↗

Variability of ethane on Jupiter

Varying stratospheric temperature profiles and C2H6 altitude distributions furnish contexts for the evaluation of ethane abundances and distributions in the Jupiter stratosphere. Substantial ethane line emission and retrieved mole fraction variability is noted near the footprint of Io's flux tube, as well as within the auroral regions. It is suggested that this and other observed phenomena are due to the modification of local stratospheric chemistry by higher-order effects, which are in turn speculated to be due to the precipitation of charged particles along magnetic field lines.

Kostiuk, Theodor↗

Lessons from an evolving rRNA: 16S and 23S rRNA structures from a comparative perspective

The 16S and 23S rRNA higher-order structures inferred from comparative analysis are now quite refined. The models presented here differ from their immediate predecessors only in minor detail. Thus, it is safe to assert that all of the standard secondary-structure elements in (prokaryotic) rRNAs have been identified, with approximately 90% of the individual base pairs in each molecule having independent comparative support, and that at least some of the tertiary interactions have been revealed. It is interesting to compare the rRNAs in this respect with tRNA, whose higher-order structure is known in detail from its crystal structure (36) (Table 2). It can be seen that rRNAs have as great a fraction of their sequence in established secondary-structure elements as does tRNA. However, the fact that the former show a much lower fraction of identified tertiary interactions and a greater fraction of unpaired nucleotides than the latter implies that many of the rRNA tertiary interactions remain to be located. (Alternatively, the ribosome might involve protein-rRNA rather than intramolecular rRNA interactions to stabilize three-dimensional structure.) Experimental studies on rRNA are consistent to a first approximation with the structures proposed here, confirming the basic assumption of comparative analysis, i.e., that bases whose compositions strictly covary are physically interacting. In the exhaustive study of Moazed et al. (45) on protection of the bases in the small-subunit rRNA against chemical modification, the vast majority of bases inferred to pair by covariation are found to be protected from chemical modification, both in isolated small-subunit rRNA and in the 30S subunit. The majority of the tertiary interactions are reflected in the chemical protection data as well (45). On the other hand, many of the bases not shown as paired in Fig. 1 are accessible to chemical attack (45). However, in this case a sizeable fraction of them are also protected against chemical modification (in the isolated rRNA), which suggests that considerable higher-order structure remains to be found (although all of it may not involve base-base interactions and so may not be detectable by comparative analysis). The agreement between the higher-order structure of the small-subunit rRNA and protection against chemical modification is not perfect, however; some bases shown to covary canonically are accessible to chemical modification (45).(ABSTRACT TRUNCATED AT 400 WORDS).

Review↗

Status of 22 GeV CEBAF Upgrade Using FFA Arcs

Jefferson Lab is exploring an upgrade to extend CEBAF’s energy to ~22 GeV within the existing tunnel by incorporating Fixed-Field Alternating-Gradient (FFA) arcs. In this scheme, the two lowest-energy electromagnetic arcs are removed, the remaining arcs are reassigned to lower energies, and a new pair of high-energy FFA arcs enables six additional recirculations using the existing SRF linacs. The non-scaling FFA lattice employs Halbach-derived permanent magnets with provisions for dipole, quadrupole, and higher-order components, offering large momentum acceptance and reduced operating cost. Implementing this concept requires updated linac optics, modifications to the remaining arcs, and a redesigned electromagnetic switchyard. Ongoing beam-dynamics studies address synchrotron-radiation–induced energy loss and emittance growth, along with strategies for their mitigation.

Ogur, S. [Thomas Jefferson National Accelerator Fa↗

Structural Reliability Analysis and Optimization: Use of Approximations

This report is intended for the demonstration of function approximation concepts and their applicability in reliability analysis and design. Particularly, approximations in the calculation of the safety index, failure probability and structural optimization (modification of design variables) are developed. With this scope in mind, extensive details on probability theory are avoided. Definitions relevant to the stated objectives have been taken from standard text books. The idea of function approximations is to minimize the repetitive use of computationally intensive calculations by replacing them with simpler closed-form equations, which could be nonlinear. Typically, the approximations provide good accuracy around the points where they are constructed, and they need to be periodically updated to extend their utility. There are approximations in calculating the failure probability of a limit state function. The first one, which is most commonly discussed, is how the limit state is approximated at the design point. Most of the time this could be a first-order Taylor series expansion, also known as the First Order Reliability Method (FORM), or a second-order Taylor series expansion (paraboloid), also known as the Second Order Reliability Method (SORM). From the computational procedure point of view, this step comes after the design point identification; however, the order of approximation for the probability of failure calculation is discussed first, and it is denoted by either FORM or SORM. The other approximation of interest is how the design point, or the most probable failure point (MPP), is identified. For iteratively finding this point, again the limit state is approximated. The accuracy and efficiency of the approximations make the search process quite practical for analysis intensive approaches such as the finite element methods; therefore, the crux of this research is to develop excellent approximations for MPP identification and also different approximations including the higher-order reliability methods (HORM) for representing the failure surface. This report is divided into several parts to emphasize different segments of the structural reliability analysis and design. Broadly, it consists of mathematical foundations, methods and applications. Chapter I discusses the fundamental definitions of the probability theory, which are mostly available in standard text books. Probability density function descriptions relevant to this work are addressed. In Chapter 2, the concept and utility of function approximation are discussed for a general application in engineering analysis. Various forms of function representations and the latest developments in nonlinear adaptive approximations are presented with comparison studies. Research work accomplished in reliability analysis is presented in Chapter 3. First, the definition of safety index and most probable point of failure are introduced. Efficient ways of computing the safety index with a fewer number of iterations is emphasized. In chapter 4, the probability of failure prediction is presented using first-order, second-order and higher-order methods. System reliability methods are discussed in chapter 5. Chapter 6 presents optimization techniques for the modification and redistribution of structural sizes for improving the structural reliability. The report also contains several appendices on probability parameters.

Grandhi, Ramana V.↗

A higher-order shear deformation theory of laminated elastic shells

A higher-order shear deformation theory of elastic shells is developed for shells laminated of orthotropic layers. The theory is a modification of the Sanders' theory and accounts for parabolic distribution of the transverse shear strains through thickness of the shell and tangential stress-free boundary conditions on the boundary surfaces of the shell. The Navier-type exact solutions for bending and natural vibration are presented for cylindrical and spherical shells under simply supported boundary conditions.

Reddy, J. N.↗

CFD Predictions of Sonic-Boom Characteristics for Unmodified and Modified SR-71 Configurations

Shaped sonic-boom signatures refer to signatures that look something other than the typical N-waves. Shaped sonic-boom signatures such as "flat-top," "ramp-type," or "hybrid-type" waveforms have been shown to reduce the subjective loudness without requiring reductions in overpressure peaks. The shaping of sonic-boom signatures requires increasing the shock rise time and changes in frequency spectra. So far, a flat-top waveform was shown to be achievable in wind tunnels; however, the influence of long propagation distance and real atmosphere on shaped signatures should be addressed using flight tests. Two different approaches have been proposed for sonic-boom minimization flight tests. The first approach, proposed by Eagle Aerospace, is for a flight test using a modified BQM-34 "FIREBEE" remotely piloted vehicle. The 30-foot long FIREBEE has a steady state flight condition at the Mach number and altitude of interest, and it can be recovered by helicopter from the water. As an alternative approach, a modified SR-71 vehicle has been proposed by the McDonnell Douglas Corporation. Benefits of the SR-71 include its variable geometry supersonic inlets, small cockpit bulge, higher Mach number capabilities, slender design, and longer length (105 foot). The present investigation addresses the sonic-boom analysis for the second vehicle.The objective of the current investigation is to assess the feasibility of a modified SR-71 configuration, with McDonnell Douglas-designed fuselage modifications, intended to produce shaped sonic-boom signatures on the ground. The present study describes the use of a higher-order computational fluid dynamics (CFD) method to predict the sonic-boom characteristics for both unmodified and modified SR-71 configurations. An Euler unstructured grid methodology is used to predict the near-field, three-dimensional pressure patterns generated by both SR-71 models. The computed near-field pressure signatures are extrapolated to specified distances below the aircraft down to impingement on the ground using the code MDBOOM. Comparisons of the near-field pressure signatures with available flight-test data are presented in the current paper.

Fouladi, Kamran↗

Analysis of Factors Affecting the Performance of RLV Thrust Cell Liners

The reusable launch vehicle (RLV) thrust cell liner, or thrust chamber, is a critical component of the Space Shuttle Main Engine (SSME). It is designed to operate in some of the most severe conditions seen in engineering practice. This requirement, in conjunction with experimentally observed 'dog-house' failure modes characterized by bulging and thinning of the cooling channel wall, provides the motivation to study the factors that influence RLV thrust cell liner performance. Factors or parameters believed to be directly related to the observed characteristic deformation modes leading to failure under in-service loading conditions are identified, and subsequently investigated using the cylindrical version of the higher-order theory for functionally graded materials in conjunction with the Robinson's unified viscoplasticity theory and the power-law creep model for modeling the response of the liner s constituents. Configurations are analyzed in which specific modifications in cooling channel wall thickness or constituent materials are made to determine the influence of these parameters on the deformations resulting in the observed failure modes in the outer walls of the cooling channel. The application of thermal barrier coatings and functional grading are also investigated within this context. Comparison of the higher-order theory results based on the Robinson and power-law creep model predictions has demonstrated that, using the available material parameters, the power-law creep model predicts more precisely the experimentally observed deformation leading to the 'dog-house' failure mode for multiple short cycles, while also providing much improved computational efficiency. However, for a single long cycle, both models predict virtually identical deformations. Increasing the power-law creep model coefficients produces appreciable deformations after just one long cycle that would normally be obtained after multiple cycles, thereby enhancing the efficiency of the analysis. This provides a basis for the development of an accelerated modeling procedure to further characterize dog-house deformation modes in RLV thrust cell liners. Additionally, the results presented herein have demonstrated that the mechanism responsible for deformation leading to 'dog-house' failure modes is driven by pressure, creep/relaxation and geometric effects.

Arnold, Steven M.↗

An Overview of the Molten Salt Thermal Properties Database--Thermophysical, Version 3.1 (MSTDB-TP v.3.1)

This report presents the current status of the Molten Salt Thermal Properties Database–Thermophysical (MSTDB-TP). Information regarding version 3.1 is provided herein, which contains 820 individual salt entries (data from 180+ independent studies); the thermophysical properties contained in the database include density, viscosity, thermal conductivity, and heat capacity. The major updates to the database include a significant expansion of pseudobinary and higher-order chloride salt mixtures, many of which bearing actinides, and an incorporation of more recent literature data (i.e., that within the past 5 years). Also, modifications have been made to the pure compound data in the database as a consequence of an external quality assessment of duplicate datasets. The user-facing API for the MSTDB-TP, Saline, has been updated to include viscosity estimation capabilities based on the Redlich–Kister formalism; this is an advancement with respect to the existing density estimation capabilities. The graphical user interface was also updated to include a density estimation capability, backed by Saline. Finally, additional preliminary efforts to include surface tension into the database, as well as an investigation on formalisms that would be appropriate for thermal conductivity estimation, are reported herein.

36 MATERIALS SCIENCE↗

Improved solution for potential flow about arbitrary axisymmetric bodies by the use of a higher-order surface source method. Part 1. Theory and results

The surface-source method of calculating potential flow is improved by refining the underlying numerical analysis. The analysis uses parabolic elements and linearly-varying source density. The result is a large increase in computing speed and accuracy. The theory is described, and the effectiveness of the modification is illustrated.

Hess, J. L.↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

Mapping strain and structural heterogeneities around bubbles in amorphous ionically conductive Bi 2 O 3

While amorphous materials are often approximated to have a statistically homogeneous atomic structure, they frequently exhibit localized structural heterogeneity that challenges simplified models. This study uses 4D scanning transmission electron microscopy to investigate the strain and structural modifications around gas bubbles in amorphous Bi 2 O 3 induced by argon irradiation. We present a method for determining strain fields surrounding bubbles that can be used to measure the internal pressure of the gas. Compressive strain is observed around the cavities, with higher-order crystalline symmetries emerging near the cavity interfaces, suggesting paracrystalline ordering as a result of bubble coarsening. This ordering, along with a compressive strain gradient, indicates that gas bubbles induce significant localized changes in atomic packing. By analyzing strain fields with maximum compressive strains of 3%, we estimate a lower bound on the internal pressure of the bubbles at 2.5 GPa. These findings provide insight into the complex structural behavior of amorphous materials under stress, particularly in systems with gas inclusions, and offer new methods for probing the local atomic structure in disordered materials. Although considering structural heterogeneity in amorphous systems is non-trivial, these features have crucial impacts on material functionalities, such as mechanical strength, ionic conductivity, and electronic mobility.

36 MATERIALS SCIENCE↗

Juncture Flow Computations using kL-Based Turbulence Models

The development and implementation of kL-based Reynolds-Average Navier-Stokes (RANS) two-equation turbulence models are reported in this paper. The kL model is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta showed that a reliable transport equation can be formed from the turbulence length scale L, and the turbulence kinetic energy k. Rotta's kL equation is well suited for a term-by-term modeling and displays useful features compared to other scale formulations. One of the important differences is the inclusion of higher-order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The k-kL scheme has been modified to include the Algebraic Reynolds Stress Model (ARSM) and the Quadratic Constitutive Relation (QCR) as nonlinear models. These models have gone through extensive validations using two-dimensional benchmark problems. The present study documents the application of the k-kL schemes to simulate flow around the Juncture Flow Model (JFM). The k-kL prediction results show generally good comparisons with measurements. The results from this formulation are similar to, or better than results using the SA-RC-QCR2000 two-equation turbulence model. The ARSM and the QCR formulations of k-kL show promise with a similar level of computational resources as basic two-equation turbulence models.

K S Abdol-Hamid↗

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗

Low and moderate x gluon contribution to exclusive Compton scattering processes

We revisit the high energy semi-classical description of the exclusive processes DVCS, TCS, and Double DVCS by explicitly keeping track of the Feynman x dependence in both the hard and the hadronic matrix elements. This is achieved by a modification of the standard shock wave approximation to derive the effective Feynman rules, which leads to a generic expression on which we then perform a partial twist expansion to get rid of quantities suppressed by the proper physical scales. We obtain a compact factorized master formula that can be used to investigate the Bjorken limit at leading twist. In particular, we recover the full one-loop result in the collinear limit for pure gluon exchange with the target. Finally, we discuss the subtleties in taking the simultaneous collinear and small x limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Flavor fragmentation function factorization

A definition of partonic jet flavor that is both theoretically well-defined and experimentally robust would have profound implications for measurements and predictions especially for heavy flavor applications. Recently, a definition of jet flavor was introduced as the net flavor flowing along the direction of the Winner-Take-All axis of a jet which is soft safe to all orders, but not collinear safe. Here, we exploit the lack of collinear safety and propose a factorization theorem of perturbative flavor fragmentation functions that resum collinear divergences and describe the evolution of flavor from the short distance of jet production to the long distance at which hadronization occurs. Collinear flavor evolution is governed by a small modification of the DGLAP equations. We present a detailed all-orders analysis and identify exact relations that must hold amongst the various anomalous dimensions by probability conservation and the existence of fixed points of the renormalization group flow. We explicitly validate the factorization theorem at one-loop order, and demonstrate its consistency at two loops in particular flavor channels. Starting at two-loops, constraints on phase space imposed by flavor measurements potentially allow for non-trivial soft contributions, but we demonstrate that they are scaleless and so explicitly vanish, ensuring that soft particles are summed inclusively and all divergences are exclusively collinear in nature. This factorization theorem opens the door to precision calculations with identified flavor in the infrared.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Flavor Fragmentation Function Factorization

A definition of partonic jet flavor that is both theoretically well-defined and experimentally robust would have profound implications for measurements and predictions especially for heavy flavor applications. Recently, a definition of jet flavor was introduced as the net flavor flowing along the direction of the Winner-Take-All axis of a jet which is soft safe to all orders, but not collinear safe. Here, we exploit the lack of collinear safety and propose a factorization theorem of perturbative flavor fragmentation functions that resum collinear divergences and describe the evolution of flavor from the short distance of jet production to the long distance at which hadronization occurs. Collinear flavor evolution is governed by a small modification of the DGLAP equations. We present a detailed all-orders analysis and identify exact relations that must hold amongst the various anomalous dimensions by probability conservation and the existence of fixed points of the renormalization group flow. We explicitly validate the factorization theorem at one-loop order, and demonstrate its consistency at two loops in particular flavor channels. Starting at two-loops, constraints on phase space imposed by flavor measurements potentially allow for non-trivial soft contributions, but we demonstrate that they are scaleless and so explicitly vanish, ensuring that soft particles are summed inclusively and all divergences are exclusively collinear in nature. This factorization theorem opens the door to precision calculations with identified flavor in the infrared.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗