On the propagation of a disturbance in a smoothly varying heterogeneous porous medium saturated with three fluid phases
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Abstract not provided.
Abstract not provided.
We recommend using the dayside Martian ionosphere as a reflector for global communication, because it has a stable density peak and usable critical frequency. This is very crucial for the future Mars ground to ground communication. The dayside ionosphere has been well modeled as a Chapman layer. We suggest performing the Martian nightside ionosphere modeling study. Because the nightside ionosphere has very little measurements available, we propose to drop a digital ionosound instrument into the Mars surface for data collection. Even though the Martian tropospheric radio refractivity had a small value, it still can cause ray bending and multipath effects. We recommend performing an accurate calculation on excess phase and group delays (range and time delays). Other effects, such as range rate errors, appearance angle deviation, defocusing loss on Mars, etc. should be estimated. Ice depolarization effects due to Martian clouds on radio waves are unknown yet, but they are expected to be small, because of lower optical depth and the thinner layer of cloud. Total Martian atmospheric gaseous attenuation is expected to be less than 1 dB on microwaves and, because the Martian atmosphere had very low concentration in uncondensed H2O and O2. An accurate calculation for zenith opacity requires the information about scale heights of H2O and O2 distribution. An accurate water vapor altitude profile at Mars is not available et. Under the normal condition, CO2 and N2 gases ddo not have electric or magnetic dipoles and do not absorb electromagnetic energy from the waves. However, they may generate the dipoles through a collision and interact with waves under a high density condition and absorb electromagnetic waves in the infrared and visible band. Dust storm is the most dominant factor on the radio wave attenuation. Large Martian dust storms can cause at least 3 dB or higher loss to Ka band wave. For a normal dust storm, the attenuation is about 1 dB. The attenuation much depends on dust mass loading, dust size distribution, etc. Most large dust storms occur in the southern hemisphere during later spring and early summer when the southern hemisphere becomes suddenly hot.
We recommend to use the dayside Martian ionosphere as a reflector for global communication, because the dayside ionosphere has stable density peak and usable critic frequency. This is very crucial for the future Mars ground to around communication. The dayside ionosphere has been well modeled as a Chapman layer. We suggest to perform the Martian nightside ionospheric modeling study. Because the nightside ionosphere has very little measurements available, we propose to drop a digital ionosond instrument into the Mars surface for data collection. Even though the Martian tropospheric radio refractivity has small value, it still can cause the ray bending and multipath effects. We recommend to perform an accurate calculation on excess phase and group delays (range and time delays). Other effects, such as range rate errors, appearance angle deviation, defocusing loss on Mars, etc. are also needed to be estimated. Ice depolarization effects due to Martian clouds on radio waves is unknown yet, which is expected to be small, because lower optical depth and thinner layer of cloud: Total Martian atmospheric gaseous attenuation is expected to be less than 1 dB on microwaves band, because the Martian atmosphere has very low concentration in uncondensed H2O and O2. An accurate calculation for zenith opacity requires the information about scale heights of H2O and O2 distribution. An accurate water vapor altitude profile at Mars is not available yet. Under the normal condition, CO2 and N2 gases do not have electric or magnetic dipoles and do not absorb electromagnetic energy from the waves. However, they may generate the dipoles through a collision and interact with waves under a high density condition and absorb electromagnetic waves in the infrared and visible band. Dust storm is most dominant factor to the radio wave attenuation. Large Martian dust storm can cause at least 3 dB or higher loss to Ka band wave. For a normal dust storm, the attenuation is about 1 dB. The attenuation much depends on dust mass loading, dust size distribution, etc. Most large dust storm occur in the southern hemisphere during later spring and early summer when the southern hemisphere become suddenly hot.
Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.
Reduced weight is a major design goal in aircraft power transmissions. Some gear designs incorporate thin rims to help meet this goal. Thin rims, however, may lead to bending fatigue cracks. These cracks may propagate through a gear tooth or into the gear rim. A crack that propagates through a tooth would probably not be catastrophic, and ample warning of a failure could be possible. On the other hand, a crack that propagates through the rim would be catastrophic. Such cracks could lead to disengagement of a rotor or propeller from an engine, loss of an aircraft, and fatalities. To help create and validate tools for the gear designer, the NASA Lewis Research Center performed in-house analytical and experimental studies to investigate the effect of rim thickness on gear-tooth crack propagation. Our goal was to determine whether cracks grew through gear teeth (benign failure mode) or through gear rims (catastrophic failure mode) for various rim thicknesses. In addition, we investigated the effect of rim thickness on crack propagation life. A finite-element-based computer program simulated gear-tooth crack propagation. The analysis used principles of linear elastic fracture mechanics, and quarter-point, triangular elements were used at the crack tip to represent the stress singularity. The program had an automated crack propagation option in which cracks were grown numerically via an automated remeshing scheme. Crack-tip stress-intensity factors were estimated to determine crack-propagation direction. Also, various fatigue crack growth models were used to estimate crack-propagation life. Experiments were performed in Lewis' Spur Gear Fatigue Rig to validate predicted crack propagation results. Gears with various backup ratios were tested to validate crack-path predictions. Also, test gears were installed with special crack-propagation gages in the tooth fillet region to measure bending-fatigue crack growth. From both predictions and tests, gears with backup ratios (rim thickness divided by tooth height) of 3.3 and 1.0 produced tooth fractures, whereas a backup ratio of 0.3 produced rim fractures. For a backup ratio of 0.5, the experiments produced rim fractures and the predictions produced both rim and tooth fractures, depending on the initial geometry of the crack. Good correlation between predicted and measured crack growth was achieved when the fatigue crack-closure concept was introduced into the analysis. As the gear rim thickness decreased, the compressive cyclic stress in the gear-tooth fillet region increased. This retarded crack growth and increased the number of crack-propagation cycles to failure.
The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.
The characterization of the rates of crack propagation by means of a Paris law is intended for steady-state fatigue conditions. While rates of crack propagation are primarily driven by 𝐽𝐽𝑚𝑚𝑚𝑚𝑚𝑚, in the presence of evolving process zones and fiber bridging they also depend on the crack propagation length: cracks with fully formed process zones propagate more slowly. Analyses of the rates of crack propagation in laminates secondarily bonded with AF163-2K film adhesive were conducted to identify the rate transients that occur during the formation of the process zone in the adhesive and during the evolution of bridging of the adhesive’s carrier knit across the crack faces. Dávila's CF22 cohesive fatigue damage model was used in finite element analyses of the double cantilever beam test to explain discrepancies in experimental rates of propagation obtained from tests conducted with either increasing or decreasing 𝐽max. The results indicate that the analyses accurately reproduce transient and steady-state rates of propagation observed in the experimental results, as well as the threshold of propagation. Transient rates of propagation always exceed those in steady state, which underscores the need to understand transients and account for them in fatigue characterization. The results also highlight that a Paris curve for steady-state propagation can be obtained most quickly using a decreasing load function, but that this curve consists of an unconservative bound to the rates of propagation. A lateral shift in the Paris curve is proposed for a conservative characterization that removes the effect of bridging and other R-curve effects.
Indirect excitons (IXs), also known as interlayer excitons, can form the medium for excitonic devices whose operation is based on controlled propagation of excitons. Here, a proof of principle for excitonic devices was demonstrated in GaAs heterostructures where the operation of excitonic devices is limited to low temperatures. IXs in van der Waals transition-metal dichalcogenide (TMD) heterostructures are characterized by high binding energies making IXs robust at room temperature and offering an opportunity to create excitonic devices operating at high temperatures suitable for applications. However, a characteristic feature of TMD heterostructures is the presence of moiré superlattice potentials, which are predicted to cause modulations of IX energy reaching tens of meV. These in-plane energy landscapes can lead to IX localization, making IX propagation fundamentally different in TMD and GaAs heterostructures and making uncertain whether long-range IX propagation, sufficiently long to allow for creating elaborate excitonic devices and circuits, can be realized in TMD heterostructures. In this work, we realize long-range IX propagation with the 1/e IX luminescence decay distances reaching 13 microns in a MoSe 2 /WSe 2 heterostructure. We trace the IX luminescence along the IX propagation path. In the presented TMD materials, the long-range IX propagation occurs up to ~50 K. This is a step toward the room temperature operation, which can be realized in TMD heterostructures due to the high IX binding energy. We also realize control of the long-range IX propagation by voltage. The IX luminescence signal in the drain of an excitonic transistor is controlled within 40 times by gate voltage. The control of the IX propagation in the MoSe 2 /WSe 2 heterostructure is governed by new mechanisms, beyond the mechanism for controlling IX transport by an energy barrier to IX propagation (or a trap for IXs) created by the gate electrode, known since the studies of GaAs heterostructures. We discuss the origin of the voltage-controlled long-range IX propagation in the MoSe 2 /WSe 2 heterostructure, in particular, the electric-field control of the moiré potential.
The NASA Propagation Program supports academic research that models various propagation phenomena in the space research frequency bands. NASA supports such research via school and institutions prominent in the field. The products of such efforts are particularly useful for researchers in the field of propagation phenomena and telecommunications systems engineers. The systems engineer usually needs a few propagation parameter values for a system design. Published literature on the subject, such as the Cunsultative Committee for International Radio (CCIR) publications, may help somewhat, but often times, the parameter values given in such publications use a particular set of conditions which may not quite include the requirements of the system design. The systems engineer must resort to programming the propagation phenomena model of interest and to obtain the parameter values to be used in the project. Furthermore, the researcher in the propagation field must then program the propagation models either to substantiate the model or to generate a new model. The researcher or the systems engineer must either be a skillful computer programmer or hire a programmer, which of course increases the cost of the effort. An increase in cost due to the inevitable programming effort may seem particularly inappropriate if the data generated by the experiment is to be used to substantiate the already well-established models, or a slight variation thereof. To help researchers and the systems engineers, it was recommended by the participants of NASA Propagation Experimenters (NAPEX) 15 held in London, Ontario, Canada on 28-29 June 1991, that propagation software should be constructed which will contain models and prediction methods of most propagation phenomenon. Moreover, the software should be flexible enough for the user to make slight changes to the models without expending a substantial effort in programming.
Analytic and numerical techniques are presented which directly address the problem of propagation of wide bandwidth signals through random ionized media. This work is applicable to the problems of satellite communication and space based radar observation through a disturbed ionospheric propagation channel that would result from a high altitude chemical release or nuclear detonation. An analytic solution is presented for the two-position, two-frequency mutual coherence function for spherical wave propagation in the strong scattering limit. This solution is used to derive simple expressions for the power impulse response function and to determine its relationship to the important parameters of decorrelation distance and coherence bandwidth which describe the disturbed propagation channel. Results for mean time delay and time delay jitter are presented and compared to direct simulation results and to other theoretical calculations. A numerical/analytical solution to the parabolic wave equation is presented in the form of a multiple phase-screen (MPS) propagation simulation. In this very general solution technique, the random medium is divided into a finite number of layers. The field fluctuation through each layer is obtained by replacing the layer by a centrally located thin phase-screen, whose statistical properties are determined from the statistics of the electron-density irregularities. The waveform then propagates from phase-screen to phase-screen via an exact solution to the Fresnel-Kirchnoff equation. For wide bandwidth waveforms, numerical solutions are obtained at a number of discrete frequencies centered about the carrier and then time-domain solutions are obtained by Fourier transform techniques. Detailed results are presented for a number of interesting cases including propagation of a 20 MHz bandwidth signal through a finite barium cloud at a carrier frequency of 100 MHz. One of the major uses of the MPS propagation simulation has been to provide realizations of the received signal after propagation through a disturbed channel. The MPS simulation obtains a general solution of the parabolic wave equation under both weak and strong scattering conditions. A second technique to directly obtain realizations of wide bandwidth waveforms is presented here. This technique is limited to the case of strong scattering but requires only a fraction of the computer resources needed for MPS signal generation. Detailed comparisons of the two signal generation techniques are presented.
Rapid growth in new satellite services incorporating very small aperture terminals (VSAT) and ultra small aperture terminals (USAT) is expected in the coming years. Small size terminals allow for widespread use of satellite services in small business and domestic applications. Due to congestion of lower frequency bands such as C and Ku, most of these services will use Ka-band (2/20 GHz) frequencies. Propagation impairments produced by the troposphere is a limiting factor for the effective use of the 20/30 GHz band and the use of smaller Earth terminals makes it difficult to provide sufficient link margins for propagation related outages. In this context, reliable prediction of propagation impairments for low margin systems becomes important. Due to the complexity of propagation phenomena propagation modeling is mainly attempted on an empirical basis. As such, the availability of reliable measured data that extend to probability levels well in excess of the traditional limit of 1 percent is of great importance in the development, validation, and refinement of propagation models. The beacon payload on the Advanced Communications Technology Satellite (ACTS) together with the propagation measurement terminals developed under the NASA ACTS propagation program provide an excellent opportunity to collect such data on a long-term basis. This paper presents the results of ACTS propagation measurements conducted in the Washington, DC metropolitan area by COMSAT Laboratories.
Gears used in current helicopters and turboprops are designed for light weight, high margins of safety, and high reliability. However, unexpected gear failures may occur even with adequate tooth design. To design an extremely safe system, the designer must ask and address the question, "What happens when a failure occurs?" With gear-tooth bending fatigue, tooth or rim fractures may occur. A crack that propagates through a rim will be catastrophic, leading to disengagement of the rotor or propeller, loss of an aircraft, and possible fatalities. This failure mode should be avoided. A crack that propagates through a tooth may or may not be catastrophic, depending on the design and operating conditions. Also, early warning of this failure mode may be possible because of advances in modern diagnostic systems. One concept proposed to address bending fatigue fracture from a safety aspect is a splittooth gear design. The prime objective of this design would be to control crack propagation in a desired direction such that at least half of the tooth would remain operational should a bending failure occur. A study at the NASA Lewis Research Center analytically validated the crack-propagation failsafe characteristics of a split-tooth gear. It used a specially developed three-dimensional crack analysis program that was based on boundary element modeling and principles of linear elastic fracture mechanics. Crack shapes as well as the crack-propagation life were predicted on the basis of the calculated stress intensity factors, mixed-mode crack-propagation trajectory theories, and fatigue crack-growth theories. The preceding figures show the effect of the location of initial cracks on crack propagation. Initial cracks in the fillet of the teeth produced stress intensity factors of greater magnitude (and thus, greater crack growth rates) than those in the root or groove areas of the teeth. Crack growth was simulated in a case study to evaluate crack-propagation paths. Tooth fracture was predicted from the crackgrowth simulation for an initial crack in the tooth fillet region. This was the desired failure mode for an ultrasafe design. Lastly, tooth loads on the uncracked mesh of the split-tooth design were up to five times greater than those on the cracked mesh if equal deflections of the cracked and uncracked teeth were considered. This effect needs to be considered in the design of a split-tooth configuration. This work was done in-house at Lewis in support of the National Rotorcraft Technology Center project, Ultra-Safe Gear Design, with the Boeing Defense and Space Group. The crack-propagation package was developed by the Cornell Fracture Group at Cornell University. The reported results, which are the initial findings of Lewis gear-crackpropagation research for the Rotorcraft Base Program, will be further investigated to develop generalized gear design guidelines.
Fracture propagation caused by hydraulic fracturing operations can be significantly influenced by adjacent fractures. In this paper, we present a detailed coupled hydro-mechanical analysis to study the effects of nearby natural fractures on hydraulically induced fracture propagation and changes in fracture permeability. Two rock domains were considered in comparison: FD1, with one fracture, and FD2, with two adjacent parallel or non-parallel fractures. It was assumed that water injection occurred in a borehole that intersected the single fracture in FD1 and one of the two fractures in FD2. Simulations were made for a time period of 3 h with an injection period of 2 h followed by 1 h of shut-in. An elastic-brittle model based on the degradation of material properties was implemented in a 2D finite-difference scheme and used for elements of the intact rock subjected to tension and shear failure. The intact rock was considered to have a low but non-negligible permeability. A verification study against analytical solutions showed that the fracture propagation and stress concentrations due to differential boundary stresses could be accurately represented by our model. Next, a base case was considered, in which the stress ratio (SR) between the magnitudes of the maximum and minimum boundary stresses, the permeability k R of the intact rock and the initial permeability k TF of the tension failure regions were fixed. In FD2, the distance d F between the two natural fractures defined by the closest distance was also fixed. The results showed that in both rock domains, the fracture started to propagate when the pore pressure was approximately 85% of the magnitude of the minimum boundary stress. The propagation of a single fracture was significantly greater than the propagation of a double fracture system because, in the latter case, the pore pressure decreased when the two fractures connected. As a result, changes in permeability in FD2 were smaller than in FD1. At the end of injection, the maximum ratios between the final and initial permeability of the natural fractures were found to be approximately 3 and 2 for rock domains FD1 and FD2, respectively. For non-parallel fractures, the controlling factor for fracture propagation was the separation between the tips of the pressurised fracture and the neighbouring non-pressurised fracture. A sensitivity study was conducted to study the influence of the key parameters d F , SR, k R and k TF on the simulation results. Lastly, fracture propagation showed more sensitivity to dF and SR than to the other parameters.
Purpose: To investigate the feasibility of using DIR to propagate the manually contoured rectum and bladder from the 1st insertion to the new CT images on subsequent insertions and evaluate the segmentation performance. Methods: Ten cervical cancer patients, who were treated by T&R brachytherapy in 3–4 insertions, were retrospectively collected. In each insertion, rectum and bladder were manually delineated on the planning CT by a physicist and verified by a radiation oncologist. Using VelocityAI (Velocity Medical Solutions, Atlanta, GA), a rigid registration was firstly employed to match the bony structures between the first insertion and each of the following insertions, then a multi-pass B-spine DIR was carried out to further map the sub volume that encompasses rectum and bladder. The resultant deformation fields propagated contours, and dice similarity coefficient (DSC) was used to quantitatively evaluate the agreement between the propagated contours and the manually-delineated organs. For the 3rd insertion, we also evaluated if the segmentation performance could be improved by propagating the contours from the most recent insertion, i.e., the 2nd insertion. Results: On average, the contour propagation took about 1 minute. The average and standard deviation of DSC over all insertions and patients was 0.67±0.10 (range: 0.44–0.81) for rectum, and 0.78±0.07 (range: 0.63–0.87) for bladder. For the 3rd insertion, propagating contours from the 2nd insertion could improve the segmentation performance in terms of DSC from 0.63±0.10 to 0.72±0.08 for rectum, and from 0.77±0.07 to 0.79±0.06 for bladder. A Wilcoxon signed rank test indicated that the improvement was statistically significant for rectum (p = 0.004). Conclusion: The preliminary results demonstrate that deformable image registration could efficiently and accurately propagate rectum and bladder contours between CT images in different T&R brachytherapy fractions. We are incorporating the propagated contours into our learning-based method to further segment these organs.