Search NASA⌕ Search

SEARCH · Search NASA

Results for “radial diffusion”

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 847 records · Page 47

Saturn Neutron Exosphere as Source for Inner and Innermost Radiation Belts

Energetic proton and electron measurements by the ongoing Cassini orbiter mission are expanding our knowledge of the highest energy components of the Saturn magnetosphere in the inner radiation belt region after the initial discoveries of these belts by the Pioneer 11 and Voyager 2 missions. Saturn has a neutron exosphere that extends throughout the magnetosphere from the cosmic ray albedo neutron source at the planetary main rings and atmosphere. The neutrons emitted from these sources at energies respectively above 4 and 8 eV escape the Saturn system, while those at lower energies are gravitationally bound. The neutrons undergo beta decay in average times of about 1000 seconds to provide distributed sources of protons and electrons throughout Saturn's magnetosphere with highest injection rates close to the Saturn and ring sources. The competing radiation belt source for energetic electrons is rapid inward diffusion and acceleration of electrons from the middle magnetosphere and beyond. Minimal losses during diffusive transport across the moon orbits, e.g. of Mimas and Enceladus, and local time asymmetries in electron intensity, suggest that drift resonance effects preferentially boost the diffusion rates of electrons from both sources. Energy dependences of longitudinal gradient-curvature drift speeds relative to the icy moons are likely responsible for hemispheric differences (e.g., Mimas, Tethys) in composition and thermal properties as at least partly produced by radiolytic processes. A continuing mystery is the similar radial profiles of lower energy (<10 MeV) protons in the inner belt region. Either the source of these lower energy protons is also neutron decay, but perhaps alternatively from atmospheric albedo, or else all protons from diverse distributed sources are similarly affected by losses at the moon' orbits, e.g. because the proton diffusion rates are extremely low. Enceladus cryovolcanism, and radiolytic processing elsewhere on the icy moon and ring surfaces, are additional sources of protons via ionization and charge exchange from breakup of water molecules. But one must then account somehow for local acceleration to the observed keV-MeV energies, since moon sweeping and E-ring absorption would remove protons diffusing inward from the middle magnetosphere. Although the main rings block further inward diffusion from the inner radiation belts, the exospheric neutron-decay source, combined with much slower diffusion of protons relative to electrons, may produce an innermost radiation belt in the gap between the upper atmosphere and the D-ring. This innermost belt will first be explored in-situ during the final proximal orbits of the Cassini mission.

Cooper, John↗

Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simu-lations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-di?usion equations using a non-negative finite element formulation for di?erent input parameters. The reactive-mixing model input parameters are: time-scale associated with flipping of velocity, spatial-scale controlling small/large vortex structures of velocity, perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. We observed that anisotropic dispersion strength/contrast is the most important feature and time-scale associated with flipping of velocity is the least important feature. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. The constructed SVR-ROMs are then used to predict scaling of QoIs. We also present estimates and inequalities on the QoIs, which inform that the species decay, mix, and produce in an exponential fashion. These inequalities also inform that a radial basis function is the most suitable kernel for the SVM/SVR models for QoIs. It is observed that R2-score for SVR-ROMs on unseen data is greater than 0.9, implying that the SVR-ROMs are able to predict the reaction-diffusion system state reasonably well. Finally, in terms of the computational cost, the proposed SVM-ROMs are O(107) times faster than running a high-fidelity finite element simulation for evaluating QoIs. This makes the proposed ML-based ROMs attractive for reactive-transport sensing and real-time monitoring applications as they are significantly faster yet reasonably accurate.

Mudunuru, Maruti K.↗

Experimentally reproduced relict enstatite in porphyritic chondrules of enstatite chondrite composition

Experiments are presented that test a model for the origin of porphyritic pyroxene (PP) chondrules in enstatite chondrites that contain phenocrysts of enstatite with blue cathodoluminescence (CL) set in a matrix of radial, dendritic enstatite with red CL. Established one-atmosphere, gas-mixing techniques were used. Relict enstatite phenocrysts with blue CL in a matrix of coarsely radial to dendritic enstatite with red CL were successfully produced. The relict crystals are preserved in runs with a melt time of 36 minutes or less at 1537 C. The relicts remain angular with smooth crystal/melt interfaces, and thus melting has occurred uniformly. Partial melting does occur along fractures produced when the blue CL enstatite was initially grown and cooled through the proto/ortho enstatite transition with the attendent volume change. There is either reaction with the melt and diffusion of Mn and Cr into the blue CL En, or there is an overgrowth of red CL En along the fractures. The bulk of the relicts remain blue. The melt enclosing the relicts crystalized to a coarsely radial to dendritic to micro porphyritic texture comprised of enstatite that has a bright red CL with decreasing melt time. The blue CL En has Mn and Cr contents at or below detection limits of the electron probe as described in earlier studies and in natural blue CL En. In the red CL En in this study, the Mn, Al2O3, and Cr are at previously observed levels and the levels change rapidly.

Lofgren, Gary E.↗

STOCHASTIC PARTICLE ACCELERATION IN TURBULENCE GENERATED BY MAGNETOROTATIONAL INSTABILITY

We investigate stochastic particle acceleration in accretion flows. It is believed that magnetorotational instability (MRI) generates turbulence inside accretion flows and that cosmic rays (CRs) are accelerated by the turbulence. We calculate equations of motion for CRs in the turbulent fields generated by MRI with the shearing box approximation and without back reaction to the field. Our results show that the CRs randomly gain or lose their energy through interaction with the turbulent fields. The CRs diffuse in the configuration space anisotropically: the diffusion coefficient in the direction of the unperturbed flow is about 20 times higher than the Bohm coefficient, while those in the other directions are only a few times higher than the Bohm. The momentum distribution is isotropic and its evolution can be described by the diffusion equation in momentum space where the diffusion coefficient is a power-law function of the CR momentum. We show that the shear acceleration works efficiently for energetic particles. We also cautiously note that in the shearing box approximation, particles that cross the simulation box many times along the radial direction undergo unphysical runaway acceleration by the Lorentz transformation, which needs to be taken into account with special care.

79 ASTRONOMY AND ASTROPHYSICS↗

Controllable transition from positive space charge to negative space charge in an inverted cylindrical magnetron

The combined effect of magnetic field (B), gas pressure (P), and the corresponding discharge voltage on the discharge properties of argon in inverted cylindrical magnetron has been investigated. In the experiment, anode is biased with continuous 10 ms sinusoidal half wave. It is observed that at a comparatively higher magnetic field (i.e., >200 gauss) and lower operating pressure (i.e., <1 × 10{sup −3} mbar), the discharge extinguishes and demands a high voltage to reignite. Discharge current increases with increase in magnetic field and starts reducing at sufficiently higher magnetic field for a particular discharge voltage due to restricted electron diffusion towards the anode. It is observed that B/P ratio plays an important role in sustaining the discharge and is constant for a discharge voltage. The discharge is transformed to negative space charge regime from positive space charge regime at certain B/P ratio and this ratio varies linearly with the discharge voltage. The space charge reversal is indicated by the radial profile of the floating potential and plasma potential in between two electrodes for different magnetic fields. At a particular higher magnetic field (beyond 100 gauss), the floating potential increases gradually with the radial distance from cathode, whereas it remains almost constant at lower magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-series forecasting using manifold learning, radial basis function interpolation, and geometric harmonics

We address a three-tier numerical framework based on nonlinear manifold learning for the forecasting of high-dimensional time series, relaxing the “curse of dimensionality” related to the training phase of surrogate/machine learning models. At the first step, we embed the high-dimensional time series into a reduced low-dimensional space using nonlinear manifold learning (local linear embedding and parsimonious diffusion maps). Then, we construct reduced-order surrogate models on the manifold (here, for our illustrations, we used multivariate autoregressive and Gaussian process regression models) to forecast the embedded dynamics. Finally, we solve the pre-image problem, thus lifting the embedded time series back to the original high-dimensional space using radial basis function interpolation and geometric harmonics. The proposed numerical data-driven scheme can also be applied as a reduced-order model procedure for the numerical solution/propagation of the (transient) dynamics of partial differential equations (PDEs). In conclusion, we assess the performance of the proposed scheme via three different families of problems: (a) the forecasting of synthetic time series generated by three simplistic linear and weakly nonlinear stochastic models resembling electroencephalography signals, (b) the prediction/propagation of the solution profiles of a linear parabolic PDE and the Brusselator model (a set of two nonlinear parabolic PDEs), and (c) the forecasting of a real-world data set containing daily time series of ten key foreign exchange rates spanning the time period 3 September 2001–29 October 2020.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Laser thermoelastic generation in metals above the melt threshold

An approach is presented for calculating thermoelastic generation of ultrasound in a metal plate exposed to nanosecond pulsed laser heating, sufficient to cause melting but not ablation. Detailed consideration is given to the spatial and temporal profiles of the laser pulse, penetration of the laser beam into the sample, the appearance and subsequent growth and then contraction of the melt pool, and the time dependent thermal conduction in the melt and surrounding solid throughout. The excitation of the ultrasound takes place during and shortly after the laser pulse and occurs predominantly within the thermal diffusion length of a micron or so beneath the surface. It is shown how, because of this, the output of the thermal simulations can be expressed as axially symmetric transient radial and normal surface force distributions. The epicentral displacement response to these force distributions is obtained by two methods, the one based on the elastodynamic Green's functions for plate geometry determined by the Cagniard generalized ray method and the other using a finite element numerical method. The two approaches are in very close agreement. Numerical simulations are reported on the epicentral displacement response of a 3.12 mm thick tungsten plate irradiated with a 4 ns pulsed laser beam with Gaussian spatial profile, at intensities below and above the melt threshold.

36 MATERIALS SCIENCE↗

Assessing high fidelity multi-component models to facilitate safeguards at Gas Centrifuge Enrichment Plants

We report that the International Atomic Energy Agency (IAEA) inspectors routinely carry out environmental sampling (ES) as a verification method. Collection of environmental swipe samples at various locations in Gas Centrifuge Enrichment Facilities (GCEPs) is an important process in detecting misuse of a declared facility and possibly the existence of undeclared nuclear material. These samples are measured for isotopic composition in uranium containing particles by Thermal Ionization Mass Spectrometry (TIMS) or Inductively Coupled Plasma Mass Spectrometry (ICP-MS). Even though ES is highly effective in detecting the absolute value of enrichments and their deviations from the declared values, it cannot explain the cause of those changes. Several potential explanations can serve as possibilities for particles detected above or other than the declared enrichment. These include normal and non-malicious events such as the design of the enrichment cascades, unintentional failures of the machines, or overshoot during startup of a cascade. It can also be the result of deliberate misuse by the facility operators. The primary objective of this work is to understand how these factors affect the enrichments produced by a cascade and quantify anticipated multi-isotopic concentrations for each case. The following methodology is employed to determine signatures at a particular facility. 1) Utilize a new two-dimensional multi-component diffusion code to obtain centrifuge performance data and use that information to design and perform cascade analysis. Compare and contrast the results with previous 1-D radially averaged solutions from the Pancake code. 2) Design a GCEP cascade with the production goal of 19.75% 235 U for each set of machine data above. Investigate two cascade scenarios that include enrichment of natural uranium (NU) feed to 19.75% 235U in a single cascade compared to a two-step process of NU to 5% and then 5% to 19.75%. 3) Simulate the intentional vs. unintentional off-normal scenarios in the cascades to assess the differences in isotopic concentrations. A non-ideal squared-off cascade model developed at the University of Virginia is used to calculate flow rates and isotopic concentrations of the process gas. The analysis is performed using the Rome machine model operated at 600 m/s rotor speed. The upper and lower bounds of normal and abnormal enrichments in a typical facility are used in conjunction with ES results to understand the root causes of such observations.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Performance and boundary-layer evaluation of a sonic inlet

Tests were conducted to determine the boundary layer characteristics and aerodynamic performance of a radial vane sonic inlet with a length/diameter ratio of 1 for several vane configurations. The sonic inlet was designed with a slight wavy wall type of diffuser geometry, which permits operation at high inlet Mach numbers (sufficiently high for good noise suppression) without boundary layer flow separation and with good total pressure recovery. A new method for evaluating the turbulent boundary layer was developed to separate the boundary layer from the inviscid core flow, which is characterized by a total pressure variation from hub to tip, and to determine the experimental boundary layer parameters.

Schmidt, J. F.↗

Impulsive solar X-ray bursts. III - Polarization, directivity, and spectrum of the reflected and total bremsstrahlung radiation from a beam of electrons directed toward the photosphere

The paper presents the spectrum, directivity, and state of polarization of the bremsstrahlung radiation expected from a beam of high-energy electrons spiraling along radial magnetic field lines toward the photosphere. A Monte Carlo method is then described for evaluation of the spectrum, directivity, and polarization of X-rays diffusely reflected from stellar photospheres. The accuracy of the technique is evaluated through comparison with analytic results. The calculated characteristics of the incident X-rays are used to evaluate the spectrum, directivity, and polarization of the reflected and total X-ray fluxes. The results are compared with observations.

Langer, S. H.↗

Payload Fairing Acoustic Trade Study: Fill Effect Modeling

A review of acoustic fill effect was conducted using vibroacoustic models of the NASA Space Launch System Block 1B 8.4-meter Payload Fairing. Fill effect is the term used to describe how sound pressure levels (SPLs) inside of a launch vehicle’s payload cavity change when a payload is present, compared to that of an unfilled cavity. Several vibroacoustic models with various payload shapes were developed and the deviation in output from an unfilled cavity model were compared to an analytical fill effect calculation. The analytical expression for the fill effect correction factor, known as the fill factor, is published in NASA-STD-7001B. The development and validation of the NASA fill factor is discussed to provide context for this study and define important variables including the volume fill percentage and the fairing-payload separation gap. The fill factor is reviewed from two perspectives; a global perspective which defines the volume ratio with respect to the entire fairing cavity, and a local perspective which defines the volume ratio with respect to a zone of interest dependent on the payload’s height. The global perspective was considered because it allows a decoupling of the volume ratio and fairing-payload separation gap variables. Both the global and local fill factor methods have difficulty capturing low frequency SPL because of a breakdown in the underlying assumptions of a diffuse acoustic field. Spatial variation in the pressure field at low frequency showed that axially oriented standing waves may be of concern for sensitive payload components; radial variation in SPL was found to be less significant. Further, a novel secondary derivation of the theoretical fill factor is posed for future work in which the absorption characteristics of the payload and fairing surfaces are considered. This work is supported by the NASA Glenn Research Center and the Space Launch System Payload Fairing Working Group.

Spacecraft environments↗

MRI Measurements and Granular Dynamics Simulation of Segregation of Granular Mixture

A counter intuitive axial segregation phenomenon in a rotating horizontal cylinder has recently captured attention of many researchers in different disciplines. There is a growing consensus that the interplay between the particle dynamics and the evolution of the internal structure during the segregation process must be carefully investigated. Magnetic resonance imaging (MRI) has been used to non-invasively obtain much needed dynamic/static information such as velocity and concentration profiles, and it has proven to be capable of depicting the evolution of segregation processes. Segregation in a rotating cylinder involves two processes: the first is to transport small particles in the radial direction to form a radial core, and the second is to transform the radial core into axially segregated bands. Percolation and/or "stopping" have been proposed as mechanisms for the radial segregation. As to mechanisms for axial band formation, much less is known. The difference in the dynamic angle of repose has been proposed to segregate different components in the axial direction. Recently, Hill and Kakalios have reported that particles mix or demix depending upon the competition between diffusion and preferential drift whose order can be determined by the dynamic angle of repose through the adjustment of the rotation rate. We claim that the dynamic angle of repose could be one of the causes, however, it fails to offer reasonable explanations for certain aspects of the axial migration. For example, we always observe that the radial segregation precedes the axial segregation and small particles migrate in the radial direction to form an axially extended radial core. It then transforms into axially segregated bands. By definition, the effects of the dynamic angle of repose are restricted near the free surface where the flowing layer is present. However, during the process of transforming from the radially segregated core to axially segregated bands, small particles located in the deep core region, which is untouched by the flowing layer, also completely disappear. Usually, the dynamics angle of repose are uniquely defined for individual species to characterize particle properties, and the dynamic angle of repose thus defined provides little information for the dynamic angle of repose of the mixture since the concentration ratio and the internal packing structure do not remain the same during the segregation processes. Under microgravity environment, the dynamics angle of repose argument does not hold since there is simply no flowing layer to influence/determine the preferred directions of segregation. We have thus designed an experiment so that the effects of the dynamic angle of repose can be minimized by filling the cylinder almost completely full. Small particles still formed a radial core and also migrated to form axial bands. As ground based experiments we have designed and conducted both 2D and 3D segregation experiments. The 2D experiments are performed using a thin cylinder (the gap between two end caps is about 5 mm) filled with different combinations of particles. The 3D experiments are conducted with a long cylinder of its length and diameter of 27cm and 7cm, respectively. Results of 2D experiments indicate that different mechanisms govern particle motion in regions near and far from the axis of rotation. Results of 3D experiments indicate that a series of collapses of microstructures of particle packing (micro-collapses) may be responsible for the creation of voids for small particles to migrate through in the axial direction. We have successfully eliminated the dynamic angle of repose as a cause for segregation, however, by almost completely filling the cylinder with the particles, we have lost an opportunity to investigate a possibility of particle "mobility" being a cause for segregation which requires a flowing surface but not the difference in the angle of repose. This is currently being investigated.

Nakagawa, M.↗

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗

A 2D/1D coupling neutron transport method based on the matrix MOC and NEM methods

A new 2D/1D coupling method based on the matrix MOC method (MMOC) and nodal expansion method (NEM) is proposed for solving the three-dimensional heterogeneous neutron transport problem. The MMOC method, used for radial two-dimensional calculation, constructs a response matrix between source and flux with only one sweep and then solves the linear system by using the restarted GMRES algorithm instead of the traditional trajectory sweeping process during within-group iteration for angular flux update. Long characteristics are generated by using the customization of commercial software AutoCAD. A one-dimensional diffusion calculation is carried out in the axial direction by employing the NEM method. The 2D and ID solutions are coupled through the transverse leakage items. The 3D CMFD method is used to ensure the global neutron balance and adjust the different convergence properties of the radial and axial solvers. A computational code is developed based on these theories. Two benchmarks are calculated to verify the coupling method and the code. It is observed that the corresponding numerical results agree well with references, which indicates that the new method is capable of solving the 3D heterogeneous neutron transport problem directly. (authors)

97 MATHEMATICS AND COMPUTING↗

SATURATION OF THE MAGNETO-ROTATIONAL INSTABILITY IN STRONGLY RADIATION-DOMINATED ACCRETION DISKS

The saturation level of the magneto-rotational instability (MRI) in a strongly radiation-dominated accretion disk is studied using a new Godunov radiation MHD code in the unstratified shearing box approximation. Since vertical gravity is neglected in this work, our focus is on how the MRI saturates in the optically thick mid-plane of the disk. We confirm that turbulence generated by the MRI is very compressible in the radiation-dominated regime, as found by previous calculations using the flux-limited diffusion approximation. We also find little difference in the saturation properties in calculations that use a larger horizontal domain (up to four times the vertical scale height in the radial direction). However, in strongly radiation pressure dominated disks (one in which the radiation energy density reaches {approx}1% of the rest mass energy density of the gas), we find that Maxwell stress from the MRI turbulence is larger than the value produced when radiation pressure is replaced with the same amount of gas pressure. At the same time, the ratio between Maxwell stress and Reynolds stress is increased by almost a factor of eight compared with the gas pressure dominated case. We suggest that this effect is caused by radiation drag, which acts like bulk viscosity and changes the effective magnetic Prandtl number of the fluid. Radiation viscosity significantly exceeds both the microscopic plasma viscosity and resistivity, ensuring that radiation-dominated systems occupy the high magnetic Prandtl number regime. Nevertheless, we find that radiative shear viscosity is negligible compared to the Maxwell stress and Reynolds stress in the flow. This may have important implications for the structure of radiation-dominated accretion disks.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A steady-state alignment front in an accretion disk subjected to Lense–Thirring torques

Using only physical mechanisms, i.e., 3D magnetohydrodynamics (MHD) with no phenomenological viscosity, we have simulated the dynamics of a moderately thin accretion disk subject to torques whose radial scaling mimics those produced by lowest-order post-Newtonian gravitomagnetism. In this simulation, we have shown how, in the presence of MHD turbulence, a time-steady transition can be achieved between an inner disk region aligned with the equatorial plane of the central mass’s spin and an outer region orbiting in a different plane. The position of the equilibrium orientation transition is determined by a balance between gravitomagnetic torque and warp-induced inward mixing of misaligned angular momentum from the outer disk. If the mixing is interpreted in terms of diffusive transport, the implied diffusion coefficient is ≃(0.6–0.8)c{sub s}{sup 2}/Ω for sound speed c{sub s} and orbital frequency Ω. This calibration permits estimation of the orientation transition’s equilibrium location given the central mass, its spin parameter, and the disk’s surface density and scaleheight profiles. However, the alignment front overshoots before settling into an equilibrium, signaling that a diffusive model does not fully represent the time-dependent properties of alignment fronts under these conditions. Because the precessional torque on the disk at the alignment front is always comparable to the rate at which misaligned angular momentum is brought inward to the front by warp-driven radial motions, no break forms between the inner and outer portions of the disk in our simulation. Our results also raise questions about the applicability to MHD warped disks of the traditional distinction between “bending wave” and “diffusive” regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

Turbine exhaust diffuser flow path with region of reduced total flow area

An exhaust diffuser system and method for a turbine engine includes an inner boundary and an outer boundary with a flow path defined therebetween. The inner boundary is defined at least in part by a hub that has an upstream end and a downstream end. The outer boundary has a region in which the outer boundary extends radially inward toward the hub. The region can begin at a point that is substantially aligned with the downstream end of the hub or, alternatively, at a point that is proximately upstream of the downstream end of the hub. The region directs at least a portion of an exhaust flow in the diffuser toward the hub. As a result, the exhaust diffuser system and method can achieve the performance of a long hub system while enjoying the costs of a short hub system.

42 ENGINEERING↗

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗