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 271 records · Page 15

Leveraging design of experiments to build chemometric models for the quantification of uranium (VI) and HNO3 by Raman spectroscopy

Partial least squares regression (PLSR) and support vector regression (SVR) models were optimized for the quantification of U(VI) (10–320 g L −1 ) and HNO 3 (0.6–6 M) by Raman spectroscopy with optimized calibration sets chosen by optimal design of experiments. The designed approach effectively minimized the number of samples in the calibration set for PLSR and SVR by selecting sample concentrations with a quadratic process model, despite complex confounding and covarying spectral features in the spectra. The top PLS2 model resulted in percent root mean square errors of prediction for U(VI), HNO 3 , and NO 3 − of 3.7%, 3.6%, and 2.9%, respectively. PLS1 models performed similarly despite modeling an analyte with a majority linear response (i.e., uranyl symmetric stretch) and another with more covarying vibrational modes (i.e., HNO 3 ). Partial least squares (PLS) model loadings and regression coefficients were evaluated to better understand the relationship between weaker Raman bands and covarying spectral features. Support vector machine models outperformed PLS1 models, resulting in percent root mean square error of prediction values for U(VI) and HNO 3 of 1.5% and 3.1%, respectively. The optimal nonlinear SVR model was trained using a similar number of samples (11) compared with the PLSR model, even though PLS is a linear modeling approach. The generic D-optimal design presented in this work provides a robust statistical framework for selecting training set samples in disparate two-factor systems. This approach reinforces Raman spectroscopy for the quantification of species relevant to the nuclear fuel cycle and provides a robust chemometric modeling approach to bolster online monitoring in challenging process environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cold dark energy constraints from the abundance of galaxy clusters

We constrain cold dark energy of negligible sound speed using galaxy cluster abundance observations. In contrast to standard quasi-homogeneous dark energy, negligible sound speed implies clustering of the dark energy fluid at all scales, allowing us to measure the effects of dark energy perturbations at cluster scales. We compare those models and set the stage for using non-linear information from semi-analytical modelling in cluster growth data analyses. For this, we recalibrate the halo mass function with non-linear characteristic quantities, the spherical collapse threshold and virial overdensity, that account for model and redshift-dependent behaviours, as well as an additional mass contribution for cold dark energy. Here in this paper, we present the first constraints from this cold dark matter plus cold dark energy mass function using our cluster abundance likelihood, which self-consistently accounts for selection effects, covariances and systematic uncertainties. We combine cluster growth data with cosmic microwave background, supernovae Ia and baryon acoustic oscillation data, and find a shift between cold versus quasi-homogeneous dark energy of up to 1σ. We make a Fisher matrix forecast of constraints attainable with cluster growth data from the ongoing Dark Energy Survey (DES). For DES, we predict ~ 50 percent tighter constraints on (Ωm, w) for cold dark energy versus wCDM models, with the same free parameters. Overall, we show that cluster abundance analyses are sensitive to cold dark energy, an alternative, viable model that should be routinely investigated alongside the standard dark energy scenario.

79 ASTRONOMY AND ASTROPHYSICS↗

On a family of nonoscillatory equations y double prime = phi(x)y

The oscillation or nonoscillation of a class of second-order linear differential equations is investigated analytically, with a focus on cases in which the functions phi(x) and y are complex-valued. Two linear transformations are introduced, and an asymptotic-decomposition procedure involving Shur triangularization is applied. The relationship of the present analysis to the nonoscillation criterion of Kneser (1896) and other more recent results is explored in two examples.

Gingold, H.↗

Small body mass determination from a flyby

This paper investigates the accuracy with which the mass of a small body can be determined from a flyby. Due to the nearly linear nature of hyperbolic flybys of small bodies, analytic solutions can be derived. Four analytic solutions are developed, each of increasing complexity with different applications. These solutions are compared to a complex computer simulation for accuracy. When these methods are applied to nearly linear trajectories, the analytic solutions agree with the computer simulation results to within 4.0 percent except in one boundary case which is discussed in the text. Also described is the linear range for which these equations are applicable. These solutions are then applied to the proposed asteroid flyby of the CRAF mission to predict the mass estimate accuracy.

Sybert, Catherine B.↗

Scalings for the Alfvén-cyclotron instability: Linear dispersion theory and hybrid particle-in-cell simulations: Scalings of Alfven-Cyclotron Instability

The Alfvén-cyclotron instability is driven by an ion temperature anisotropy such that T ⊥ /T || >1 where ⊥ and || denote directions perpendicular and parallel to a uniform background magnetic field B o , respectively. The computations presented here consider a model of a magnetized, homogeneous, collisionless plasma. Two representations of the proton velocity distribution are considered: a single bi-Maxwellian and a magnetospheric-like configuration of two components, a more dense, relatively cool, isotropic component and a less dense, relatively hot, bi-Maxwellian component which drives the instability. Only wave propagation parallel to B o is considered. Here, using numerical solutions of the full kinetic linear dispersion equation, concise analytic expressions for the scaling of the dimensionless maximum instability growth rate and the corresponding dimensionless real frequency are derived as functions of three dimensionless variables: the hot proton temperature anisotropy, the relative hot proton density, and the hot proton β || . Furthermore, using one-dimensional hybrid particle-in-cell simulations of this same instability, a third relation for the scaling of the maximum amplitude of the dimensionless fluctuating magnetic field energy density is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal smoothing length scale for actuator line models of wind turbine blades based on Gaussian body force distribution: Wind energy, actuator line model

The actuator line model (ALM) is a commonly used method to represent lifting surfaces such as wind turbine blades within large-eddy simulations (LES). In the ALM, the lift and drag forces are replaced by an imposed body force that is typically smoothed over several grid points using a Gaussian kernel with some prescribed smoothing width e. To date, the choice of e has most often been based on numerical considerations related to the grid spacing used in LES. However, especially for finely resolved LES with grid spacings on the order of or smaller than the chord length of the blade, the best choice of e is not known. In this work, a theoretical approach is followed to determine the most suitable value of e, based on an analytical solution to the linearized inviscid flow response to a Gaussian force. We find that the optimal smoothing width eopt is on the order of 14%-25% of the chord length of the blade, and the center of force is located at about 13%-26% downstream of the leading edge of the blade for the cases considered. These optimal values do not depend on angle of attack and depend only weakly on the type of lifting surface. It is then shown that an even more realistic velocity field can be induced by a 2-D elliptical Gaussian lift-force kernel. Some results are also provided regarding drag force representation.

17 WIND ENERGY↗

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗

Fast numerical treatment of nonlinear wave equations by spectral methods

A method is presented that accelerates spectral methods for numerical solution of a broad class of nonlinear partial differential wave equations that are first order in time and that arise in plasma wave theory. The approach involves exact analytical treatment of the linear part of the wave evolution including growth and damping as well as dispersion. After introducing the method for general scalar and vector equations, we discuss and illustrate it in more detail in the context of the coupling of high- and low-frequency plasma wave modes, as modeled by the electrostatic and electromagnetic Zakharov equations in multiple dimensions. For computational efficiency, the method uses eigenvector decomposition, which is particularly advantageous when the wave damping is mode-dependent and anisotropic in wavenumber space. In this context, it is shown that the method can significantly speed up numerical integration relative to standard spectral or finite difference methods by allowing much longer time steps, especially in the limit in which the nonlinear Schroedinger equation applies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Photon Impact Factor for DIS at NLO: analytic result

Using the Operator Product Expansion for high-energy scattering processes, we compute the photon impact factor at next-to-leading order accuracy. We obtain an analytic expression as a linear combination of five independent conformal tensor structures.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Boussinesq approximation for Rayleigh-Taylor and Richtmyer-Meshkov instabilities

We apply numerical and analytic techniques to study the Boussinesq approximation in Rayleigh-Taylor and Richtmyer-Meshkov instabilities. In this approximation, one sets the Atwood number A equal to zero except where it multiplies the acceleration g or velocity-jump Δv. While this approximation is generally applied to low-A systems, we show that it can be applied to high-A systems also in certain regimes and to the “bubble” part of the instability, i.e., the penetration depth of the lighter fluid into the heavier fluid. It cannot be applied to the spike. We extend the Boussinesq approximation for incompressible fluids and show that it always overestimates the penetration depth but the error is never more than about 41%. The effect of compressibility is studied by analytic techniques in the linear regime which indicate that compressibility has the opposite effect and the Boussinesq approximation underestimates bubbles by about 14%. We also present direct numerical simulations of two compressible systems which have approximately the same AΔv: a low-A air/CO2 system shocked at Ms = 1.57, and a high-A air/SF6 system shocked at Ms = 1.24. While the bubbles are approximately equal, the lower-A system has a shorter (less penetrating) spike; however, because its mushrooms are more tightly wound, the low-A system has the larger interface area.

Mikaelian, Karnig O.↗

Poloidal asymmetry in the narrow heat flux feature in the TCV scrape-off layer

Heat flux profiles inferred from a reciprocating probe at the outer midplane of the TCV tokamak during inner wall limited discharges feature radial fall-off lengths that shorten near the last closed flux surface (LCFS) consistent with the so-called narrow feature. The narrow feature is significantly wider on the outboard side compared with that measured on the inner wall by infrared thermography, so it is difficult to discern from the main scrape-off layer feature. After small shifts were applied for alignment, the fraction of the power contained in the narrow feature matches between inboard and outboard measurements, and they scale together with plasma current Ip, suggesting that we are observing the same phenomenon. The outboard side fall-off length within the narrow feature is found to scale closely with the radial correlation length of the edge turbulence as expected if the narrow feature arises due to radially sheared E x B flows. This is found to hold true even for cases where the narrow feature is weak and the fall-off lengths are approaching that of the far scrape-off layer. After the small shifts for alignment, non-zero floating potential profiles were found to match between inboard and outboard sides. A simple model of polarization and diamagnetic cross-field currents is described, which is consistent with the shape of these floating potential profiles. The model predicts that the floating potential at the LCFS must be negative, which supports the argument to shift the upstream measurements. The predicted currents are also consistent with the E x B flows believed to cause the narrow feature. The model is used to predict the magnitude of the floating potential of the LCFS, and the results are found to match measurements for all values of I p . This paper therefore demonstrates consistency between the measurements of the narrow feature on the inboard and outboard sides of the plasma, as well as consistency between the measurements, non-linear turbulence simulations, and analytical models of the narrow feature arising from sheared E x B flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DIRECT INTEGRATION OF THE COLLISIONLESS BOLTZMANN EQUATION IN SIX-DIMENSIONAL PHASE SPACE: SELF-GRAVITATING SYSTEMS

We present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov-Poisson equations in six-dimensional phase space. Using the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations, including the stability of King spheres, the gravitational instability, and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64{sup 6} grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of non-linearities on magnetic field generation

Magnetic fields are present on all scales in the Universe. While we understand the processes which amplify the fields fairly well, we do not have a ''natural'' mechanism to generate the small initial seed fields. By using fully relativistic cosmological perturbation theory and going beyond the usual confines of linear theory we show analytically how magnetic fields are generated. This is the first analytical calculation of the magnetic field at second order, using gauge-invariant cosmological perturbation theory, and including all the source terms. To this end, we have rederived the full set of governing equations independently. Our results suggest that magnetic fields of the order of 10{sup -30}- 10{sup -27} G can be generated (although this depends on the small scale cut-off of the integral), which is largely in agreement with previous results that relied upon numerical calculations. These fields are likely too small to act as the primordial seed fields for dynamo mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗

Elastically mediated interactions between grain boundaries and precipitates in two-phase coherent solids

In this work, we investigate analytically and numerically the interaction between grain boundaries and second phase precipitates in two-phase coherent solids in the presence of misfit strain. Our numerical study uses amplitude equations that describe the interaction of composition and stress [R. Spatschek and A. Karma, Phys. Rev. B 81, 214201 (2010)] and free-energies corresponding to two-dimensional hexagonal and three-dimensional BCC crystal structures that exhibit isotropic and anisotropic elastic properties, respectively. We consider two experimentally motivated geometries where (i) a lamellar precipitate nucleates along a planar grain boundary that is centered inside the precipitate, and (ii) a circular precipitate nucleates inside a grain at a finite distance to an initially planar grain boundary. For the first geometry, we find that the grain boundary becomes morphologically unstable due to the combination of long-range elastic interaction between the grain boundary and compositional domain boundaries, and shear-coupled grain boundary motion. We characterize this instability analytically by extending the linear stability analysis carried out recently [P.-A. Geslin, Y.-C. Xu, and A. Karma, Phys. Rev. Lett. 114, 105501 (2015)] to the more general case of elastic anisotropy. The analysis predicts that elastic anisotropy hinders but does not suppress the instability. Simulations also reveal that, in a well-developed nonlinear regime, this instability can lead to the break-up of low-angle grain boundaries when the misfit strain exceeds a threshold that depends on the grain boundary misorientation. For the second geometry, simulations show that the elastic interaction between an initially planar grain boundary and an adjacent circular precipitate causes the precipitate to migrate to and anchor at the grain boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytical solutions for transient and steady state beam loading in arbitrary traveling wave accelerating structures

Analytical solutions are derived for both transient and steady state gradient distributions in the traveling wave (TW) accelerating structures with arbitrary variation of parameters over the structure length. The results of the unloaded and beam loaded cases are presented. Finally, the exact analytical shape of the rf pulse waveform was found in order to apply the transient beam loading compensation scheme during the structure filling time. The obtained theoretical formulas were cross-checked by direct numerical simulations on the CLIC main linac accelerating structure and demonstrated a good agreement. The proposed methods provide a fast and reliable tool for the initial stage of the TW structure analysis.

43 PARTICLE ACCELERATORS↗