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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Coupling MOOSE-Wrapped MPACT to BISON

As part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort is being made to leverage codes like Virtual Environment for Reactor Analysis (VERA) and Michigan Parallel Characteristics Transport (MPACT) by coupling them with other NEAMS codes. To facilitate coupling with other Multiphysics Objected Oriented Simulation Environment (MOOSE) applications, a MOOSE-wrapped MPACT app is created. The capability of this app, named Trogdor, is demonstrated by coupling it with another MOOSE app, BISON. Several single pin problems were run to test the Trogdor and BISON coupling.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2010

Utilizing In Situ Directional Hyperspectral Measurements to Validate Bio-Indicator Simulations for a Corn Crop Canopy

Two radiative transfer canopy models, SAIL and the two-layer Markov-Chain Canopy Reflectance Model (MCRM), were coupled with in situ leaf optical properties to simulate canopy-level spectral band ratio vegetation indices with the focus on the photochemical reflectance index in a cornfield. In situ hyperspectral measurements were made at both leaf and canopy levels. Leaf optical properties were obtained from both sunlit and shaded leaves. Canopy reflectance was acquired for eight different relative azimuth angles (psi) at three different view zenith angles (Theta (sub v)), and later used to validate model outputs. Field observations of photochemical reflectance index (PRI) for sunlit leaves exhibited lower values than shaded leaves, indicating higher light stress. Canopy PRI expressed obvious sensitivity to viewing geometry, as a function of both Theta (sub v) and psi . Overall, simulations from MCRM exhibited better agreements with in situ values than SAIL. When using only sunlit leaves as input, the MCRM-simulated PRI values showed satisfactory correlation and RMSE, as compared to in situ values. However, the performance of the MCRM model was significantly improved after defining a lower canopy layer comprised of shaded leaves beneath the upper sunlit leaf layer. Four other widely used band ratio vegetation indices were also studied and compared with the PRI results. MCRM simulations were able to generate satisfactory simulations for these other four indices when using only sunlit leaves as input; but unlike PRI, adding shaded leaves did not improve the performance of MCRM. These results support the hypothesis that the PRI is sensitive to physiological dynamics while the others detect static factors related to canopy structure. Sensitivity analysis was performed on MCRM in order to better understand the effects of structure related parameters on the PRI simulations. Leaf area index (LAI) showed the most significant impact on MCRM-simulated PRI among the parameters studied. This research shows the importance of hyperspectral and narrow band sensor studies, and especially the necessity of including the green wavelengths (e.g., 531 nm) on satellites proposing to monitor carbon dynamics of terrestrial ecosystems.

Cheng, Yen-Ben↗
NASA NTRSDate not supplied

Parametric Mechanism Design Through Numerical Optimization and Physics Simulation

Design-Build-Test approaches for developing spaceflight hardware are prohibitively time and cost intensive and often lead to suboptimal mechanism designs. Approaches that couple machine learning and high-fidelity physics simulation could eliminate the need for hardware prototyping and dramatically accelerate the engineering design cycle, ultimately reducing cost. This work presents a modular NASA-developed toolchain to optimize hardware mechanisms in a virtual environment using numerical optimization and multi-body physics simulation. The toolchain enables multi-objective optimization, generates parametric CAD files that can be further post-processed by an end user, and can be expanded to optimize full systems and non-mechanical parameters such as feedback control variables. We demonstrate the toolchain through an independently verifiable design problem that optimizes wheel radius to achieve a desired linear velocity in a rigid-body physics environment when the wheel rotates at a constant angular speed, and then post-process the parametric CAD file of the optimal design generated by the tool before ultimately manufacturing it via 3D printing. We end with a discussion of how the toolchain can incorporate other analysis tools, including finite element analysis, computational fluid dynamics, and granular media simulations.

Optimization↗
DOE OSTI2020

Generic characterization of electrical test benches for AC- and HVDC-connected wind power plants

The electrical test and assessment of wind turbines go hand in hand with standards and network connection requirements. In this paper, the generic structure of advanced electrical test benches, including grid emulator or controllable grid interface, wind torque emulator, and device under test, is proposed to harmonize state-of-the-art test sites. On the other hand, modern wind turbines are under development towards new features, concerning grid-forming, black-start, and frequency support capabilities as well as harmonic stability and control interaction considerations, to secure the robustness and stability of renewable-energy-based power systems. Therefore, it is necessary to develop new and revised test standards and methodologies to address the new features of wind turbines. This paper proposes a generic test structure within two main groups, including open-loop and closed-loop tests. The open-loop tests include the IEC 61400-21-1 standard tests as well as the additional proposed test options for the new capabilities of wind turbines, which replicate grid connection compliance tests using open-loop references for the grid emulator. In addition, the closed-loop tests evaluate the device under test as part of a virtual wind power plant and perform real-time simulations considering the grid dynamics. The closed-loop tests concern grid connection topologies consisting of AC and HVDC, as well as different electrical characteristics, including impedance, short-circuit ratio, inertia, and background harmonics. The proposed tests can be implemented using available advanced test benches by adjusting their control systems. The characteristics of a real power system can be emulated by a grid emulator coupled with real-time digital simulator systems through a high-bandwidth power-hardware-in-the-loop interface.

17 WIND ENERGY↗
DOE OSTI2021

A variational phase-field model For ductile fracture with coalescence dissipation

A novel phase-field model for ductile fracture is presented. The model is developed within a consistent variational framework in the context of finite-deformation kinematics. A novel coalescence dissipation introduces a new coupling mechanism between plasticity and fracture by degrading the fracture toughness as the equivalent plastic strain increases. The proposed model is compared with a recent alternative where plasticity and fracture are strongly coupled. Several representative numerical examples motivate specific modeling choices. In particular, a linear crack geometric function provides an “unperturbed” ductile response prior to crack initiation, and Lorentz-type degradation functions ensure that the critical fracture strength remains independent of the phase-field regularization length. In addition, the response of the model is demonstrated to converge with a vanishing phase-field regularization length. Additionally, the model is then applied to calibrate and simulate a three-point bending experiment of an aluminum alloy specimen with a complex geometry. The effect of the proposed coalescence dissipation coupling on simulations of the experiment is first investigated in a two-dimensional plane strain setting. The calibrated model is then applied to a three-dimensional calculation, where the calculated load-deflection curves and the crack trajectory show excellent agreement with experimental observations. Finally, the model is applied to simulate crack nucleation and growth in a specimen from a recent Sandia Fracture Challenge.

42 ENGINEERING↗
DOE OSTI2022

Coherently coupled quantum oscillators for quantum reservoir computing

Here, we analyze the properties of a quantum system composed of two coherently coupled quantum oscillators and show through simulations that it fulfills the two properties required for reservoir computing: non-linearity and fading memory. We first show that the basis states of this system apply a set of nonlinear transformations on the input signals and thus can implement neurons. We then show that the system exhibits a fading memory that can be controlled by its dissipation rates. Finally we show that a strong coupling between the oscillators is important in order to ensure complex dynamics and to populate a number of basis state neurons that is exponential in the number of physical devices.

Dudas, Julien↗
DOE OSTI2021

High spatial and energy resolution electron energy loss spectroscopy of the magnetic and electric excitations in plasmonic nanorod oligomers

We leverage the high spatial and energy resolution of monochromated aberration-corrected scanning transmission electron microscopy to study the hybridization of cyclic assemblies of plasmonic gold nanorods. Detailed experiments and simulations elucidate the hybridization of the coupled long-axis dipole modes into collective magnetic and electric dipole plasmon resonances. We resolve the magnetic dipole mode in these closed loop oligomers with electron energy loss spectroscopy and confirm the mode assignment with its characteristic spectrum image. The energy splitting of the magnetic mode and antibonding modes increases with the number of polygon edges (n). For the n=3-6 oligomers studied, optical simulations using normal incidence and s-polarized oblique incidence show the respective electric and magnetic modes’ extinction efficiencies are maximized in the n=4 arrangement.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Predictive dynamic wetting, fluid–structure interaction simulations for braze run-out

Brazing and soldering are metallurgical joining techniques that use a wetting molten metal to create a joint between two faying surfaces. Here, the quality of the brazing process depends strongly on the wetting properties of the molten filler metal, namely the surface tension and contact angle, and the resulting joint can be susceptible to various defects, such as run-out and underfill, if the material properties or joining conditions are not suitable. In this work, we implement a finite element simulation to predict the formation of such defects in braze processes. This model incorporates both fluid–structure interaction through an arbitrary Eulerian–Lagrangian technique and free surface wetting through conformal decomposition finite element modeling. Upon validating our numerical simulations against experimental run-out studies on a silver-Kovar system, we then use the model to predict run-out and underfill in systems with variable surface tension, contact angles, and applied pressure. Finally, we consider variable joint/surface geometries and show how different geometrical configurations can help to mitigate run-out. This work aims to understand how brazing defects arise and validate a coupled wetting and fluid–structure interaction simulation that can be used for other industrial problems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Lithographic quantum dots for quantum computation and quantum simulation

Lithographic quantum dots (QDs) are highly controllable few-level quantum systems created in semiconductor nanoelectronic devices, with a variety of scientific applications. These include technologically-driven applications like quantum computing and more fundamental applications in which they serve as a platform for exploring basic many-body physics. This document is a brief summary of my Ph.D. research so far and the directions with which I intend to continue it. Highlights include theoretical efforts to understand and design qubits in germanium hole QDs, as well as explorations of the possibility of using QDs coupled to nearby baths for analog simulation of quantum impurity models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Fresh Water Content Variability in the Arctic Ocean

Arctic Ocean model simulations have revealed that the Arctic Ocean has a basin wide oscillation with cyclonic and anticyclonic circulation anomalies (Arctic Ocean Oscillation; AOO) which has a prominent decadal variability. This study explores how the simulated AOO affects the Arctic Ocean stratification and its relationship to the sea ice cover variations. The simulation uses the Princeton Ocean Model coupled to sea ice. The surface forcing is based on NCEP-NCAR Reanalysis and its climatology, of which the latter is used to force the model spin-up phase. Our focus is to investigate the competition between ocean dynamics and ice formation/melt on the Arctic basin-wide fresh water balance. We find that changes in the Atlantic water inflow can explain almost all of the simulated fresh water anomalies in the main Arctic basin. The Atlantic water inflow anomalies are an essential part of AOO, which is the wind driven barotropic response to the Arctic Oscillation (AO). The baroclinic response to AO, such as Ekman pumping in the Beaufort Gyre, and ice meldfreeze anomalies in response to AO are less significant considering the whole Arctic fresh water balance.

Hakkinen, Sirpa↗
DOE OSTI2022

Access to Ru(IV)–Ru(V) and Ru(V)–Ru(VI) Redox in Layered Li 7 RuO 6 via Intercalation Reactions

The push to increase the capacity of oxides to store charge as Li-ion battery cathodes has challenged our knowledge of the boundaries of redox chemistry and its relationship to the atomic structure. Layered phases with Ru are now classical models in the quest for refined theories. Current understanding dictates that the highest formal state that can be reached in them is Ru(V), after which activity centered at the oxide anions is triggered. Yet, this picture is challenged by this investigation of the intercalation chemistry of Li 7 RuO 6 . Coupling structural and spectroscopic analysis with computational simulations, we show that Li 7 RuO 6 is capable of highly reversible lithiation to Li 8 RuO 6 through the Ru(V)/Ru(IV) redox couple. Further, Li 7 RuO 6 can also undergo anodic Li deintercalation, showing clear evidence of the reversible formation of Ru(VI) with octahedral coordination, an unusual state that is not accessible in other layered oxides. The results highlight the versatility of the Ru-O bond to undergo distinct redox transitions depending on the specific layered arrangements. They enrich our understanding of redox chemistry in solids while underscoring the need for descriptions of charge compensation that reflect the nuance of covalent interactions.

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

Sea surface warming and ocean-to-atmosphere feedback driven by large-scale offshore wind farms under seasonally stratified conditions

Offshore wind farms may induce changes in the upper ocean and near-surface atmosphere through coupled ocean-atmosphere feedbacks. Yet, the role of air-sea interactions mediated by offshore wind farms remains poorly understood. Using fully coupled ocean-atmosphere-wave model simulations for seasonally stratified conditions along the US East Coast, we show that simulated cumulative reductions in wind stress due to large-scale wind farm clusters lead to sea surface warming of 0.3° to 0.4°C and a shallower mixed layer. This warming drives upward heat fluxes, destabilizing the atmospheric boundary layer and enhancing wind stress, which partially offsets wake-induced wind deficits. These wake-ocean interactions influence near-surface meteorology and air-sea fluxes, suggesting that a coupled modeling approach may be necessary for assessing potential oceanographic impacts of offshore wind developments. However, ocean coupling exerts limited influence on winds at turbine-relevant heights or within downstream wakes, resulting in minimal impact on long-term energy. These findings suggest that models without ocean coupling may be adequate for wind energy applications.

17 WIND ENERGY↗
DOE OSTI2022

Simulation of Peach Bottom turbine trip transient with RETRAN-3D using BXFGEN

The Peach Bottom 2 Turbine Trip transient is used to evaluate the capabilities and accuracy of coupled 3-D neutronics/thermal hydraulics best estimate calculations for support of safety analysis. Selected results, obtained with RETRAN-3D, are compared to available measurements and benchmark results. The comparisons are analyzed and the application of RETRAN-3D for a best estimate analysis approach within the safety analysis methodology is discussed. Overall, both the thermal-hydraulic and neutronic results show a good comparison to the experimental data, and the model is considered to accurately predict test data. In particular, the initial acoustic response of the main steam line and the subsequent main pressurization stages are well-predicted. During the main pressurization stage, the timing and magnitude of the in-core local pressurization and associated void compression are properly predicted based upon the agreement between local power range monitors (LPRM) predictions upwards, but the shift is slightly underpredicted. This indicates the methodology is well-posed for simulating this class of transient response.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Partnership Center for High-fidelity Boundary Plasma Simulation

The University of Colorado Boulder contributed to the DOE Partnership Center for High‑fidelity Boundary Plasma Simulation (award DE‑SC0018271) by advancing gyrokinetic modeling of core–edge plasma interactions in tokamaks. Using the GEM and XGC codes, the project investigated turbulence spreading from the plasma edge to the core, finding that while edge‑driven trapped electron mode turbulence enhances core ion heat flux in DIII‑D L‑mode plasmas, it does not fully explain experimentally observed transport shortfalls, indicating a need for flux‑driven models. The team also demonstrated that edge micro‑instabilities are highly sensitive to small variations in magnetic equilibrium, particularly near the pedestal and edge, highlighting the importance of accurately characterizing equilibrium uncertainties. Significant progress was made in spatially coupling the core GEM and edge XGC simulations, with successful validation against standalone XGC results in realistic tokamak geometries, laying the groundwork for future extensions to kinetic electrons and electromagnetic effects.

42 ENGINEERING↗
DOE OSTI2022

Coupling Subsurface and Above-Surface Models for Optimizing the Design of Borefields and District Heating and Cooling Systems

Accurate dynamic energy simulation is important for the design and sizing of district heating and cooling systems with geothermal heat exchange for seasonal energy storage. Current modeling approaches in building and district energy simulation tools typically consider heat conduction through the ground between boreholes without flowing groundwater. While detailed simulation tools for subsurface heat and mass transfer exist, these fall short in simulating above-surface energy systems. To support the design and operation of such systems, the study developed a coupled model including a software package for building and district energy simulation, and software for detailed heat and mass transfer in the subsurface. For the first, it uses the open-source Modelica Buildings Library, which includes dynamic simulation models for building and district energy and control systems. For the heat and mass transfer in the soil, it uses the TOUGH simulator. The TOUGH family of codes can model heat and multi-phase, multi-component mass transport for a variety of fluid systems, as well as chemical reactions, in fractured porous media. The study validated the coupled modeling approach by comparing the simulation results with one from the g-function based ground response model. It then looked into effects when the water table and the regional groundwater flow are considered in the ground, from the perspective of heat exchange between borehole and ground, and the electrical consumption of the district heating and cooling systems. To access the simulation models, please find the links in the submission: -- For coupled approach validation: see model Buildings.Fluid.Geothermal.Borefields.Examples.BorefieldsWithTough and Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from the "Modelica Building Library" resource, branch issue1495_tough_interface, commit a2667c0. -- For the study of the effect of water table: see model Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO, commit 760de49. -- For the study of the effect of regional groundwater flow: see Buildings.Examples.DistrictReservoirNetworks.Examples.Reservoir3Variable_TOUGH from he "Modelica Building Library" resource, branch issue1495_tough_interface_moreIO_3D, commit c2a2d2a. The coupling interface script "GrounResponse.py" can be found from the above links in the folder Buildings/Resources/Python-Sources. Also, the needed files for TOUGH simulation are in the folder Buildings/Resources/Python-Sources/ToughFiles that can be accessed through the above links. A brief description of these files is given below; detailed specifications for the first three files may be found in the TOUGH3 Users Guide (Jung et al., 2018) https://tough.lbl.gov/documentation/tough-manuals/. (1) INCON - initial conditions for each grid block (2) INFILE - main input file with material properties and control parameters (3) MESH - description of the computational grid (4) readsave - Modelica/TOUGH interface program: read the final output of TOUGH simulation after TOUGH time step and prepare for transfer to Modelica for next Modelica time step (5) readsave.inp - input parameters for program readsave (6) writeincon - Modelica/TOUGH interface program: write the output of Modelica after Modelica time step and prepare for transfer to TOUGH as initial conditions for the next TOUGH step (7) writeincon.inp - input parameters for program writeincon

15 GEOTHERMAL ENERGY↗
NASA NTRSDate not supplied

Fluid-Structure Interaction Simulations of the ASPIRE SR01 Supersonic Parachute

High-fidelity computational fluid dynamics (CFD) simulations have so far only played exploratory and supporting roles in the study and qualification of parachutes for planetary entry. The last few years have seen the maturation of coupled computational methods that are capable of modeling the complex fluid-structure interactions between a parachute canopy and the supersonic flow in the wake of an entry vehicle in flight conditions. One of the primary goals of these methods is to predict the peak opening load experienced by the parachute during inflation. The Launch, Ascent, and Vehicle Aerodynamics (LAVA)team is developing efficient, high-fidelity numerical methods to perform such challenging fluid-structure interaction simulations. A loose coupling approach is used to advance the solutions of a Cartesian ghost cell immersed boundary method CFD solver and a finite element computational structural dynamics (CSD) solver in space and time. The coupled solver is employed to simulate the ASPIRE SR01 flight test, where a build-to-print version of the Mars Science Laboratory parachute was inflated in supersonic conditions in the upper terrestrial atmosphere. The simulations conducted in the current paper predict a peak opening load that is within 10% of that from the flight test. Grid convergence with respect to the volume and structural domains is demonstrated, and less than 1% variation in the peak opening load is predicted between all grid resolutions.

ARMD↗
NASA NTRS2011

Influence of Aerosol Heating on the Stratospheric Transport of the Mt. Pinatubo Eruption

On June 15th, 1991 the eruption of Mt. Pinatubo (15.1 deg. N, 120.3 Deg. E) in the Philippines injected about 20 Tg of sulfur dioxide in the stratosphere, which was transformed into sulfuric acid aerosol. The large perturbation of the background aerosol caused an increase in temperature in the lower stratosphere of 2-3 K. Even though stratospheric winds climatological]y tend to hinder the air mixing between the two hemispheres, observations have shown that a large part of the SO2 emitted by Mt. Pinatubo have been transported from the Northern to the Southern Hemisphere. We simulate the eruption of Mt. Pinatubo with the Goddard Earth Observing System (GEOS) version 5 global climate model, coupled to the aerosol module GOCART and the stratospheric chemistry module StratChem, to investigate the influence of the eruption of Mt. Pinatubo on the stratospheric transport pattern. We perform two ensembles of simulations: the first ensemble consists of runs without coupling between aerosol and radiation. In these simulations the plume of aerosols is treated as a passive tracer and the atmosphere is unperturbed. In the second ensemble of simulations aerosols and radiation are coupled. We show that the set of runs with interactive aerosol produces a larger cross-equatorial transport of the Pinatubo cloud. In our simulations the local heating perturbation caused by the sudden injection of volcanic aerosol changes the pattern of the stratospheric winds causing more intrusion of air from the Northern into the Southern Hemisphere. Furthermore, we perform simulations changing the injection height of the cloud, and study the transport of the plume resulting from the different scenarios. Comparisons of model results with SAGE II and AVHRR satellite observations will be shown.

Aquila, Valentina↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗