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,009 records · Page 56

Translation and rotation equivariant normalizing flow (TRENF) for optimal cosmological analysis

ABSTRACT Our Universe is homogeneous and isotropic, and its perturbations obey translation and rotation symmetry. In this work, we develop translation and rotation equivariant normalizing flow (TRENF), a generative normalizing flow (NF) model which explicitly incorporates these symmetries, defining the data likelihood via a sequence of Fourier space-based convolutions and pixel-wise non-linear transforms. TRENF gives direct access to the high dimensional data likelihood p(x|y) as a function of the labels y, such as cosmological parameters. In contrast to traditional analyses based on summary statistics, the NF approach has no loss of information since it preserves the full dimensionality of the data. On Gaussian random fields, the TRENF likelihood agrees well with the analytical expression and saturates the Fisher information content in the labels y. On non-linear cosmological overdensity fields from N-body simulations, TRENF leads to significant improvements in constraining power over the standard power spectrum summary statistic. TRENF is also a generative model of the data, and we show that TRENF samples agree well with the N-body simulations it trained on, and that the inverse mapping of the data agrees well with a Gaussian white noise both visually and on various summary statistics: when this is perfectly achieved the resulting p(x|y) likelihood analysis becomes optimal. Finally, we develop a generalization of this model that can handle effects that break the symmetry of the data, such as the survey mask, which enables likelihood analysis on data without periodic boundaries.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistics of errors in fibre communication lines with a phase-modulation format and optical phase conjugation

Analytical formulas are derived to approximate the probability density functions of 'zero' and 'one' bits in a linear communication channel with a binary format of optical signal phase modulation. Direct numerical simulation of the propagation of optical pulses in a communication line with optical phase conjugation is performed. The results of the numerical simulation are in good agreement with the analytical approximation. (fibreoptic communication lines)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical evaluation of the three-point scalar-tensor cross-correlations and the tensor bi-spectrum

Utilizing the Maldacena formalism and extending the earlier efforts to compute the scalar bi-spectrum, we construct a numerical procedure to evaluate the three-point scalar-tensor cross-correlations as well as the tensor bi-spectrum in single field inflationary models involving the canonical scalar field. We illustrate the accuracy of the adopted procedure by comparing the numerical results with the analytical results that can be obtained in the simpler cases of power law and slow roll inflation. We also carry out such a comparison in the case of the Starobinsky model described by a linear potential with a sudden change in the slope, which provides a non-trivial and interesting (but, nevertheless, analytically tractable) scenario involving a brief period of deviation from slow roll. We then utilize the code we have developed to evaluate the three-point correlation functions of interest (and the corresponding non-Gaussianity parameters that we introduce) for an arbitrary triangular configuration of the wavenumbers in three different classes of inflationary models which lead to features in the scalar power spectrum, as have been recently considered by the Planck team. We also discuss the contributions to the three-point functions during preheating in inflationary models with a quadratic minimum. We conclude with a summary of the main results we have obtained.

79 ASTRONOMY AND ASTROPHYSICS↗

Electrophoretic mobility of cells in a vertical Ficoll gradiant

The upward migration of living cells and test particles under the influence of a constant electric field in a low conductivity Ficoll gradient occurs at nearly constant velocity. Viscosity and neutral polymer concentration affect migration rate. Decreasing viscosity speeds up the particle migration, decreasing neutral polymer (Ficoll) concentration, slows particle migration, since electrophoretic mobility increases approximately linearly with neutral polymer concentration. Neutral polymers interact with the cell surface to effectively raise its zeta potential. An analytic function was developed from the known dependence of these physical variables on migration distance; the analysis expresses migration velocity as an explicit function of position in the density gradient. It predicts an almost linear increase in velocity of about 12 to 16% over the working region of the gradient. It was numerically integrated and correctly predicts cell migration distance vs time curves without the use of any fitted parameters. The resulting migration curves follow the expected slowly varying exponential form that closely resembles a straight line. The ability to determine standard electrophoretic mobilities from such curves depends on knowledge of the effect of Ficoll on the zeta potential of the cell type that is separated.

Plank, L. D.↗

Rapid distortion analysis of high speed homogeneous turbulence subject to periodic shear

The effect of unsteady shear forcing on small perturbation growth in compressible flow is investigated. In particular, flow-thermodynamic field interaction and the resulting effect on the phase-lag between applied shear and Reynolds stress are examined. Simplified linear analysis of the perturbation pressure equation reveals crucial differences between steady and unsteady shear effects. The analytical findings are validated with numerical simulations of inviscid rapid distortion theory (RDT) equations. In contrast to steadily sheared compressible flows, perturbations in the unsteady (periodic) forcing case do not experience an asymptotic growth phase. Further, the resonance growth phenomenon found in incompressible unsteady shear turbulence is absent in the compressible case. Overall, the stabilizing influence of both unsteadiness and compressibility is compounded leading to suppression of all small perturbations. As a result, the underlying mechanisms are explained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Surface analysis insight note: Illustrating the effect of adventitious contamination on Pt photoemission peak intensities

Adventitious carbon contaminations are not only omnipresent and used for charge referencing of XPS spectra but also can alter the apparent presence of the element peaks that span over the large spectral window of binding energies. This Insight note describes the effect of an adventitious contamination layer on Pt and presents, in brief, the approach whereby the component spectra are derived for ion beam cleaned Pt samples that can then utilize linear mathematics to peak fit said spectra thus quantifying the amount of each component including that assigned to the contamination itself of Pt metal.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sensitivity analysis for aeroacoustic and aeroelastic design of turbomachinery blades

A new method for computing the effect that small changes in the airfoil shape and cascade geometry have on the aeroacoustic and aeroelastic behavior of turbomachinery cascades is presented. The nonlinear unsteady flow is assumed to be composed of a nonlinear steady flow plus a small perturbation unsteady flow that is harmonic in time. First, the full potential equation is used to describe the behavior of the nonlinear mean (steady) flow through a two-dimensional cascade. The small disturbance unsteady flow through the cascade is described by the linearized Euler equations. Using rapid distortion theory, the unsteady velocity is split into a rotational part that contains the vorticity and an irrotational part described by a scalar potential. The unsteady vorticity transport is described analytically in terms of the drift and stream functions computed from the steady flow. Hence, the solution of the linearized Euler equations may be reduced to a single inhomogeneous equation for the unsteady potential. The steady flow and small disturbance unsteady flow equations are discretized using bilinear quadrilateral isoparametric finite elements. The nonlinear mean flow solution and streamline computational grid are computed simultaneously using Newton iteration. At each step of the Newton iteration, LU decomposition is used to solve the resulting set of linear equations. The unsteady flow problem is linear, and is also solved using LU decomposition. Next, a sensitivity analysis is performed to determine the effect small changes in cascade and airfoil geometry have on the mean and unsteady flow fields. The sensitivity analysis makes use of the nominal steady and unsteady flow LU decompositions so that no additional matrices need to be factored. Hence, the present method is computationally very efficient. To demonstrate how the sensitivity analysis may be used to redesign cascades, a compressor is redesigned for improved aeroelastic stability and two different fan exit guide vanes are redesigned for reduced downstream radiated noise. In addition, a framework detailing how the two-dimensional version of the method may be used to redesign three-dimensional geometries is presented.

Lorence, Christopher B.↗

Analytic in a Sector Resolving Families of Operators for Degenerate Evolution Fractional Equations

We introduce a class of pairs of operators defining a linear homogeneous degenerate evolution fractional differential equation in a Banach space. Reflexive Banach spaces are represented as the direct sums of the phase space of the equation and the kernel of the operator at the fractional derivative. In a sector of the complex plane containing the positive half-axis, we construct an analytic family of resolving operators that degenerate only on the kernel. The results are used in the study of the solvability of initial-boundary value problems for partial differential equations containing fractional time-derivatives and polynomials in the Laplace operator with respect to the spatial variable.

97 MATHEMATICS AND COMPUTING↗

An optimally weighted estimator of the linear power spectrum disentangling the growth of density perturbations across galaxy surveys

Measuring the clustering of galaxies from surveys allows us to estimate the power spectrum of matter density fluctuations, thus constraining cosmological models. This requires careful modelling of observational effects to avoid misinterpretation of data. In particular, signals coming from different distances encode information from different epochs. This is known as ''light-cone effect'' and is going to have a higher impact as upcoming galaxy surveys probe larger redshift ranges. Generalising the method by Feldman, Kaiser and Peacock (1994) [1], I define a minimum-variance estimator of the linear power spectrum at a fixed time, properly taking into account the light-cone effect. An analytic expression for the estimator is provided, and that is consistent with the findings of previous works in the literature. I test the method within the context of the Halofit model, assuming Planck 2014 cosmological parameters [2]. I show that the estimator presented recovers the fiducial linear power spectrum at present time within 5% accuracy up to k ∼ 0.80 h Mpc{sup −1} and within 10% up to k ∼ 0.94 h Mpc{sup −1}, well into the non-linear regime of the growth of density perturbations. As such, the method could be useful in the analysis of the data from future large-scale surveys, like Euclid.

79 ASTRONOMY AND ASTROPHYSICS↗

Variational design principles for nonequilibrium colloidal assembly

Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. Furthermore, by studying 21 unique clusters made of six, seven, or eight particles, we uncover basic design principles for targeted assembly out of equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the application of a linearized unsteady potential-flow analysis to fan-tip cascades

A linearized potential flow analysis, which accounts for the effects of nonuniform steady flow phenomena on the unsteady response to prescribed blade motions, has been applied to five two-dimensional cascade configurations. These include a flat-plate cascade and three cascades which are representative of the tip sections of current fan designs. Here the blades are closely spaced, highly staggered, and operate at low mean incidence. The fifth configuration is a NASA Lewis cascade of symmetric biconvex airfoils for which experimental measurements are available. Numerical solutions are presented that clearly illustrate the effects and importance of blade geometry and mean blade loading on the linearized unsteady response at high subsonic inlet Mach number and high blade-vibrational frequency. In addition, a good qualitative agreement is shown between the analytical predictions and experimental measurements for the cascade of symmetric biconvex airfoils. Finally, recommendations on the research needed to extend the range of application of linearized unsteady aerodynamic analyses are provided.

Usab, W. J., Jr.↗

Stability study of aqueous foams under high-temperature and high-pressure conditions relevant to Enhanced Geothermal Systems (EGS)

In Enhanced Geothermal System (EGS) hydraulic fracturing is carried out by injecting cold water into deep Hot Dry Rocks (HDR) under carefully controlled conditions to create new or reopen existing fractures. Water-based fracturing fluids demonstrate some challenges including immense quantity of water usage, water sensitivity of the formations, water blocking, and lack of proppant carrying capacity and transportation. Thus, an alternative is to use foam-based fracturing fluid which offers potential advantage over conventional water-based fracturing fluid such as minimum water usage, reduced wellbore damage, high proppant carrying capacity, and less environmental damage. However, foams are complex mixture of gaseous phase and liquid phase which are thermodynamically unstable at downhole conditions, and their stability can decrease over time due to liquid drainage, bubble coarsening, and coalescence. This paper shows laboratory experiments executed to study foam stability at high temperature (200°C) and high pressure (6.9MPa) conditions which simulates the geothermal environment. Foam stability was characterized by half-life of foam, which is defined as the time taken by the foam to decreases by 50% of its original height due to drainage. In this paper, two types of gaseous phases, nitrogen (N 2 ) and carbon dioxide (CO 2 ) were investigated. Also, based on successful practice of foam-based fracturing fluid in oil and gas industries, four surfactants, including Alpha olefin sulfonate (AOS), Sodium dodecyl sulfonate (SDS), Tergitol™ (NP – 40), and Cetyltrimethylammonium chloride (CTAC) at optimum concentration of 1wt.% were tested for best stability performance. In addition, different stabilizing agents including guar gum, bentonite clay, crosslinker, silicon dioxide nanoparticles (SiO 2 ), graphene oxide (GO) were also studied. Experimental results showed that N 2 foams were more stable than CO 2 foams. It was observed that foam half-life decreased with the increase in temperature. Among all the surfactants, AOS foams showed the most promising thermal stability at high temperatures. Moreover, with the addition of stabilizing agents, foam's half-life was enhanced. Stabilizing agents such as crosslinker and GO dispersion showed the most stable foams with half-life recorded at 20min and 17min, respectively, at 200°C and 6.9MPa. Finally, pressure also showed a positive effect on foam stability; with increased pressure, foam half-life was increased. Based on the experimental data, analytical models for the effect of temperature and pressure were developed, considering foam degradation is a first-order kinetic reaction that linearly depends on the foam drainage mechanism. The effect of temperature on foam half-life was studied as an exponential decay model. In this model, foam half-life is a function of drainage rate constant (D A ) and activation energy (E a ) of the foam system. The effect of pressure on foam half-life was found to obey a power-law model where an increase in pressure showed an increase in foam half-life. Furthermore, a linear relation was studied for the effect of pressure on foam activation energy and drainage rate. Then the combined effects of temperature and pressure were studied, which yielded an analytical model to predict the foam stabilities in terms of half-life for different foam compositions. In conclusion, this research indicates that with an appropriate selection of surfactants and stabilizing agents, it is possible to obtain stable foams, which could replace conventional water fracturing fluid under EGS conditions.

15 GEOTHERMAL ENERGY↗

Consistent Fluid Property Evaluation in RELAP5-3D

Previous investigations have developed and tested a methodology to improve the consistency between used to determine the fluid property derivatives used in the numerical solution scheme and the solution scheme itself. Within a bounding rectangle of the fluid property table for a fluid state point, linear interpolants for both fluid property values and their derivatives are calculated. A weighted average between these interpolants the analytical approximations of the fluid properties and derivatives is formed with a weighting factor between 0.0 and 1.0. This average has shown potential to reduce mass error relating to the fluid properties in RELAP5-3D. This linear interpolation and weighted averaging methodology has been implemented in RELAP5-3D with two additional options that adjust the weighting factor in the average according to the location within the bounding rectangle. The first calculates the weighting factor based only on the position between the left and right internal energy endpoints. The second also factors the position between the upper and lower pressure grid points into the calculation of the weighting factor. The new coding passes verification testing and expanded testing of the various options. Graphs are presented from three kinds of test cases, a two volume basic model, Edward’s O’Brien blowdown, and typical PWR. Some trends of the effect of crossing fluid property table pressure and internal energy grid line are identified and discussed for the simple and Edwards test cases. Preliminary results show that it the consistent fluid methodology can reduce mass error in simple cases, but in more complex input models, there is no clear best.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

TH-C-BRD-02: Analytical Modeling and Dose Calculation Method for Asymmetric Proton Pencil Beams

Purpose: A dynamic collimation system (DCS), which consists of two pairs of orthogonal trimmer blades driven by linear motors has been proposed to decrease the lateral penumbra in pencil beam scanning proton therapy. The DCS reduces lateral penumbra by intercepting the proton pencil beam near the lateral boundary of the target in the beam's eye view. The resultant trimmed pencil beams are asymmetric and laterally shifted, and therefore existing pencil beam dose calculation algorithms are not capable of trimmed beam dose calculations. This work develops a method to model and compute dose from trimmed pencil beams when using the DCS. Methods: MCNPX simulations were used to determine the dose distributions expected from various trimmer configurations using the DCS. Using these data, the lateral distribution for individual beamlets was modeled with a 2D asymmetric Gaussian function. The integral depth dose (IDD) of each configuration was also modeled by combining the IDD of an untrimmed pencil beam with a linear correction factor. The convolution of these two terms, along with the Highland approximation to account for lateral growth of the beam along the depth direction, allows a trimmed pencil beam dose distribution to be analytically generated. The algorithm was validated by computing dose for a single energy layer 5×5 cm{sup 2} treatment field, defined by the trimmers, using both the proposed method and MCNPX beamlets. Results: The Gaussian modeled asymmetric lateral profiles along the principal axes match the MCNPX data very well (R{sup 2}≥0.95 at the depth of the Bragg peak). For the 5×5 cm{sup 2} treatment plan created with both the modeled and MCNPX pencil beams, the passing rate of the 3D gamma test was 98% using a standard threshold of 3%/3 mm. Conclusion: An analytical method capable of accurately computing asymmetric pencil beam dose when using the DCS has been developed.

07 ISOTOPE AND RADIATION SOURCES↗

Non‐Linear Optical Properties of the ( RE ) 3 CuGeS 7 Family of Compounds

Abstract Non‐linear optical materials must possess a balanced combination of laser‐induced damage threshold (LDT) and second‐harmonic generation (SHG) and be phase matchable. In our previous work, chiral and polar La 3 CuGeS 7 was identified as a promising non‐linear optical material. Herein, we report the optimization of non‐linear optical properties through replacement of La with smaller lanthanides. It is determined that Gd 3 CuGeS 7 exhibits the best combination of SHG (1.6× AgGaS 2 at 88–105 μm particle size) and LDT (3× AgGaS 2 , 89 MW/cm 2 ) and is phase matchable. Based on changes in metal‐sulfur bond lengths and angles, we further propose structural optimization through solid‐solution formation and doping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LINEAR THEORY OF WEAKLY AMPLIFIED, PARALLEL PROPAGATING, TRANSVERSE TEMPERATURE-ANISOTROPY INSTABILITIES IN MAGNETIZED THERMAL PLASMAS

A rigorous analytical study of the dispersion relations of weakly amplified transverse fluctuations with wavevectors ( k-vector {parallel} B-vector ) parallel to the uniform background magnetic field B-vector in an anisotropic bi-Maxwellian magnetized electron-proton plasma is presented. A general analytical instability condition is derived that holds for different values of the electron (A{sub e}) and proton (A{sub p}) temperature anisotropies. We determine the conditions for which the weakly amplified left-handed (LH) polarized Alfven proton cyclotron and right-handed (RH) polarized Alfven Whistler electron cyclotron branches can be excited. For different regimes of the electron plasma frequency phase speed w = {omega}{sub p,e}/(kc) these branches reduce to the RH- and LH-polarized Alfven waves, RH-polarized high and low phase speed Whistler, RH-polarized proton, and LH-polarized electron cyclotron modes. Analytic instability threshold conditions are derived in terms of the combined temperature anisotropy A = T{sub perpendicular}/T{sub {parallel}}, the parallel plasma beta {beta}{sub {parallel}} = 8 {pi} n{sub e}k{sub B}T{sub {parallel}}/B {sup 2} and the electron plasma frequency phase speed w = {omega}{sub p,e}/(kc) for each mode. The results of our instability study are applied to the observed solar wind magnetic turbulence at values of 90 {<=} w {<=} 330. According to the existence conditions of the different instabilities, only the LH- and RH-polarized Alfven wave instabilities can operate here. Besides the electron-proton mass ratio {mu} = 1836, the Alfvenic instability threshold conditions are controlled by the single observed plasma parameter w. The Alfvenic instability diagram explains well the observed confinement limits at small parallel plasma beta values in the solar wind.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytic theory for the determination of velocity and stability of bubbles in a Hele-Shaw cell. Part 1: Velocity selection

An asymptotic theory is presented for the determination of velocity and linear stability of a steady symmetric bubble in a Hele-Shaw cell for small surface tension. In the first part, the bubble velocity U relative to the fluid velocity at infinity is determined for small surface tension T by determining transcendentally small correction to the asymptotic series solution. It is found that for any relative bubble velocity U in the interval (U(c),2), solutions exist at a countably infinite set of values of T (which has zero as its limit point) corresponding to the different branches of bubble solutions. U(c) decreases monotonically from 2 to 1 as the bubble area increases from 0 to infinity. However, for a bubble of arbitrarily given size, as T approaches 0, solution exists on any given branch with relative bubble velocity U satisfying the relation 2-U = cT to the 2/3 power, where c depends on the branch but is independent of the bubble area. The analytical evidence further suggests that there are no solutions for U greater than 2. These results are in agreement with earlier analytical results for a finger. In Part 2, an analytic theory is presented for the determination of the linear stability of the bubble in the limit of zero surface tension. Only the solution branch corresponding to the largest possible U for given surface tension is found to be stable, while all the others are unstable, in accordance with earlier numerical results.

Tanveer, Saleh↗

Two-level system coupled to phonons: Full analytical solution

In this study, we propose an analytical procedure to fully solve a two-level system coupled to phonons. Instead of using the common formulation in terms of linear and quadratic system-phonon couplings, we introduce different phonons depending on the system electronic level. We use this approach to recover known results for the linear-coupling limit in a simple way. More importantly, we derive results for the quadratic coupling induced by a phonon frequency change, a problem considered up to now as not analytically solvable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗