Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation modeling”

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 271 records · Page 15

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗

Titan's Ionospheric Chemistry, Fullerenes, Oxygen, Galactic Cosmic Rays and the Formation of Exobiological Molecules on and within its Surfaces and Lakes

We discuss the formation of aerosols within Titan's thermosphere-ionosphere and the different chemical pathways. Negative ion measurements by the Cassini Plasma Spectrometer (CAPS) Electron Spectrometer (ELS) give evidence for formation of unsaturated anion carbon chains, while positive ion measurements of the Cassini Ion Neutral Mass Spectrometer (INMS) indicate formation of more aromatic cation hydrocarbons. There is presently no direct observational evidence for large neutral molecule growth in Titan's thermosphere-ionosphere. The hydrocarbon cations are expected to form Polycyclic Aromatic Hydrocarbons (PAH), those with the addition of nitrogen being called PAHNs. We theorize anion carbon chains can eventually become long enough to fold into fullerene C(60,70) carbon shells, of various charge states. Based on laboratory data the fullerenes can trap incoming O(sup +) magnetospheric ions that have relatively high energy collisions with the fullerenes and, once trapped, protect the oxygen atom from Titan's reducing thermosphere-ionosphere. The fullerenes can form into larger onion fullerenes and condense into larger embryo aerosols (i.e., m/q > 10,000 amu/q anions as observed by CAPS/ELS) eventually falling onto Titan's surface and precipitating to the bottom of its hydrocarbon lakes. Molecule production composed of H, C, N is known to occur in Titan's atmosphere with energy input from the magnetosphere, solar UV, and deeppenetrating irradiation from galactic cosmic rays (GCR). Space radiation effects by GCR irradiation of Titan's surface and lakes can lead to the manufacture of exobiological molecules with oxygen as the new ingredient. We have developed a model of galactic cosmic ray irradiation of Titan's atmosphere, surface, subsurface and bottoms of Titan lakes. GCR would provide further energy for processing of the aerosols into more complex organic forms such as tholins and precursor molecules for amino acids. A second process called hydrolysis then converts the precursor molecules into amino acids. Hydrolysis is provided via meteor impacts with size > 10 km and cryovolcanism both which can produce liquid water on Titan's surface for episodic periods > several 100 to 1000 years. Our model shows that GCR secondary particles can penetrate ~ 100 m below the ice surface (including the bottom of Titan's less dense hydrocarbon lakes ~ 150 m depths) and produce chemically significant dosages over very long timescales ~ 450 Myrs. The GCR model is combined with laboratory data from experiments in which dry methyl ices were irradiated to doses producing prebiotic amino acids such as glycine. The model calculations show glycine can form to ~ 2.5 ppb levels near the surface after ~450 Myrs of GCR proton irradiation and potentially to 5 ppb if heavy-ion GCRs up through Fe are included. If such molecules were detected, this would not only confirm this model but indicate that life forms different from ours may not be required.

Edward C. Sittler Jr↗

Understanding the Impacts of Soil Moisture Initial Conditions on NWP in the Context of Land-Atmosphere Coupling

The role of soil moisture in NWP has gained more attention in recent years, as studies have demonstrated impacts of land surface states on ambient weather from diurnal to seasonal scales. However, soil moisture initialization approaches in coupled models remain quite diverse in terms of their complexity and observational roots, while assessment using bulk forecast statistics can be simplistic and misleading. In this study, a suite of soil moisture initialization approaches is used to generate short-term coupled forecasts over the U.S. Southern Great Plains using NASA’s Land Information System (LIS) and NASA Unified WRF (NU-WRF) modeling systems. This includes a wide range of currently used initialization approaches, including soil moisture derived from ‘off-the-shelf’ products such as atmospheric models and land data assimilation systems, high-resolution land surface model spinups, and satellite-based soil moisture products from SMAP. Results indicate that the spread across initialization approaches can be quite large in terms of soil moisture conditions and spatial resolution, and that SMAP performs well in terms of heterogeneity and temporal dynamics when compared against high resolution land surface model and in-situ soil moisture estimates. Case studies are analyzed using the local land-atmosphere coupling (LoCo) framework that relies on integrated assessment of soil moisture, surface flux, boundary layer, and ambient weather, with results highlighting the critical role of inherent model background biases. In addition, simultaneous assessment of land vs. atmospheric initial conditions in an integrated, process-level fashion can help address the question of whether improvements in traditional NWP verification statistics are achieved for the right reasons.

Southern Great Plains↗

Titan's Ionospheric Chemistry, Fullerenes, Oxygen, Galactic Cosmic Rays and the Formation of Exobiological Molecules on and Within Its Surfaces and Lakes

We discuss the formation of aerosols within Titan's thermosphere-ionosphere and the different chemical pathways. Negative ion measurements by the Cassini Plasma Spectrometer (CAPS) Electron Spectrometer (ELS) give evidence for formation of unsaturated anion carbon chains, while positive ion measurements of the Cassini Ion Neutral Mass Spectrometer (INMS) indicate formation of more aromatic cation hydrocarbons. There is presently no direct observational evidence for large neutral molecule growth in Titan's thermosphere-ionosphere. The hydrocarbon cations are expected to form Polycyclic Aromatic Hydrocarbons (PAHs), those with the addition of nitrogen being called PAHNs. We theorize anion carbon chains can eventually become long enough to fold into fullerene C60,70 carbon shells, of various charge states. Based on laboratory data the fullerenes can trap incoming O+ magnetospheric ions that have relatively high energy collisions with the fullerenes and, once trapped, protect the oxygen atom from Titan's reducing thermosphere-ionosphere. The fullerenes can form into larger onion fullerenes and condense into larger embryo aerosols (i.e., m/q (mass-to-charge ratio) greater than 10,000 amu/q (atomic mass unit-to-charge ratio) anions as observed by CAPS/ELS) eventually falling onto Titan's surface and precipitating to the bottom of its hydrocarbon lakes. Molecule production composed of H, C, N is known to occur in Titan's atmosphere with energy input from the magnetosphere, solar UV (ultraviolet), and deep-penetrating irradiation from galactic cosmic rays (GCR). Space radiation effects by GCR irradiation of Titan's surface and lakes can lead to the manufacture of exobiological molecules with oxygen as the new ingredient. We have developed a model of galactic cosmic ray irradiation of Titan's atmosphere, surface, subsurface and bottoms of Titan lakes. GCR would provide further energy for processing of the aerosols into more complex organic forms such as tholins and precursor molecules for amino acids. A second process called hydrolysis then converts the precursor molecules into amino acids. Hydrolysis is provided via meteor impacts with size greater than 10 kilometers and cryovolcanism both which can produce liquid water on Titan's surface for episodic periods more than several 100 to 1000 years. Our model shows that GCR secondary particles can penetrate ~100 meters below the ice surface (including the bottom of Titan's less dense hydrocarbon lakes ~150-meter depths) and produce chemically significant dosages over very long timescales ~450 million years. The GCR model is combined with laboratory data from experiments in which dry methyl ices were irradiated to doses producing prebiotic amino acids such as glycine. The model calculations show glycine can form to ~2.5 parts-per-billion levels near the surface after ~450 million years of GCR proton irradiation and potentially to 5 parts-per-billion if heavy-ion GCRs up through Fe are included. If such molecules were detected, this would not only confirm this model but indicate that life forms different from ours may not be required.

Edward C. Sittler Jr↗

The production curve for agglutinates in planetary regoliths

Models for the production of agglutinates are developed that can be applied to the lunar surface or to any planetary or asteroidal body lacking an atmosphere. Models are developed using rate equations for progressively more complex situations and range from Model 1, which is a simple linear increase of agglutinate content with time, to Model 4, which includes provision for recycling of existing agglutinates and replenishment and burial of exposed soil. Model 4 has some aspects of a steady state because, depending on the rate constants, agglutinate content may be limited to an intermediate value, even for long exposure times. In an extreme case, agglutinate content may be limited to a value near zero. These models predict that agglutinates should be low in abundance in areas of thin regolith, such as the Lunokhod-2 site on the moon, and on asteroids. The models may also help explain the apparent low agglutinate abundances of lunar regolith breccias and meteorite regolith breccias.

Mckay, D. S.↗

Development of Temperature Sensitive Paints for the Detection of Small Temperature Differences

Temperature sensitive paints (TSP s) have recently been used to detect small temperature differences on aerodynamic model surfaces. These types of applications impose stringent performance requirements on a paint system. The TSP s must operate over a broad temperature range, must be physically robust (cannot chip or peel), must be polishable to at least the smoothness of the model surface, and must have sufficient sensitivity to detect small temperature differences. TSP coatings based on the use of metal complexes in polymer binders were developed at NASA Langley Research Center which meet most of the requirements for detection of small temperature differences under severe environmental conditions.

Oglesby, Donald M.↗

Numerical Modeling of Ablation Heat Transfer

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Ewing, Mark E.↗

Improvements to Wire Bundle Thermal Modeling for Ampacity Determination

Determining current carrying capacity (ampacity) of wire bundles in aerospace vehicles is critical not only to safety but also to efficient design. Published standards provide guidance on determining wire bundle ampacity but offer little flexibility for configurations where wire bundles of mixed gauges and currents are employed with varying external insulation jacket surface properties. Thermal modeling has been employed in an attempt to develop techniques to assist in ampacity determination for these complex configurations. Previous developments allowed analysis of wire bundle configurations but was constrained to configurations comprised of less than 50 elements. Additionally, for vacuum analyses, configurations with very low emittance external jackets suffered from numerical instability in the solution. A new thermal modeler is presented allowing for larger configurations and is not constrained for low bundle infrared emissivity calculations. Formulation of key internal radiation and interface conductance parameters is discussed including the effects of temperature and air pressure on wire to wire thermal conductance. Test cases comparing model-predicted ampacity and that calculated from standards documents are presented.

Rickman, Steve L.↗

Implementation of a Very Large Atmospheric Correction Lookup Table for ASTER Using a Relational Database Management System

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER)is designed to provide a high resolution map of the Earth in both visible, near-infrared, and thermal spectral regions of the electromagnetic spectrum. The ASTER Science Team has developed several standard data product algorithms, but the most complex aod computing-intensive of these is the estimation of surface radiance and reflectance values, which is done by modeling and correcting for the effects of the atmosphere. The algorithm for atmospheric correction in the visible bands sensed by ASTER calls for use of a very large Atmospheric Correction Look Up Table (ACLUT).

ASTER↗

Empirical Requirements Analysis for Mars Surface Operations Using the Flashline Mars Arctic Research Station

Living and working on Mars will require model-based computer systems for maintaining and controlling complex life support, communication, transportation, and power systems. This technology must work properly on the first three-year mission, augmenting human autonomy, without adding-yet more complexity to be diagnosed and repaired. One design method is to work with scientists in analog (mars-like) setting to understand how they prefer to work, what constrains will be imposed by the Mars environment, and how to ameliorate difficulties. We describe how we are using empirical requirements analysis to prototype model-based tools at a research station in the High Canadian Arctic.

Clancey, William J.↗

Modeling Composite Laminate Crushing for Crash Analysis

Crash modeling of composite structures remains limited in application and has not been effectively demonstrated as a predictive tool. While the global response of composite structures may be well modeled, when composite structures act as energy-absorbing members through direct laminate crushing the modeling accuracy is greatly reduced. The most efficient composite energy absorbing structures, in terms of energy absorbed per unit mass, are those that absorb energy through a complex progressive crushing response in which fiber and matrix fractures on a small scale dominate the behavior. Such failure modes simultaneously include delamination of plies, failure of the matrix to produce fiber bundles, and subsequent failure of fiber bundles either in bending or in shear. In addition, the response may include the significant action of friction, both internally (between delaminated plies or fiber bundles) or externally (between the laminate and the crushing surface). A figure shows the crushing damage observed in a fiberglass composite tube specimen, illustrating the complexity of the response. To achieve a finite element model of such complex behavior is an extremely challenging problem. A practical crushing model based on detailed modeling of the physical mechanisms of crushing behavior is not expected in the foreseeable future. The present research describes attempts to model composite crushing behavior using a novel hybrid modeling procedure. Experimental testing is done is support of the modeling efforts, and a test specimen is developed to provide data for validating laminate crushing models.

Fleming, David C.↗

Synthesis Analysis of Mulit-Dimensional Ozone Lidar Measurements in Coastal Environments Toward Improving Ozone Simulation

Large surface ozone variations are frequently observed in coastal regions which makes the modeling of air quality in these regions challenging. The ability to properly model coastal regions relies heavily on the proper detailed observations. There have been three air quality studies (over Chesapeake Bay and the Long Island Sound) that focused on thorough measurements in these complex coastal regions; OWLETS 1 & 2 and LISTOS. Over the course of the campaigns, there were surface ozone observations taken from LIDAR (LIght Detection And Ranging) instruments for a total of 89 days at a high temporal rate. We developed a clustering method that characterizes the multiple timestamps of measurements into a range of ozone events that are specific to these coastal regions. We simulated coastal ozone for these cluster groups using WRF-GC and WRF-Chem to identify gaps within the model. Using the characterized ozone clusters, we were able to identify multiple coastal mechanisms that led to incorrect model simulations.

Claudia Bernier↗

Deriving Scaling Factors Using a Global Hydrological Model to Restore GRACE Total Water Storage Changes for China's Yangtze River Basin

This study used a global hydrological model (GHM), PCR-GLOBWB, which simulates surface water storage changes, natural and human induced groundwater storage changes, and the interactions between surface water and subsurface water, to generate scaling factors by mimicking low-pass filtering of GRACE signals. Signal losses in GRACE data were subsequently restored by the scaling factors from PCR-GLOBWB. Results indicate greater spatial heterogeneity in scaling factor from PCR-GLOBWB and CLM4.0 than that from GLDAS-1 Noah due to comprehensive simulation of surface and subsurface water storage changes for PCR-GLOBWB and CLM4.0. Filtered GRACE total water storage (TWS) changes applied with PCR-GLOBWB scaling factors show closer agreement with water budget estimates of TWS changes than those with scaling factors from other land surface models (LSMs) in China's Yangtze River basin. Results of this study develop a further understanding of the behavior of scaling factors from different LSMs or GHMs over hydrologically complex basins, and could be valuable in providing more accurate TWS changes for hydrological applications (e.g., monitoring drought and groundwater storage depletion) over regions where human-induced interactions between surface water and subsurface water are intensive.

water↗

Future Improvements to Impact for Long Duration Exploration Spaceflight

BACKGROUND: While NASA currently uses the Integrated Medical Model (IMM) to model spaceflight medical risk for Low Earth Orbit missions, IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) was created to improve the fidelity of medical risk analysis for exploration missions involving the Lunar surface and beyond. OVERVIEW: While completion work continues on the initial version of the IMPACT tool, work has already begun to investigate ways to further enhance the fidelity of the model and its utility to stakeholders. DISCUSSION: Future additions leveraging the more robust and flexible IMPACT architecture include analyzing mission segments to account for segment-specific environments (ie. Lunar surface or microgravity), clinically tracible outcomes based on partial treatment of conditions, the ability to model mission outcomes based on provider Knowledge, Skills, and Abilities (KSAs), and the ability to affect the incidence of conditions based on the occurrence of related conditions (e.g. UTI progressing to sepsis). IMPACT also enables the addition of new conditions to its database as required to meet future risk assessment needs. Ultimately, the desire is to broaden IMPACT from a tool that models/trades medical risk to one that does so for all crew health and performance relevant systems (e.g., food, exercise, etc.).

E Stratton↗

Properties of Rubble-Pile Asteroid (101955) Bennu from OSIRIS-REx Imaging and Thermal Analysis

Establishing the abundance and physical properties of regolith and boulders on asteroids is crucial for understanding the formation and degradation mechanisms at work on their surfaces. Using images and thermal data from NASA's Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer (OSIRIS-REx) spacecraft, we show that asteroid (101955) Bennu's surface is globally rough, dense with boulders, and low in albedo. The number of boulders is surprising given Bennu's moderate thermal inertia, suggesting that simple models linking thermal inertia to particle size do not adequately capture the complexity relating these properties. At the same time, we find evidence for a wide range of particle sizes with distinct albedo characteristics. Our findings imply that ages of Bennu's surface particles span from the disruption of the asteroid's parent body (boulders) to recent in situ production (micrometre-scale particles).

DellaGiustina, D. N.↗

Empirical Modeling of Plant Gas Fluxes in Controlled Environments

As humans extend their reach beyond the earth, bioregenerative life support systems must replace the resupply and physical/chemical systems now used. The Controlled Ecological Life Support System (CELSS) will utilize plants to recycle the carbon dioxide (CO2) and excrement produced by humans and return oxygen (O2), purified water and food. CELSS design requires knowledge of gas flux levels for net photosynthesis (PS(sub n)), dark respiration (R(sub d)) and evapotranspiration (ET). Full season gas flux data regarding these processes for wheat (Triticum aestivum), soybean (Glycine max) and rice (Oryza sativa) from published sources were used to develop empirical models. Univariate models relating crop age (days after planting) and gas flux were fit by simple regression. Models are either high order (5th to 8th) or more complex polynomials whose curves describe crop development characteristics. The models provide good estimates of gas flux maxima, but are of limited utility. To broaden the applicability, data were transformed to dimensionless or correlation formats and, again, fit by regression. Polynomials, similar to those in the initial effort, were selected as the most appropriate models. These models indicate that, within a cultivar, gas flux patterns appear remarkably similar prior to maximum flux, but exhibit considerable variation beyond this point. This suggests that more broadly applicable models of plant gas flux are feasible, but univariate models defining gas flux as a function of crop age are too simplistic. Multivariate models using CO2 and crop age were fit for PS(sub n), and R(sub d) by multiple regression. In each case, the selected model is a subset of a full third order model with all possible interactions. These models are improvements over the univariate models because they incorporate more than the single factor, crop age, as the primary variable governing gas flux. They are still limited, however, by their reliance on the other environmental conditions under which the original data were collected. Three-dimensional plots representing the response surface of each model are included. Suitability of using empirical models to generate engineering design estimates is discussed. Recommendations for the use of more complex multivariate models to increase versatility are included.

Cornett, Jessie David↗

SIBatt-3D: In-Space/On-Surface 3D Printing of Sodium Ion Batteries from ISRU Materials

Constructed more than 20 years ago, the International Space Station’s primary power system originally used nickel-hydrogen batteries with a lifetime of 6.5 years, until NASA began the process of replacing them in 2016 with lithium-ion batteries with a lifetime of 10 years. The demanding and costly process was accomplished after four flights of the Japanese H-II Transfer Vehicle cargo spacecraft (with a cost of about $10,000 per pound of payload), and 13 different astronauts conducting 14 spacewalks. Besides utilization in the ISS, rechargeable batteries are present in many space applications: they are installed in exploration robots, life support systems and in portable communication devices, to mention some. In this context, this project is focused on the in-space manufacturing of shape-conformable batteries using in-situ resources, and aims to address the NASA’s gaps related to the development of next generation of energy storage devices (TX03), as well as in-space manufacturing and in-situ resource utilization (TX07). The proposed work also tackles the HEOMD’s objectives targeting the in-space additive manufacturing (AM) from Lunar/Martian materials (regolith as AM feedstock) to reinvigorate America’s Human Space Exploration Program (SPD-1). This project is in direct alignment with the STMD’s objectives to demonstrate in-space autonomous manufacturing and assembly of complete systems by 2030, and to enable humans to live and explore in space and on planetary surfaces by 2040 thanks to in-space habitation, infrastructure development and in-situ resource utilization (ST1 and ST5). Manufacturing of shape conformable batteries directly in-space and using in-situ resources would also contribute to reducing the payload weight and volume (TX12) for future missions, thus reducing risk for long term Mars missions where rapid resupply is logistically infeasible. Nowadays, commercial batteries consist of stacked two-dimensional (2D) sheets, which are only manufactured in restricted geometries (cylindrical and coin cell). Evolving from conventional 2D, complex 3D battery architectures have been proven to increase the electrochemical active surface area and ion diffusion path, leading to improved areal energy density and power performance. This tendency was illustrated in our recent in-depth modeling studies by simulating a classical Ragone plot exhibiting the energy-power relationship. Our team demonstrated through modeling that a complex gyroidal 3D printed battery architecture exhibits significantly improved power performances (>150% at the current density of 6C; full discharge in 10 minutes) in comparison to a traditional 3D printed planar geometry. Motivated by these results and as the fabrication of intricate 3D battery design is only possible experimentally thanks to the geometric freedom offered by additive manufacturing (AM), our team has already initiated leveraging thermoplastic material extrusion at the laboratory scale. While 3D printing of batteries is relatively recent (2013), it has witnessed a growing interest during the last recent years, as next-generation shape-conformable 3D batteries can be co-designed with the system. Consequently, dead-volume and mass brought from Earth are minimized, in addition to improved battery performance, in alignment with the aforementioned NASA’s objectives. Further, while this project is specifically dedicated to batteries, it lends itself towards the maturation of in-space manufacturing via 3D printing using in-situ resources, stated in HEOMD and STMD goals.

In-Space Manufacturing↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗