Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

A numerical study of the effects of ambipolar diffusion on the collapse of magnetic gas clouds

The gravitational collapse of isothermal, nonrotating magnetic gas clouds have been calculated numerically, including the effects of ambipolar diffusion. The fractional ionization in the clouds is approximated by a power-law function of the gas density, f = K/n to the q-power, where K and q are adjustable parameters. Eleven numerical experiments were run, and the results indicate that the asymptotic character of collapse is determined mainly by the value of q and is largely independent of the other parameters characterizing a cloud (e.g., K, cloud mass). In particular, there is nearly a one-to-one correspondence between q and the slope, x, of the central magnetic field strength-gas density relationship. If q is no more than 0.8, a cloud collapses asymptotically, as though the magnetic field were 'frozen' to the neutral matter. The magnetic field strength at the center of a collapsing cloud is strongly amplified during collapse even for values of q of about 1, despite extremely low values of fractional ionization. A discussion of the theoretical basis for this unexpected behavior is given. Possible implications of our results for the problems of magnetic braking of rotating protostars and star formation in general are also presented.

Black, D. C.↗
NASA NTRS2000

Laboratory Astrochemistry: Interstellar PAHs

Polycyclic aromatic hydrocarbons (PAHs) are now considered to be an important and ubiquitous component of the organic material in space. PAHs are found in a large variety of extraterrestrial materials such as interplanetary dust particles (IDPs) and meteoritic materials. PAHs are also good candidates to account for the infrared emission bands (UIRs) and the diffuse interstellar optical absorption bands (DIBs) detected in various regions of the interstellar medium. The recent observations made with the Infrared Space Observatory (ISO) have confirmed the ubiquitous nature of the UIR bands and their carriers. PAHs are thought to form through chemical reactions in the outflow from carbon-rich stars in a process similar to soot formation. Once injected in the interstellar medium, PAHs are further processed by the interstellar radiation field, interstellar shocks and energetic particles. A major, dedicated, laboratory effort has been undertaken to measure the physical and chemical characteristics of these complex molecules and their ions under experimental conditions that mimic the interstellar conditions. These measurements require collision-free conditions where the molecules and ions are cold and chemically isolated. The spectroscopy of PAHs under controlled conditions represents an essential diagnostic tool to study the evolution of extraterrestrial PAHs. The Astrochemistry Laboratory program will be discussed through its multiple aspects: (1) objectives, (2) approach and techniques adopted, (3) adaptability to the nature of the problem(s), and (4) results and implications for astronomy as well as for molecular spectroscopy. A review of the data generated through laboratory simulations of space environments and the role these data have played in our current understanding of the properties of interstellar PAHs will be presented. The discussion will also introduce the newest generation of laboratory experiments that are currently being developed in order to provide a closer simulation of space environments and a better support to space missions.

Salama, Farid↗
NASA NTRS2000

Laboratory Astrochemistry: Interstellar PAH Analogs

Polycyclic aromatic hydrocarbons (PAHs) are now considered to be an important and ubiquitous component of the organic material in space. PAHs are found in a large variety of extraterrestrial materials such as interplanetary dust particles (IDPs) and meteoritic materials. PAHs are also good candidates to account for the infrared emission bands (UIRs) and the diffuse interstellar optical absorption bands (DIBs) detected in various regions of the interstellar medium. The recent observations made with the Infrared Space Observatory (ISO) have confirmed the ubiquitous nature of the UIR bands and their carriers. PAHs are though to form through chemical reactions in the outflow from carbon-rich stars in a process similar to soot formation. Once injected in the interstellar medium, PAHs are further processed by the interstellar radiation field, interstellar shocks and energetic particles. A major, dedicated, laboratory effort has been undertaken over the past years to measure the physical and chemical characteristics of these complex molecules and their ions under experimental conditions that mimic the interstellar conditions. These measurements require collision-free conditions where the molecules and ions are cold and chemically isolated. The spectroscopy of PAHs under controlled conditions represents an essential diagnostic tool to study the evolution of extraterrestrial PAHs. The Astrochemistry Laboratory program will be discussed through its multiple aspects: objectives, approach and techniques adopted, adaptability to the nature of the problem(s), results and implications for astronomy as well as for molecular spectroscopy. A review of the data generated through laboratory simulations of space environments and the role these data have played in our current understanding of the properties of interstellar PAHs will be presented. The discussion will also introduce the newest generation of laboratory experiments that are currently being developed in order to provide a closer simulation of space environments and a better support to space missions.

Salama, Farid↗
NASA NTRS2000

Laboratory Studies of Interstellar PAH Analogs

Polycyclic aromatic hydrocarbons (PAHs) are now considered to be an important and ubiquitous component of the organic material in space. PAHs are found in a large variety of extraterrestrial materials such as interplanetary dust particles (IDPs) and meteoritic materials. PAHs are also good candidates to account for the infrared emission bands (UIRs) and the diffuse interstellar optical absorption bands (DIBs) detected in various regions of the interstellar medium. The recent observations made with the Infrared Space Observatory (ISO) have confirmed the ubiquitous nature of the UIR bands and their carriers. PAHs are though to form through chemical reactions in the outflow from carbon-rich stars in a process similar to soot formation. Once injected in the interstellar medium, PAHs are further processed by the interstellar radiation field, interstellar shocks and energetic particles. A major, dedicated, laboratory effort has been undertaken over the past years to measure the physical and chemical characteristics of these complex molecules and their ions under experimental conditions that mimic the interstellar conditions. These measurements require collision-free conditions where the molecules and ions are cold and chemically isolated. The spectroscopy of PAHs under controlled conditions represents an essential diagnostic tool to study the evolution of extraterrestrial PAHs. The Astrochemistry Laboratory program will be discussed through its multiple aspects: objectives, approach and techniques adopted, adaptability to the nature of the problem(s), results and implications for astronomy as well as for molecular spectroscopy. A review of the data generated through laboratory simulations of space environments and the role these data have played in our current understanding of the properties of interstellar PAHs will be presented. The discussion will also introduce the newest generation of laboratory experiments that are currently being developed in order to provide a closer simulation of space environments and a better support to space missions.

Salama, Farid↗
NASA NTRS1991

Leading-edge receptivity for blunt-nose bodies

This research program investigates boundary-layer receptivity in the leading-edge region for bodies with blunt leading edges. Receptivity theory provides the link between the unsteady distrubance environment in the free stream and the initial amplitudes of the instability waves in the boundary layer. This is a critical problem which must be addressed in order to develop more accurate prediction methods for boundary-layer transition. The first phase of this project examines the effects of leading-edge bluntness and aerodynamic loading for low Mach number flows. In the second phase of the project, the investigation is extended to supersonic Mach numbers. Singular perturbation techniques are utilized to develop an asymptotic theory for high Reynolds numbers. In the first year, the asymptotic theory was developed for leading-edge receptivity in low Mach number flows. The case of a parabolic nose is considered. Substantial progress was made on the Navier-Sotkes computations. Analytical solutions for the steady and unsteady potential flow fields were incorporated into the code, greatly expanding the types of free-stream disturbances that can be considered while also significantly reducing the the computational requirements. The time-stepping algorithm was modified so that the potential flow perturbations induced by the unsteady pressure field are directly introduced throughout the computational domain, avoiding an artificial 'numerical diffusion' of these from the outer boundary. In addition, the start-up process was modified by introducing the transient Stokes wave solution into the downstream boundary conditions.

Kerschen, Edward J.↗
NASA NTRS2007

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗
NASA NTRS1992

A status of the activities of the NASA/MSFC pump stage technology team

The Consortium for Computational Fluid Dynamics (CFD) Application in Propulsion Technology was established to aid the transfer of CFD related advancements among academia, government agencies, and industry. The specific goals of the Consortium are to develop CFD methodologies necessary to solve propulsion problems, to validate these methodologies, and to apply these methodologies in the design process. To accomplish these goals, a team of experts in various related fields was formed, a schedule of activities necessary to meet the goals was generated, and funding for the activities was obtained from NASA. During the past year (Mar. 1991 - Mar. 1992) the team's activities have focused on preliminary code validation and on the design of an advanced impeller. Six codes were used to calculate the flow in a Rocketdyne 0.3 flow coefficient inducer, and the results were compared to L2F data available for the inducer. This activity identified shortcomings in the experimental data sets and in the analytical solutions which must be surmounted in any future team activity. The design of the advanced impeller relied heavily on CFD results to obtain an optimized geometry. The optimized geometry was analyzed using four different codes, at design and off-design conditions. Activities for the next year include the optimization of a tandem blade impeller design, benchmark of CFD codes for diffuser and volute flows, the collection of L2F data for 'state-of-the-art' impeller and inducer, and the verification of the advanced pump team impeller design in a water rig.

Garcia, R.↗
NASA NTRS1992

Computational fluid dynamics study of the variable-pitch split-blade fan concept

A computational fluid dynamics study was conducted to evaluate the feasibility of the variable-pitch split-blade supersonic fan concept. This fan configuration was conceived as a means to enable a supersonic fan to switch from the supersonic through-flow type of operation at high speeds to a conventional fan with subsonic inflow and outflow at low speeds. During this off-design, low-speed mode of operation, the fan would operate with a substantial static pressure rise across the blade row like a conventional transonic fan; the front (variable-pitch) blade would be aligned with the incoming flow, and the aft blade would remain fixed in the position set by the supersonic design conditions. Because of these geometrical features, this low speed configuration would inherently have a large amount of turning and, thereby, would have the potential for a large total pressure increase in a single stage. Such a high-turning blade configuration is prone to flow separation; it was hoped that the channeling of the flow between the blades would act like a slotted wing and help alleviate this problem. A total of 20 blade configurations representing various supersonic and transonic configurations were evaluated using a Navier Stokes CFD program called ADAPTNS because of its adaptive grid features. The flow fields generated by this computational procedure were processed by another data reduction program which calculated average flow properties and simulated fan performance. These results were employed to make quantitative comparisons and evaluations of blade performance. The supersonic split-blade configurations generated performance comparable to a single-blade supersonic, through-flow fan configuration. Simulated rotor total pressure ratios of the order of 2.5 or better were achieved for Mach 2.0 inflow conditions. The corresponding fan efficiencies were approximately 75 percent or better. The transonic split-blade configurations having large amounts of turning were able to generate large amounts of total turning and achieve simulated total pressure ratios of 3.0 or better with subsonic inflow conditions. These configurations had large losses and low fan efficiencies in the 70's percent. They had large separated regions and low velocity wakes. Additional turning and diffusion of this flow in a subsequent stator row would probably be very inefficient. The high total pressure ratios indicated by the rotor performance would be substantially reduced by the stators, and the stage efficiency would be substantially lower. Such performance leaves this dual-mode fan concept less attractive than originally postulated.

Kepler, C. E.↗
NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗
NASA NTRS2001

Spherical Ethylene/Air Diffusion Flames Subject to Concentric DC Electric Field in Microgravity

It is well known that microgravity conditions, by eliminating buoyant flow, enable many combustion phenomena to be observed that are not possible to observe at normal gravity. One example is the spherical diffusion flame surrounding a porous spherical burner. The present paper demonstrates that by superimposing a spherical electrical field on such a flame, the flame remains spherical so that we can study the interaction between the electric field and flame in a one-dimensional fashion. Flames are susceptible to electric fields that are much weaker than the breakdown field of the flame gases owing to the presence of ions generated in the high temperature flame reaction zone. These ions and the electric current of the moving ions, in turn, significantly change the distribution of the electric field. Thus, to understand the interplay between the electric field and the flame is challenging. Numerous experimental studies of the effect of electric fields on flames have been reported. Unfortunately, they were all involved in complex geometries of both the flow field and the electric field, which hinders detailed study of the phenomena. In a one-dimensional domain, however, the electric field, the flow field, the thermal field and the chemical species field are all co-linear. Thus the problem is greatly simplified and becomes more tractable.

Yuan, Z. -G.↗
NASA NTRS1989

Report of the 1988 2-D Intercomparison Workshop, chapter 3

Several factors contribute to the errors encountered. With the exception of the line-by-line model, all of the models employ simplifying assumptions that place fundamental limits on their accuracy and range of validity. For example, all 2-D modeling groups use the diffusivity factor approximation. This approximation produces little error in tropospheric H2O and CO2 cooling rates, but can produce significant errors in CO2 and O3 cooling rates at the stratopause. All models suffer from fundamental uncertainties in shapes and strengths of spectral lines. Thermal flux algorithms being used in 2-D tracer tranport models produce cooling rates that differ by as much as 40 percent for the same input model atmosphere. Disagreements of this magnitude are important since the thermal cooling rates must be subtracted from the almost-equal solar heating rates to derive the net radiative heating rates and the 2-D model diabatic circulation. For much of the annual cycle, the net radiative heating rates are comparable in magnitude to the cooling rate differences described. Many of the models underestimate the cooling rates in the middle and lower stratosphere. The consequences of these errors for the net heating rates and the diabatic circulation will depend on their meridional structure, which was not tested here. Other models underestimate the cooling near 1 mbar. Suchs errors pose potential problems for future interactive ozone assessment studies, since they could produce artificially-high temperatures and increased O3 destruction at these levels. These concerns suggest that a great deal of work is needed to improve the performance of thermal cooling rate algorithms used in the 2-D tracer transport models.

Jackman, Charles H.↗
NASA NTRS2010

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗
NASA NTRS1995

A new approach to the linear theory of single-species tearing in two-dimensional quasi-neutral sheets

We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Brittnacher, M.↗
NASA NTRS2001

Magnetized Target Fusion: Prospects for Low-Cost Fusion Energy

Magnetized Target Fusion (MTF) has attracted renewed interest in recent years because it has the potential to resolve one of the major problems with conventional fusion energy research - the high cost of facilities to do experiments and in general develop practical fusion energy. The requirement for costly facilities can be traced to fundamental constraints. The Lawson condition implies large system size in the case of conventional magnetic confinement, or large heating power in the case of conventional inertial confinement. The MTF approach is to use much higher fuel density than with conventional magnetic confinement (corresponding to megabar pressures), which results in a much-reduced system size to achieve Lawson conditions. Intrinsically the system must be pulsed because the pressures exceed the strength of any known material. To facilitate heating the fuel (or "target") to thermonuclear conditions with a high-power high-intensity source of energy, magnetic fields are used to insulate the high-pressure fuel from material surroundings (thus "magnetized target"). Because of magnetic insulation, the required heating power intensity is reduced by many orders of magnitude compared to conventional inertial fusion, even with relatively poor energy confinement in the magnetic field, such as that characterized by Bohm diffusion. In this paper we show semi-quantitatively why MTF-should allow fusion energy production without costly facilities within the same generally accepted physical constraints used for conventional magnetic and inertial fusion. We also briefly discuss potential applications of this technology ranging from nuclear rockets for space propulsion to a practical commercial energy system. Finally, we report on the exploratory research underway, and the interesting physics issues that arise in the MTF regime of parameters. Experiments at Los Alamos are focused on formation of a suitable plasma target for compression, utilizing the knowledge base for compact toroids called Field-Reversed Configurations. As reported earlier, it appears that the existing pulsed-power Shiva Star facility at the Air Force Research Laboratory in Albuquerque, NM can satisfy the heating requirements by means of imploding a thin metal cylinder (called a "liner") surrounding an FRC of the type presently being developed. The proposed next step is an integrated liner-on-plasma experiment in which an FRC would be heated to 10 keV by the imploding liner.

Siemon, Richard E.↗
NASA NTRS2001

Additive Runge-Kutta Schemes for Convection-Diffusion-Reaction Equations

Additive Runge-Kutta (ARK) methods are investigated for application to the spatially discretized one-dimensional convection-diffusion-reaction (CDR) equations. First, accuracy, stability, conservation, and dense output are considered for the general case when N different Runge-Kutta methods are grouped into a single composite method. Then, implicit-explicit, N = 2, additive Runge-Kutta ARK2 methods from third- to fifth-order are presented that allow for integration of stiff terms by an L-stable, stiffly-accurate explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method while the nonstiff terms are integrated with a traditional explicit Runge-Kutta method (ERK). Coupling error terms are of equal order to those of the elemental methods. Derived ARK2 methods have vanishing stability functions for very large values of the stiff scaled eigenvalue, z(exp [I]) goes to infinity, and retain high stability efficiency in the absence of stiffness, z(exp [I]) goes to zero. Extrapolation-type stage-value predictors are provided based on dense-output formulae. Optimized methods minimize both leading order ARK2 error terms and Butcher coefficient magnitudes as well as maximize conservation properties. Numerical tests of the new schemes on a CDR problem show negligible stiffness leakage and near classical order convergence rates. However, tests on three simple singular-perturbation problems reveal generally predictable order reduction. Error control is best managed with a PID-controller. While results for the fifth-order method are disappointing, both the new third- and fourth-order methods are at least as efficient as existing ARK2 methods while offering error control and stage-value predictors.

Kennedy, Christopher A.↗
NASA NTRS1996

Non-linear resonance of fluids in a crystal growth cavity

In the microgravity environment, the effect of gravity on fluid motion is much reduced. Hence, secondary effects such as vibrations, jitters, surface tension, capillary effects, and electromagnetic forces become the dominant mechanism of fluid convection. Numerous studies have been conducted to investigate fluid behavior in microgravity with the ultimate goal of developing processes with minimal influence from convection. Industrial applications such as crystal growth from solidification of melt and protein growth for pharmatheutical application are just a few examples of the vast potential benefit that can be reaped from material processing in space. However, a space laboratory is not immune from all undesirable disturbances and it is imperative that such disturbances be well understood, quantifiable, and controlled. Non-uniform and transient accelerations such as vibrations, jitters, and impulsive accelerations exist as a result of crew activities, space vehicle maneuvering, and the operations of on-board equipment. Measurements conducted on-board a U.S. Spacelab showed the existence of vibrations in the frequency range of 1 to 100 Hz with a dominant mode of 17 Hz and harmonics of 54 Hz. The observed vibration is not limited to any coordinate plane but exists in all directions. Similar situation exists on-board the Russian MIR space station. Due to the large structure of its design, the future International Space Station will have its own characteristic vibration spectrum. It is well known that vibration can exert substantial influence on heat and mass transfer processes, thus hindering any attempts to achieve a diffusion-limited process. Experiments on vibration convection for a liquid-filled enclosure under one-g environment showed the existence of different flow regimes as vibration frequency and intensity changes. Results showed the existence of a resonant frequency, near which the enhancement is the strongest, and the existence of a high frequency asymptote. Numerical simulations of vibration convection have been conducted by Yurkov, Fu and Shieh, and by Wang. These analyses considered a two-dimensional air-filled cell under weightlessness condition and showed results similar to those of the experiments. It is not yet known whether resonance convection can be triggered by jitter alone or whether it requires the interaction of jitter with other convective forces in low gravity. An order of magnitude analysis, however, can be used to show the dependence of the resonance frequency on the fluid Prandtl number. Even though the onset of resonance convection may depend on other factors, results indicates that fluids with low Prandtl numbers are more susceptible to resonance than those with high Prandtl numbers. The current study is aimed at gaining additional insights to this problem using germanium as working fluid. Germanium was chosen for this analysis because of its common usage in solidification process and its relatively low Prandtl number (Pr = 0.02).

Wang, Francis C.↗
NASA NTRS2002

Fatigue Crack and Porosity Measurement in Composite Materials by Thermographic and Ultrasonic Methods

Many nondestructive methods exist for the detection of localized material anomalies in an otherwise good composite structure. The problem arises when the material system as a whole has degraded during service or was improperly manufactured. Porosity and intra-ply microcracking are two such conditions that in unlined composite pressure vessels can be very troublesome to detect and when linked through the thickness can be critical to mission success. These leak paths may lead to loss of pressure/propellant, increased risk of explosion and possible cryo-pumping. Research sought nondestructive methods for quantifying porosity and microcracking in composite tankage. Both thermographic and resonance ultrasound methods have been utilized with artificial neural network and statistical approaches to analyze the data. Resonant ultrasound spectroscopy provides measurements, which are sensitive to fine details in the materials character, such as micro-cracking and porosity. Here, the higher frequency (shorter wavelength) components of the signal train provide more significant interaction with the defects causing the spectral characteristics to shift toward lower amplitudes at the higher frequencies. As the density of the defects increases more interactions occur and more drastic amplitude changes are observed. From a thermal perspective, the higher the defect density the lower the through thickness thermal diffusivity will be. Utilizing a point heat source, and thermographically recording the heat profile with time, diffusivity calculations can be made which in turn can be related to the relative quality of the material. Preliminary experiments to verify the measurable effect on the resonance spectrum of the ultrasonic data to detect microcracking and for porosity detection thermographically are presented. Methods involving supervised and unsupervised artificial neural networks as well as other clustering algorithms are developed for signal identification.

Walker, James L.↗
NASA NTRS2015

Free-Energy Calculations. A Mathematical Perspective

Ion channels are pore-forming assemblies of transmembrane proteins that mediate and regulate ion transport through cell walls. They are ubiquitous to all life forms. In humans and other higher organisms they play the central role in conducting nerve impulses. They are also essential to cardiac processes, muscle contraction and epithelial transport. Ion channels from lower organisms can act as toxins or antimicrobial agents, and in a number of cases are involved in infectious diseases. Because of their important and diverse biological functions they are frequent targets of drug action. Also, simple natural or synthetic channels find numerous applications in biotechnology. For these reasons, studies of ion channels are at the forefront of biophysics, structural biology and cellular biology. In the last decade, the increased availability of X-ray structures has greatly advanced our understanding of ion channels. However, their mechanism of action remains elusive. This is because, in order to assist controlled ion transport, ion channels are dynamic by nature, but X-ray crystallography captures the channel in a single, sometimes non-native state. To explain how ion channels work, X-ray structures have to be supplemented with dynamic information. In principle, molecular dynamics (MD) simulations can aid in providing this information, as this is precisely what MD has been designed to do. However, MD simulations suffer from their own problems, such as inability to access sufficiently long time scales or limited accuracy of force fields. To assess the reliability of MD simulations it is only natural to turn to the main function of channels - conducting ions - and compare calculated ionic conductance with electrophysiological data, mainly single channel recordings, obtained under similar conditions. If this comparison is satisfactory it would greatly increase our confidence that both the structures and our computational methodologies are sufficiently accurate. Channel conductance, defined as the ratio of ionic current through the channel to applied voltage, can be calculated in MD simulations by way of applying an external electric field to the system and counting the number of ions that traverse the channel per unit time. If the current is small, a voltage significantly higher than the experimental one needs to be applied to collect sufficient statistics of ion crossing events. Then, the calculated conductance has to be extrapolated to the experimental voltage using procedures of unknown accuracy. Instead, we propose an alternative approach that applies if ion transport through channels can be described with sufficient accuracy by the one-dimensional diffusion equation in the potential given by the free energy profile and applied voltage. Then, it is possible to test the assumptions of the equation, recover the full voltage/current dependence, determine the reliability of the calculated conductance and reconstruct the underlying (equilibrium) free energy profile, all from MD simulations at a single voltage. We will present the underlying theory, model calculations that test this theory and simulations on ion conductance through a channel that has been extensively studied experimentally. To our knowledge this is the first case in which the complete, experimentally measured dependence of the current on applied voltage has been reconstructed from MD simulations.

free energy↗