SEARCH · Engineering Papers

Results for “numerical modelling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Activating a Natural Fault Zone in the Swiss Alps

One major hurdle for understanding earthquake mechanics are observational limitations. Important phenomena like strain localisation, fault dilation, and fault healing are readily studied in rock mechanical laboratory experiments and with numerical models. At the scale of natural earthquakes, however, these phenomena are often unresolvable, even by state-of-the-art observatories. To overcome this limitation, we are currently building the Earthquake Physics Testbed at the Bedretto Underground Laboratory for Geosciences and Geoenergies (BedrettoLab), an experimental testbed where we can activate an extensively instrumented natural fault zone via hydraulic stimulation. The goal of the Fault Activation and Earthquake Rupture (FEAR) project is to induce earthquakes of up to Mw~1.0 on this exceptionally well characterised and instrumented fault zone. Here we summarize the main scientific goals and current FEAR project status, and present first results from conducted experiments. We discuss how this large-scale experimental approach may allow us to tackle both fundamental science as well as practical questions on earthquake physics, induced seismicity and seismic hazard.

Meier, Men-Andrin↗

DOE OSTI2024

A Novel Approach to Investigate Thermal Protection Systems Materials

The Koo Research Group (KRG) at The University of Texas at Austin (UT) and KAI has specialized in “Ablation Research” for more than fifteen years. Recently, the group has developed several incredibly unique capabilities that can advance “Thermal Protection Systems (TPS) Materials Research & Development” using an integrated experimental and numerical approach. The paper aims to introduce the methodology KRG has developed to solve this challenging problem. It will discuss how the KRG develops “Process-Properties-Performance” relationships of novel TPS materials in a systematical approach using (a) processing and fabrication, (b) thermal characterization of properties, (c) aerothermal testing, (d) microstructures characterization and analysis, and (e) numerical modeling. Progress and challenges of this research will also be discussed.

Engineering↗

DOE OSTI2023

Editorial: Sorption processes in nuclear waste management: data knowledge management and new methodologies for data acquisition/prediction

A fundamental approach to Nuclear Waste Repository research involves the collection of experimental data in a laboratory setting, development of empirical and/or mechanistic numerical models representing those observations, and application of these models (or Reduced Order Models) into reactive transport and performance assessment models as predictive tools for informing society of impacts and risks associated with nuclear waste repository scenarios (Stevens et al., 2020). Therefore, the assimilation and interpretation of experimental data must take advantage of both new data and the rich historical data available in the literature and apply novel modeling approaches to improve predictive tools, particularly from the standpoint of uncertainty quantification, for nuclear waste repository performance assessment (Zavarin et al., 2022).

sorption↗

DOE OSTI2021

OC6 Phase Ib: Floating Wind Component Experiment for Difference-Frequency Hydrodynamic Load Validation

A new validation campaign was conducted at the W2 Harold Alfond Ocean Engineering Laboratory at the University of Maine to investigate the hydrodynamic loading on floating offshore wind substructures, with a focus on the low-frequency contributions that tend to drive extreme and fatigue loading in semisubmersible designs. A component-level approach was taken to examine the hydrodynamic loads on individual parts of the semisubmersible in isolation and then in the presence of other members to assess the change in hydrodynamic loading. A variety of wave conditions were investigated, including bichromatic waves, to provide a direct assessment of difference-frequency wave loading. An assessment of the impact of wave uncertainty on the loading was performed, with the goal of enabling validation with this dataset of numerical models with different levels of fidelity. The dataset is openly available for public use and can be downloaded from the U.S. Department of Energy Data Archive and Portal.

17 WIND ENERGY↗

DOE OSTI2022

Numerical and Experimental Study of a Novel Additively Manufactured Metal-Polymer Composite Heat-Exchanger for Liquid Cooling Electronics

In order to meet increasing power-dissipation requirements of the electronics industry, compact, low-cost, and lightweight heat exchangers (HXs) are desired. With proper design, materials, and manufacture, polymer composite heat exchangers could meet these requirements. This paper presents a novel crossflow air-to-water, low-cost, and lightweight metal-polymer composite HX. This HX, which is entirely additively manufactured, utilizes a novel cross-media approach that provides direct heat exchange between air and liquid sides by using connecting fins. A robust numerical model was developed, which includes the dimensional effects of additive manufacturing. The study consists of a simplified 3D CFD model based on ellipsoidal-shaped staggered tube banks for the laminar range. It then uses an analytical approach to compute entire HX performance. The model is validated experimentally within 8% for thermal performance, 12% for air-side impedance, and 18% for water-side impedance. Finally, HX is compared with a conventional CPU radiator and performs within 10% of the conventional unit for reasonable flow rates and pressure-drop ranges. Moreover, HX also provides added design and cost advantages over the conventional unit, which makes the HX a potential candidate for electronic cooling applications.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2023

Impact of Limited Degree of Freedom Drag Coefficients on a Floating Offshore Wind Turbine Simulation

The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.

17 WIND ENERGY↗

DOE OSTI2024

Assessment of Cementitious Composites for High-Temperature Geothermal Wells

High-temperature (HT) geothermal wells can provide green power 24 hours a day, 7 days a week. Under harsh environmental and operational conditions, the long-term durability requirements of such wells require special cementitious composites for well construction. This paper reports a comprehensive assessment of geothermal cement composites in cyclic pressure function laboratory tests and field exposures in an HT geothermal well (300–350 °C), as well as a numerical model to complement the experimental results. Performances of calcium–aluminate cement (CAC)-based composites and calcium-free cement were compared against the reference ordinary Portland cement (OPC)/silica blend. The stability and degradation of the tested materials were characterized by crystalline composition, thermo-gravimetric and elemental analyses, morphological studies, water-fillable porosity, and mechanical property measurements. All CAC-based formulations outperformed the reference blend both in the function and exposure tests. The reference OPC/silica lost its mechanical properties during the 9-month well exposure through extensive HT carbonation, while the properties of the CAC-based blends improved over that period. The Modified Cam-Clay (MCC) plasticity parameters of several HT cement formulations were extracted from triaxial and Brazilian tests and verified against the experimental results of function cyclic tests. These parameters can be used in well integrity models to predict the field-scale behavior of the cement sheath under geothermal well conditions.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Shear Band Formation in Thin-Film Multilayer Columns Under Compressive Loading: A Mechanistic Study

Micro-pillar compression is a popular experimental technique used for characterizing the mechanical behavior of nano- and micro-laminates. The compressive stress–strain response of the column-shaped thin-film composite can be measured, and the deformation and damage features can be revealed by post-test cross-section microscopy. The development of plastic instability in the form of localized strain concentration (shear bands), leading to eventual failure, is frequently observed. In the present study, a computational approach is used to illustrate the commonality of shear band formation from a continuum standpoint. Systematic finite element analyses are conducted, showing that the strain field tends to become localized once plastic yielding commences. Distinct shear offsets of the layered structure can be revealed from the numerical model, which is similar to those observed in experiments. The actual appearance of shear bands depends on the materials’ constitutive behavior and precise geometries. Post-yield strain hardening reduces the propensity of shear band formation, while strain softening enhances it. Imperfections such as the undulated layer geometry, as well as the frictional characteristics between the specimen and test apparatus, can also influence the shear band morphology and overall stress–strain response.

finite element modeling↗

DOE OSTI2020

Velocity Measurements in Channel Gas Flows in the Slip Regime by means of Molecular Tagging Velocimetry

Direct measurements of the slip velocity in rarefied gas flows produced by local thermodynamic non-equilibrium at the wall represent crucial information for the validation of existing theoretical and numerical models. In this work, molecular tagging velocimetry (MTV) by direct phosphorescence is applied to argon and helium flows at low pressures in a 1-mm deep channel. MTV has provided accurate measurements of the molecular displacement of the gas at average pressures of the order of 1 kPa. To the best of our knowledge, this work reports the very first flow visualizations of a gas in a confined domain and in the slip flow regime, with Knudsen numbers up to 0.014. MTV is cross-validated with mass flowrate measurements by the constant volume technique. The two diagnostic methods are applied simultaneously, and the measurements in terms of average velocity at the test section are in good agreement. Moreover, preliminary results of the slip velocity at the wall are computed from the MTV data by means of a reconstruction method

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

Ultra-Wideband Microwave Imaging System for Root Phenotyping

The roots are a vital organ for plant growth and health. The opaque surrounding environment of the roots and the complicated growth process means that in situ and non-destructive root phenotyping face great challenges, which thus spur great research interests. The existing methods for root phenotyping are either unable to provide high-precision and high accuracy in situ detection, or they change the surrounding root environment and are destructive to root growth and health. Thus,we propose and develop an ultra-wideband microwave scanning method that uses time reversal to achieve in situ root phenotyping nondestructively. To verify the method’s feasibility, we studied an electromagnetic numerical model that simulates the transmission signal of two ultra-wideband microwave antennas. The simulated signal of roots with different shapes shows the proposed system’s capability to measure the root size in the soil. Experimental validations were conducted considering three sets of measurements with different sizes, numbers and locations, and the experimental results indicate that the developed imaging system was able to differentiate root sizes and numbers with high contrast. The reconstruction from both simulations and experimental measurements provided accurate size estimation of the carrots in the soil, which indicates the system’s potential for root imaging.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Can Simple Machine Learning Tools Extend and Improve Temperature-Based Methods to Infer Streambed Flux?

Temperature-based methods have been developed to infer 1D vertical exchange flux between a stream and the subsurface. Current analyses rely on fitting physically based analytical and numerical models to temperature time series measured at multiple depths to infer daily average flux. These methods have seen wide use in hydrologic science despite strong simplifying assumptions including a lack of consideration of model structural error or the impacts of multidimensional flow or the impacts of transient streambed hydraulic properties. We performed a “perfect-model experiment” investigation to examine whether regression trees, with and without gradient boosting, can extract sufficient information from model-generated subsurface temperature time series, with and without added measurement error, to infer the corresponding exchange flux time series at the streambed surface. Using model-generated, synthetic data allowed us to assess the basic limitations to the use of machine learning; further examination of real data is only warranted if the method can be shown to perform well under these ideal conditions. We also examined whether the inherent feature importance analyses of tree-based machine learning methods can be used to optimize monitoring networks for exchange flux inference.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Onset of Plasmoid Reconnection during Magnetorotational Instability

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two- and three-dimensional numerical modeling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical (z) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal (rz) plane, driving magnetic reconnection when the Lundquist number S > 3 × 10 2 , making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at S = 5 × 10 3 , self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When nonaxisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets more unstable and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Double-power-law Feature of Energetic Particles Accelerated at Coronal Shocks

Recent observations have shown that in many large solar energetic particle (SEP) events the event-integrated differential spectra resemble double power laws. We perform numerical modeling of particle acceleration at coronal shocks propagating through a streamer-like magnetic field by solving the Parker transport equation, including protons and heavier ions. We find that for all ion species the energy spectra integrated over the simulation domain can be described by a double power law, and the break energy depends on the ion charge-to-mass ratio as E B ~ (Q/A) α , with α varying from 0.16 to 1.2 by considering different turbulence spectral indices. We suggest that the double-power-law distribution may emerge as a result of the superposition of energetic particles from different source regions where the acceleration rates differ significantly due to particle diffusion. The diffusion and mixing of energetic particles could also provide an explanation for the increase of Fe/O at high energies as observed in some SEP events. Although further mixing processes may occur, our simulations indicate that either a power-law break or rollover can occur near the Sun and predict that the spectral forms vary significantly along the shock front, which may be examined by upcoming near-Sun SEP measurements from the Parker Solar Probe and Solar Orbiter.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Large-scale Linear Magnetic Holes with Magnetic Mirror Properties in Hybrid Simulations of Solar Wind Turbulence

Magnetic holes (MHs) are coherent magnetic field dips whose size ranges from fluid to kinetic scale, ubiquitously observed in the heliosphere and in planetary environments. Despite the long-standing effort in interpreting the abundance of observations, the origin and properties of MHs are still debated. In this Letter, we investigate the interplay between plasma turbulence and MHs, using a 2D hybrid simulation initialized with solar wind parameters. We show that fully developed turbulence exhibits localized elongated magnetic depressions, whose properties are consistent with linear MHs frequently encountered in space. The observed MHs develop self-consistently from the initial magnetic field perturbations by trapping hot ions with large pitch angles. Ion trapping produces an enhanced perpendicular temperature anisotropy that makes MHs stable for hundreds of ion gyroperiods, despite the surrounding turbulence. We introduce a new quantity, based on local magnetic field and ion temperature values, to measure the efficiency of ion trapping, with potential applications to the detection of MHs in satellite measurements. We complement this method by analyzing the ion velocity distribution functions inside MHs. Our diagnostics reveal the presence of trapped gyrotropic ion populations, whose velocity distribution is consistent with a loss cone, as expected for the motion of particles inside a magnetic mirror. Our results have potential implications for the theoretical and numerical modeling of MHs.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Self-potential tomography preconditioned by particle swarm optimization—Self-potential monitoring and streamflow data acquired March 26–September 14, 2023 at East Fork Poplar Creek near Oak Ridge Tennessee

This data release contains self-potential (SP) monitoring data measured on the flood plain of East Fork Poplar Creek (East Fork) in Oak Ridge, Tennessee and streamflow data measured at streamgage EFK5.4 about 310 meters upstream from the SP monitoring site. Additionally, forward and inverse numerical modeling scripts used to model the electrical-potential field on the East Fork flood plain are provided. SP monitoring data included in this data release were measured at 39 different data-collection points on the east flood plain; 30 points were spaced 3-m apart along an 87-m profile parallel to the edge of the streambank, and 9 points were spaced 5-m apart along a 40-m profile approximately perpendicular to the streambank. The two profiles of SP data-collection points intersected at the approximate midpoint of the profile parallel to the streambank. Transient voltages were measured at each data-collection point every 60 seconds between 16:13 Eastern Standard Time (EST) on March 26, 2023, and 11:41 EST on September 14, 2023. Streamflow data included in this data release overlap the time-period of self-potential monitoring and were measured every 900 seconds between 16:23 on March 26, 2023, and 23:53 on September 14, 2023.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Long-range prediction and the stratosphere

Abstract. Over recent years there have been concomitant advances in the development of stratosphere-resolving numerical models, our understanding of stratosphere–troposphere interaction, and the extension of long-range forecasts to explicitly include the stratosphere. These advances are now allowing for new and improved capability in long-range prediction. We present an overview of this development and show how the inclusion of the stratosphere in forecast systems aids monthly, seasonal, and annual-to-decadal climate predictions and multidecadal projections. We end with an outlook towards the future and identify areas of improvement that could further benefit these rapidly evolving predictions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Impact of phase state and non-ideal mixing on equilibration timescales of secondary organic aerosol partitioning

Abstract. Evidence has accumulated that secondary organic aerosols (SOAs) exhibit complex morphologies with multiple phases that can adopt amorphous semisolid or glassy phase states. However, experimental analysis and numerical modeling on the formation and evolution of SOA still often employ equilibrium partitioning with an ideal mixing assumption in the particle phase. Here we apply the kinetic multilayer model of gas–particle partitioning (KM-GAP) to simulate condensation of semi-volatile species into a core–shell phase-separated particle to evaluate equilibration timescales of SOA partitioning. By varying bulk diffusivity and the activity coefficient of the condensing species in the shell, we probe the complex interplay of mass transfer kinetics and the thermodynamics of partitioning. We found that the interplay of non-ideality and phase state can impact SOA partitioning kinetics significantly. The effect of non-ideality on SOA partitioning is slight for liquid particles but becomes prominent in semisolid or solid particles. If the condensing species is miscible with a low activity coefficient in the viscous shell phase, the particle can reach equilibrium with the gas phase long before the dissolution of concentration gradients in the particle bulk. For the condensation of immiscible species with a high activity coefficient in the semisolid shell, the mass concentration in the shell may become higher or overshoot its equilibrium concentration due to slow bulk diffusion through the viscous shell for excess mass to be transferred to the core phase. Equilibration timescales are shorter for the condensation of lower-volatility species into semisolid shell; as the volatility increases, re-evaporation becomes significant as desorption is faster for volatile species than bulk diffusion in a semisolid matrix, leading to an increase in equilibration timescale. We also show that the equilibration timescale is longer in an open system relative to a closed system especially for partitioning of miscible species; hence, caution should be exercised when interpreting and extrapolating closed-system chamber experimental results to atmosphere conditions. Our results provide a possible explanation for discrepancies between experimental observations of fast particle–particle mixing and predictions of long mixing timescales in viscous particles and provide useful insights into description and treatment of SOA in aerosol models.

54 ENVIRONMENTAL SCIENCES↗