Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

COG User's Manual: A Multiparticle Monte Carlo Transport Code (Sixth Edition)

COG is a high-resolution code for the Monte Carlo simulation of coupled particle transport in arbitrary 3-D geometry. COG will transport neutrons, protons, deuterons, alpha particles with energies up to hundreds of GeV, and photons with energy ranges limited by the available cross section sets and physics models. Electrons can be transported via the EGS5 electron transport kernel, electrons can also be transported. The COG code is a significant upgrade from earlier Monte Carlo transport codes and has been written specifically to make it more versatile, accurate, and easy to use. COG has provisions for calculating deep penetration (shielding) problems, criticality problems, and neutron activation problems while retains all of the standard capabilities found in other Monte Carlo transport codes. COG uses high-resolution pointwise cross-section databases and makes no compromises in the transport physics, so that the results of a COG run are limited only by the accuracy of the databases used. COG runs primarily on Linux Operating System workstations with MPICH software installed – currently, Red Hat 7 & 8, Windows 10 (Windows Subsystem for Linux –WSL), Ubuntu 16, 18 & 20, OpenSUSE Leap 15.2, Fedora 32, Apple Power Mac with Intel CPU (with MacPorts installed) workstations, and LLNL LC supercomputer CTS-1 cluster with TOSS 3 are supported.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

COG11.3 Abstract

COG is a high-resolution code for the Monte Carlo simulation of coupled particle transport in arbitrary3-D geometry. COG will transport neutrons, protons, deuterons, alpha particles with energies up to hundreds of GeV, and photons with energy ranges limited by the available cross section sets and physics models. Electrons can be transported via the EGS5 electron transport kernel, electrons can also be transported. The COG code is a significant upgrade from earlier Monte Carlo transport codes and has been written specifically to makes it more versatile, accurate, and easy to use. COG has provisions for calculating deep penetration (shielding) problems, criticality problems, and neutron activation problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗
DOE OSTI2025

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗
DOE OSTI2025

Toward Higher-order Accuracy in Self-gravitating Hydrodynamics

High-order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution Δ x ) with only a moderate increase to computational expense. Significant effort has been placed in the development of finite-volume algorithms for (magneto)hydrodynamics; however, state-of-the-art astrophysical simulations tightly couple a plenitude of physics, additionally including gravity, photon transport, cosmic-ray transport, chemistry, and/or diffusion, to name a few. Algorithms frequently operator-split this additional physics (often a first-order error in time) and/or adopt a model wherein their evaluation is limited to second-order accuracy in space. In this work, we present a fourth-order-accurate finite-volume scheme for self-gravitating hydrodynamics on a uniform Cartesian grid. The method supplies source terms for the gravitational acceleration ( ρ g ) and gravitational energy release ( ρ v · g ) associated with fourth-order-accurate solutions to the Poisson equation. Our scheme (1) guarantees the conservation of total linear momentum while (2) decreasing (in proportion to Δ x 4 ) the effects of spurious heating and/or cooling associated with truncation error in the gravity. We demonstrate expected convergence rates for the algorithm by measuring errors in test problems evolving self-gravity modified linear waves and 3D polytropic equilibria. We test robustness of the algorithm by integrating an induced “inside-out” adiabatic collapse. We also discuss a method to smoothly downgrade the solution to second-order spatial accuracy to avoid spurious overshoots near steep density and/or pressure gradients.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Can we safely go to 200 °C? An integrated approach to assessing impacts to the engineered barrier system in a high-temperature repository

Abstract. This presentation gives on overview of the complex thermo-hydro-mechanical-chemical (THMC) processes occurring during the disposal of heat-producing high-level radioactive waste in geologic repositories. A specific focus is on the role of compacted bentonite, which is commonly used as an engineered backfill material for emplacement tunnels because of its low permeability, high swelling pressure, and radionuclide retention capacity. Laboratory and field tests integrated with THMC modeling have provided an effective way to deepen our understanding of temperature-related perturbations in the engineered barrier system; however, most of this work has been conducted for maximum temperatures around 100 ∘C. In contrast, some international disposal programs have recently started investigations to understand whether local temperatures in the bentonite of up to 200 ∘C could be tolerated with no significant changes to safety relevant properties. Raising the maximum temperature is attractive for economical and safety reasons but faces the challenge of exposing the bentonite to significant temperature increases. Strong thermal gradients may induce complex moisture transport processes while geochemical processes, such as cementation and perhaps also illitization effects may occur, all of which could strongly affect the bentonite and near-field rock properties. Here, we present initial investigations of repository behavior exposed to strongly elevated temperatures. We will start discussing our current knowledge base for temperature effects in repositories exposed to a maximum temperature of 100 ∘C, based on data and related modeling analysis from a large heater experiment conducted for over 18 years in the Grimsel Test Site in Switzerland. We then show results from coupled THMC simulations of a nuclear waste repository in a clay formation exposed to a maximum temperature of 200 ∘C. We also explore preliminary data from a bench-scale laboratory mock-up experiment, which was designed to represent the strong THMC gradients occurring in a “hot” repository, and we finally touch on a full-scale field heater test to be conducted soon in the Grimsel Test Site underground research laboratory in Switzerland (referred to as HotBENT).

Birkholzer, Jens T. (ORCID:0000000279891912)↗
DOE OSTI2025

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗
NASA NTRS1973

A roll-pitch interaction simulator and a control position command encoder for remote piloting of spin-entry research models

The Langley Research Center uses radio-controlled, scaled aircraft models to study the spin-entry characteristics of aircraft. Recent spin-entry studies required the use of an electronic proportional-control system for manipulating model control surfaces. In order to meet control system requirements, a special-purpose analog computer was designed to simulate the coupling between roll and pitch controls. A digital encoder was designed to encode the voltage analogs of control-surface position into a special pulse format for transmission to the model. This paper describes the two special developments and their relationship to the functions of the overall control system.

Meissner, C. W., Jr.↗
NASA NTRS1973

Convection in the earth's mantle

The computer is used to solve for thermal convection within the earth's mantle. A review of the knowledge of surface displacements and of the present understanding of the mantle and its relevant physical and chemical properties is contained in the paper. Applicable equations assume a Newtonian fluid layer heated from below and within, with gravity acting downward. The numerical method employs finite differences and was constructed with a view toward the faithful simulation of coupling mechanisms. It enables surveying the effect of a parameter using a relatively coarse computing mesh. Some of the results obtained are presented.

Turcotte, D. L.↗
NASA NTRS1977

Meteorological radar facility (MRF) slot conductance investigations

A preliminary meteorological radar facility (MRF) array design was completed in support of the slot conductance measurement program. Three different slot measurement techniques were evaluated. The selection of the probe comparison measurement technique was selected as the principal experimental method with the impedance measurement technique chosen to measure a few higher conductance slots to be used as reference slots. The impedance of 43 slots in 0.9 x 0.4 inch standard waveguide and of 40 slots in 0.835 x 0.4 inch waveguide was measured. Also, impedance measurements were made of a few slots using image planes to simulate mutual coupling effects. The measured and theoretical conductance, susceptance, and radiation phase data are presented in graphic form as a function of slot displacement for constant slot length, and of slot length for constant slot displacement. It is concluded that the proposed MRF array design approach is a feasible one.

Ratkevich, A. E.↗
NASA NTRS1985

Transonic separated solutions for an augmentor-wing

The viscous transonic flow about a multielement airfoil (augmentor-wing) is simulated by coupling full-potential and direct/inverse differential boundary-layer algorithms. Solutions have been obtained for a variety of conditions and are in fair agreement with available experimental data. Typical results from this transonic augmentor-wing code (TAUG-V) require approximately three minutes of CRAY-XMP CPU time. Since this viscous transonic code accounts for most of the important flow physics, yet is still economical, it is a practical tool for the design aerodynamicist.

Flores, J.↗
NASA NTRS1990

Numerical Analysis Of Flows With FIDAP

Report presents an evaluation of accuracy of Fluid Dynamics Package (FIDAP) computer program. Finite-element code for analysis of flows of incompressible fluids and transfers of heat in multidimensional domains. Includes both available methods for treatment of spurious numerical coupling between simulated velocity and simulated pressure; namely, penalty method and mixed-interpolation method with variable choices of interpolation polynomials for velocity and pressure. Streamwise upwind (STU) method included as option for flows dominated by convection.

Sohn, Jeong L.↗
NASA NTRS1990

The effect of interplanetary trajectory options on a manned Mars aerobrake configuration

Manned Mars missions originating in low Earth orbit (LEO) in the time frame 2010 to 2025 were analyzed to identify preferred mission opportunities and their associated vehicle and trajectory characteristics. Interplanetary and Mars atmospheric trajectory options were examined under the constraints of an initial manned exploration scenario. Two chemically propelled vehicle options were considered: (1) an all propulsive configuration, and (2) a configuration which employs aerobraking at Earth and Mars with low lift/drag (L/D) shapes. Both the interplanetary trajectory options as well as the Mars atmospheric passage are addressed to provide a coupled trajectory simulation. Direct and Venus swingby interplanetary transfers with a 60 day Mars stopover are considered. The range and variation in both Earth and Mars entry velocity are also defined. Two promising mission strategies emerged from the study: (1) a 1.0 to 2.0 year Venus swingby mission, and (2) a 2.0 to 2.5 year direct mission. Through careful trajectory selection, 11 mission opportunities are identified in which the Mars entry velocity is between 6 and 10 km/sec and Earth entry velocity ranges from 11.5 to 12.5 km/sec. Simulation of the Earth return aerobraking maneuver is not performed. It is shown that a low L/D configuration is not feasible for Mars aerobraking without substantial improvements in the interplanetary navigation system. However, even with an advanced navigation system, entry corridor and aerothermal requirements restrict the number of potential mission opportunities. It is also shown that for a large blunt Mars aerobrake configuration, the effects of radiative heating can be significant at entry velocities as low as 6.2 km/sec and will grow to dominate the aerothermal environment at entry velocities above 8.5 km/sec. Despite the additional system complexity associated with an aerobraking vehicle, the use of aerobraking was shown to significantly lower the required initial LEO weight. In comparison with an all propulsive mission, savings between 19 and 59 percent were obtained depending upon launch date.

Braun, Robert D.↗
NASA NTRS1990

Aeroelastic stability analysis of a high-energy turbine blade

The dynamic analysis for the SSME HPOTP first stage turbine blade is presented wherein the rotor aeroelastic stability is assessed. The method employs normal modes analysis to simulate the coupled blade/fluid system. A three-dimensional finite element model of the blade is used in conjunction with a two-dimensional linearized unsteady aerodynamic theory which accounts for steady aerodynamic loading effects. This unsteady aerodynamic model is applied in stacked axisymmetric strips along the airfoil span. The blade dynamic and aerodynamic behaviors are coupled within modal space by expressing the unsteady aerodynamic forces in the frequency domain. A complex eigenvalue problem is solved to determine the stability of the rotor assuming tuned blades. The present analysis indicates that the HPOTP rotor experiences very low aerodynamic damping in the first four vibrational modes. The edgewise mode was found to be dynamically unstable. This mode of the blade became stable when the effect of mechanical damping was considered.

Smith, Todd E.↗
NASA NTRS1991

Navier-Stokes simulation of a close-coupled canard-wing-body configuration

The thin-layer Navier-Stokes equations are solved for the flow about a coplanar close-coupled canard-wing-body configuration at a transonic Mach number of 0.90 and at angles of attack ranging from 0 to 12 degrees. The influence of the canard on the wing flowfield, including canard-wing vortex interaction and wing vortex breakdown, is investigated. A study of canard downwash and canard leading-edge vortex effects, which are the primary mechanisms of the canard-wing interaction, is emphasized. Comparisons between the computations and experimental measurements of surface pressure coefficients, lift, drag and pitching moment data are favorable. A grid refinement study for configurations with and without canard shows that accurate results are obtained using a refined grid for angles of attack where vortex burst is present. At an angle of attack of approximately 12 deg, favorable canard-wing interaction which delays wing vortex breakdown is indicated by the computations and is in good agreement with experimental findings.

Tu, Eugene L.↗
NASA NTRS1991

Electromechanical simulation and testing of actively controlled rotordynamic systems with piezoelectric actuators

A method is presented for simulating the coupled 'electromechanical' system to predict rotordynamic stability and unbalance response along with control system stability. The piezoelectric actuators and their amplifiers are represented as equivalent linear electrical circuits. The electromechanical system modeling approach is utilized to correlate test results from a double overhung rotor rig. The test results also show the effectiveness of the control system for suppressing the unbalance response of two modes using active stiffness and active damping.

Lin, Reng Rong↗