Results for “Theoretical description and modelling” — Page 25

NASA NTRS1975

Contribution of the National Aeronautics and Space Administration Langley Research Center

As part of a special international effort, three nozzles were designed and tested on single nacelle models in wind tunnels of several nations belonging to the North Atlantic Treaty Organization. All three of these nozzles were investigated in the Langley 16-foot transonic wind tunnel at the National Aeronautics and Space Administration's Langley Research Center. Langley Research Center also contributed theoretical calculations of the jet plume boundary and afterbody pressures. The calculations were obtained using an iterative solution which combined the inviscid Douglas Neumann method for the external flow with the method of characteristics for the flow in the jet plume. For the investigation, the nozzles were mounted on a single nacelle model 15.24 centimeters in diameter and 162.56 centimeters long. Tests were made at free stream Mach number from 0.4 to 1.2, and at Reynolds numbers per meter from 7.38 million to 13.78 million depending on the Mach number. Four types of data were recorded: afterbody pressure data, afterbody force data, model boundary layer data, and tunnel wall pressure data. The ratio of jet total pressure to free stream static pressure ranged up to 8.5. A description of the wind tunnel, model, and test procedure is included.

Compton, W. B., III↗

DOE OSTI2023

Electronic structure of Co 3 d states in the Kitaev material candidate honeycomb cobaltate Na 3 Co 2 SbO 6 probed with x-ray dichroism

The recent prediction that honeycomb lattices of Co 2+ (3d 7 ) ions could host dominant Kitaev interactions provides an exciting direction for exploration of new routes to stabilizing Kitaev’s quantum spin liquid in real materials. Na 3 Co 2 SbO 6 has been singled out as a potential material candidate provided that spin and orbital moments couple into a J eff = $\frac{1}{2}$ ground state, and that the relative strength of trigonal crystal field and spin-orbit coupling acting on Co ions can be tailored. Using x-ray linear dichroism (XLD) and x-ray magnetic circular dichroism (XMCD) experiments, alongside configuration interaction calculations, we confirm the counterintuitive positive sign of the trigonal crystal field acting on Co 2+ ions and test the validity of the J eff = $\frac{1}{2}$ description of the electronic ground state. Furthermore, the results lend experimental support to recent theoretical predictions that a compression (elongation) of CoO 6 octahedra along (perpendicular to) the trigonal axis would drive this cobaltate toward the Kitaev limit, assuming the J eff = $\frac{1}{2}$ character of the electronic ground state is preserved.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Superconducting diode effect in quasi-one-dimensional systems

The recent observations of the superconducting diode effect pose the challenge to fully understand the necessary ingredients for nonreciprocal phenomena in superconductors. In this theoretical work, we focus on the nonreciprocity of the critical current in a quasi-one-dimensional superconductor. We define the critical current as the value of the supercurrent at which the quasiparticle excitation gap closes (depairing). Once the critical current is exceeded, the quasiparticles can exchange energy with the superconducting condensate, giving rise to dissipation. Our minimal model can be microscopically derived as a low-energy limit of a Rashba spin-orbit coupled superconductor in a Zeeman field. Within the proposed model, we explore the nature of the nonreciprocal effects of the critical current both analytically and numerically. Further, our results quantify how system parameters such as spin-orbit coupling and quantum confinement affect the strength of the superconducting diode effect. Our theory provides a complementary description to Ginzburg-Landau theories of the effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Toward a unified treatment of Δ S = 0 parity violation in low-energy nuclear processes

We revisit the unified treatment of low-energy hadronic parity violation espoused by Desplanques, Donoghue, and Holstein to the end of an ab initio treatment of parity violation in low-energy nuclear processes within the standard model. Here, we use our improved effective Hamiltonian and precise nonperturbative assessments of the quark charges of the nucleon within lattice quantum chromodynamics (QCD) to make new assessments of the parity-violating meson-nucleon coupling constants. Comparing with recent, precise measurements of hadronic parity violation in few-body nuclear reactions, we find improved agreement with these experimental results, though some tensions remain. We thus note the broader problem of comparing low-energy constants from nuclear and few-nucleon systems, considering, too, unresolved theoretical issues in connecting an ab initio, effective Hamiltonian approach to chiral effective theories. We note how future experiments and lattice QCD studies could sharpen the emerging picture, promoting the study of hadronic parity violation as a laboratory for testing “end-to-end” theoretical descriptions of weak processes in hadrons and nuclei at low energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Quantum Many-Body Scars: A Quasiparticle Perspective

Weakly interacting quasiparticles play a central role in the low-energy description of many phases of quantum matter. At higher energies, however, quasiparticles cease to be well defined in generic many-body systems owing to a proliferation of decay channels. In this review, we discuss the phenomenon of quantum many-body scars, which can give rise to certain species of stable quasiparticles throughout the energy spectrum. This goes along with a set of unusual nonequilibrium phenomena including many-body revivals and nonthermal stationary states. We provide a pedagogical exposition of this physics via a simple yet comprehensive example, that of a spin-1 XY model. We place our discussion in the broader context of symmetry-based constructions of many-body scar states, projector embeddings, and Hilbert space fragmentation. We conclude with a summary of experimental progress and theoretical puzzles.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Building surrogate models of nuclear density functional theory with Gaussian processes and autoencoders

From the lightest Hydrogen isotopes up to the recently synthesized Oganesson (Z = 118), it is estimated that as many as about 8,000 atomic nuclei could exist in nature. Most of these nuclei are too short-lived to be occurring on Earth, but they play an essential role in astrophysical events such as supernova explosions or neutron star mergers that are presumed to be at the origin of most heavy elements in the Universe. Understanding the structure, reactions, and decays of nuclei across the entire chart of nuclides is an enormous challenge because of the experimental difficulties in measuring properties of interest in such fleeting objects and the theoretical and computational issues of simulating strongly-interacting quantum many-body systems. Nuclear density functional theory (DFT) is a fully microscopic theoretical framework which has the potential of providing such a quantitatively accurate description of nuclear properties for every nucleus in the chart of nuclides. Thanks to high-performance computing facilities, it has already been successfully applied to predict nuclear masses, global patterns of radioactive decay like β or γ decay, and several aspects of the nuclear fission process such as, e.g., spontaneous fission half-lives. Yet, predictive simulations of nuclear spectroscopy—the low-lying excited states and transitions between them—or of nuclear fission, or the quantification of theoretical uncertainties and their propagation to basic or applied nuclear science applications, would require several orders of magnitude more calculations than currently possible. However, most of this computational effort would be spent into generating a suitable basis of DFT wavefunctions. Such a task could potentially be considerably accelerated by borrowing tools from the field of machine learning and artificial intelligence. In this paper, we review different approaches to applying supervised and unsupervised learning techniques to nuclear DFT.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1990

Modeling the pressure-strain correlation of turbulence: An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗

NASA NTRS1991

Modelling the pressure-strain correlation of turbulence - An invariant dynamical systems approach

The modeling of the pressure-strain correlation of turbulence is examined from a basic theoretical standpoint with a view toward developing improved second-order closure models. Invariance considerations along with elementary dynamical systems theory are used in the analysis of the standard hierarchy of closure models. In these commonly used models, the pressure-strain correlation is assumed to be a linear function of the mean velocity gradients with coefficients that depend algebraically on the anisotropy tensor. It is proven that for plane homogeneous turbulent flows the equilibrium structure of this hierarchy of models is encapsulated by a relatively simple model which is only quadratically nonlinear in the anisotropy tensor. This new quadratic model - the SSG model - is shown to outperform the Launder, Reece, and Rodi model (as well as more recent models that have a considerably more complex nonlinear structure) in a variety of homogeneous turbulent flows. Some deficiencies still remain for the description of rotating turbulent shear flows that are intrinsic to this general hierarchy of models and, hence, cannot be overcome by the mere introduction of more complex nonlinearities. It is thus argued that the recent trend of adding substantially more complex nonlinear terms containing the anisotropy tensor may be of questionable value in the modeling of the pressure-strain correlation. Possible alternative approaches are discussed briefly.

Speziale, Charles G.↗

DOE OSTI2021

Addressing Human and Organizational Factors In Nuclear Industry Modernization: A Sociotechnically-Based Strategic Framework

The modernization of nuclear power plants will require an advanced concept of operations, involving an integrated set of tightly coupled systems in which all stakeholders act in a coordinated manner. For this modernization effort to be enabled, we developed a human and organizational factors approach, based on a broad sociotechnical framework. Starting from core human factors principles, we conducted a literature review of the methods and approaches relevant to the modernization problem. These included core disciplines such as cognitive systems engineering, systems theoretic accident modeling and processes (STAMP), human systems integration, resilience engineering, and macroergonomics but also related topics of safety culture and organizational change. From this literature, we developed a conceptual framework centered around the work system with its four interacting components: people, technology, process, and governance. In an effective work system, these four components are jointly optimized according to three systems criteria: efficiency, effectiveness, and safety. System failure may result from excessive emphasis on any one criterion. The actual work of attaining joint optimization in a given work system can be accomplished by utilizing three high level functions: knowledge elicitation, knowledge representation, and cross-functional integration. We illustrated the utility of this approach by applying it to practical problems and case studies. The complete report is available in [1]; this paper is a description of our approach and report contents.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2025

Practical Introduction to Benchmarking and Characterization of Quantum Computers

Rapid progress in quantum technology has transformed quantum computing and quantum information science from theoretical possibilities into tangible engineering challenges. Breakthroughs in quantum algorithms, quantum simulations, and quantum error correction are bringing useful quantum computation closer to fruition. These remarkable achievements have been facilitated by advances in quantum characterization, verification, and validation (QCVV). QCVV methods and protocols enable scientists and engineers to scrutinize, understand, and enhance the performance of quantum information-processing devices. In this tutorial, we review the fundamental principles underpinning QCVV, and introduce a diverse array of QCVV tools used by quantum researchers. We define and explain QCVV’s core models and concepts—quantum states, measurements, and processes—and illustrate how these building blocks are leveraged to examine a target system or operation. We survey and introduce protocols ranging from simple qubit characterization to advanced benchmarking methods. Along the way, we provide illustrated examples and detailed descriptions of the protocols, highlight the advantages and disadvantages of each, and discuss their potential scalability to future large-scale quantum computers. This tutorial serves as a guidebook for researchers unfamiliar with the benchmarking and characterization of quantum computers, and also as a detailed reference for experienced practitioners.

open quantum systems & decoherence↗

NASA NTRS1993

Comprehensive investigation of HgCdTe metalorganic chemical vapor deposition

The principal objective of this experimental and theoretical research program was to explore the possibility of depositing high quality epitaxial CdTe and HgCdTe at very low pressures through metalorganic chemical vapor deposition (MOCVD). We explored two important aspects of this potential process: (1) the interaction of molecular flow transport and deposition in an MOCVD reactor with a commercial configuration, and (2) the kinetics of metal alkyl source gas adsorption, decomposition and desorption from the growing film surface using ultra high vacuum surface science reaction techniques. To explore the transport-reaction issue, we have developed a reaction engineering analysis of a multiple wafer-in-tube ultrahigh vacuum chemical vapor deposition (UHV/CVD) reactor which allows an estimate of wafer or substrate throughput for a reactor of fixed geometry and a given deposition chemistry with specified film thickness uniformity constraints. The model employs a description of ballistic transport and reaction based on the pseudo-steady approximation to the Boltzmann equation in the limit of pure molecular flow. The model representation takes the form of an integral equation for the flux of each reactant or intermediate species to the wafer surfaces. Expressions for the reactive sticking coefficients (RSC) for each species must be incorporated in the term which represents reemission from a wafer surface. The interactions of MOCVD precursors with Si and CdTe were investigated using temperature programmed desorption (TPD) in ultra high vacuum combined with Auger electron spectroscopy (AES). These studies revealed that diethyltellurium (DETe) and dimethylcadmium (DMCd) adsorb weakly on clean Si(100) and desorb upon heating without decomposing. These precursors adsorb both weakly and strongly on CdTe(111)A, with DMCd exhibiting the stronger interaction with the surface than DETe.

Raupp, Gregory B.↗

DOE OSTI2022

The impact of neutrino-nucleus interaction modeling on new physics searches

Accurate neutrino-nucleus interaction modeling is an essential requirement for the success of the accelerator-based neutrino program. As no satisfactory description of cross sections exists, experiments tune neutrino-nucleus interactions to data to mitigate mis-modeling. In this work, we study how the interplay between near detector tuning and cross section mis-modeling affects new physics searches. We perform a realistic simulation of neutrino events and closely follow NOvA’s tuning, the first published of such procedures in a neutrino experiment. We analyze two illustrative new physics scenarios, sterile neutrinos and light neutrinophilic scalars, presenting the relevant experimental signatures and the sensitivity regions with and without tuning. While the tuning does not wash out sterile neutrino oscillation patterns, cross section mis-modeling can bias the experimental sensitivity. In the case of light neutrinophilic scalars, variations in cross section models completely dominate the sensitivity regardless of any tuning. Our findings reveal the critical need to improve our theoretical understanding of neutrino-nucleus interactions, and to estimate the impact of tuning on new physics searches. We urge neutrino experiments to follow NOvA’s example and publish the details of their tuning procedure, and to develop strategies to more robustly account for cross section uncertainties, which will expand the scope of their physics program.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1983

A comparison of wind observations of the upper thermosphere from the Dynamics Explorer satellite with the predictions of a global time-dependent model

Experimental and theoretical investigations of high-latitude circulation in the upper thermosphere are reported. Vector-wind plots constructed by combining remotely sensed meridional-wind and in situ zonal-wind data obtained during four south-polar and three north-polar passes of the NASA Dynamics Explorer satellite in October and December, 1981, are presented and analyzed, taking the geomagnetic activity level and local solar time into account. The results are then compared with the predictions of 3D time-dependent global models of thermospheric neutral winds (Fuller-Rowell and Rees, 1980, 1981, 1983) adjusted to account for both solar-UV/EUV heating and ionization effects and quiet (Kp from 1 to 2) or moderately disturbed (Kp from 3 to 4) geomagnetic conditions (models Q and MD). An MD model incorporating a self-consistent description of the high-latitude ionosphere and a Q model excluding high-latitude effects are found to give the most accurate predictions for the respective geomagnetic conditions.

Rees, D.↗

NASA NTRS2020

Designing and Training for Appropriate Trust in Increasingly Autonomous Advanced Air Mobility Operations: A Mental Model Approach: Version 1

To enable effective human-autonomy teaming (HAT) in Advanced Air Mobility (AAM) operations, the current paper presents a theoretical framework to design and train for appropriate trust in automation. The novel contribution of this work resides in connecting the construct of trust to mental models and showing how this method could be used to enable emerging HAT concepts such as Adaptive Trust Calibration. To contextualize this framework, in section 2 we discuss simplified vehicle operations (SVO) and remote vehicle operations (RVO), which are leading operational concepts within AAM. In section 3 we describe our perspective on automation and increasingly autonomous systems and present a brief discussion on human-automation interaction and human-autonomy teaming. In section 4 we provide a detailed discussion on the construct of trust in automation. In section 5 we present a framework that associates mental models with trust through principles of transparent design. Finally, in section 6 we present three descriptive models for designing and training for appropriate trust in increasingly autonomous systems.

Human-Autonomy Teaming↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

NASA NTRS1991

A guide to the use of theoretical models of the solar nebula for the interpretation of the meteoritic record

Attempts to derive a theoretical framework for the interpretation of the meteoritic record have been frustrated by our incomplete understanding of the fundamental processes that controlled the evolution of the primitive solar nebula. Nevertheless, it is possible to develop qualitative models of the nebula that illuminate its dynamic character, as well as the roles of some key parameters. These models draw on the growing body of observational data on the properties of disks around young, solar-type stars, and are constructed by applying the results of known solutions of protostellar collapse problems; making simple assumptions about the radial variations of nebular variables; and imposing the integral constraints demanded by conservation of mass, angular momentum, and energy. The models so constructed are heuristic, rather than predictive; they are intended to help us think about the nebula in realistic ways, but they cannot provide a definitive description of conditions in the nebula.

Cassen, Pat↗

NASA NTRS1999

Potential Flow Theory and Operation Guide for the Panel Code PMARC

The theoretical basis for PMARC, a low-order panel code for modeling complex three-dimensional bodies, in potential flow, is outlined. PMARC can be run on a wide variety of computer platforms, including desktop machines, workstations, and supercomputers. Execution times for PMARC vary tremendously depending on the computer resources used, but typically range from several minutes for simple or moderately complex cases to several hours for very large complex cases. Several of the advanced features currently included in the code, such as internal flow modeling, boundary layer analysis, and time-dependent flow analysis, including problems involving relative motion, are discussed in some detail. The code is written in Fortran77, using adjustable-size arrays so that it can be easily redimensioned to match problem requirements and computer hardware constraints. An overview of the program input is presented. A detailed description of the input parameters is provided in the appendices. PMARC results for several test cases are presented along with analytic or experimental data, where available. The input files for these test cases are given in the appendices. PMARC currently supports plotfile output formats for several commercially available graphics packages. The supported graphics packages are Plot3D, Tecplot, and PmarcViewer.

Ashby, Dale L.↗