Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential geometry”

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 469 records · Page 26

Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations

One of the major achievements in engineering science has been the development of computer algorithms for solving nonlinear differential equations such as the Navier-Stokes equations. In the past, limited computer resources have motivated the development of efficient numerical schemes in computational fluid dynamics (CFD) utilizing structured meshes. The use of structured meshes greatly simplifies the implementation of CFD algorithms on conventional computers. Unstructured grids on the other hand offer an alternative to modeling complex geometries. Unstructured meshes have irregular connectivity and usually contain combinations of triangles, quadrilaterals, tetrahedra, and hexahedra. The generation and use of unstructured grids poses new challenges in CFD. The purpose of this note is to present recent developments in the unstructured grid generation and flow solution technology.

Barth, Timothy J.↗

Automated fiber placement composite manufacturing: The mission at MSFC's Productivity Enhancement Complex

Automated fiber placement is a manufacturing process used for producing complex composite structures. It is a notable leap to the state-of-the-art in technology for automated composite manufacturing. The fiber placement capability was established at the Marshall Space Flight Center's (MSFC) Productivity Enhancement Complex in 1992 in collaboration with Thiokol Corporation to provide materials and processes research and development, and to fabricate components for many of the Center's Programs. The Fiber Placement System (FPX) was developed as a distinct solution to problems inherent to other automated composite manufacturing systems. This equipment provides unique capabilities to build composite parts in complex 3-D shapes with concave and other asymmetrical configurations. Components with complex geometries and localized reinforcements usually require labor intensive efforts resulting in expensive, less reproducible components; the fiber placement system has the features necessary to overcome these conditions. The mechanical systems of the equipment have the motion characteristics of a filament winder and the fiber lay-up attributes of a tape laying machine, with the additional capabilities of differential tow payout speeds, compaction and cut-restart to selectively place the correct number of fibers where the design dictates. This capability will produce a repeatable process resulting in lower cost and improved quality and reliability.

Vickers, John H.↗

Fiber Optic Sensors for Smart Materials and Structures

In this paper we describe recently developed fiber sensors which are capable of monitoring the health of smart-structures. The unobstrusive geometry of these sensors make them an excellent choice for embedding the sensor in composite materials to measure internal states of strain in structures and materials. Some of these sensors have gage lengths that can be tailored from tens of microns to many meters. We will describe various demodulation schemes (Pseudo-Heterodyne, Synthetic-Heterodyne, Homodyne, Differential-Cross Multiplier, and Single Channel Phase-Tracker) to obtain high bandwidth measurements, enabling measurement of static to high frequency impact generated strains with a dynamic response exceeding tens of thousands of microstrains. In addition, we will show that we can tailor the fiber sensor to either measure only strain and reject temperature response or measure only the temperature, or measure both temperature and strain simultaneously. We will also demonstrate the ability to measure multiple strain components inside a host simultaneously using a single fiber sensor embedded in the host using a certain sensor type and transverse strain immunity using another sensor type. Additionally we will show the ability to measure temperature up to 100 C using fiber optic sensors.

Singh, H.↗

hIPPYlib-MUQ: A Bayesian Inference Software Framework for Integration of Data with Complex Predictive Models under Uncertainty

Bayesian inference provides a systematic framework for integration of data with mathematical models to quantify the uncertainty in the solution of the inverse problem. However, the solution of Bayesian inverse problems governed by complex forward models described by partial differential equations (PDEs) remains prohibitive with black-box Markov chain Monte Carlo (MCMC) methods. We present hIPPYlib-MUQ, an extensible and scalable software framework that contains implementations of state-of-the art algorithms aimed to overcome the challenges of high-dimensional, PDE-constrained Bayesian inverse problems. These algorithms accelerate MCMC sampling by exploiting the geometry and intrinsic low-dimensionality of parameter space via derivative information and low rank approximation. The software integrates two complementary open-source software packages, hIPPYlib and MUQ. hIPPYlib solves PDE-constrained inverse problems using automatically-generated adjoint-based derivatives, but it lacks full Bayesian capabilities. MUQ provides a spectrum of powerful Bayesian inversion models and algorithms, but expects forward models to come equipped with gradients and Hessians to permit large-scale solution. By combining these two complementary libraries, we created a robust, scalable, and efficient software framework that realizes the benefits of each and allows us to tackle complex large-scale Bayesian inverse problems across a broad spectrum of scientific and engineering disciplines. To illustrate the capabilities of hIPPYlib-MUQ, we present a comparison of a number of MCMC methods available in the integrated software on several high-dimensional Bayesian inverse problems. These include problems characterized by both linear and nonlinear PDEs, various noise models, and different parameter dimensions. The results demonstrate that large (~ 50×) speedups over conventional black box and gradient-based MCMC algorithms can be obtained by exploiting Hessian information (from the log-posterior), underscoring the power of the integrated hIPPYlib-MUQ framework.

97 MATHEMATICS AND COMPUTING↗

First Temperature Profile of a Stellar Flare Using Differential Chromatic Refraction

We present the first derivation of a stellar flare temperature profile from single-band photometry. Stellar flare DWF 030225.574−545707.45129 was detected in 2015 by the Dark Energy Camera as part of the Deeper, Wider, Faster program. The brightness (Δm g = −6.12) of this flare, combined with the high air mass (1.45 ≲ X ≲ 1.75) and blue filter (DES g, 398–548 nm) in which it was observed, provided ideal conditions to measure the zenithward apparent motion of the source due to differential chromatic refraction (DCR) and, from that, infer the effective temperature of the event. We model the flare’s spectral energy distribution as a blackbody to produce the constraints on flare temperature and geometric properties derived from single-band photometry. We additionally demonstrate how simplistic assumptions on the flaring spectrum, as well as on the evolution of flare geometry, can result in solutions that overestimate the effective temperature. Exploiting DCR enables studying chromatic phenomena with ground-based astrophysical surveys, and stellar flares on M dwarfs are a particularly enticing target for such studies due to their ubiquity across the sky and the heightened color contrast between their red quiescent photospheres and the blue flare emission. Our novel method will enable similar temperature constraints for a large sample of objects in upcoming photometric surveys like the Vera C. Rubin Legacy Survey of Space and Time.

Clarke, Riley W. [University of Delaware, Newark, ↗

A viscous-inviscid interaction algorithm for three-dimensional turbulent subsonic aerodynamic juncture region flow

An order-of magnitude analysis of the subsonic, three-dimensional, steady time-averaged Navier-Stokes equations for semibounded aerodynamic juncture geometries yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three-dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement first-order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. The algorithm is employed to predict the three-dimensional turbulent flowfield in the juncture region formed by two intersecting parabolic arcs at a free-stream Mach-number of 0.08, and Re/C = 600,000/m.

Baker, A. J.↗

Towards Robust Designs Via Multiple-Objective Optimization Methods

Fabricating and operating complex systems involves dealing with uncertainty in the relevant variables. In the case of aircraft, flow conditions are subject to change during operation. Efficiency and engine noise may be different from the expected values because of manufacturing tolerances and normal wear and tear. Engine components may have a shorter life than expected because of manufacturing tolerances. In spite of the important effect of operating- and manufacturing-uncertainty on the performance and expected life of the component or system, traditional aerodynamic shape optimization has focused on obtaining the best design given a set of deterministic flow conditions. Clearly it is important to both maintain near-optimal performance levels at off-design operating conditions, and, ensure that performance does not degrade appreciably when the component shape differs from the optimal shape due to manufacturing tolerances and normal wear and tear. These requirements naturally lead to the idea of robust optimal design wherein the concept of robustness to various perturbations is built into the design optimization procedure. The basic ideas involved in robust optimal design will be included in this lecture. The imposition of the additional requirement of robustness results in a multiple-objective optimization problem requiring appropriate solution procedures. Typically the costs associated with multiple-objective optimization are substantial. Therefore efficient multiple-objective optimization procedures are crucial to the rapid deployment of the principles of robust design in industry. Hence the companion set of lecture notes (Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks ) deals with methodology for solving multiple-objective Optimization problems efficiently, reliably and with little user intervention. Applications of the methodologies presented in the companion lecture to robust design will be included here. The evolutionary method (DE) is first used to solve a relatively difficult problem in extended surface heat transfer wherein optimal fin geometries are obtained for different safe operating base temperatures. The objective of maximizing the safe operating base temperature range is in direct conflict with the objective of maximizing fin heat transfer. This problem is a good example of achieving robustness in the context of changing operating conditions. The evolutionary method is then used to design a turbine airfoil; the two objectives being reduced sensitivity of the pressure distribution to small changes in the airfoil shape and the maximization of the trailing edge wedge angle with the consequent increase in airfoil thickness and strength. This is a relevant example of achieving robustness to manufacturing tolerances and wear and tear in the presence of other objectives.

Man Mohan, Rai↗

LET spectra measurements on LDEF: variations with shielding and location

LET spectra measurements made with passive plastic nuclear track detectors (PNTDs) were found to depend on detector orientation, shielding and experiment location. LET spectra were measured at several locations on LDEF as part of the P0006 LETSME experiment (Benton and Parnell, 1984), the P0004 Seeds in Space experiment (Parks and Alston, 1984), the A00l5 Free Flyer Biostacks and the M0004 Fiber Optics Data Link experiment (Taylor, 1984). Locations included the east, west and Earth sides of the LDEF satellite. The LET spectra measured with PNTDs deviated significantly from calculations, especially for high LET particles (LET infinity H2O > or = 100 keV/micrometer). At high LETs, short-range inelastic secondary particles produced by trapped proton interactions with the nuclei of the detector were found to be the principal contributor to LET spectra. At lower LETs, the spectra appeared to be due to short-range, inelastic and stopping primary protons, with primary GCR particles making a smaller contribution. The dependence of LET spectra on detector orientation and shielding was studied using the four orthogonal stacks in the P0006 experiment. Both measurements of total track density and LET spectra showed a greater number of particles arriving from the direction of space than from Earth. Measurements of LET spectra in CR-39 PNTD on the east (leading) and west (trailing) sides of LDEF showed a higher rate of production at the west side. This was caused by a larger flux of trapped protons on the west side as predicted by the east/west trapped proton anisotropy in the South Atlantic Anomaly (SAA). Track density measured in CR-39 PNTDs increased as a function of shielding depth in the detector stack. A similar measurement made in a thick stack of CR-39 interspersed with layers of Al and exposed to 154 MeV protons at a ground-based accelerator showed a similar result, indicating that a significant fraction of the particle events counted were from secondaries and that the total cross-section for production of proton-induced secondaries increased as the energy of primary protons attenuated. Little change was seen in either total differential or integral LET spectra as a function of shielding depth, indicating that the increase in cross section with decreasing proton energy affected mostly the shorter range secondary components. Similarity in the slopes of LET spectra from ground-based proton exposures and the A00l5 LET spectra showed that modeling of a monoenergetic proton beam transported through a 1-D geometry was a useful first step in modeling the production of secondary particles by trapped protons in the SAA.

long duration↗

Influence of Humidity on the Aerosol Scattering Coefficient and Its Effect on the Upwelling Radiance During ACE-2

Aerosol scattering coefficients (sigma(sub sp)) have been measured over the ocean at different relative humidities (RH) as a function of attitude in the region surrounding the Canary Islands during the Second Aerosol Characterization Experiment (ACE-2) in June and July 1997. The data were collected by the University of Washington passive humidigraph (UWPH) mounted on the Pelican research aircraft. Concurrently, particle size distributions absorption coefficients and aerosol optical depth were measured throughout 17 flights. A parameterization of sigma(sub sp) as a function of RH was utilized to assess the impact of aerosol hydration on the upwelling radiance (normalized to the solar constant and cosine of zenith angle). The top of the atmosphere radiance signal was simulated at wavelengths corresponding to visible and near-infrared bands or the EOS-AM ("Terra") detectors, MODIS and MISR. The UWPH measured (sigma(sub sp)) at 2 RHs, one below and the other above ambient conditions. Ambient (sigma(sub sp)) was obtained by interpolation of these 2 measurements. The data were stratified in terms of 3 types of aerosols: Saharan dust, clean marine (marine boundary layer background) and polluted marine aerosols (i.e., 2- or 1-day old polluted aerosols advected from Europe). An empirical relation for the dependence of (sigma(sub sp)) on RH, defined by (sigma(sub sp))(RH) = k. ((1 - RH/100)(exp -gamma), was used with the hygroscopic exponent gamma derived from the data. The following gamma values were obtained for the 3 aerosol types: gamma(dust) = 0.23 +/- 0.05, gamma(clean marine) = 0.69 +/- 0.06 and gamma(polluted marine) = 0.57 + 0.06. Based on the measured (gamma)(s), the above equation was utilized to derive aerosol models with different hygroscopicities. The satellite simulation signal code 6S was used to compute the upwelling radiance corresponding to each of those aerosol models at several ambient humidities. For the pre-launch estimated precision of the sensors and the assumed viewing geometry of the instrument, the simulations suggest that the spectral and angular dependence of the reflectance measured by MISR is not sufficient to distinguish aerosol models with various different combinations of values ror dry composition. y and ambient RH. A similar behavior is observed for MODIS at visible wavelengths. However, the 2100 nm band of MODIS appears to be able to differentiate between at least some aerosol models with different aerosol hygroscopicity given the MODIS calibration error requirements. This result suggests the possibility of retrieval of aerosol hygroscopicity by MODIS.

Gasso, B. S.↗

Influence of Humidity On the Aerosol Scattering Coefficient and Its Effect on the Upwelling Radiance During ACE-2

Aerosol scattering coefficients (sigma(sub sp)) have been measured over the ocean at different relative humidities (RH) as a function of altitude in the region surrounding the Canary Islands during the Second Aerosol Characterization Experiment (ACE-2) in June and July 1997. The data were collected by the University of Washington passive humidigraph (UWPH) mounted on the Pelican research aircraft. Concurrently, particle size distributions, absorption coefficients and aerosol optical depth were measured throughout 17 flights. A parameterization of sigma(sub sp) as a function of RH was utilized to assess the impact of aerosol hydration on the upwelling radiance (normalized to the solar constant and cosine of zenith angle). The top of the atmosphere radiance signal was simulated at wavelengths corresponding to visible and near-infrared bands of the EOS (Earth Observing System) AM-1 (Terra) detectors, MODIS (Moderate Resolution Imaging Spectroradiometer) and MISR (Multi-angle Imaging Spectroradiometer). The UWPH measured sigma(sub sp) at two RHs, one below and the other above ambient conditions. Ambient sigma(sub sp) was obtained by interpolation of these two measurements. The data were stratified in terms of three types of aerosols: Saharan dust, clean marine (marine boundary layer background) and polluted marine aerosols (i.e., two- or one-day old polluted aerosols advected from Europe). An empirical relation for the dependence of sigma(sub sp) on RH, defined by sigma(sub sp)(RH) = k.(1 - RH/100)(sup gamma), was used with the hygroscopic exponent gamma derived from the data. The following gamma values were obtained for the 3 aerosol types: gamma(dust) = 0.23 +/- 0.05, gamma(clean marine) = 0.69 +/- 0.06 and gamma(polluted marine) = 0.57 +/- 0.06. Based on the measured gammas, the above equation was utilized to derive aerosol models with different hygroscopicities. The satellite simulation signal code 6S was used to compute the upwelling radiance corresponding to each of those aerosol models at several ambient humidities. For the prelaunch estimated precision of the sensors and the assumed viewing geometry of the instrument, the simulations suggest that the spectral and angular dependence of the reflectance measured by MISR is not sufficient to distinguish aerosol models with various different combinations of values for dry composition, gamma and ambient RH. A similar behavior is observed for MODIS at visible wavelengths. However, the 2100 nm band of MODIS appears to be able to differentiate between at least same aerosol models with different aerosol hygroscopicity given the MODIS calibration error requirements. This result suggests the possibility of retrieval of aerosol hygroscopicity by MODIS.

Gasso, S.↗

Characterization of Pinhole Collimators for High-Resolution Gamma Imaging of Irradiated Fuel

Post-irradiation examination (PIE) of nuclear fuels requires imaging tools capable of resolving isotopic and spatial features with high throughput. This project contributes to a proof-of-concept effort aimed at advancing gamma emission tomography (GET) by evaluating novel fine-aperture pinhole collimators. Two Rose’s metal collimators, 100 µm (20° acceptance angle) and 350 µm (30° acceptance angle), were prototyped and characterized for their effectiveness in transporting gamma rays through the pinhole aperture. To support data collection, a Python-based data acquisition system was developed to coordinate a rotation stage, linear stage, and CZT detector, reducing latency in high-rate gamma event logging to one second per acquisition. Queue-based file writing enabled seamless real-time data capture for count rates up to 35,000 counts per second (cps). List-mode parsing algorithms were implemented to differentiate single and simultaneous gamma interactions for future tomographic reconstruction. Detector response was evaluated in both spectroscopy and list mode acquisition methods across varying source distances to confirm absolute and collimator efficiencies. Preliminary efficiency figures suggest effective collimation of gamma-rays with energies below 700 keV, with ~4% residual intensity through the aperture for Cs-137. The impact of collimator geometry on image quality is currently being evaluated. This groundwork supports the ongoing development of a sub mm resolution cone-beam CT system for imaging fuel phantoms, an essential step toward improving GET efficiency and accelerating nuclear fuel qualification efforts.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

The Omics of Stem Cell Mediated Regeneration: A Pilot Single Cell RNA-Seq Study of Mechanotransduction

Mechanical forces are potent modulators of stem cell based tissue regenerative mechanisms, inducing cell fate decisions and tissue specific commitment. A unique platform for investigating mechanotransduction is spaceflight, where microgravity and altered fluid mechanics provide a loading-null experimental condition. Seminal investigations of regenerative capacity in a wholly regenerative species, the newt model, and in a variety of totipotent and adult stem cell populations have demonstrated the detrimental effects of unloading on maintenance of stem cell based regeneration. Of particular interest is the observation that unloading interferes with the transition of stem cell pools from proliferative state to differentiation commitment. In this work we sought to test the hypothesis that gravity mechanotransduction regulates stem cell tissue regenerative processes by modulating stem cell proliferation and differentiation fates at specific cell cycle stages. To do this, clonally-derived ESCs were plated on a collagen matrix and expanded for 36 hours before re-plating on a non-adherent culture dish in the absence of leukemia inhibitory factor (LIF) to form spheroid aggregate EBs. After formation, the EBs were transferred to a collagen matrix coated culture dishes and given 4 days to allow implantation and outgrowth. In parallel, totipotent ESCs were plated 24 hours before mechanical stimulation on collagen matrix culture dishes in the presence of LIF to maintain totipotency and serve as un-differentiation committed controls. The EBs and ESCs were then subjected to either a 60 minute pulse of gravity (static loading) or 60 minutes of cyclic stretch (dynamic loading) mechanotransduction. Six hours post-stimulation, we used a 10X Genomics Single Cell controller to generate bar-coded single cell Illumina libraries and sequenced expressomes for 5,000 static loaded cells, representative of a change in gravity mechanotransduction, 5,000 dynamic loaded cells, representative of tissue loading associate with physiologic function, and 5,000 unstimulated 1g control cells. The comparison of these 3 libraries by cluster assignment based on like gene expression patterns show substantial alteration in cluster geometry due to mechanical loading. Specifically the mechanically loaded EB outgrowth cells to retain potency markers (PAX6, SOX2, CD34) and suppress early commitment markers (Dhh, VCAN, Igf1). Whereas the EBs cultured under the non-stimulated conditions display clear departure from the ESC expressome with lineage commitment markers upregulated and several tissue specific markers being expressed (BMP "early musculoskeletal development, Mesp1" early cardiovascular cell lineage). These markers are not seen in the mechano-stimulated cultures or the totipotent ESC cultures. Comparison of like clusters between our experimental conditions revealed an array of regenerative and stem cell genes are significantly mechano-regulated. Of particular importance CDKN1a/p21, a gene shown by previous investigation of our research team to be significantly upregulated in unloading, was suppressed in the static and dynamic loaded EBS. In addition to CDKN1a/p21 many genes related to cell cycle and transitory differentiation markers had elevated expression in the mechano-stimulated EBs, but surprisingly these trends were not observed in the ESC cultures. This study is the first of its kind investigating for mechano-signaling and mechano-regulated pathways, and has alre

Single cell RNA-sequecing↗

A Generalized Stress Inversion Approach With Application to Residual Stress Estimation

In this work develop a generalized stress inversion technique (or the generalized inversion method) capable of recovering stresses in linear elastic bodies subjected to arbitrary cuts. Specifically, given a set of displacement measurements found experimentally from digital image correlation (DIC), we formulate a stress estimation inverse problem as a partial differential equation-constrained optimization problem. We use gradient-based optimization methods, and we accordingly derive the necessary gradient and Hessian information in a matrix-free form to allow for parallel, large-scale operations. By using a combination of finite elements, DIC, and a matrix-free optimization framework, the generalized inversion method can be used on any arbitrary geometry, provided that the DIC camera can view a sufficient part of the surface. We present numerical simulations and experiments, and we demonstrate that the generalized inversion method can be applied to estimate residual stress.

42 ENGINEERING↗

Hover performance tests of full scale variable geometry rotors

Full scale whirl tests were conducted to determine the effects of interblade spatial relationships and pitch variations on the hover performance and acoustic signature of a 6-blade main rotor system. The variable geometry rotor (VGR) variations from the conventional baseline were accomplished by: (1) shifting the axial position of alternate blades by one chord-length to form two tip path planes; and (2) varying the relative azimuthal spacing from the upper rotor to the lagging hover rotor in four increments from 25.2 degrees to 62.1 degrees. For each of these four configurations, the differential collective pitch between upper and lower rotors was set at + or - 1 deg, 0 deg and -1 deg. Hover performance data for all configurations were acquired at blade tip Mach numbers of 0.523 and 0.45. Acoustic data were recorded at all test conditions, but analyzed only at 0 deg differential pitch at the higher rotor speed. The VGR configurations tested demonstrated improvements in thrust at constant power as high as 6 percent. Reductions of 3 PNdb in perceived noise level and of 4 db in blade passage frequency noise level were achieved at the higher thrust levels. Consistent correlation exists between performance and acoustic improvements. For any given azimuth spacing, performance was consistently better for the differential pitch condition of + or - 1 degree, i.e. with the upper rotor pitch one degree higher than the lower rotor.

Rorke, J. B.↗

Quantification of Aerosol Transmission through Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using more significant backfill pressures in the canister, up to approximately 800 kPa. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO 2 , were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Tensor Network Space-Time Spectral Collocation Method for Time-Dependent Convection-Diffusion-Reaction Equations

Emerging tensor network techniques for solutions of partial differential equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultra-fast numerical solutions of high-dimensional problems. Here, we introduce a Tensor Train (TT) Chebyshev spectral collocation method, in both space and time, for the solution of the time-dependent convection-diffusion-reaction (CDR) equation with inhomogeneous boundary conditions, in Cartesian geometry. Previous methods for numerical solution of time-dependent PDEs often used finite difference for time, and a spectral scheme for the spatial dimensions, which led to a slow linear convergence. Spectral collocation space-time methods show exponential convergence; however, for realistic problems they need to solve large four-dimensional systems. We overcome this difficulty by using a TT approach, as its complexity only grows linearly with the number of dimensions. We show that our TT space-time Chebyshev spectral collocation method converges exponentially, when the solution of the CDR is smooth, and demonstrate that it leads to a very high compression of linear operators from terabytes to kilobytes in TT-format, and a speedup of tens of thousands of times when compared to a full-grid space-time spectral method. These advantages allow us to obtain the solutions at much higher resolutions.

97 MATHEMATICS AND COMPUTING↗

Axially grooved heat pipe study

A technology evaluation study on axially grooved heat pipes is presented. The state-of-the-art is reviewed and present and future requirements are identified. Analytical models, the Groove Analysis Program (GAP) and a closed form solution, were developed to facilitate parametric performance evaluations. GAP provides a numerical solution of the differential equations which govern the hydrodynamic flow. The model accounts for liquid recession, liquid/vapor shear interaction, puddle flow as well as laminar and turbulent vapor flow conditions. The closed form solution was developed to reduce computation time and complexity in parametric evaluations. It is applicable to laminar and ideal charge conditions, liquid/vapor shear interaction, and an empirical liquid flow factor which accounts for groove geometry and liquid recession effects. The validity of the closed form solution is verified by comparison with GAP predictions and measured data.

Source record↗

A statistical rain attenuation prediction model with application to the advanced communication technology satellite project. Part 2: Theoretical development of a dynamic model and application to rain fade durations and tolerable control delays for fade countermeasures

A dynamic rain attenuation prediction model is developed for use in obtaining the temporal characteristics, on time scales of minutes or hours, of satellite communication link availability. Analagous to the associated static rain attenuation model, which yields yearly attenuation predictions, this dynamic model is applicable at any location in the world that is characterized by the static rain attenuation statistics peculiar to the geometry of the satellite link and the rain statistics of the location. Such statistics are calculated by employing the formalism of Part I of this report. In fact, the dynamic model presented here is an extension of the static model and reduces to the static model in the appropriate limit. By assuming that rain attenuation is dynamically described by a first-order stochastic differential equation in time and that this random attenuation process is a Markov process, an expression for the associated transition probability is obtained by solving the related forward Kolmogorov equation. This transition probability is then used to obtain such temporal rain attenuation statistics as attenuation durations and allowable attenuation margins versus control system delay.

Manning, Robert M.↗