Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural model decomposition”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mechanism and kinetic study of rare earth extraction from allanite by direct acid leaching

Direct acid leaching was conducted to recover rare earth elements (REEs) from an allanite ore assaying 9,861.4 ppm of REEs. The effects of operational variables such as acid type, sulfuric acid (H 2 SO 4 ) concentration, temperature, solid/liquid (S/L) ratio, and particle size on the total REE (TREE) recovery were systematically examined. The REE mineralogy and the involved leaching mechanisms were investigated by microscopic characterization and residue re-leaching tests. It was found that approximately 80% of REEs was extracted using 1 M H 2 SO 4 at 75 °C for 2 h. Increasing temperature significantly improved the REE extraction, while changing acid type, H 2 SO 4 concentration, and particle size did not. The TREE recovery also gradually decreased at higher S/L ratios. The low fraction of Si leached (<4.0%) suggested that the extraction did not significantly alter the crystal structure of the dominant mineral. The leaching kinetics of REEs were fast within the first 10 min, followed by a much slower stage (i.e., 10-120 min). The kinetic modeling results showed that each stage may be controlled by both chemical reaction and diffusion (i.e., mix control). The activation energies for REE leaching using 1 M H 2 SO 4 in the 0-10 min and 10-120 min were determined to be 20.3 and 10.8 kJ/mol, respectively. Finally, the residue re-leaching results suggested that the fast leaching kinetics at the beginning may be explained by the preferential decomposition of metamict allanite rather than the formation of silica gel, while the following slower leaching stage was due to the dissolution of difficult-to-leach REE-bearing minerals such as well-crystallized allanite.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

High Pressure X-ray Diffraction and Equation of State of Hydrazine

Synchrotron X-ray diffraction has been used to investigate the structure and equation of state (EOS) of hydrazine (N 2 H 4 ) up to 54.3 GPa at 298 K. The diffraction patterns could be fit to a monoclinic unit-cell structure and put strong constraints on previously reported phase transitions documented by vibrational spectroscopy over this pressure range. Pressure–volume ( P–V ) data were fit using a Vinet EOS, yielding parameters: V 0 = 45.2 Å 3 /molecule (fixed), K 0 = 11.8(7) GPa, and K 0 ′ = 6.5(2). Previously measured high-pressure vibrational frequency shifts were used to estimate the vibrational free energy and model P–V–T isotherms from 0 to 1200 K. The results of the P–V–T isotherms are compared to existing shock Hugoniot data on hydrazine and 298 K isotherms for assemblages of possible decomposition products. This comparison suggests dissociation at high density under shock loading. Good correspondence was found between the static lattice EOS as calculated by the model and the previously reported EOS as calculated by density functional theory. Finally, these results resolve existing uncertainties about the EOS and crystal symmetry of hydrazine at high pressure and provide valuable baseline information on this important energetic material.

diffraction↗
NASA NTRS2017

Aeroelastic Uncertainty Quantification Studies Using the S4T Wind Tunnel Model

This paper originates from the joint efforts of an aeroelastic study team in the Applied Vehicle Technology Panel from NATO Science and Technology Organization, with the Task Group number AVT-191, titled "Application of Sensitivity Analysis and Uncertainty Quantification to Military Vehicle Design." We present aeroelastic uncertainty quantification studies using the SemiSpan Supersonic Transport wind tunnel model at the NASA Langley Research Center. The aeroelastic study team decided treat both structural and aerodynamic input parameters as uncertain and represent them as samples drawn from statistical distributions, propagating them through aeroelastic analysis frameworks. Uncertainty quantification processes require many function evaluations to asses the impact of variations in numerous parameters on the vehicle characteristics, rapidly increasing the computational time requirement relative to that required to assess a system deterministically. The increased computational time is particularly prohibitive if high-fidelity analyses are employed. As a remedy, the Istanbul Technical University team employed an Euler solver in an aeroelastic analysis framework, and implemented reduced order modeling with Polynomial Chaos Expansion and Proper Orthogonal Decomposition to perform the uncertainty propagation. The NASA team chose to reduce the prohibitive computational time by employing linear solution processes. The NASA team also focused on determining input sample distributions.

Nikbay, Melike↗
NASA NTRS2018

Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds

The organic compounds found in carbonaceous chondrite meteorites provide insight into primordial solar system chemistry. Evaluating the formation and decomposition mechanisms of meteoritic amino acids may aid our understanding of the origins of life and homochirality on Earth. The amino acid glycine is widespread in meteorites and other extraterrestrial environments; other amino acids, such as isovaline, are found with enantiomeric excesses in some meteorites. The relationship between meteoritic amino acids and other compounds with similar molecular structures, such as aliphatic monoamines and monocarboxylic acids is unclear; experimental results evaluating the decomposition of amino acids have produced inconclusive results about the preferred pathways, reaction intermediates, and if the conditions applied may be compatible with those occurring inside meteoritic parent bodies. In this work, we performed extensive tandem metadynamics, umbrella sampling, and committor analysis to simulate the neutral mild hydrothermal decomposition mechanisms of glycine and isovaline and put them into context for the origins of meteoritic organic compounds. Our ab initio simulations aimed to determine free energy profiles and decomposition pathways for glycine and isovaline. We found that under our modeled conditions, methylammonium, glycolic acid, and sec-butylamine are the most likely decomposition products. These results suggest that meteoritic aliphatic monocarboxylic acids are not produced from decomposition of meteoritic amino acids. Our results also indicate that the decomposition of L-isovaline prefers an enantioselective pathway resulting in the production of (S)-sec-butylamine.

Enhanced sampling↗
DOE OSTI2013

A new approach to detailed structural decomposition from the splash and phat surveys: Kicked-up disk stars in the Andromeda galaxy?

We characterize the bulge, disk, and halo subcomponents in the Andromeda galaxy (M31) over the radial range 4 kpc < R {sub proj} < 225 kpc. The cospatial nature of these subcomponents renders them difficult to disentangle using surface brightness (SB) information alone, especially interior to ∼20 kpc. Our new decomposition technique combines information from the luminosity function (LF) of over 1.5 million bright (20 < m {sub 814W} < 22) stars from the Panchromatic Hubble Andromeda Treasury survey, radial velocities of over 5000 red giant branch stars in the same magnitude range from the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo survey, and integrated I-band SB profiles from various sources. We use an affine-invariant Markov chain Monte Carlo algorithm to fit an appropriate toy model to these three data sets. The bulge, disk, and halo SB profiles are modeled as a Sérsic, exponential, and cored power law, respectively, and the LFs are modeled as broken power laws. We present probability distributions for each of 32 parameters describing the SB profiles and LFs of the three subcomponents. We find that the number of stars with a disk-like LF is 5.2% ± 2.1% larger than the number with disk-like (dynamically cold) kinematics, suggesting that some stars born in the disk have been dynamically heated to the point that they are kinematically indistinguishable from halo members. This is the first kinematical evidence for a 'kicked-up disk' halo population in M31. The fraction of kicked-up disk stars is consistent with that found in simulations. We also find evidence for a radially varying disk LF, consistent with a negative metallicity gradient in the stellar disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Quantified Gamow shell model interaction for p s d -shell nuclei

Background: The structure of weakly bound and unbound nuclei close to particle drip lines is one of the major science drivers of nuclear physics. A comprehensive understanding of these systems goes beyond the traditional configuration interaction approach formulated in the Hilbert space of localized states (nuclear shell model) and requires an open quantum system description. The complex-energy Gamow shell model (GSM) provides such a framework as it is capable of describing resonant and nonresonant many-body states on equal footing. Purpose: To make reliable predictions, quality input is needed that allows for the full uncertainty quantification of theoretical results. In this study, we carry out the optimization of an effective GSM (one-body and two-body) interaction in the psdf-shell-model space. The resulting interaction is expected to describe nuclei with 5 ≤ A ≲ 12 at the p-sd-shell interface. Method: The one-body potential of the 4 He core is modeled by a Woods-Saxon + spin-orbit + Coulomb potential, and the finite-range nucleon-nucleon interaction between the valence nucleons consists of central, spin-orbit, tensor, and Coulomb terms. The GSM is used to compute key fit observables. The χ 2 optimization is performed using the Gauss-Newton algorithm augmented by the singular value decomposition technique. The resulting covariance matrix enables quantification of statistical errors within the linear regression approach. Results: The optimized one-body potential reproduces nucleon- 4 He scattering phase shifts up to an excitation energy of 20 MeV. The two-body interaction built on top of the optimized one-body field is adjusted to the bound and unbound ground-state binding energies and selected excited states of the helium, lithium, and beryllium isotopes up to A = 9 . A very good agreement with experimental results was obtained for binding energies. First applications of the optimized interaction include predictions for two-nucleon correlation densities and excitation spectra of light nuclei with quantified uncertainties. In conclusion: The new interaction will enable comprehensive and fully quantified studies of structure and reactions aspects of nuclei from the psd region of the nuclear chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Insights into Acetic Acid Binding and Ketene Formation on Anatase TiO 2 (101)

Understanding the adsorption and reactivity of carboxylic acids on oxide surfaces is of great interest in catalysis for biomass upgrading via ketonization, a carbon–carbon coupling reaction. Herein, we investigate the adsorption and reaction of acetic acid on anatase TiO 2 (101) using scanning tunneling microscopy, infrared spectroscopy, temperature programmed reaction, and density functional theory calculations. We demonstrate the adsorption of acetic acid can form two intermediates: (1) dissociated, bidentate acetate with an associated bridging hydroxyl, and (2) molecular, monodentate acetic acid. The coexistence of ordered phases with increasing monolayer (ML) saturation coverages consisting of (1) pure acetate (0.5 ML), (2) mixed acetate/acetic acid (0.67 ML), (3) mixed acetate/acetic acid (1.0 ML) and (4) pure acetic acid demonstrates similar energetics for both acetate and acetic acid species. Under ultra-high vacuum conditions, the presence of both monodentate acetic acid and bidentate acetate was observed below room temperature, while solely bidentate acetate was observed up to 575 K. The deprotonation of acetic acid produces water at 280 K, while the thermal decomposition of bidentate acetate produces ketene and acetic acid at 645 K. In conclusion, this model study provides detailed insight into the stability and reactivity of carboxylic acid surface-bound intermediates, which could participate during ketonization reactions for biomass upgrading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Towards space from Hilbert space: finding lattice structure in finite-dimensional quantum systems

We report that field theories place one or more degrees of freedom at every point in space. Hilbert spaces describing quantum field theories, or their finite-dimensional discretizations on lattices, therefore have large amounts of structure: they are isomorphic to the tensor product of a smaller Hilbert space for each lattice site or point in space. Local field theories respecting this structure have interactions which preferentially couple nearby points. The emergence of classicality through decoherence relies on this framework of tensor-product decomposition and local interactions. We explore the emergence of such lattice structure from Hilbert-space considerations alone. We point out that the vast majority of finite-dimensional Hilbert spaces cannot be isomorphic to the tensor product of Hilbert-space subfactors that describes a lattice theory. A generic Hilbert space can only be split into a direct sum corresponding to a basis of state vectors spanning the Hilbert space, we consider setups in which the direct sum is naturally decomposed into two pieces. We define a notion of direct-sum locality which characterizes states and decompositions compatible with Hamiltonian time evolution. Finally, we illustrate these notions for a toy model that is the finite-dimensional discretization of the quantum-mechanical double-well potential. We discuss their relevance in cosmology and field theory, especially for theories which describe a landscape of vacua with different spacetime geometries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Hyper-Differential Sensitivity Analysis of Uncertain Parameters in PDE-Constrained Optimization

Many problems in engineering and sciences require the solution of large scale optimization constrained by partial differential equations (PDEs). Though PDE-constrained optimization is itself challenging, most applications pose additional complexity, namely, uncertain parameters in the PDEs. Uncertainty quantification (UQ) is necessary to characterize, prioritize, and study the influence of these uncertain parameters. Sensitivity analysis, a classical tool in UQ, is frequently used to study the sensitivity of a model to uncertain parameters. In this article, we introduce "hyper-differential sensitivity analysis" which considers the sensitivity of the solution of a PDE-constrained optimization problem to uncertain parameters. Our approach is a goal-oriented analysis which may be viewed as a tool to complement other UQ methods in the service of decision making and robust design. We formally define hyper-differential sensitivity indices and highlight their relationship to the existing optimization and sensitivity analysis literatures. Assuming the presence of low rank structure in the parameter space, computational efficiency is achieved by leveraging a generalized singular value decomposition in conjunction with a randomized solver which converts the computational bottleneck of the algorithm into an embarrassingly parallel loop. Two multi-physics examples, consisting of nonlinear steady state control and transient linear inversion, demonstrate efficient identification of the uncertain parameters which have the greatest influence on the optimal solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Evaluation of Advanced Signal Processing Techniques to Improve Detection and Identification of Embedded Defects

By the end of 1996, 109 Nuclear Power Plants were operating in the United States, producing 22% of the Nation’s electricity [1]. At present, more than two thirds of these power plants are more than 40 years old. The purpose of the U.S. Department of Energy Office of Nuclear Energy’s Light Water Reactor Sustainability (LWRS) Program is to develop technologies and other solutions that can improve the reliability, sustain the safety, and extend the operating lifetimes of nuclear power plants (NPPs) beyond 60 years [2]. The most important safety structures in an NPP are constructed of concrete. The structures generally do not allow for destructive evaluation and access is limited to one side of the concrete element. Therefore, there is a need for techniques and technologies that can assess the internal health of complex, reinforced concrete structures nondestructively. Previously, we documented the challenges associated with Non-Destructive Evaluation (NDE) of thick, reinforced concrete sections and prioritized conceptual designs of specimens that could be fabricated to represent NPP concrete structures [3]. Consequently, a 7 feet tall, by 7 feet wide, by 3 feet and 4-inch-thick concrete specimen was constructed with 2.257-inch-and 1-inch-diameter rebar every 6 to 12 inches. In addition, defects were embedded the specimen to assess the performance of existing and future NDE techniques. The defects were designed to give a mix of realistic and controlled defects for assessment of the necessary measures needed to overcome the challenges with more heavily reinforced concrete structures. Information on the embedded defects is documented in [4]. We also documented the superiority of Frequency Banded Decomposition (FBD) Synthetic Aperture Focusing Technique (SAFT) over conventional SAFT when probing defects under deep concrete cover. Improvements include seeing an intensity corresponding to a defect that is either not visible at all in regular, full frequency content SAFT, or an improvement in contrast over conventional SAFT reconstructed images. This report documents our efforts in four fronts: 1) Comparative study between traditional SAFT and FBD SAFT for concrete specimen with and without Alkali-Silica Reaction (ASR) damage, 2) improvement of our Model-Based Iterative Reconstruction (MBIR) for thick reinforced concrete [5], 3) development of a universal framework for sharing, reconstruction, and visualization of ultrasound NDE datasets, and 4) application of machine learning techniques for automated detection of ASR inside concrete. Our comparative study between FBD and traditional SAFT reconstruction images shows a clear difference between images of ASR and non-ASR specimens. In particular, the left first harmonic shows an increased contrast and sensitivity to ASR damage. For MBIR, we show the superiority of model-based techniques over delay and sum techniques such as SAFT. Improvements include elimination of artifacts caused by direct arrival signals, and increased contrast and Signal to Noise Ratio. For the universal framework, we document a format for data storage based on the HDF5 file format, and also propose a modular Graphic User Interface (GUI) for easy customization of data conversion, reconstruction, and visualization routines. Finally, two techniques for ASR automated detection are presented. The first technique is based on an analysis of the frequency content using Hilbert Transform Indicator (HTI) and the second technique employees Artificial Neural Network (ANN) techniques for training and classification of ultrasound data as ASR or non-ASR damaged classes. The ANN technique shows great potential with classification accuracy above 95%. These approaches are extensible to the detection of additional reinforced, thick concrete defects and damage.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Photocatalytic degradation of cooking fume on a TiO2-coated carbon nanotubes composite filter

Highlights: • The TiO{sub 2}-coated carbon nanotube composite filter was fabricated for treating cooking fume particles. • MWCNTs formed a network-like structure surface layer of composite filter. • The composite filter showed high surface area. • The composite filter was tested in oleic acid decomposition reaction. • Parallel reaction mechanism was attempted to describe oleic acid degradation. A TiO{sub 2}/MWCNTs/Al{sub 2}O{sub 3}/NF (TiO{sub 2} nanoparticles/multi-walled carbon nanotubes /Al{sub 2}O{sub 3}-modified nickel foams) composite filter was prepared through grown of MWCNTs upon Al{sub 2}O{sub 3}-modified nickel foams using chemical vapor deposition (CVD) method and then coated TiO{sub 2} nanoparticles onto MWCNTs using dip-coating method. The MWCNTs forms a network-like structure on the surface layer of composite filter. The BET surface area of TiO{sub 2}/MWCNTs/Al{sub 2}O{sub 3}/NF composite filter is more than 375 times higher than that of the pristine nickel foam. The degradation performance of oil particles was evaluated. Under optimized reaction conditions, 90% of the oleic acid on the filter was decomposed after 25 h. In addition, the effect of reaction temperature and light intensity were studied. Increased reaction temperature and light intensity were found to significantly improve the degradation rate and restrain the formation of refractory product. Furthermore, a parallel reaction kinetic model was adopted for degradation of oleic acid on the composite filter. The composite filter can reduces the accumulation of oil particles on the filter and the energy consumption of the filter. This approach proposes a new way to treat the cooking fume emission at residential building roofs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Electrochemical Stability of Metastable Materials

We present a first-principles-based formalism to provide a quantitative measure of the thermodynamic instability and propensity for electrochemical stabilization, passivation, or corrosion of metastable materials in aqueous media. Here, we demonstrate that this formalism can assess the relative Gibbs free energy of candidate materials in aqueous media as well as their decomposition products, combining solid and aqueous phases, as a function of pH and potential. On the basis of benchmarking against 20 stable as well as metastable materials reported in the literature and also our experimental characterization of metastable triclinic-FeVO 4 , we present quantitative estimates for the relative Gibbs free energy and corresponding aqueous regimes where these materials are most likely to be stable, form inert passivating films, or steadily corrode to aqueous species. Furthermore, we show that the structure and composition of the passivating films formed on triclinic-FeVO 4 are also in excellent agreement with the Point Defect Model, as proposed by the corrosion community.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Assembled Solute Networks in Crowded Electrolyte Solutions and Nanoconfinement of Charged Redoxmer Molecules

Redoxmers are electrochemically active organic molecules storing charge and energy in electrolyte fluids circulating through redox flow batteries (RFBs). Such molecules typically have solvent-repelling cores and solvent-attracting pendant groups introduced to increase solubility in liquid electrolytes. These two features can facilitate nanoscale aggregation of the redoxmer molecules in crowded solutions. In some cases, this aggregation leads to the emergence of continuous networks of solute molecules in contact, and the solution becomes microscopically heterogeneous. In this study, we used small-angle X-ray scattering (SAXS) and molecular dynamics modeling to demonstrate formation of such networks and examine structural factors controlling this self-assembly. We also show that salt ions become excluded from these solute aggregates into small pockets of electrolytes, where these ions strongly associate. This confinement by exclusion is also likely to occur to charged redoxmer molecules in a “sea” of neutral precursors coexisting in the same solution. Here, we demonstrate that the decay lifetime of the confined charged molecules in such solutions can increase several fold compared to dilute solutions. We attribute this behavior to a “microreactor effect” on reverse reactions of the confined species during their decomposition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Bond Order Conservation Strategies in Catalysis Applied to the NH 3 Decomposition Reaction

On the basis of an extensive set of density functional theory calculations, it is shown that a simple scheme provides a fundamental understanding of variations in the transition state energies and structures of reaction intermediates on transition metal surfaces across the periodic table. The scheme is built on the bond order conservation principle and requires a limited set of input data, still achieving transition state energies as a function of simple descriptors with an error smaller than those of approaches based on linear fits to a set of calculated transition state energies. Here, we have applied this approach together with linear scaling of adsorption energies to obtain the energetics of the NH 3 decomposition reaction on a series of stepped fcc(211) transition metal surfaces. Moreover, this information is used to establish a microkinetic model for the formation of N 2 and H 2 , thus providing insight into the components of the reaction that determines the activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1992

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗
DOE OSTI2016

A FOUR-FLUID MHD MODEL OF THE SOLAR WIND/INTERSTELLAR MEDIUM INTERACTION WITH TURBULENCE TRANSPORT AND PICKUP PROTONS AS SEPARATE FLUID

We have developed a four-fluid, three-dimensional magnetohydrodynamic model of the solar wind interaction with the local interstellar medium. The unique features of the model are: (a) a three-fluid description for the charged components of the solar wind and interstellar plasmas (thermal protons, electrons, and pickup protons), (b) the built-in turbulence transport equations based on Reynolds decomposition and coupled with the mean-flow Reynolds-averaged equations, and (c) a solar corona/solar wind model that supplies inner boundary conditions at 40 au by computing solar wind and magnetic field parameters outward from the coronal base. The three charged species are described by separate energy equations and are assumed to move with the same velocity. The fourth fluid in the model is the interstellar hydrogen which is treated by separate continuity, momentum, and energy equations and is coupled with the charged components through photoionization and charge exchange. We evaluate the effects of turbulence transport and pickup protons on the global heliospheric structure and compute the distribution of plasma, magnetic field, and turbulence parameters throughout the heliosphere for representative solar minimum and maximum conditions. We compare our results with Voyager 1 observations in the outer heliosheath and show that the relative amplitude of magnetic fluctuations just outside the heliopause is in close agreement with the value inferred from Voyager 1 measurements by Burlaga et al. The simulated profiles of magnetic field parameters in the outer heliosheath are in qualitative agreement with the Voyager 1 observations and with the analytical model of magnetic field draping around the heliopause of Isenberg et al.

79 ASTRONOMY AND ASTROPHYSICS↗