Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical homogenization”

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 253 records · Page 14

Influence of rotation on the maximum mass of pulsationally stable stars

The critical mass for stability against radial pulsations in rotating, homogeneous main-sequence stars is found to be greater than in the case of no rotation. Analytic and detailed numerical models show that the critical mass rises steeply with increasing concentration of angular momentum to the center of the star. For uniform rotation near breakup velocity at the star's equator the critical mass is about 850 solar masses if an electron-scattering opacity is used, or about 5000 solar masses if the opacities of Cox and Stewart are used. For nonuniform rotation with a constant ratio of centrifugal force to gravity in the star, the critical mass becomes 'infinite' long before breakup velocity is attained. The relevance of the present results to several observational problems is noted.

Stothers, R.↗

Electromagnetic diffraction efficiencies for plane reflection diffraction gratings

The theory and computer programs, based on electromagnetic theory, for the analysis and design of echelle gratings were developed. The gratings are designed for instruments that operate in the ultraviolet portion of the spectrum. The theory was developed so that the resulting computer programs will be able to analyze deep (up to 30 wavelengths) gratings by including as many as 100 real or homogeneous diffraction orders. The program calculates the complex amplitude coefficient for each of the diffracted orders. A check on the numerical method used to solve the integral equations is provided by a conservation of energy calculation.

Marathay, A. S.↗

Laser amplification of incoherent radiation

The amplification of noise in a laser amplifier is treated theoretically. The model for the active medium and its description using density-matrix techniques, are taken from the theory of laser operation. The spectral behavior of the radiation in the nonlinear regime is studied and the formalism is written from the onset in the frequency domain. The statistics of the light are gradually modified by the nonlinear amplification process, and expressions are derived for the rate of change of fluctuations in intensity as a measure of statistical changes. In addition, the range of validity of Litvak's Gaussian-statistics approximation is discussed. In the homogeneous-broadening case, the evolution of initially broadband Gaussian radiation toward quasimonochromatic oscillations with laserlike statistics is explored in several numerical examples. The connections of this study with the time-domain work on self-pulsing in a ring-laser configuration, are established. Finally, spectral-narrowing and -rebroadening effects in Doppler-broadened media are discussed both analytically and with numerical examples. These examples show the distinct contribution of pulsations in the population ('Raman-type terms'), and saturation phenomena.

Menegozzi, L. N.↗

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗

Laboratory scaled simulation of lidar cloud sounding experiments

The results of lidar measurements carried out on laboratory scale models of clouds are presented. Measurements on laboratory scale models are important since one has the knowledge of the relevant parameters of the diffusing medium, such as: scattering and absorption coefficients, phase function, homogeneity, shape, etc. Knowledge of these parameters enables one to use the results to test the reliability of theoretical and numerical investigations. To obtain a laboratory scaled model of a lidar system sounding a cloud, it is necessary to scale down all the geometrical quantities by the same factor to reduce distances of the order of kilometers to the order of meters, keeping the size and the optical depth of the diffusers unchanged. If a time resolution of the order of nanoseconds is necessary for a lidar sounding actual clouds, the corresponding time resolution for the laboratory model should be of the order of picoseconds. It is possible to obtain this resolution by using picosecond laser systems and fast electrooptical detectors like the streak camera. The results of the laboratory measurements showed that the multiple scattering effect strongly depends on the size of the diffusers, as well as on the concentration. The experimental results were compared with the numerical results of a Monte Carlo code. A generally good agreement was obtained.

Zaccanti, G.↗

Nonlinear decay of magnetic helicity in magnetohydrodynamic turbulence with a mean magnetic field

We show that the magnetic helicity associated with fluctuations in homogeneous incompressible magnetohydrodynamic (MHD) turbulence with a mean magnetic field decays in time because of nonlinear processes. Evidence is obtained numerically, by use of both dissipative and nondissipative spectral method simulations. The described effect stands in contrast to expectations based on studies of MHD turbulence without an applied mean field, in which magnetic helicity is transferred nonlinearly to long wavelengths and is preserved in time because of selective decay when dissipation is present. The process of nonlinear decay is described in terms of a generalized ideal helicity invariant and is characterized by the transient production of a mean induced electric field aligned with the applied magnetic field, an effect reminiscent of the alpha effect in dynamo theory. A simple phenomenological model for the decay process is proposed.

Stribling, Troy↗

Numerical modeling of coronal mass ejections based on various pre-event model atmospheres

We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution at lower boundary (i.e. solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) Constant boundary conditions and a polytropic gas with gamma = 1.05; (2) Non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) Constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) Non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using MHD characteristic theory. The simulations show how different CME's, including loop-like transients, clouds, and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.

Suess, S. T.↗

Numerical Modeling of Coronal Mass Ejections Based on Various Pre-event Model Atmospheres

We examine how the initial state (pre-event corona) affects the numerical MHD simulation for a coronal mass ejection (CME). Earlier simulations based on a pre-event corona with a homogeneous density and temperature distribution, at the lower boundary (i.e., solar surface) have been used to analyze the role of streamer properties in determining the characteristics of loop-like transients. The present paper extends these studies to show how a broader class of global coronal properties leads not only to different types of CME's, but also modifies the adjacent quiet corona and/or coronal holes. We consider four pre-event coronal cases: (1) constant boundary conditions and a polytropic gas with gamma = 1.05; (2) non-constant (latitude dependent) boundary conditions and a polytropic gas with gamma = 1.05; (3) constant boundary conditions with a volumetric energy source and gamma = 1.67; (4) non-constant (latitude dependent) boundary conditions with a volumetric energy source and gamma = 1.67. In all models, the pre-event magnetic fields separate the corona into closed field regions (streamers) and open field regions. The CME's initiation is simulated by introducing at the base of the corona, within the streamer region, a standard pressure pulse and velocity change. Boundary values are determined using magnetohydrodynamic (MHD) characteristic theory. The simulations show how different CME's, including loop-like transients, clouds and bright rays, might occur. There are significant new features in comparison to published results. We conclude that the pre-event corona is a crucial factor in dictating CME's properties.

Wang, A. H.↗

Broken Ergodicity in Two-Dimensional Homogeneous Magnetohydrodynamic Turbulence

Two-dimensional (2-D) homogeneous magnetohydrodynamic (MHD) turbulence has many of the same qualitative features as three-dimensional (3-D) homogeneous MHD turbulence.The se features include several ideal invariants, along with the phenomenon of broken ergodicity. Broken ergodicity appears when certain modes act like random variables with mean values that are large compared to their standard deviations, indicating a coherent structure or dynamo.Recently, the origin of broken ergodicity in 3-D MHD turbulence that is manifest in the lowest wavenumbers was explained. Here, a detailed description of the origins of broken ergodicity in 2-D MHD turbulence is presented. It will be seen that broken ergodicity in ideal 2-D MHD turbulence can be manifest in the lowest wavenumbers of a finite numerical model for certain initial conditions or in the highest wavenumbers for another set of initial conditions.T he origins of broken ergodicity in ideal 2-D homogeneous MHD turbulence are found through an eigen analysis of the covariance matrices of the modal probability density functions.It will also be shown that when the lowest wavenumber magnetic field becomes quasi-stationary, the higher wavenumber modes can propagate as Alfven waves on these almost static large-scale magnetic structures

Shebalin, John V.↗

The growth of aspherical structure in the universe - Is the Local Supercluster an unusual system

The growth and subsequent collapse of homogeneous ellipsoidal perturbations in a uniform expanding background is considered as a simple model for the formation of large-scale aspherical structures in the observed universe. Numerical calculations of the evolution of such perturbations turn out to be well described by an approximate analytic solution of the equations of motion, and simple relationships are found between the initial shape of a perturbation and its shape and kinematic properties at the time of collapse. Perturbations do not change their shape significantly until they reach a density contrast of order unity. As a result, structures with the kinematic properties of the Local Supercluster should form much more commonly in a low-density universe than in a flat universe. The homogeneity of the local Hubble flow, the motion of the Milky Way with respect to the microwave background, and the flattening of the Local Supercluster can be successfully accounted for by these models, provided that the initial perturbation is sufficiently flattened. Viable models are obtained only if the ratio of the lengths of the two smaller axes of the initial perturbation is at least 3:1 in an Einstein-de Sitter universe or at least 1.8:1 in a universe for which the density parameter (Omega) is of order 0.1, when the protocluster pancakes.

White, S. D. M.↗

Asteroids as meteorite parent-bodies - The astronomical perspective

A review of astronomical evidence suggests that asteroids are the parent-bodies for most meteorites. The surfaces of most asteroids are like carbonaceous chondrites while a significant minority are of stony-iron composition. Other meteorite types are recognized in the belt but are rare. The size distribution of stony-iron asteroids implies they are remnant cores of about 100 differentiated bodies subjected to collisional fragmentation by carbonaceous objects initially about 300 times more numerous than now. Incomplete evidence on parent-body collisions exists in data on Hirayama families, asteroid lightcurves, and the compositional homogeneity of individual asteroids. Modern-day asteroid regoliths are thin and cannot have been environments for formation of most brecciated, gas-rich meteorites; such meteorites formed during early accretion of the asteroids. A scenario for the origin and evolution of meteorite parent-bodies is presented which includes: (1) interruption of planet-formation by processes due to Jupiter; (2) substantial asteroidal collisions during the first 0.5 b.y.; and (3) formation of most meteorite types within the differentiated bodies.

Chapman, C. R.↗

The ion-ion acoustic instability

This paper considers the linear theory of ion-acoustic-like instabilities in a homogeneous Vlasov plasma with two ion components, a less dense beam and a more dense core, with a relative drift velocity. Numerical solutions of the full electrostatic dispersion equation are presented, and the properties of the ion-ion acoustic instability are studied in detail. At relatively cold beam temperatures, the instability is fluid-like, but it becomes a beam resonant kinetic instability as the beam temperature becomes of the order of the core temperature; if the mode is unstable, its threshold lies well below the threshold of the electron-ion acoustic instability. An electron temperature anisotropy enhances the instability, and at sufficiently large beam-core relative drift speeds, electron magnetization can either detract from or enhance the instability.

Gary, S. Peter↗

The Stickney impact of Phobos - A dynamical model

The hypervelocity impact that excavated the Stickney crater on Phobos is numerically modeled in order to understand the effects this collision had on the interior and surface of Phobos. The model assumes homogeneity prior to impact based on Fujiwara's (1991) results. Constitutive relations for water ice and basalt are applied in tandem models. Fracture energetics is used to show that the impact itself is sufficient to cause the cracks observed on Phobos. It is further shown that the Stickney impact itself could not have created large voids; thus, the low density of Phobos must either be compositional in nature or else the result of porosity at a scale sufficiently small to avoid scattering the impact energy substantially.

Asphaug, E.↗

Rapid distortion analysis and direct simulation of compressible homogeneous turbulence at finite Mach number

The effect of rapid mean compression on compressible turbulence at a range of turbulent Mach numbers is investigated. Rapid distortion theory (RDT) and direct numerical simulation results for the case of axial (one-dimensional) compression are used to illustrate the existence of two distinct rapid compression regimes. These regimes are set by the relationships between the timescales of the mean distortion, the turbulence, and the speed of sound. A general RDT formulation is developed and is proposed as a means of improving turbulence models for compressible flows.

Cambon, C.↗

A Review of the Ginzburg-Syrovatskii's Galactic Cosmic-Ray Propagation Model and its Leaky-Box Limit

Phenomenological models of galactic cosmic-ray propagation are based on a diffusion equation known as the Ginzburg-Syrovatskii s equation, or variants (or limits) of this equation. Its one-dimensional limit in a homogeneous volume, known as the leaky-box limit or model, is sketched here. The justification, utility, limitations, and a typical numerical implementation of the leaky-box model are examined in some detail.

Barghouty, A. F.↗

Spacecraft detumbling using movable telescoping appendages

The dynamics of detumbling a randomly spinning spacecraft using externally mounted, movable telescoping appendages were studied both analytically and numerically. Two types of telescoping appendages are considered: where an end mass is mounted at the end of an (assumed) massless boom; and where the appendage is assumed to consist of a uniformly distributed homogeneous mass throughout its length. From an application of Liapunov's second method, boom extension maneuvers were determined to approach either of two desired final states: close to a zero inertial angular velocity state, and a final spin rate about only one of the principal axes. Recovery dynamics are evaluated analytically for the case of symmetrical deployment. Numerical examination of other asymmetrical cases verifies the practicality of using movable appendages to recover a randomly tumbling spacecraft.

Bainum, P. M.↗

Thermal convection at infinite Prandtl number

Numerical solutions are developed for the full scaled three dimensional thermal convection problem, at infinite Prandtl number, and for rigid boundaries. The procedure is designed to give an approximate account of horizontally homogeneous and isotropic flow situations. Results, which included a maximum of 36 horizontal wave number vectors and 20 vertical wave numbers, appear to adequately describe the flow fill up to R = 100,000; beyond this R the results appear to show horizontal wave number truncation error. This error seems to affect the boundary slope of the mean temperature field more than other mean quantities. Despite some numerical uncertainities, certain of the qualitative features of the flow fill are predicted with reasonable confidence.

Herring, J. R.↗

Rapid Distortion Analysis and Direct Simulation of Compressible Homogeneous Turbulence at Finite Mach Number

The effect of rapid mean compression on compressible turbulence at a range of turbulent Mach numbers is investigated. Rapid dist'ortion theory (RDT) and direct numerical simulation results for the case of axial (one-dimensional) compression are used to illustrate the existence of two distinct rapid compression regimes. These regimes - the nearly solenoidal and the 'pressure-released' - are defined by a single parameter involving the timescales of the mean distortion, the turbulence, and the speed of sound. A general RDT formulation is developed and is proposed as a means of improving turbulence models for compressible flows. In contrast to the well-documented observation that 'compressibility' (measured, for example, by the turbulent Mach number) is often associated with a decrease in the growth rate of turbulent kinetic energy, we find that under rapid distortion compressibility can produce an amplification of the kinetic energy growth rate. We also find that as the compressibility increases, the magnitude of the pressure-dilation correlation increases, in absolute terms, but its relative importance decreases compared to the magnitude of the kinetic energy production.

Cambon, C.↗