SEARCH · Engineering Papers

Results for “Conditional diffusion models”

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.

SourceResearch connectionFrom yearThrough year

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

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Regional-scale analysis of karst underground flow deduced from tracing experiments: examples from carbonate aquifers in Malaga province, southern Spain

Tracer concentration data from field experiments conducted in several carbonate aquifers (Malaga province, southern Spain) were analyzed following a dual approach based on the graphical evaluation method (GEM) and solute transport modeling to decipher flow mechanisms in karst systems at regional scale. The results show that conduit system geometry and flow conditions are the principal factors influencing tracer migration through the examined karst flow routes. Solute transport is mainly controlled by longitudinal advection and dispersion throughout the conduit length, but also by flow partitioning between mobile and immobile fluid phases, while the matrix diffusion process appears to be less relevant. The simulation of tracer breakthrough curves (BTCs) suggests that diffuse and concentrated flow through the unsaturated zone can have equivalent transport properties under extreme recharge, with high flow velocities and efficient mixing due to the high hydraulic gradients generated. Tracer mobilization within the saturated zone under low flow conditions mainly depends on the hydrodynamics (rather than on the karst conduit development), which promote a lower longitudinal advection and retardation in the tracer migration, resulting in a marked tailing effect of BTCs. The analytical advection-dispersion equation better approximates the effective flow velocity and longitudinal dispersion estimations provided by the GEM, while the non-equilibrium transport model achieves a better adjustment of most asymmetric and long-tailed BTCs. The assessment of karst underground flow properties from tracing tests at regional scale can aid design of groundwater management and protection strategies, particularly in large hydrogeological systems (i.e. transboundary carbonate aquifers) and/or in poorly investigated ones.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Evaluation of Spectral Band Adjustment Factors for Cross-Calibration of Visible Imagers

The CERES EBAF dataset provides TOA SW and LW fluxes for long-term monitoring of the Earth’s energy balance and to validate climate models. The EBAF products, based on the Terra and Aqua CERES instrument observed radiances, rely on coincident measurements from the onboard MODIS imagers to determine cloud properties that are used in the angular distribution model scene selection required to convert the CERES observed radiances into flux values. Once the Terra and Aqua orbits start drifting outside of their 15-minute window, the CERES project will rely solely on the NOAA-20 CERES observations. A seamless transition of fluxes and clouds can only occur if the analogous MODIS and VIIRS channels are properly inter-calibrated. The visible band (0.65 µm) in MODIS (band 1) and VIIRS (M5 or I1) is critical for retrieving the cloud mask and optical depth. The spectral response functions (SRF) of these bands differ noticeably and will require scene dependent spectral band adjustment factors (SBAF) for proper radiometric scaling between them. Fortunately, the VIIRS I1 and M5 bands are calibrated using the same solar diffuser as a reference, which implies that band reflectance differences must be due to the spectral disparity. The coincident I1 and M5 TOA reflectance measurements provide the optimal opportunity to validate SBAFs over many surface and cloud conditions. The CERES project maintains SCIAMACHY, GOME-2, Hyperion, and radiative transfer model-based scene-stratified hyper-spectral reflectance measurements that can be convolved with sensor pair SRFs to compute the corresponding SBAF. This study will evaluate the SCIAMACHY, GOME-2, Hyperion, and RTM based SBAFs over tropical ocean targets used in SNO inter-calibration of MODIS and VIIRS, including clear-sky ocean, deep convective clouds (DCC), and liquid cloud targets in terms of their applicability to absolute intercomparison of visible imagers. The large SCIAMACHY and GOME-2 footprints and the lack of Hyperion cloudy scenes will impact the resulting SBAFs and necessitates the need for future CLARREO hyper-spectral reflectances to decrease SBAF uncertainties Preliminary results indicate that the SCIAMACHY and GOME-2 based SBAFs for the VIIRS I1 and M5 band pairs may differ by 1.5% for some ATO scene types. Because most of the modern GEOs visible band SRFs encompass the VIIRS I1 band SRF, these evaluations are critical to ensure that the MODIS, VIIRS, and GEO clouds and fluxes are consistent. SBAFs for Libya-4, Dome-C and other Earth invariant target will also be evaluated.

D. R. Doelling↗

DOE OSTI2025

Preliminary Plan to Inform Testing of a Heat Exchanger Test Article

This report presents a preliminary plan to guide the qualification testing of advanced heat exchanger (HX) components for nuclear-to-industrial heat transfer applications. The objective is to establish a defensible, physics-based methodology that integrates computational modeling, targeted experimentation, and in-service inspection considerations to demonstrate component performance and reliability under representative reactor conditions. The analysis identifies Sodium-cooled Fast Reactor (SFR) and High-Temperature Gas-cooled Reactor (HTGR) systems as reference configurations in terms of temperature, pressure, and chemical environment. Within these operating envelopes, dominant degradation mechanisms— including creep–fatigue interaction, flow-induced vibration, corrosion, and diffusion-bond deterioration—were evaluated to define test requirements. A comprehensive computationalexperimental framework is proposed to support life prediction and qualification activities. The framework couples high-fidelity structural-mechanics, thermal-hydraulic, and fluid-structure interaction models with accelerated degradation testing to produce a traceable linkage between microstructural evolution, mechanical performance, and remaining useful life (RUL). The approach adheres to established Verification, Validation, and Uncertainty Quantification (VVUQ) standards (ASME V&V 10/20; NUREG-2152) and incorporates a digital-twin architecture for continuous model refinement through data assimilation. The plan further outlines testing methodologies, including pre-test analyses, test-loop design parameters, and sensor placement strategies that maximize information yield while maintaining mechanistic fidelity. Complementary sections describe in-service inspection (ISI), on-line monitoring (OLM), and structural-health-monitoring (SHM) techniques applicable to compact HX geometries typical of advanced reactors. Collectively, these activities establish the technical foundation for demonstrating 40-60-year equivalent service life of advanced heat exchangers in support of the U.S. Department of Energy’s Advanced Reactor and Integrated Energy Systems programs. The forthcoming phase will execute the defined pre-test analyses, initiate hardware fabrication, and implement the integrated testing campaign to validate the proposed qualification methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1993

Hypersonic panel flutter in a rarefied atmosphere

Panel flutter is a form of dynamic aeroelastic instability resulting from the interaction between motion of an aircraft structural panel and the aerodynamic loads exerted on that panel by air flowing past one of the faces. It differs from lifting surface flutter in the sense that it is not usually catastrophic, the panel's motion being limited by nonlinear membrane stresses produced by the transverse displacement. Above some critical airflow condition, the linear instability grows to a limit cycle . The present investigation studies panel flutter in an aerodynamic regime known as 'free molecule flow', wherein intermolecular collisions can be neglected and loads are caused by interactions between individual molecules and the bounding surface. After collision with the panel, molecules may be reflected specularly or reemitted in diffuse fashion. Two parameters characterize this process: the 'momentum accommodation coefficient', which is the fraction of the specularly reflected molecules; and the ratio between the panel temperature and that of the free airstream. This model is relevant to the case of hypersonic flight vehicles traveling at very high altitudes and especially for panels oriented parallel to the airstream or in the vehicle's lee. Under these conditions the aerodynamic shear stress turns out to be considerably larger than the surface pressures, and shear effects must be included in the model. This is accomplished by means of distributed longitudinal and bending loads. The former can cause the panel to buckle. In the example of a simply-supported panel, it turns out that the second mode of free vibration tends to dominate the flutter solution, which is carried out by a Galerkin analysis. Several parametric studies are presented. They include the effects of (1) temperature ratio; (2) momentum accommodation coefficient; (3) spring parameters, which are associated with how the panel is connected to adjacent structures; (4) a parameter which relates compressive end load to its value which would cause classical column buckling; (5) a parameter proportional to the pressure differential between the front and back faces; and (6) initial curvature. The research is completed by an investigation into the possibility of accounting for molecular collisions, which proves to be infeasible given the speeds of current mainframe supercomputers.

Resende, Hugo B.↗

DOE OSTI2026

SOLPS-ITER modeling of a dedicated divertor for negative triangularity operation on DIII-D

The design of a new dedicated divertor for negative triangularity (NT) operation on DIII-D with neutral baffles and pumping is informed by SOLPS-ITER transport modeling. This dedicated NT divertor is the latest step in a progression of NT shapes with various divertor characteristics explored on DIII-D, including NT shapes at reduced triangularity and a campaign with stronger shaping that included new armored components on the outboard side. SOLPS simulations played a key role in these divertor designs. Interpretive simulations, using cross-field diffusivities constrained by experimental data in the NT Shelf shape were used to inform the design of the 2023 armor campaign components. A similar procedure used armor campaign data to predict conditions for the dedicated NT divertor. The predictive simulations were used to assess the divertor fluxes, detachment threshold, pumped flux, and neutral leakage. For the dedicated NT divertor, SOLPS simulations and two-point-modeling were used to show the relative impact of magnetic topology (mainly longer connection length) and divertor closure on the divertor conditions relative to the armor campaign. It is predicted that the dedicated NT divertor reaches detachment (measured by target ion flux rollover) at a lower upstream density (≈(1.75−−2.4)×1⁢019m−3) as compared to the armor campaign shape. For the preliminary design geometry, divertor closure reduces the neutral leakage by ≈10%. Parametric optimization indicating further ≈20%–60% improvement in the leakage flux and recycled flux crossing the pump entrance is possible for relatively minor changes to the divertor and baffle layout.

Lore, Jeremy [ORNL] (ORCID:000000029192465X)↗

NASA NTRS1984

A one-dimensional photochemical model of the troposphere with planetary boundary-layer parameterization

A time-dependent, one-dimensional photochemical model of the troposphere is used to describe the vertical distribution of atmospheric trace constituents for summer-time conditions at midlatitudes in the Northern Hemisphere. The model incorporates a planetary boundary layer (PBL) parametrization and a detailed chemical mechanism that includes the photochemistry of important nonmethane hydrocarbon species formed during the oxidation process. One result of the parametrized PBL is that the concentrations of some trace species in the free troposphere are 20-30 percent higher than when mixing processes are described by a vertical eddy diffusion coefficient which is held constant with respect to height and time. The lifetime of the oxides of nitrogen against photochemical conversion to nitric acid during summertime conditions is on the order of six hours. This lifetime is short enough to deplete most of the NO(x) in the PBL so that other reactive nitrogen species are more abundant than NO(x) throughout the free troposphere.

Fishman, J.↗

DOE OSTI2017

Effects of alpha stopping power modelling on the ignition threshold in a directly-driven inertial confinement fusion capsule

The alpha-particle energy deposition mechanism modifies the ignition conditions of the thermonuclear Deuterium-Tritium fusion reactions, and constitutes a key issue in achieving high gain in Inertial Confinement Fusion implosions. One-dimensional hydrodynamic calculations have been performed with the code Multi-IFE to simulate the implosion of a capsule directly irradiated by a laser beam. The diffusion approximation for the alpha energy deposition has been used to optimize three laser profiles corresponding to different implosion velocities. A Monte-Carlo package has been included in Multi-IFE to calculate the alpha energy transport, and in this case the energy deposition uses both the LP and the BPS stopping power models. Homothetic transformations that maintain a constant implosion velocity have been used to map out the transition region between marginally-igniting and high-gain configurations. Furthermore, the results provided by the two models have been compared and it is found that – close to the ignition threshold – in order to produce the same fusion energy, the calculations performed with the BPS model require about 10% more invested energy with respect to the LP model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Implications of a two-fluid heat diffusion model for RF condensation effects in magnetic islands

Theoretical calculations have predicted an “RF condensation” effect through which the temperature perturbation and current drive may be nonlinearly enhanced in electron cyclotron (EC) heated magnetic islands. Previously, two-fluid effects had been neglected in the OCCAMI code used to model RF condensation. For realistic experimental conditions, we show that two-fluid effects in an EC heated plasma may have a major impact on RF condensation. Less collisional plasmas, where electron and ion temperatures are decoupled, exhibit an enhanced linear temperature perturbation in the island, and an even stronger nonlinear enhancement of the perturbed temperature. This can occur even when nonlinear effects are small in the single-fluid limit. We also calculate corresponding hysteresis effects. Finally, we demonstrate that other physical effects, including the rotation of the magnetic island and the Gaussian spread of the beam from a gyrotron, can significantly decrease the predicted observable RF condensation effect.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

1-D coupled surface flow and transport equations revisited via the physics-informed neural network approach

The de Saint-Venant equation (SVE) and advection–diffusion equation (ADE) are commonly employed to solve solute transport problems in surface water. In this work, we propose a mesh-free method based on the physics-informed neural network (PINN) to solve the one dimensional (1-D) SVE, ADE, and the coupled SVE and ADE (SVE-ADE) under various initial and boundary conditions. The PINN model extends the architecture of deep neural networks (DNNs) with implementation of loss function, which are additionally subject to constraints imposed by the physical laws of SVE and ADE, along with their initial and boundary conditions. In such a manner, PINNs can be quickly steered to the true solution while obeying the physical laws. The results of PINN model are compared with the analytical and/or numerical solutions under various conditions to investigate its accuracy and efficiency in solving the SVE, ADE, and SVE-ADE. Our results indicate PINN can accurately simulate the shock wave morphology and avoid numerical dissipation in unsteady flow condition. The PINN method outweighed traditional numerical methods in several aspects, including its ability to function with small amounts of data, no grid discretization, and random selection of sampling points, etc. Additionally, the PINN method is also suitable for solving inverse problems with sparse and noisy data. With 1% noise and 2000 initial and boundary condition points (N u ), the errors of the estimated flow rate (v) and diffusion coefficient (D) are 0.003% and 0.105%, respectively, which indicate the accuracy and robustness of the proposed method. Finally, our results indicate the capability and robustness of the proposed PINN methodology for solving multi-physics problems, irrespective of the presence of sparse and noisy data in the training dataset.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

Counterflow ignition of n-butanol at atmospheric and elevated pressures

Vital to the development of predictive combustion models is a robust understanding of the coupled effects of chemical kinetics and convective–diffusive transport at both atmospheric and elevated pressures. The current study describes a new variable-pressure non-premixed counterflow ignition experiment designed to address the need for well-characterized reference data to validate such models under conditions sensitive to both chemical and transport processes. A comprehensive characterization of system behavior is provided to demonstrate boundary condition and ignition quality as well as adherence to the assumption of quasi-one-dimensionality, and suggest limitations and best practices for counterflow ignition experiments. This effort reveals that the counterflow ignition experiment requires special attention to ignition location in order to ensure that the assumption of quasi-one-dimensionality is valid, particularly at elevated pressures. This experimental tool is then applied to the investigation of n-butanol for pressures of 1–4 atm, pressure-weighted strain rates of 200–400 s -1 , and fuel mole fractions of 0.05–0.25. Results are simulated using two n-butanol models available in the literature and used to validate and assess model performance. Comparison of experimental and numerical ignition results for n-butanol demonstrates that while existing models largely capture the trends observed with varying pressure, strain rate, and fuel loading, the models universally over-predict experimental ignition temperatures. While several transport coefficients are found to exhibit order-of-magnitude or greater sensitivities relative to reaction rates, variation of transport parameters is not able to account for the large deviations observed between experimental and numerical results. Additional comparison of ignition kernel structure and fuel breakdown pathways between two literature models suggests that an under-prediction in the radical pool growth with respect to temperature variation may be responsible for both the deviation from the experimental results and the discrepancy in ignition temperature results observed between models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Experimental study of ethanol oxidation behind reflected shock waves: Ignition delay time and H 2 O laser-absorption measurements

The oxidation of ethanol was investigated in shock tubes by measuring ignition delay times and water time history profiles. A laser absorption technique was used in the latter case. Ignition delay times were investigated for three equivalence ratios, 0.5, 1.0, and 2.0, and for pressures ranging from 1.3 to 53 atm. The fuel concentration and the type of diluent used were also varied. The ignition delay times were in very good agreement with results from the literature for the very few overlapping conditions, and most of the experimental conditions have never been studied before. A high level of accuracy was observed for the predictions of the ignition delay times from detailed kinetics mechanisms. Inhomogeneous ignition events were observed in some specific conditions, validating the notions of thermal diffusivity and flame thickness developed in the literature for the reasons behind such inhomogeneous ignition. The water time-history profiles showed the formation of water in two stages. These profiles were also modeled using well-known detailed kinetic mechanisms from the literature. Severe discrepancies were often observed for laser-absorption measurements, and this was the case with most of the models. A numerical analysis (sensitivity and rate-of-production) was conducted using the two best-performing literature mechanisms (namely, AramcoMech3.0 and CRECK).

42 ENGINEERING↗

DOE OSTI2017

Condensed-phase biogenic-anthropogenic interactions with implications for cold cloud formation

Anthropogenic and biogenic gas emissions contribute to the formation of secondary organic aerosol (SOA). When present, soot particles from fossil-fuel combustion can acquire a coating of SOA. We investigate SOA-soot biogenic-anthropogenic interactions and their impact on ice nucleation in relation to the particles’ organic phase state. SOA particles were generated from the OH oxidation of naphthalene, α-pinene, longifolene, or isoprene, with or without presence of sulfate or soot particles. Corresponding particle glass transition (T g ) and full deliquescence relative humidity (FDRH) were estimated by a numerical diffusion model. Longifolene SOA particles are solid-like and all biogenic SOA sulfate mixtures exhibit a core-shell configuration (i.e. a sulfate-rich core coated with SOA). Biogenic SOA with or without sulfate formed ice at conditions expected for homogeneous ice nucleation in agreement with respective T g and FDRH. α-pinene SOA coated soot particles nucleated ice above the homogeneous freezing temperature with soot acting as ice nuclei (IN). At lower temperatures the α-pinene SOA coating can be semisolid inducing ice nucleation. Naphthalene SOA coated soot particles acted as IN above and below the homogeneous freezing limit, which can be explained by the presence of a highly viscous SOA phase. Our results suggest that biogenic SOA does not play a significant role in mixed-phase cloud formation and the presence of sulfate further renders this even less likely. Furthermore, anthropogenic SOA may have an enhancing effect on cloud glaciation under mixed-phase and cirrus cloud conditions compared to biogenic SOA that dominate during preindustrial times or in pristine areas.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Mesoscopic analyses of the impact of morphology and operating conditions on the transport resistances in a proton-exchange-membrane fuel-cell catalyst layer

Exploring the origins of local transport resistance and characterizing the oxygen transport resistances in the catalyst layer (RCL) are critical for cost reduction. In this paper, a comprehensive mesoscopic model for simulating coupled transport processes of oxygen and water vapor for different structural parameters under different operating conditions in reconstructed microstructures is proposed. The local transport resistance is calculated after achieving the limiting current density and the effective diffusivity of oxygen. The results demonstrate that RCL increases greatly with decreasing platinum loading (LPt) and the transport resistances in other components of the cell dominate for high-loadings. Both the reduced oxygen permeation coefficient in the ionomer thin-film and the adsorption resistance account for the origins of local transport resistance. The local transport resistance increases with the bare carbon ratio for a constant LPt and Pt/C ratio due to the decreased effective ionomer surface, and increases with the I/C ratio due to the increased ionomer thickness and decreased Knudsen diffusivity. Due to the presence of liquid water, a slight decrease followed by an increase of the local transport resistance versus relative humidity is obtained. The contribution of ionomer thin-films to RCL is more sensitive to liquid saturation compared with that of pores. This journal is

Mu, YT↗

DOE OSTI2021

Energetics of boron near tungsten surfaces: A first-principles study

Interest exists in utilizing boron (B) wall conditioning of fusion tokamaks containing tungsten (W) plasma facing components, in order to improve plasma confinement. To understand the interactions of B with W surfaces, first-principles density functional theory calculations have been performed to model the adsorption, diffusion, and solution of B near the W(100), W(110), and W(111) surfaces. The results show that B within a distance of 0.6 nm above the surfaces is adsorbed to the surfaces without activation barriers. B atoms are strongly adsorbed on the W(100) surface with an adsorption energy of 7.80 eV, which is 1.22 and 1.35 eV larger than on the W(110) and W(111) surfaces. B diffusion on the W(100), W(110), and W(111) surfaces has an activation energies of 2.08, 1.12, and 1.47 eV, respectively; while, diffusion from the adsorption sites into the bulk requires 2.2–2.3 eV. The B solution energy below a clean W(100) surface is the lowest, followed by the W(111) and W(110) surfaces. B clustering and B-induced surface deformation as a function of B coverage has been investigated. B on the W(100) surface occupy epitaxial sites at coverages of 0–1.25 ML, but form clusters at higher coverages. B clustering on the W(110) and W(111) surfaces is expected throughout the adsorption process. Compared to a clean surface, B atoms on the W(100) reduce the surface effect on the B solution energy below the surface, while the presence of B on the W(110) and W(111) surfaces generally decreases or increases the B solution energy below surfaces, respectively.

Yang, L. (ORCID:0000000322166071)↗

DOE OSTI2020

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Study of The Compression of an Imploding Plasma and the Magnetic- Field Flux Using Advanced, Non-Intrusive, Spectroscopic Techniques

The main motivation of the project was to measure the magnetic field in pulsed power driven plasma compression using advanced non-intrusive spectroscopic techniques. The strength of the project is that it encompasses a combination of laboratory experiments with diagnostics for both the magnetic field and plasma properties and computational modeling. The following issues were addressed: (1) Plasma and magnetic field dynamics, (2) Field compression and diffusion, (3) Energy balance in the system, (4) Onset of instabilities, and (5) Correlation of the plasma and field evolutions with the initial conditions. The unique feature of the research was the use of advanced spectroscopy to measure both the magnetic field and various plasma properties in multiple locations using both side-on and end-on views, simultaneously. The experimental system allowed for imploding different gas mixtures (which facilitated local measurements) and for flexibility in modifying the initial conditions. Two current generators were used with significantly different current rise times (200 ns and 1.6 µs) but at the same current level of 400 kA. The slower pulsed power generator (at rise time 1.6 µs) was based at the Weizmann Institute of Science (WIS) and imploded plasma into which an axial magnetic field (up to 1 T) was initially embedded. The faster current rise time generator (at rise time 200 ns) based at University of California San Diego (UCSD) was used for plasma parameter studies initially, with collaborative efforts we successfully implemented the advanced spectroscopy techniques at UC San Diego. One graduate student (Nick Aybar) and a postdoctoral scholar (Maylis Dozeries) were trained. They both are based at the Lawrence Livermore National Laboratory now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1973

Preliminary numerical analysis of improved gas chromatograph model

A mathematical model for the gas chromatograph was developed which incorporates the heretofore neglected transport mechanisms of intraparticle diffusion and rates of adsorption. Because a closed-form analytical solution to the model does not appear realizable, techniques for the numerical solution of the model equations are being investigated. Criteria were developed for using a finite terminal boundary condition in place of an infinite boundary condition used in analytical solution techniques. The class of weighted residual methods known as orthogonal collocation is presently being investigated and appears promising.

Woodrow, P. T.↗

NASA NTRS1998

Long Term Degradation of Polyimide Composites

The durability of polymer matrix composites exposed to harsh environments is a major concern. Surface degradation and damage are observed in polyimide composites used in air at 125-300 C. It is believed that diffusion of oxygen into the material and oxidative chemical reactions are responsible. Previous work has characterized and modeled diffusion behavior, and thermogravimetric analyses (TGAS) have been carried out in nitrogen, air, and oxygen to provide quantitative information on thermal and oxidative reactions. However, the model developed using these data did not successfully extrapolate TGA data down to conditions seen in service. A test program that focuses on lower temperatures and makes use of isothermal tests was undertaken to achieve a better understanding of the degradation reactions under use conditions. A new, low-cost technique was developed to collect chemical degradation data for isothermal tests lasting over 200 hours in the temperature range 125-300 C. Results indicate complex behavior not captured by the previous model, including the presence of weight-adding reactions. Weight gain reactions dominated in the 125-225 C temperature range, while weight loss reactions dominated beyond 225 C. The data obtained from isothermal tests and earlier TGAs is used to develop an advanced model of the material behavior.

Patekar, Kaustubh A.↗