Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic linearization”

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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

An analytical method for ditching analysis of an airborne vehicle

A simple analytical method has been introduced for aerohydrodynamic load analysis of an airborne configuration during water ditching. The method employs an aerodynamic panel code, based on linear potential flow theory, to simulate the flow of air and water around an aircraft configuration. The free surface separating the air and water region is represented by doublet sheet singularities. Although all the theoretical load distributions are computed for air, provisions are made to correct the pressure coefficients obtained on the configuration wetted surfaces to account for the water density. As an analytical tool, the Vortex Separation Aerodynamic (VSAERO) code is chosen to carry out the present investigation. After assessing the validity of the method, its first application is to analyze the water ditching of the Space Shuttle configuration at a 12 degree attitude.

Linear potential flow↗
DOE OSTI2017

A Simple Analytical Model for Rocky Planet Interiors

This work aims to explore the scaling relations among rocky exoplanets. The assumption that the internal gravity increases linearly in the core and stays constant in the mantle is tested against numerical simulations, and a simple model is constructed, applicable to rocky exoplanets of CMF (core mass fraction) = 0.2–0.35 and mass = 0.1–30 M ⊕. Various scaling relations are derived: (1) CMF ≈ CRF 2 (core radius fraction squared), (2) ${P}_{\mathrm{typical}}\sim {g}_{s}^{2}$ (typical interior pressure scales as surface gravity squared), (3) energy released in core formation is $\sim \tfrac{1}{10}$ the total gravitational energy, (4) effective heat capacity of the mantle $\approx \left(\tfrac{{M}_{p}}{{M}_{\oplus }}\right)\times 7.5\times {10}^{27}$ J K -1 , (5) moment of inertia $\approx \tfrac{1}{3}{M}_{p}{R}_{p}^{2}$. These approximate relations are handy for quick use owing to their simplicity and lucidity, and provide insights into the interior structures of those exoplanets.

planets and satellites: interiors↗
DOE OSTI2015

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Symmetric and asymmetric solitons trapped in H-shaped potentials

We report results of numerical and analytical studies of the spontaneous symmetry breaking in solitons, both two- and one-dimensional (2D and 1D), which are trapped in H-shaped potential profiles built of two parallel potential troughs linked by a narrow rung in the transverse direction. This system can be implemented in self-attractive Bose-Einstein condensates (BECs), as well as in a nonlinear bulk optical waveguide. We demonstrate that the introduction of the transverse link changes the character of the symmetry-breaking bifurcation (SBB) in the system from subcritical to supercritical. (In terms of the corresponding phase transition, it is a change between the first and second kinds.) A noteworthy feature of the SBB in this setting is a nonmonotonous dependence of the soliton's norm at the bifurcation point on the strength of the transverse link. In the full 2D system, the results are obtained in a numerical form. An exact analytical solution is found for the bifurcation in the 1D version of the model, with the transverse rung modeled by the local linear coupling between the parallel troughs with the {delta}-functional longitudinal profile. By replacing the {delta} function by its finite-width Gaussian counterpart, similar results are obtained by means of the variational approximation (VA). The VA is also applied to the 1D system with a mixed linear and nonlinear transverse localized coupling. Comparison of the results produced by the different varieties of the system clearly reveals basic features of its the symmetry-breaking transition.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Simulations of x-ray absorption spectra for CO desorbing from Ru(0001) with transition-potential and time-dependent density functional theory approaches

The desorption of a carbon monoxide molecule from a Ru(0001) surface was studied by means of X-ray Absorption Spectra (XAS) computed with Transition Potential (TP-DFT) and Time Dependent (TD-DFT) DFT methods. By unraveling the evolution of the CO electronic structure upon desorption, we observed that at 2.3 Å from the surface, the CO molecule has already predominantly gas-phase character. While C 1s XAS is quite insensitive to changes in the C–O bond length, the O 1s excitation is very sensitive with the π* coming down in energy upon CO bond stretching, which competes with the increase in orbital energy due to the repulsive interaction with the metallic surface. We show in a systematic way that the TP-DFT method can describe the XAS rather well at the endpoints (chemisorbed and gas phase) but is affected by artificial charge transfer and/or incorrect spin treatment in the transition region in cases like CO, where there are low-lying π* orbitals and large exchange interactions between the core 1s and valence-acceptor π* orbitals. As an alternative, we demonstrate by comparing with experimental data that a linear response approach using TD-DFT employing common exchange-correlation functionals and finite-size clusters can yield a good description of the spectral evolution of the 1s → π* transition with correct spin and gas-to-chemisorbed chemical shifts in good agreement with experiment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1988

Helicopter stability and control modeling improvements and verification on two helicopters

A linearized model of helicopter flight dynamics is developed which includes the flapping, lead-lag, and dynamic inflow degrees of freedom (DOF). The model is a combination of analytical terms and numerically determined stability derivatives, and is used to investigate the importance of the rotor DOF to stability and control modeling. The results show that the rotor DOF can have a significant impact on some of the natural modes in a linear model. The flap and dynamic inflow DOF show the greatest influence. Flapping exhibits strong coupling to the body, dynamic inflow, and to lead-lag to a lesser extent. Dynamic inflow tends to damp the high-frequency flapping modes, and reduces the damping on coupled body-flap motion. Dynamic inflow also couples to the flapping motion to produce complex roots. With body-flap and lag regressing modes as exceptions, the results show essentially similar behavior for most modes of articulated and hingeless rotor helicopters.

Schrage, D. P.↗
DOE OSTI2021

Closed bipolar electrode-enabled electrochromic detector for chemical sensing

The described BPE-enabled device includes two separated chambers which perform detection and reporting independently. Analytical reaction of a target molecule in the analytical cell is coupled to and monitored by an electrochromic reaction in the reporting cell. The color change in the reporting cell can be determined spectrophotometrically by RGB analysis of a CCD image acquired via smartphone. This detection method provides a linear response and a low limit of detection due to the redox cycling behavior in both chambers. The BPE based electrochromic detector can be modified for sensing of multiple analytes by integrating three or more sets of detection chemistries into one single device. Multiple analytes with different concentrations can be detected within this device simultaneously. The BPE based electrochromic device can be used for metabolite detection, wherein a redox mediator can be combined with specific oxidases to form an electrochemical mediator-electrocatalyst pair that completes redox cycling reactions.

Bohn, Paul W.↗
DOE OSTI2011

Explicit expressions for three-dimensional boundary integrals in linear elasticity

On employing isoparametric, piecewise linear shape functions over a flat triangle, exact formulae are derived for all surface potentials involved in the numerical treatment of three-dimensional singular and hyper-singular boundary integral equations in linear elasticity. These formulae are valid for an arbitrary source point in space and are represented as analytical expressions along the edges of the integration triangle. They can be employed to solve integral equations defined on triangulated surfaces via a collocation method or may be utilized as analytical expressions for the inner integrals in a Galerkin technique. A numerical example involving a unit triangle and a source point located at various distances above it, as well as sample problems solved by a collocation boundary element method for the Lame equation are included to validate the proposed formulae.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

A scalable exponential-DG approach for nonlinear conservation laws: With application to Burger and Euler equations

In this work, we propose an Exponential DG framework for partial differential equations. We decompose 7 governing equations into linear and nonlinear parts to which we apply the discontinuous Galerkin 8 (DG) spatial discretization. In particular, we construct the linear part using Jacobian that effectively 9 capture stiff characteristics in the system. The former is integrated analytically, whereas the latter 10 is approximated. This approach i) is stable with a large Courant number (Cr > 1); ii) supports 11 high-order solutions both in time and space; iii) is computationally favorable compared to IMEX 12 DG methods with no preconditioner; iv) becomes comparable to explicit RKDG methods on uniform 13 mesh and beneficial on non-uniform grid for Euler equations; v) is scalable in a modern massively 14 parallel computing architecture due to its explicit nature of exponential time integrators and com15 pact communication stencil of DG method. Numerical results demonstrate the performance of our 16 proposed methods through various examples. We also discuss the stability and convergence analysis 17 for our exponential DG scheme in the context of Burgers equation.

42 ENGINEERING↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2013

Effects of collisions on energetic particle-driven chirping bursts

In the presence of an energetic particle population in a dissipative plasma, self-trapped structures in phase-space (holes and clumps) emerge from nonlinear wave-particle interactions. Their dynamics can lead to a nonlinear continuous shifting of the wave frequency (chirping). The effects of collisions on chirping characteristics are investigated, with a one-dimensional kinetic model. Existing analytic theory is extended to account for Krook-like collisions, which quantitatively explains a significant departure from widely accepted square-root time dependency. Relaxation oscillations, associated with chirping bursts, are investigated in the presence of dynamical friction and velocity-diffusion. The period increases with decreasing drag and weakly increases with decreasing diffusion. The mechanism is clarified with a simple semi-analytic model of hole/clump pair, which satisfies a Fokker-Planck equation. The model shows that the linear growth rate cannot be obtained simply by fitting an exponential to the amplitude time-series.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Kinetic transverse dispersion relation for relativistic magnetized electron-positron plasmas with Maxwell-Jüttner velocity distribution functions

We use a kinetic treatment to study the linear transverse dispersion relation for a magnetized isotropic relativistic electron-positron plasma with finite relativistic temperature. The explicit linear dispersion relation for electromagnetic waves propagating along a constant background magnetic field is presented, including an analytical continuation to the whole complex frequency plane for the case of Maxwell-Jüttner velocity distribution functions. This dispersion relation is studied numerically for various temperatures. For left-handed solutions, the system presents two branches, the electromagnetic ordinary mode and the Alfvén mode. In the low frequency regime, the Alfvén branch has two dispersive zones, the normal zone (where ∂ω/∂k > 0) and an anomalous zone (where ∂ω/∂k < 0). We find that in the anomalous zone of the Alfvén branch, the electromagnetic waves are damped, and there is a maximum wave number for which the Alfvén branch is suppressed. We also study the dependence of the Alfvén velocity and effective plasma frequency with the temperature. We complemented the analytical and numerical approaches with relativistic full particle simulations, which consistently agree with the analytical results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Expansion of a cold non-neutral plasma slab

Expansion of the ion and electron fronts of a cold non-neutral plasma slab with a quasi-neutral core bounded by layers containing only ions is investigated analytically and exact solutions are obtained. It is found that on average, the plasma expansion time scales linearly with the initial inverse ion plasma frequency as well as the degree of charge imbalance, and no expansion occurs if the cold plasma slab is stationary and overall neutral. However, in both cases, there can exist prominent oscillations on the electron front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A numerically efficient damping model for acoustic resonances in microfluidic cavities

Bulk acoustic wave devices are typically operated in a resonant state to achieve enhanced acoustic amplitudes and high acoustofluidic forces for the manipulation of microparticles. Among other loss mechanisms related to the structural parts of acoustofluidic devices, damping in the fluidic cavity is a crucial factor that limits the attainable acoustic amplitudes. In the analytical part of this study, we quantify all relevant loss mechanisms related to the fluid inside acoustofluidic micro-devices. Subsequently, a numerical analysis of the time-harmonic visco-acoustic and thermo-visco-acoustic equations is carried out to verify the analytical results for 2D and 3D examples. The damping results are fitted into the framework of classical linear acoustics to set up a numerically efficient device model. For this purpose, all damping effects are combined into an acoustofluidic loss factor. Since some components of the acoustofluidic loss factor depend on the acoustic mode shape in the fluid cavity, we propose a two-step simulation procedure. In the first step, the loss factors are deduced from the simulated mode shape. Subsequently, a second simulation is invoked, taking all losses into account. Owing to its computational efficiency, the presented numerical device model is of great relevance for the simulation of acoustofluidic particle manipulation by means of acoustic radiation forces or acoustic streaming. For the first time, accurate 3D simulations of realistic micro-devices for the quantitative prediction of pressure amplitudes and the related acoustofluidic forces become feasible.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Study of the behaviour of the light field orbital angular momentum upon astigmatic mode conversion

We consider the behaviour of the orbital angular momentum of the light field as a superposition of Hermite – Gaussian beams upon astigmatic mode conversion. An analytical expression is obtained for the orbital angular momentum of such fields. Expressions are found for various linear combinations of Hermite – Gaussian beams. The astigmatic mode conversion of the initial light fields that are equivalent from the point of view of the orbital angular momentum is shown to lead to significantly different final results. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Multiple lensing of the cosmic microwave background anisotropies

We study the gravitational lensing effect on the Cosmic Microwave Background (CMB) anisotropies performing a ray-tracing of the primordial CMB photons through intervening large-scale structures (LSS) distribution predicted by N-Body numerical simulations with a particular focus on the precise recovery of the lens-induced polarized counterpart of the source plane. We apply both a multiple plane ray-tracing and an effective deflection approach based on the Born approximation to deflect the CMB photons trajectories through the simulated lightcone. We discuss the results obtained with both these methods together with the impact of LSS non-linear evolution on the CMB temperature and polarization power spectra. We compare our results with semi-analytical approximations implemented in Boltzmann codes like, e.g., CAMB. We show that, with our current N-body setup, the predicted lensing power is recovered with good accuracy in a wide range of multipoles while excess power with respect to semi-analytic prescriptions is observed in the lensing potential on scales ℓ ∼> 3000. We quantify the impact of the numerical effects connected to the resolution in the N-Body simulation together with the resolution and band-limit chosen to synthesise the CMB source plane. We found these quantities to be particularly important for the simulation of B-mode polarization power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗