Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Study”

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 757 records · Page 42

Development of a hydraulic model of the human systemic circulation

Physical and numeric models of the human circulation are constructed for a number of objectives, including studies and training in physiologic control, interpretation of clinical observations, and testing of prosthetic cardiovascular devices. For many of these purposes it is important to quantitatively validate the dynamic response of the models in terms of the input impedance (Z = oscillatory pressure/oscillatory flow). To address this need, the authors developed an improved physical model. Using a computer study, the authors first identified the configuration of lumped parameter elements in a model of the systemic circulation; the result was a good match with human aortic input impedance with a minimum number of elements. Design, construction, and testing of a hydraulic model analogous to the computer model followed. Numeric results showed that a three element model with two resistors and one compliance produced reasonable matching without undue complication. The subsequent analogous hydraulic model included adjustable resistors incorporating a sliding plate to vary the flow area through a porous material and an adjustable compliance consisting of a variable-volume air chamber. The response of the hydraulic model compared favorably with other circulation models.

NASA Program Biomedical Research and Countermeasur↗

Supercritical flow about a thick circular-arc airfoil

The supercritical flow about a biconvex circular-arc airfoil is being thoroughly documented at Ames Research Center in order to provide experimental test cases suitable for guiding and evaluating current and future computer codes. The effects of angle of attack, effects of leading and trailing-edge splitter plates, additional unsteady pressure fluctuation (buffeting) measurements and glow-field shadowgraphs, and application of an oil-film technique to display separated-wake streamlines were studied. Computed and measured pressure distributions for steady and unsteady flows, using a recent computer code representative of current methodology, are compared. It was found that the numerical solutions are often fundamentally incorrect in that only strong (shock-polar terminology) shocks are captured, whereas experimentally, both strong and weak shock waves appear.

Mcdevitt, J. B.↗

Reynolds Number Effects on Leading Edge Radius Variations of a Supersonic Transport at Transonic Conditions

A computational study focused on leading-edge radius effects and associated Reynolds number sensitivity for a High Speed Civil Transport configuration at transonic conditions was conducted as part of NASA's High Speed Research Program. The primary purposes were to assess the capabilities of computational fluid dynamics to predict Reynolds number effects for a range of leading-edge radius distributions on a second-generation supersonic transport configuration, and to evaluate the potential performance benefits of each at the transonic cruise condition. Five leading-edge radius distributions are described, and the potential performance benefit including the Reynolds number sensitivity for each is presented. Computational results for two leading-edge radius distributions are compared with experimental results acquired in the National Transonic Facility over a broad Reynolds number range.

Rivers, S. M. B.↗

An Experimental and Numerical Study of Icing Effects on the Performance and Controllability of a Twin Engine Aircraft

In September 1997 the National Transportation Safety Board (NTSB) requested assistance from the NASA Lewis Research Center (LeRC) Icing Branch in the investigation of an aircraft accident that was suspected of being caused by ice contamination. In response to the request NASA agreed to perform an experimental and computational study. The main activities that NASA performed were LERC Icing Research Tunnel (IRT) testing to define ice shapes and 2-D Navier-Stokes analysis to determine the performance degradation that those ice shapes would have caused. An IRT test was conducted in January 1998. Most conditions for the test were based upon raw and derived data from the Flight Data Recorder (FDR) recovered from the accident and upon the current understanding of the Meteorological conditions near the accident. Using a two-dimensional Navier-Stokes code, the flow field and resultant lift and drag were calculated for the wing section with various ice shapes accreted in the IRT test. Before the final calculations could be performed extensive examinations of geometry smoothing and turbulence were conducted. The most significant finding of this effort is that several of the five-minute ice accretions generated in the IRT were found by the Navier-Stokes analysis to produce severe lift and drag degradation. The information generated by this study suggests a possible scenario for the kind of control upset recorded in the accident. Secondary findings were that the ice shapes accreted in the IRT were mostly limited to the protected pneumatic boot region of the wing and that during testing, activation of the pneumatic boots cleared most of the ice.

Reehorst, A.↗

Deciphering the Olefin Isomerization-Polymerization Paradox of Palladium(II) Diimine Catalysts: Discovery of Simultaneous and Independent Pathways of Olefin Isomerization and Living Polymerization

This work elucidates a long-standing unexplained paradox commonly observed within the polymerization of α-olefin using palladium (Pd)(II)–diimine catalysts, in which isomerization and living polymerization of α-olefins are both observed. With a classical mechanistic understanding of these complexes, this behavior is often dismissed and interpreted as experimental error. Herein, we present a comprehensive mechanistic investigation into this phenomenon that supports the existence of a novel mechanistic pathway for Pd(II)–diimine complexes. Part one of the mechanistic study lays the foundation of the proposed mechanism, in which neutral Pd(II)–diimine complexes were found to exhibit a moderate to good catalytic activity for olefin isomerization of α-olefins despite the established notion that catalyst activation is required. Extensive experimental and computational studies reveal the possibility of a partial dissociation of the diimine ligand, which frees up one coordination site and enables coordination–insertion. This finding is significant as the coexistence of two reactive coordination sites at the palladium center becomes a valid proposal for the activated cationic Pd(II)–diimine complexes. In part two, we examined and validated the simultaneously observed α-olefin isomerization and living polymerization using the cationic Pd(II)–diimine catalyst, which supports the presence of two independent reaction pathways of isomerization and polymerization, respectively. Furthermore, the addition of a strong Lewis acid, such as AlCl 3 , accelerates the ligand dissociation and the consequential isomerization as it weakens the palladium–nitrogen bond through competitive binding. In part three, Lewis acid-triggered olefin isomerization-polymerization is employed to prepare living olefinic block copolymers and further synthesize novel polyolefin-polar block copolymers with unique architectures, distinct levels of branching, crystallinity, and polar functionality in a one-pot manner.

Catalysts↗

Computational Analysis of Hydraulic Efficiency of Michigan DOT Covers J and K

Drainage structures are used in urban street and highway systems to capture stormwater runoff. These structures typically consist of catch basins fitted with grates, inlets, or combination grate/inlet configurations that collect runoff and convey it through buried drainage systems. They are strategically placed within curb-and-gutter systems to enhance public safety by efficiently removing water from roadways and thereby reducing the risk of hydroplaning. The performance of drainage structures is commonly evaluated in terms of hydraulic efficiency, defined as the percentage of flow captured by the basin relative to the total flow reaching the structure. Understanding the hydraulic performance of these structures allows designers to properly space inlets, resulting in cost-effective designs that also help ensure the safety of the traveling public. MDOT uses a variety of drainage structures for runoff capture, as documented in Michigan Department of Transportation (MDOT) Drainage Manual. Many of these structures incorporate sinusoidal-type grates that are not addressed in HEC-22. Because physical modeling of these structures has been limited, further analysis is needed to verify their capture efficiency. Under current MDOT practice, the capture efficiency of these grates is estimated by assuming performance similar to that of a comparably sized reticuline grate described in HEC22. The first phase of this effort, titled Computational Analysis of Hydraulic Efficiency of Michigan DOT Cover C, focused on evaluating the hydraulic performance of MDOT’s Cover C grate. Cover C was selected as the initial test candidate because its sinusoidal pattern is representative of other MDOT grates, while it is typically used in high-volume, higher-speed applications. A similar version, Cover CX, is used on interstate highways but does not include transverse bars for bicycle safety. The current second phase of the study expands this work to evaluate MDOT’s Covers J and K. These grates were selected for additional analysis to further assess the hydraulic performance of MDOT drainage structures that are not directly represented by grate configurations in HEC-22. The results of this phase will build on the findings from the Cover C analysis and support improved understanding of the capture efficiency of MDOT’s standard drainage grates. This report is intended to serve as a companion document to the earlier study, Computational Analysis of Hydraulic Efficiency of Michigan DOT Cover C [3]. The present work applies the same overall CFD-based evaluation approach to MDOT Covers J and K and compares the resulting performance trends with those previously identified for Cover C. In particular, both studies assess on-grade interception efficiency, sag-location hydraulic capacity, and the effects of partial obstruction relative to HEC-22-based design estimates.

42 ENGINEERING↗

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of single-defect memristor in an MoS2 atomic sheet

Non-volatile resistive switching, also known as memristor effect, where an electric field switches the resistance states of a two-terminal device, has emerged as an important concept in the development of high-density information storage, computing and reconfigurable systems. The past decade has witnessed substantial advances in non-volatile resistive switching materials such as metal oxides and solid electrolytes. It was long believed that leakage currents would prevent the observation of this phenomenon for nanometre-thin insulating layers. However, the recent discovery of non-volatile resistive switching in two-dimensional monolayers of transition metal dichalcogenide and hexagonal boron nitride sandwich structures (also known as atomristors) has refuted this belief and added a new materials dimension owing to the benefits of size scaling. Here we elucidate the origin of the switching mechanism in atomic sheets using monolayer MoS 2 as a model system. Atomistic imaging and spectroscopy reveal that metal substitution into a sulfur vacancy results in a non-volatile change in the resistance, which is corroborated by computational studies of defect structures and electronic states. These findings provide an atomistic understanding of non-volatile switching and open a new direction in precision defect engineering, down to a single defect, towards achieving the smallest memristor for applications in ultra-dense memory, neuromorphic computing and radio-frequency communication systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Singlet Fission

This report summarizes key findings from DOE Award DE‑SC0007004 – Singlet Fission, documenting more than a decade of research led by Professor Josef Michl. The project advanced foundational understanding of singlet fission (SF) through integrated synthetic, spectroscopic, crystallographic, and computational studies. Across a diverse set of chromophores—including tetracenes, cibalackrot derivatives, BODIPY dimers, TDPPs, and fluorinated isobenzofurans—the team identified structural and packing features that enable or inhibit SF, clarified the roles of charge‑transfer intermediates and excimer formation, and developed simplified theoretical models and computational screening tools for predicting optimal geometries and excited‑state-energetics. The work significantly broadened design principles for SF‑active materials and provided insights relevant to next‑generation solar energy conversion and molecular photophysics.

14 SOLAR ENERGY↗

Estimation of unsteady lift on a pitching airfoil from wake velocity surveys

The results of a joint experimental and computational study on the flowfield over a periodically pitched NACA0012 airfoil, and the resultant lift variation, are reported in this paper. The lift variation over a cycle of oscillation, and hence the lift hysteresis loop, is estimated from the velocity distribution in the wake measured or computed for successive phases of the cycle. Experimentally, the estimated lift hysteresis loops are compared with available data from the literature as well as with limited force balance measurements. Computationally, the estimated lift variations are compared with the corresponding variation obtained from the surface pressure distribution. Four analytical formulations for the lift estimation from wake surveys are considered and relative successes of the four are discussed.

Zaman, K. B. M. Q.↗

Estimation of unsteady lift on a pitching airfoil from wake velocity surveys

The results of a joint experimental and computational study on the flowfield over a periodically pitched NACA0012 airfoil, and the resultant lift variation, are reported in this paper. The lift variation over a cycle of oscillation, and hence the lift hysteresis loop, is estimated from the velocity distribution in the wake measured or computed for successive phases of the cycle. Experimentally, the estimated lift hysteresis loops are compared with available data from the literature as well as with limited force balance measurements. Computationally, the estimated lift variations are compared with the corresponding variation obtained from the surface pressure distribution. Four analytical formulations for the lift estimation from wake surveys are considered and relative successes of the four are discussed.

Zaman, K. B. M. Q.↗

Molecular Simulations in Astrobiology

One of the main goals of astrobiology is to understand the origin of cellular life. The most direct approach to this problem is to construct laboratory models of protocells. Such efforts, currently underway in the NASA Astrobiology Program, are accompanied by computational studies aimed at explaining self-organization of simple molecules into ordered structures that are capable of performing protocellular functions. Many of these functions, such as importing nutrients, capturing energy and responding to changes in the environment, are carried out by proteins bound to membranes. We use computer simulations to address the following questions about these proteins: (1) How do small proteins self-organize into ordered structures at water-membrane interfaces and insert into membranes? (2) How do peptides form membrane-spanning structures (e.g. channels)? (3) By what mechanisms do such structures perform their functions? The simulations are performed using the molecular dynamics method. In this method, Newton's equations of motion for each atom in the system are solved iteratively. At each time step, the forces exerted on each atom by the remaining atoms are evaluated by dividing them into two parts. Short-range forces are calculated in real space while long-range forces are evaluated in reciprocal space, using a particle-mesh algorithm which is of order O(NInN). With a time step of 2 femtoseconds, problems occurring on multi-nanosecond time scales (10(exp 6)-10(exp 8) time steps) are accessible. To address a broader range of problems, simulations need to be extended by three orders of magnitude, which requires algorithmic improvements and codes scalable to a large number of processors. Work in this direction is in progress. Two series of simulations are discussed. In one series, it is shown that nonpolar peptides, disordered in water, translocate to the nonpolar interior of the membrane and fold into helical structures (see Figure). Once in the membrane, the peptides exhibit orientational flexibility with changing conditions, which may have provided a mechanism of transmitting signals between the protocell and its environment. In another series of simulations, the mechanism by which a simple protein channel efficiently mediates proton transport across membranes was investigated. This process is a key step in cellular bioenergetics. In the channel under study, proton transport is gated by four histidines that occlude the channel pore. The simulations identify the mechanisms by which protons move through the gate.

Pohorille, Andrew↗

Substructure analysis using NICE/SPAR and applications of force to linear and nonlinear structures

Parallel computing studies are presented for a variety of structural analysis problems. Included are the substructure planar analysis of rectangular panels with and without a hole, the static analysis of space mast, using NICE/SPAR and FORCE, and substructure analysis of plane rigid-jointed frames using FORCE. The computations are carried out on the Flex/32 MultiComputer using one to eighteen processors. The NICE/SPAR runstream samples are documented for the panel problem. For the substructure analysis of plane frames, a computer program is developed to demonstrate the effectiveness of a substructuring technique when FORCE is enforced. Ongoing research activities for an elasto-plastic stability analysis problem using FORCE, and stability analysis of the focus problem using NICE/SPAR are briefly summarized. Speedup curves for the panel, the mast, and the frame problems provide a basic understanding of the effectiveness of parallel computing procedures utilized or developed, within the domain of the parameters considered. Although the speedup curves obtained exhibit various levels of computational efficiency, they clearly demonstrate the excellent promise which parallel computing holds for the structural analysis problem. Source code is given for the elasto-plastic stability problem and the FORCE program.

Razzaq, Zia↗

On the three families of instability waves of high-speed jets

An analytical and computational study of the normal-mode small-amplitude waves of high-speed jets is presented. Three families of instability waves have been identified: (1) the familiar Kelvin-Helmholtz instability waves; (2) supersonic instability waves; and (3) subsonic waves. It is demonstrated that the computed wave patterns and propagation characteristics of these three wave types are consistent with the findings of Oetel (1979, 1980, 1982). The subsonic waves are shown to be unstable only for jets with mixing layers of finite thickness.

Tam, Christopher K. W.↗

Three-Dimensional Imaging and Numerical Reconstruction of Graphite/Epoxy Composite Microstructure Based on Ultra-High Resolution X-Ray Computed Tomography

A combined experimental and computational study aimed at high-resolution 3D imaging, visualization, and numerical reconstruction of fiber-reinforced polymer microstructures at the fiber length scale is presented. To this end, a sample of graphite/epoxy composite was imaged at sub-micron resolution using a 3D X-ray computed tomography microscope. Next, a novel segmentation algorithm was developed, based on concepts adopted from computer vision and multi-target tracking, to detect and estimate, with high accuracy, the position of individual fibers in a volume of the imaged composite. In the current implementation, the segmentation algorithm was based on Global Nearest Neighbor data-association architecture, a Kalman filter estimator, and several novel algorithms for virtualfiber stitching, smoothing, and overlap removal. The segmentation algorithm was used on a sub-volume of the imaged composite, detecting 508 individual fibers. The segmentation data were qualitatively compared to the tomographic data, demonstrating high accuracy of the numerical reconstruction. Moreover, the data were used to quantify a) the relative distribution of individual-fiber cross sections within the imaged sub-volume, and b) the local fiber misorientation relative to the global fiber axis. Finally, the segmentation data were converted using commercially available finite element (FE) software to generate a detailed FE mesh of the composite volume. The methodology described herein demonstrates the feasibility of realizing an FE-based, virtual-testing framework for graphite/fiber composites at the constituent level.

Czabaj, M. W.↗

Study of a Satellite Attitude Control System Using Integrating Gyros as Torque Sources

This report considers the use of single-degree-of-freedom integrating gyros as torque sources for precise control of satellite attitude. Some general design criteria are derived and applied to the specific example of the Orbiting Astronomical Observatory. The results of the analytical design are compared with the results of an analog computer study and also with experimental results from a low-friction platform. The steady-state and transient behavior of the system, as determined by the analysis, by the analog study, and by the experimental platform agreed quite well. The results of this study show that systems using integrating gyros for precise satellite attitude control can be designed to have a reasonably rapid and well-damped transient response, as well as very small steady-state errors. Furthermore, it is shown that the gyros act as rate sensors, as well as torque sources, so that no rate stabilization networks are required, and when no error sensor is available, the vehicle is still rate stabilized. Hence, it is shown that a major advantage of a gyro control system is that when the target is occulted, an alternate reference is not required.

White, John S.↗

Space station Simulation Computer System (SCS) study for NASA/MSFC. Volume 3: Refined conceptual design report

The results of the refined conceptual design phase (task 5) of the Simulation Computer System (SCS) study are reported. The SCS is the computational portion of the Payload Training Complex (PTC) providing simulation based training on payload operations of the Space Station Freedom (SSF). In task 4 of the SCS study, the range of architectures suitable for the SCS was explored. Identified system architectures, along with their relative advantages and disadvantages for SCS, were presented in the Conceptual Design Report. Six integrated designs-combining the most promising features from the architectural formulations-were additionally identified in the report. The six integrated designs were evaluated further to distinguish the more viable designs to be refined as conceptual designs. The three designs that were selected represent distinct approaches to achieving a capable and cost effective SCS configuration for the PTC. Here, the results of task 4 (input to this task) are briefly reviewed. Then, prior to describing individual conceptual designs, the PTC facility configuration and the SSF systems architecture that must be supported by the SCS are reviewed. Next, basic features of SCS implementation that have been incorporated into all selected SCS designs are considered. The details of the individual SCS designs are then presented before making a final comparison of the three designs.

Source record↗