Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear expansion”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Roses in the nonperturbative current response of artificial crystals

In two-dimensional artificial crystals with large real-space periodicity, the nonlinear current response to a large applied electric field can feature a strong angular dependence, which encodes information about the band dispersion and Berry curvature of isolated electronic Bloch minibands. Within the relaxation-time approximation, we obtain analytic expressions up to infinite order in the driving field for the current in a band-projected theory with time-reversal and trigonal symmetry. For a fixed field strength, the dependence of the current on the direction of the applied field is given by rose curves whose petal structure is symmetry constrained and is obtained from an expansion in real-space translation vectors. We illustrate our theory with calculations on periodically buckled graphene and twisted double bilayer graphene, wherein the discussed physics can be accessed at experimentally relevant field strengths.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗
DOE OSTI2020

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗
DOE OSTI2026

Dynamic evolution of ultracold electron beams from nanostructured cathodes

Dynamic Coulomb expansion of dense, nonuniform particle beamlets can lead to causticlike formations in the charge density profile of the resultant beam. Such features only occur for cold beams where thermal diffusion of the particles is minimized. These phenomena have been observed experimentally for select cases with ion beams but not yet for electron beams as the electron dynamics occur over far faster timescales. We simulate the interaction of electron bunches while varying the initial transverse temperature and density profiles to determine the thresholds that characterize this pattern formation. Furthermore, we extend the simulations to geometric arrangements that are more pertinent to the nanotip cathode application and consider the effects of asymmetries and the impact of a low-density halo on the overall process. The simulations are conducted using a previously documented high-fidelity algorithm for collisional particle dynamics. These studies identify the temperature and density thresholds required for caustic formation and provide guidance for optimizing the design of nanoscale structured cathodes. We also demonstrate that the correlation between the final density patterns and the underlying beam parameters enables a virtual diagnostic tool for measuring ultralow emittance beams solely based on a transverse charge density profile.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Power counting to saturation

We present a description of saturation in small 𝑥 deep inelastic scattering from power counting in a top-down effective theory derived from quantum chromodynamics (QCD). A factorization formula isolates the universal physics of the nucleus at leading power in 𝑥. The onset of saturation is then understood as a breakdown in the expansion in an emergent power counting parameter, which is defined by the matrix element of a gauge invariant operator. We identify a new radiation mode, which enables us to extend previous literature by distinguishing the appearance of the saturation scale from the transition to nonlinear evolution.

Astronomy & Astrophysics↗
NASA NTRS1972

Nonlinear axisymmetric flexural vibration of spherical shells

Axisymmetric responses are presented of a nonshallow thin-walled spherical shell on the basis of nonlinear bending theory. An ordinary differential equation with nonlinearity of quadratic as well as cubic terms associated with variable time is derived. The derivation is based on the assumption that the deflection mode is the sum of four Legendre polynomials, and the Galerkin procedure is applied. The equation is solved by asymptotic expansion, and a first approximate solution is adopted. Unstable regions of this solution are discussed.

Kunieda, H.↗
NASA NTRS1973

Noniterative estimation of a nonlinear parameter

An algorithm is described which solves the parameters X = (x1,x2,...,xm) and p in an approximation problem Ax nearly equal to y(p), where the parameter p occurs nonlinearly in y. Instead of linearization methods, which require an approximate value of p to be supplied as a priori information, and which may lead to the finding of local minima, the proposed algorithm finds the global minimum by permitting the use of series expansions of arbitrary order, exploiting an a priori knowledge that the addition of a particular function, corresponding to a new column in A, will not improve the goodness of the approximation.

Bergstroem, A.↗
NASA NTRS1991

Computational study of three dimensional viscous flow through a turbine cascade using a multi-domain spectral technique

The three dimensional viscous flow through a planar turbine cascade is numerically simulated by direct solution of the incompressible Navier-Stokes equations. Flow dependence in the spanwise direction is represented by direct expansion in Chebyshev polynomials, while the discretization on planes parallel to the endwalls is accomplished using the spectral element method. Elemental mapping from the physical to the computational space uses an algebraic mapping technique. A fractional time stepping method that consists of an explicit nonlinear convective step, an implicit pressure correction step, and an implicit viscous step is used to advance the Navier-Stokes equations forward in time. Results computed at moderate Reynolds numbers show a three dimensional endwall flow separation, a midspan separation of the blade suction surface boundary layer, and other three-dimensional features such as the presence of a saddle point flow in the endwall region. In addition, the computed skin friction lines are shown to be orthogonal to the surface vorticity lines, demonstrating the accuracy achievable in the present method.

Renaud, Earl W.↗
NASA NTRS1972

Nonlinear behavior of shells of revolution under cyclic loading

A large deflection elastic-plastic analysis is presented, applicable to orthotropic axisymmetric plates and shells of revolution subjected to monotonic and cyclic loading conditions. The analysis is based on the finite-element method. It employs a new higher order, fully compatible, doubly curved orthotropic shell-of-revolution element using cubic Hermitian expansions for both meridional and normal displacements. Both perfectly plastic and strain hardening behavior are considered. Strain hardening is incorporated through use of the Prager-Ziegler kinematic hardening theory, which predicts an ideal Bauschinger effect. Numerous sample problems involving monotonic and cyclic loading conditions are analyzed. The monotonic results are compared with other theoretical solutions.

Levine, H. S.↗
NASA NTRS1973

Nonlinear behavior of shells of revolution under cyclic loading.

A large deflection elastic-plastic analysis is presented applicable to orthotropic axisymmetric plates and shells of revolution subjected to monotonic and cyclic loading conditions. The analysis is based on the finite-element method. It employs a new higher order, fully compatible, doubly curved orthotropic shell-of-revolution element using cubic Hermitian expansions for both meridional and normal displacements. Both perfectly plastic and strain hardening behavior are considered. Strain hardening is incorporated through use of the Prager-Ziegler kinematic hardening theory, which predicts an ideal Bauschinger effect. Numerous sample problems involving monotonic and cyclic loading conditions are analyzed.

Levine, H. S.↗
NASA NTRS1976

Nonlinear motion analysis of flexible satellites

The article reviews research on pointing accuracies of flexible satellites with regard to derived equations of motion of individual elements and composite bodies. The solution of these equations may be achieved through the formulation of a set of singularly perturbed equations by eliminating the generalized position vector. These equations are generally solved by a method of asymptotic expansions if certain convergence criteria are satisfied. The analysis may be applied to dual-spin flexible satellites.

Huang, T. C.↗
DOE OSTI2025

Acoustic misalignment mechanism for axion dark matter

A rotation in the field space of a complex scalar field corresponds to a Bose-Einstein condensation of U(1) charges. We point out that fluctuations in this rotating condensate exhibit sound-wave modes, which can be excited by cosmic perturbations and identified with axion fluctuations once the U(1) charge condensate has been sufficiently diluted by cosmic expansion. We consider the possibility that these axion fluctuations constitute dark matter and develop a formalism to compute its abundance. We carefully account for the growth of fluctuations during the epoch where the complex scalar field rotates on the body of the potential and possible nonlinear evolution when the fluctuations become non-relativistic. We find that the resultant dark matter abundance can exceed the conventional and kinetic misalignment contributions if the radial direction of the complex scalar field is sufficiently heavy. The axion dark matter may also be warm enough to leave imprints on structure formation. We discuss the implications of this novel dark matter production mechanism — acoustic misalignment mechanism — for the axion rotation cosmology, including kination domination and baryogenesis from axion rotation, as well as for axion searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1973

A study of selected radiation and propagation problems related to antennas and probes in magneto-ionic media

Research consisted of computations toward the solution of the problem of the current distribution on a cylindrical antenna in a magnetoplasma. The case of an antenna parallel to the applied magnetic field was investigated. A systematic method of asymptotic expansion was found which simplifies the solution in the general case by giving the field of a dipole even at relatively short range. Some useful properties of the dispersion surfaces in a lossy medium have also been found. A laboratory experiment was directed toward evaluating nonlinear effects, such as those due to power level, bias voltage and electron heating. The problem of reflection and transmission of waves in an electron heated plasma was treated theoretically. The profile inversion problem has been pursued. Some results are very encouraging, however, the general question of stability of the solution remains unsolved.

Source record↗
NASA NTRS1982

A second-order Budkyo-type parameterization of landsurface hydrology

A simple, second order parameterization of the water fluxes at a land surface for use as the appropriate boundary condition in general circulation models of the global atmosphere was developed. The derived parameterization incorporates the high nonlinearities in the relationship between the near surface soil moisture and the evaporation, runoff and percolation fluxes. Based on the one dimensional statistical dynamic derivation of the annual water balance, it makes the transition to short term prediction of the moisture fluxes, through a Taylor expansion around the average annual soil moisture. A comparison of the suggested parameterization is made with other existing techniques and available measurements. A thermodynamic coupling is applied in order to obtain estimations of the surface ground temperature.

Andreou, S. A.↗
NASA NTRS2015

Kodiak: An Implementation Framework for Branch and Bound Algorithms

Recursive branch and bound algorithms are often used to refine and isolate solutions to several classes of global optimization problems. A rigorous computation framework for the solution of systems of equations and inequalities involving nonlinear real arithmetic over hyper-rectangular variable and parameter domains is presented. It is derived from a generic branch and bound algorithm that has been formally verified, and utilizes self-validating enclosure methods, namely interval arithmetic and, for polynomials and rational functions, Bernstein expansion. Since bounds computed by these enclosure methods are sound, this approach may be used reliably in software verification tools. Advantage is taken of the partial derivatives of the constraint functions involved in the system, firstly to reduce the branching factor by the use of bisection heuristics and secondly to permit the computation of bifurcation sets for systems of ordinary differential equations. The associated software development, Kodiak, is presented, along with examples of three different branch and bound problem types it implements.

Smith, Andrew P.↗
DOE OSTI2024

Accelerating uncertainty quantification in incremental dynamic analysis using dimension reduction-based surrogate modeling

We propose a surrogate modeling framework based on dimension reduction to facilitate the quantification of seismic risk of structural systems in performance-based earthquake engineering. The framework adopts incremental dynamic analysis (IDA) for addressing hazard variability, and promotes significant computational efficiency improvement for propagating epistemic uncertainties associated with the structural models. It utilizes both linear and nonlinear dimension reduction approaches, equipped with inverse mappings, to learn a functional between the input parameter space (e.g., the epistemic uncertainties of the structure) to the high-dimensional output space created through the IDA implementation across different ground motions and seismic intensity levels. Polynomial chaos expansion is adopted as the surrogate model to learn this functional in the reduced space. A nine-story steel moment-resisting frame with uncertain structural properties is used as a testbed. Furthermore, we select the seismic fragility curves as a measure of the structure’s seismic performance, since it provides an estimate of the probability of entering specified damage states for given levels of ground shaking.

42 ENGINEERING↗
DOE OSTI2024

Polyconvex neural network models of thermoelasticity

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. Here, in this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

Averaged initial Cartesian coordinates for long lifetime satellite studies

A set of initial Cartesian coordinates, which are free of ambiguities and resonance singularities, is developed to study satellite mission requirements and dispersions over long lifetimes. The method outlined herein possesses two distinct advantages over most other averaging procedures. First, the averaging is carried out numerically using Gaussian quadratures, thus avoiding tedious expansions and the resulting resonances for critical inclinations, etc. Secondly, by using the initial rectangular Cartesian coordinates, conventional, existing acceleration perturbation routines can be absorbed into the program without further modifications, thus making the method easily adaptable to the addition of new perturbation effects. The averaged nonlinear differential equations are integrated by means of a Runge Kutta method. A typical step size of several orbits permits rapid integration of long lifetime orbits in a short computing time.

Pines, S.↗