Results for “MACROSCOPIC EQUATION” — Page 6

NASA NTRS2015

Polarized Bidirectional Reflectance of Optically Thick Sparse Particulate Layers: an Efficient Numerically Exact Radiative-Transfer Solution

We describe a simple yet efficient numerical algorithm for computing polarized bidirectional reflectance of an optically thick (semi-infinite), macroscopically flat layer composed of statistically isotropic and mirror symmetric random particles. The spatial distribution of the particles is assumed to be sparse, random, and statistically uniform. The 44 Stokes reflection matrix is calculated by iterating the Ambartsumian's vector nonlinear integral equation. The result is a numerically exact solution of the vector radiative transfer equation and as such fully satisfies the energy conservation law and the fundamental reciprocity relation. Since this technique bypasses the computation of the internal radiation field, it is very fast and highly accurate. The FORTRAN implementation of the technique is publicly available on the World Wide Web at http://www.giss.nasa.gov/staff/ mmishchenko/brf. It can be combined with several existing computer programs providing the requisite single-scattering properties of spherical or morphologically complex particles and applied to a wide range of optical characterization problems. Benchmark results obtained with this program can be used for testing alternative solvers of the vector radiative transfer equation.

radiative transfer↗

NASA NTRS1987

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗

NASA NTRS1983

A Potential Function Derivation of a Constitutive Equation for Inelastic Material Response

Physical and thermodynamic concepts are used to develop a potential function for application to high temperature polycrystalline material response. Inherent in the formulation is a differential relationship between the potential function and constitutive equation in terms of the state variables. Integration of the differential relationship produces a state variable evolution equation that requires specification of the initial value of the state variable and its time derivative. It is shown that the initial loading rate, which is directly related to the initial hardening rate, can significantly influence subsequent material response. This effect is consistent with observed material behavior on the macroscopic and microscopic levels, and may explain the wide scatter in response often found in creep testing.

Stouffer, D. C.↗

NASA NTRS1989

Nonlinear mesomechanics of composites with periodic microstructure

This work is concerned with modeling the mechanical deformation or constitutive behavior of composites comprised of a periodic microstructure under small displacement conditions at elevated temperature. A mesomechanics approach is adopted which relates the microimechanical behavior of the heterogeneous composite with its in-service macroscopic behavior. Two different methods, one based on a Fourier series approach and the other on a Green's function approach, are used in modeling the micromechanical behavior of the composite material. Although the constitutive formulations are based on a micromechanical approach, it should be stressed that the resulting equations are volume averaged to produce overall effective constitutive relations which relate the bulk, volume averaged, stress increment to the bulk, volume averaged, strain increment. As such, they are macromodels which can be used directly in nonlinear finite element programs such as MARC, ANSYS and ABAQUS or in boundary element programs such as BEST3D. In developing the volume averaged or efective macromodels from the micromechanical models, both approaches will require the evaluation of volume integrals containing the spatially varying strain distributions throughout the composite material. By assuming that the strain distributions are spatially constant within each constituent phase-or within a given subvolume within each constituent phase-of the composite material, the volume integrals can be obtained in closed form. This simplified micromodel can then be volume averaged to obtain an effective macromodel suitable for use in the MARC, ANSYS and ABAQUS nonlinear finite element programs via user constitutive subroutines such as HYPELA and CMUSER. This effective macromodel can be used in a nonlinear finite element structural analysis to obtain the strain-temperature history at those points in the structure where thermomechanical cracking and damage are expected to occur, the so called damage critical points of the structure.

Walker, Kevin P.↗

NASA NTRS2013

Numerical Predictions of Damage and Failure in Carbon Fiber Reinforced Laminates Using a Thermodynamically-Based Work Potential Theory

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, referred to as enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Consistent characteristic lengths are introduced into the formulation to govern the evolution of the failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs are derived. The theory is implemented into a commercial finite element code. The model is verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared against the experimental results.

composite materials↗

NASA NTRS1979

Can a macroscopic gyroscope feel torsion

We demonstrate that for a large class of Lagrangian-based torsion theories a macroscopic gyroscope is insensitive to the torsion field: there can be no coupling of the torsion to the gyroscope's angular momentum of rotation. To detect torsion a polarized system with a net elementary particle spin is needed. These conclusions are evident from the conservation laws, which form the basis for deriving the equations of motion.

Stoeger, W. R.↗

NASA NTRS1988

Bistable direction switching in an off-axis pumped continuous wave ruby laser

A report is presented of the observation of hysteretic bistable direction switching in a single-mode CW ruby laser system. This effect is only observed when the pump beam which is focused into the ruby rod is misaligned with respect to the rod end faces. At low pump powers, the ruby lases in a mode nearly collinear with the pump axis. At a higher pump power the ruby switches to a mode that is collinear with the rod end faces and preserves the original polarization. The effect is large enough to switch the beam by an angle equal to twice the diffraction angle. The observations show that under steady-state pumping, a CW ruby laser can exhibit bistable operation in its output direction and power. A calculation using the heat equation with two concentric cylinders with one as a heat source (pump laser) and the outer wall of the other held at 77 K, gives an increase in core temperature of about 0.01 K. Therefore, the increase in temperature is not large enough to change the index of refraction to account for such large macroscopic effects.

Afzal, R. Sohrab↗

NASA NTRS1991

Solutions to bi-Maxwellian transport equations for radial solar wind beyond 28 R(S)

This paper presents solar wind solutions for radial flow between 28 solar radii and 1 AU using the bi-Maxwellian-based 16-moment set of transport equations. In addition to the number density, drift velocity, and parallel and perpendicular temperatures, the 16-moment equations account for the transport of both longitudinal and transverse thermal energies as well as stress. Also, using the 16-moment approximation for the distribution function and assuming plasma parameter values characteristic of the solar wind, contour plots are generated for the proton velocity distribution function. It is shown how the shape of these plots depends on various macroscopic plasma parameters.

Demars, H. G.↗

NASA NTRS1982

A nonlocal theory of an electrostatic sinusoidal density drift instability

The stability of space plasmas in which the macroscopic variables vary in the direction x perpendicular to the magnetic field is usually studied by means of the local approximation, in which the fluctuating potential is assumed to be independent of x. To remove this approximation, a nonlocal problem is studied in which the background density, and hence the fluctuating potentials, are periodic in x. From the linear Vlasov/Poisson equations, a set of coupled, linear, homogeneous, algebraic equations relating the Fourier amplitudes of the eigensolution is derived. The equations are solved numerically and a hierarchy of exact eigenmodes characterized by different growth rates and spatial structures is found. At points where the density gradient is locally zero, the mode amplitude is generally several orders of magnitude lower than the peak amplitude. The dependence of these solutions on the parameters which define the background density perturbation is studied and a correspondence between these nonlocal results and those obtained with the local approximation is demonstrated.

Thomsen, M. F.↗

NASA NTRS1990

Equilibrium Statistical Thermodynamics of a Many-Particle System Coupled to an External Scalar Field

The equilibrium thermodynamics of a many-particle assembly in the presence of an external scalar field is examined. Two types of scalar coupling are considered: an external field coupled to the particle density and an external scalar field coupled to the energy density. It is shown that the broken translational and rotational invariance of the system due to the external field is reflected in the macroscopic physics by loss of the usual extensivity property of the system and by means of anisotropy in the response of the system to changes in the system lengths or to the system shape. In addition, the assumptions used in local equilibrium analyses are shown to be incorrect in principle. Nonlocal effects due to the external field must be included in the determination of the equation of state. Simple model calculations for a system in an external gravitational field and an externally imposed temperature field are presented as illustrations.

Salvino, R. E.↗

NASA NTRS1990

A mathematical model of the structure and evolution of small scale discrete auroral arcs

A three dimensional fluid model which includes the dispersive effect of electron inertia is used to study the nonlinear macroscopic plasma dynamics of small scale discrete auroral arcs within the auroral acceleration zone and ionosphere. The motion of the Alfven wave source relative to the magnetospheric and ionospheric plasma forms an oblique Alfven wave which is reflected from the topside ionosphere by the negative density gradient. The superposition of the incident and reflected wave can be described by a steady state analytical solution of the model equations with the appropriate boundary conditions. This two dimensional discrete auroral arc equilibrium provides a simple explanation of auroral acceleration associated with the parallel electric field. Three dimensional fully nonlinear numerical simulations indicate that the equilibrium arc configuration evolves three dimensionally through collisionless tearing and reconnection of the current layer. The interaction of the perturbed flow and the transverse magnetic field produces complex transverse structure that may be the origin of the folds and curls observed to be associated with small scale discrete arcs.

Seyler, C. E.↗

NASA NTRS2011

Numerical Implementation of a Multiple-ISV Thermodynamically-Based Work Potential Theory for Modeling Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗

NASA NTRS2012

A Thermodynamically-Based Mesh Objective Work Potential Theory for Predicting Intralaminar Progressive Damage and Failure in Fiber-Reinforced Laminates

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Damage is considered to be the effect of any structural changes in a material that manifest as pre-peak non-linearity in the stress versus strain response. Conversely, failure is taken to be the effect of the evolution of any mechanisms that results in post-peak strain softening. It is assumed that matrix microdamage is the dominant damage mechanism in continuous fiber-reinforced polymer matrix laminates, and its evolution is controlled with a single ISV. Three additional ISVs are introduced to account for failure due to mode I transverse cracking, mode II transverse cracking, and mode I axial failure. Typically, failure evolution (i.e., post-peak strain softening) results in pathologically mesh dependent solutions within a finite element method (FEM) setting. Therefore, consistent character element lengths are introduced into the formulation of the evolution of the three failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs is derived. The theory is implemented into commercial FEM software. Objectivity of total energy dissipated during the failure process, with regards to refinements in the FEM mesh, is demonstrated. The model is also verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared to the experiments.

Pineda, Evan J.↗

NASA NTRS1991

Investigation of the structure of the electromagnetic field and related phenomena, generated by the Active Satellite

The altitude dependencies of the moduli of the electric field E in the VLF and LF frequency bands (f sub B much less than F less than f sub B) and in the altitude range of the ionosphere Z equals (400 to 2500) km up to Z equals 6000 km (the bottom of the magnetosphere) were calculated by the linear theory. The amplitudes of the field have large maxima in four regions: the axis field (E sub o) close to the direction of the Earth's magnetic field line B sub o, beta approximately 0 degrees, the fields (E sub St), (E sub RevSt) and (E sub Res) in the Storey, Reversed Story and Resonance cones, beta approximately (0 approaches 20) degrees. Their maxima are very pronounced close to the low hybrid frequency F sub L. The nonlinear heating of the magnetoplasma under the action of an electric field Ee (sup iwt) was recently expanded by the macroscopic theory by the author. The velocities, collision frequencies and temperatures of all the constituents of a magnetoplasma - electrons, ions, and neutral particles - are taken into account. Formulae and numerical results are presented for the ionosphere in the frequency band F equals (1 to 10 exp 4) kHz and altitude range Z approximately (100 - 1000) km. Some results of calculations by the self consistent solution of the basis system of equations are also discussed.

Alpert, Yakov L.↗

NASA NTRS1995

Modeling the heliolatitudinal gradient of the solar wind parameters with exact MHD solutions

The heliolatitudinal dependence of observations of the solar wind macroscopic quantities such as the averaged proton speed, density and the mass and momentum flux are modeled. The published observations covering the last two and a half solar cycles, are obtained either via the technique of interplanetary scintillations for the last 2 solar cycles (1970-1990), or, from the plasma experiment aboard the ULYSSES spacecraft for the recent period 1990-1994. Exact, two dimensional solutions of the full set of the steady MHD equations are used which are obtained through a nonlinear separation of the variables in the MHD equations. The three parameters emerging from the solutions are fixed from these observations, as well as from observations of the solar rotation. It is found that near solar maximum the solar wind speed is uniformly low, around the 400 km/s over a wide range of latitudes. On the other hand, during solar minimum and the declining phase of the solar activity cycle, there is a strong heliolatitudinal gradient in proton speed between 400-800 from equator to pole. This modeling also agrees with previous findings that the gradient in wind speed with the latitude is offset by a gradient in density such that the mass and momentum flux vary relatively little.

Lima, J. J. G.↗

NASA NTRS1996

From Individuals to Epidemics

Heterogeneous mixing fundamentally changes the dynamics of infectious diseases; finding ways to incorporate it into models represents a critical challenge. Phenomenological approaches are deficient in their lack of attention to underlying processes; individual-based models, on the other hand, may obscure the essential interactions in a sea of detail. The challenge then is to find ways to bridge these levels of description, starting from individual-based models and deriving macroscopic descriptions from them that retain essential detail, and filter out the rest. In this paper, attempts to achieve this transformation are described for a class of models where non-random mixing arises from the spatial localization of interactions. In general, the epidemic threshold is found to be larger owing to spatial localization than for a homogeneous mixing population. An improved estimate of the dynamics is developed by the use of moment equations, and a simple estimate of the threshold in terms of a 'dyad heuristic'. For more general models in which local infection is not described by mass action, the connection with related partial differential equations is investigated.

Levin, Simon A.↗

NASA NTRS1976

High-beta turbulence in two-dimensional magnetohydrodynamics

Equations of ideal magnetohydrodynamics are used to study incompressible turbulent flows in a specified geometry where all the field quantities vary with only two spatial dimensions. The procedures adopted are basically those of Kraichnan (1967), in which classical equilibrium ensembles are built around constants of the motion identified from the Fourier-transformed equations of motion. Once the constants of the motion are identified, the statistical formulation of the problem is presented in a phase space whose coordinates are the real and imaginary parts of the Fourier coefficients. Canonical ensembles are constructed in this phase space by classical arguments. The theory developed permits isolation of some qualitatively new gross physical effects which have so far not been calculated. One of the more novel physical effects is the appearance of macroscopic structures involving long-wavelength, self-generated, magnetic fields for a wide range of initial parameters.

Fyfe, D.↗

NASA NTRS1978

Theory of flux anisotropies in a guiding center plasma

The one particle distribution function f on the scale of the bounce motion of particles in a magnetic field B is considered. The Vlasov equation is expanded through O(epsilon) in the adiabatic parameter which is the ratio of particle gyroradius to scale length of the magnetic field. Because f is directly proportional to particle flux differential in kinetic energy and solid angle, f is in principle measurable in space experiments, and the analysis is tailored to be explicitly applicable to space problems. To O(1), f is gyrotropic; its first velocity moment is (if non-vanishing) parallel to B, and hence macroscopic parallel flow is included in this term. The O(epsilon) contribution is non-gyrotropic and macroscopic flow parallel to B plus additional parallel flow results from these terms. The degree of non-gyrotropy and the amount of cross-field macroscopic flow depend on the perpendicular component of the electric field, on curvature and shear in the magnetic field, and on the spatial gradient, pitch angle derivative, and speed derivative of the lowest order distribution function.

Birmingham, T. J.↗