Search NASA⌕ Search

SEARCH · Search NASA

Results for “PERTURBATION”

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 199 records · Page 11

Early Time Evolution of Circumferential Perturbation of Initial Particle Volume Fraction in Explosive Cylindrical Multiphase Dispersion

When an annular bed of solid particles that surrounds a cylindrical high-energy explosive core gets radially dispersed after detonation, the expanding front of particles undergoes instabilities. One of the possible causes of the instabilities is an inhomogeneous initial distribution of particles. This study explores this possibility by introducing two-dimensional perturbations to the initial distribution of particles within the annular bed and quantifying the growth of these perturbations over time using two-dimensional simulations. The initial perturbations are in the form of superposition of up to three sinusoidal azimuthal modal variations in the initial particle volume fraction (PVF, ratio of particle to cell volume). These are observed to impact the particle distribution at later times through a channeling instability whose effects are: (i) to decrease the velocity in regions of larger particle volume (PV) and (ii) to facilitate circumferential particle migration into the slow moving high PV sectors. These departures from axisymmetry are quantified by introducing two metrics. The effect of varying the number of azimuthal modes contained in the initial PVF perturbation, along with their amplitudes, wavelengths, and relative phases is investigated. The proposed metrics do not vary substantially with the relative phases; however, there is a strong variation in the metrics due to changes in the wavenumber. Unimodal perturbations were found to amplify both metrics the most.

Engineering↗

Monte Carlo simulations of patient dose perturbations in rotational-type radiotherapy due to a transverse magnetic field: A tomotherapy investigation

Purpose: Several groups are exploring the integration of magnetic resonance (MR) image guidance with radiotherapy to reduce tumor position uncertainty during photon radiotherapy. The therapeutic gain from reducing tumor position uncertainty using intrafraction MR imaging during radiotherapy could be partially offset if the negative effects of magnetic field-induced dose perturbations are not appreciated or accounted for. The authors hypothesize that a more rotationally symmetric modality such as helical tomotherapy will permit a systematic mediation of these dose perturbations. This investigation offers a unique look at the dose perturbations due to homogeneous transverse magnetic field during the delivery of Tomotherapy{sup ®} Treatment System plans under varying degrees of rotational beamlet symmetry. Methods: The authors accurately reproduced treatment plan beamlet and patient configurations using the Monte Carlo code GEANT4. This code has a thoroughly benchmarked electromagnetic particle transport physics package well-suited for the radiotherapy energy regime. The three approved clinical treatment plans for this study were for a prostate, head and neck, and lung treatment. The dose heterogeneity index metric was used to quantify the effect of the dose perturbations to the target volumes. Results: The authors demonstrate the ability to reproduce the clinical dose–volume histograms (DVH) to within 4% dose agreement at each DVH point for the target volumes and most planning structures, and therefore, are able to confidently examine the effects of transverse magnetic fields on the plans. The authors investigated field strengths of 0.35, 0.7, 1, 1.5, and 3 T. Changes to the dose heterogeneity index of 0.1% were seen in the prostate and head and neck case, reflecting negligible dose perturbations to the target volumes, a change from 5.5% to 20.1% was observed with the lung case. Conclusions: This study demonstrated that the effect of external magnetic fields can be mitigated by exploiting a more rotationally symmetric treatment modality.

60 APPLIED LIFE SCIENCES↗

Exploration of magnetic perturbation effects on advanced divertor configurations in NSTX-U

The control of divertor heat loads - both steady state and transient - remains a key challenge for the successful operation of ITER and FNSF. Magnetic perturbations provide a promising technique to control ELMs (transients), but understanding their detailed impact is difficult due to their symmetry breaking nature. One approach for reducing steady state heat loads are so called 'advanced divertors' which aim at optimizing the magnetic field configuration: the snowflake and the (super-)X-divertor. It is likely that both concepts - magnetic perturbations and advanced divertors - will have to work together, and we explore their interaction based on the NSTX-U setup. An overview of different divertor configurations under the impact of magnetic perturbations is presented, and the resulting impact on plasma edge transport is investigated with the EMC3-EIRENE code. Variations in size of the magnetic footprint of the perturbed separatrix are found, which is related to the level of flux expansion on the divertor target. Non-axisymmetric peaking of the heat flux related to the perturbed separatrix is found at the outer strike point, but only in locations where flux expansion is not too large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Orographic Precipitation Response to Microphysical Parameter Perturbations for Idealized Moist Nearly Neutral Flow

Abstract This study explores the sensitivity of clouds and precipitation to microphysical parameter perturbations using idealized simulations of moist, nearly neutral flow over a bell-shaped mountain. Numerous parameters are perturbed within the Morrison microphysics scheme. The parameters that most affect cloud and precipitation characteristics are the snow fall speed coefficient As, snow particle density ρs, rain accretion (WRA), and ice–cloud water collection efficiency (ECI). Surface precipitation rates are affected by As and ρs through changes to the precipitation efficiency caused by direct and indirect impacts on snow fall speed, respectively. WRA and ECI both affect the amount of cloud water removed, but the precipitation sensitivity differs. Large WRA results in increased precipitation efficiency and cloud water removal below the freezing level, indirectly decreasing cloud condensation rates; the net result is little precipitation sensitivity. Large ECI removes cloud water above the freezing level but with little influence on overall condensation rates. Two environmental experiments are performed to test the robustness of the results: 1) reduction of the wind speed profile by 30% (LowU) and 2) decreasing the surface potential temperature to induce a freezing level below the mountain top (LowFL). Parameter perturbations within LowU result in similar mechanisms acting on precipitation, but a much weaker sensitivity compared to the control. The LowFL case shows ρs is no longer a dominant parameter and As now induces changes to cloud condensation, since more of the cloud depth is present above the freezing level. In general, perturbations to microphysical parameters affect the location of peak precipitation, while the total amount of precipitation is more sensitive to environmental parameter perturbations.

54 ENVIRONMENTAL SCIENCES↗

Assessing the Influence of Microphysical and Environmental Parameter Perturbations on Orographic Precipitation

Microphysical (MP) schemes contain parameters whose values can impact the amount and location of forecasted precipitation, and sensitivity is typically explored by perturbing one parameter at a time while holding the rest constant. Although much can be learned from these “one-at-a-time” studies, the results are limited as these methods do not allow for nonlinear interactions of multiple perturbed parameters. This study applies the Morris one-at-a-time (MOAT) method, a robust statistical tool allowing for simultaneous perturbation of numerous parameters, to explore orographic precipitation sensitivity to changes in microphysical and environmental parameters within an environment characteristic of an atmospheric river. Results show parameters associated with snow fall speed coefficient A s and density ρ s , ice-cloud water collection efficiency (ECI), rain accretion (WRA), relative humidity, zonal wind speed, and surface potential temperature cause the largest influence on simulated precipitation. MP and environmental parameter perturbations can cause precipitation changes of similar magnitude, but results vary by location on the mountain. Different environments are also tested, with A s being the most influential MP parameter regardless of environment. Fewer MP parameters influence precipitation in a faster-wind-speed environment, possibly due to the stronger dynamical forcing upwind and different wave dynamics downwind, compared to a slower-wind-speed environment. Finally, perturbing MP parameters within a single scheme can result in precipitation variations of similar magnitude compared to using entirely different microphysics schemes. MOAT results presented in this study have implications for Bayesian parameter estimation methods and stochastic parameterization within ensemble forecasting.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Role of Perturbing Ocean Initial Condition in Simulated Regional Sea Level Change

Multiple lines of observational evidence indicate that the global climate has been getting warmer since the early 20th century. This warmer climate has led to a global mean sea level rise of about 18 cm during the 20th century, and over 6 cm for the first 15 years of the 21st century. Regionally the sea level rise is not uniform due in large part to internal climate variability. To better serve the community, the uncertainties of predicting/projecting regional sea level changes associated with internal climate variability need to be quantified. Previous research on this topic has used single-model large ensembles with perturbed atmospheric initial conditions (ICs). Here we compare uncertainties associated with perturbing ICs in just the atmosphere and just the ocean using a state-of-the-art coupled climate model. We find that by perturbing the oceanic ICs, the uncertainties in regional sea level changes increase compared to those with perturbed atmospheric ICs. In order for us to better assess the full spectrum of the impacts of such internal climate variability on regional and global sea level rise, approaches that involve perturbing both atmospheric and oceanic initial conditions are thus necessary.

54 ENVIRONMENTAL SCIENCES↗

Perturbation of a Schwarzschild Black Hole Due to a Rotating Thin Disk

Will, in 1974, treated the perturbation of a Schwarzschild black hole due to a slowly rotating, light, concentric thin ring by solving the perturbation equations in terms of a multipole expansion of the mass-and-rotation perturbation series. In the Schwarzschild background, his approach can be generalized to perturbation by a thin disk (which is more relevant astrophysically), but, due to rather bad convergence properties, the resulting expansions are not suitable for specific (numerical) computations. However, we show that Green’s functions, represented by Will’s result, can be expressed in closed form (without multipole expansion), which is more useful. In particular, they can be integrated out over the source (a thin disk in our case) to yield good converging series both for the gravitational potential and for the dragging angular velocity. The procedure is demonstrated, in the first perturbation order, on the simplest case of a constant-density disk, including the physical interpretation of the results in terms of a one-component perfect fluid or a two-component dust in a circular orbit about the central black hole. Free parameters are chosen in such a way that the resulting black hole has zero angular momentum but non-zero angular velocity, as it is just carried along by the dragging effect of the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Shielding application of perturbation theory to determine changes in neutron and gamma doses due to changes in shield layers

Perturbation theory formulas were derived and applied to determine changes in neutron and gamma-ray doses due to changes in various radiation shield layers for fixed sources. For a given source and detector position, the perturbation method enables dose derivatives with respect to density, or equivalently thickness, for every layer to be determined from one forward and one inhomogeneous adjoint calculation. A direct determination without the perturbation approach would require two forward calculations to evaluate the dose derivative due to a change in a single layer. Hence, the perturbation method for obtaining dose derivatives requires fewer computations for design studies of multilayer shields. For an illustrative problem, a comparison was made of the fractional change in the dose per unit change in the thickness of each shield layer in a two-layer spherical configuration as calculated by perturbation theory and by successive direct calculations; excellent agreement was obtained between the two methods.

Fieno, D.↗

The spectrum of density perturbations in an expanding universe

The basic dynamic equations that govern the evolution of perturbations in a Friedmann-Lemaitre universe are derived. General solutions describing the evolution of adiabatic perturbations in the density of matter are obtained, and the choice of the appropriate initial conditions is examined. The various perturbation modes are compared, and the effects of decoupling on the perturbation spectrum are studied. The scheme used to follow the evolution of density perturbations through decoupling is based on an extension of the Eddington approximation to the radiative transfer equation, and is strictly valid in both optically thick and thin limits.

Silk, J.↗

Perturbations of a close-earth satellite due to sunlight diffusely reflected from the earth. I - Uniform albedo

A semianalytic method has been developed to calculate the radiation-pressure perturbations of a close-earth satellite due to sunlight reflected from the earth. The assumptions made are that the satellite is spherically symmetric and that the solar radiation is reflected from the earth according to Lambert's Law with uniform albedo. By using expressions for the components of the radiation-pressure force due to Lochry, the expressions for the perturbations of the elements were developed into series in the true anomaly. The perturbations within a given revolution can be obtained analytically by integrating with respect to v while holding all slowly varying quantities constant. The long-range perturbations are then obtained by accumulating the net perturbations at the end of each revolution.

Lautman, D. A.↗

The correlation function for density perturbations in an expanding universe. I - Linear theory

The evolution of the two-point correlation function for adiabatic density perturbations in the early universe is studied. Analytical solutions are obtained for the evolution of linearized spherically symmetric adiabatic density perturbations and the two-point correlation function for these perturbations in the radiation-dominated portion of the early universe. The results are then extended to the regime after decoupling. It is found that: (1) adiabatic spherically symmetric perturbations comparable in scale with the maximum Jeans length would survive the radiation-dominated regime; (2) irregular fluctuations are smoothed out up to the scale of the maximum Jeans length in the radiation era, but regular fluctuations might survive on smaller scales; (3) in general, the only surviving structures for irregularly shaped adiabatic density perturbations of arbitrary but finite scale in the radiation regime are the size of or larger than the maximum Jeans length in that regime; (4) infinite plane waves with a wavelength smaller than the maximum Jeans length but larger than the critical dissipative damping scale could survive the radiation regime; and (5) black holes would also survive the radiation regime and might accrete sufficient mass after decoupling to nucleate the formation of galaxies.

Mcclelland, J.↗

Application of functional analysis to perturbation theory of differential equations

The deviation of the solution of the differential equation y' = f(t, y), y(O) = y sub O from the solution of the perturbed system z' = f(t, z) + g(t, z), z(O) = z sub O was investigated for the case where f and g are continuous functions on I x R sup n into R sup n, where I = (o, a) or I = (o, infinity). These functions are assumed to satisfy the Lipschitz condition in the variable z. The space Lip(I) of all such functions with suitable norms forms a Banach space. By introducing a suitable norm in the space of continuous functions C(I), introducing the problem can be reduced to an equivalent problem in terminology of operators in such spaces. A theorem on existence and uniqueness of the solution is presented by means of Banach space technique. Norm estimates on the rate of growth of such solutions are found. As a consequence, estimates of deviation of a solution due to perturbation are obtained. Continuity of the solution on the initial data and on the perturbation is established. A nonlinear perturbation of the harmonic oscillator is considered a perturbation of equations of the restricted three body problem linearized at libration point.

Bogdan, V. M.↗

Comparison of uniform perturbation solutions and numerical solutions for some potential flows past slender bodies

Approximate solutions for potential flow past an axisymmetric slender body and past a thin airfoil are calculated using a uniform perturbation method and then compared with either the exact analytical solution or the solution obtained using a purely numerical method. The perturbation method is based upon a representation of the disturbance flow as the superposition of singularities distributed entirely within the body, while the numerical (panel) method is based upon a distribution of singularities on the surface of the body. It is found that the perturbation method provides very good results for small values of the slenderness ratio and for small angles of attack. Moreover, for comparable accuracy, the perturbation method is simpler to implement, requires less computer memory, and generally uses less computation time than the panel method. In particular, the uniform perturbation method yields good resolution near the regions of the leading and trailing edges where other methods fail or require special attention.

Wong, T. C.↗

Instabilities in line-driven stellar winds. I - Dependence on perturbation wavelength

An analytical study is presented of the stability of absorption line-driven flows, such as found in stellar winds, in the presence of small-amplitude disturbances. A generalized calculation of the perturbed direct extinction force is performed and the evolution of the perturbation into the nonlinear regime and the dynamical results of perturbation in the scattered radiation field are examined. An expression is derived for the wavenumber variation of the perturbed line force from the milieu of nonoverlapping lines which have a power-law distribution in opacity. A linear dispersion analysis is carried out to model the growth and propagation of radiative-acoustic waves in absorption line-driven flows, which are found unstable to perturbations that may be smaller than the Sobolev length. No damping mechanism was found that would eliminate the absorption line-driven flows.

Owocki, S. P.↗

Comparison of uniform perturbation and numerical solutions for some potential flows past slender bodies

Approximate solutions for potential flow past an axisymmetric slender body and past a thin airfoil are calculated using a uniform perturbation method and then compared with either the exact analytical solution or the solution obtained using a purely numerical method. The perturbation method is based upon a representation of the disturbance flow as the superposition of singularities distributed entirely within the body, while the numerical (panel) method is based upon a distribution of singularities on the surface of the body. It is found that the perturbation method provides very good results for small values of the slenderness ratio and for small angles of attack. Moreover, for comparable accuracy, the perturbation method is simpler to implement, requires less computer memory, and generally uses less computation time than the panel method. In particular, the uniform perturbation method yields good resolution near the regions of the leading and trailing edges where other methods fail or require special attention.

Wong, T.-C.↗

A hybrid perturbation Galerkin technique with applications to slender body theory

A two step hybrid perturbation-Galerkin method to solve a variety of applied mathematics problems which involve a small parameter is presented. The method consists of: (1) the use of a regular or singular perturbation method to determine the asymptotic expansion of the solution in terms of the small parameter; (2) construction of an approximate solution in the form of a sum of the perturbation coefficient functions multiplied by (unknown) amplitudes (gauge functions); and (3) the use of the classical Bubnov-Galerkin method to determine these amplitudes. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is applied to some singular perturbation problems in slender body theory. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the degree of applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗

Improvements in the perturbation simulations of the global reference atmospheric model

The Global Reference Atmospheric Model (GRAM) program includes the capability for simulating pseudo-random perturbations in density, temperature, pressure, or wind components along a simulated reentry trajectory or other path through the atmosphere. Some concerns were expressed by GRAM users, however, that the mean-square perturbation gradients may be too large for small values of the vertical separation Delta z. The present GRAM perturbation simulations, based on a one-step autoregressive model, yield a power spectrum versus wavenumber k which is proportional to k sup -2 at high wavenumbers. This feature also produces mean-square perturbation differences which are directly proportional to Delta z, and mean-square perturbation gradients which are inversely proportional to Delta z. Thus, root-mean-square gradients, (Delta f/Delta z) sub rms, increase with decreasing Delta a as Delta z sup -1/2. A simple modification to GRAM is suggested which overcomes this problem, i.e., which produce root-mean-square gradient that remain bound as Delta z approaches zero. Possible applications of more sophisticated simulation approaches, e.g., second order autoregressive models, or fractal model techniques, were also explored briefly but found to yield improvements which appear too small to justify their considerable added complexity for use in the GRAM programs.

Justus, C. G.↗

A hybrid perturbation Galerkin technique with applications to slender body theory

A two-step hybrid perturbation-Galerkin method to solve a variety of applied mathematics problems which involve a small parameter is presented. The method consists of: (1) the use of a regular or singular perturbation method to determine the asymptotic expansion of the solution in terms of the small parameter; (2) construction of an approximate solution in the form of a sum of the perturbation coefficient functions multiplied by (unknown) amplitudes (gauge functions); and (3) the use of the classical Bubnov-Galerkin method to determine these amplitudes. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is applied to some singular perturbation problems in slender body theory. The results obtained from the hybrid method are compared with approximate solutions obtained by other methods, and the degree of applicability of the hybrid method to broader problem areas is discussed.

Geer, James F.↗