Engineering Papers⌕ Search

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.

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

At least 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Hybrid simulation of Alfvén wave parametric decay instability in a laboratory relevant plasma

Large-amplitude Alfvén waves are subject to parametric decays, which can have important consequences in space, astrophysical, and fusion plasmas. Although the Alfvén wave parametric decay instability (PDI) was predicted decades ago, observational evidence is limited, stimulating considerable interest in laboratory demonstration of the instability and associated numerical modeling. Here, we report an important step toward direct hybrid simulation of the Alfvén wave PDI in a laboratory plasma, using antenna-like wave injection of a circularly polarized wave and realistic wave-plasma parameters. Considering collisionless damping, we identify the threshold Alfvén wave amplitudes and frequencies required for triggering the instability. These threshold behaviors are corroborated by simple theoretical analysis, considering the balance between PDI growth and ion Landau damping and PDI development in a bounded plasma. Other effects not included in the present model such as finite transverse wave scale and ion–neutral collisions are briefly discussed. These hybrid simulations promise to be a useful tool for investigating laboratory Alfvén wave dynamics and may provide guidance for future laboratory demonstration of the PDI process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗
DOE OSTI2021

Absolute contrast estimation for soft X-ray photon fluctuation spectroscopy using a variational droplet model

Abstract X-ray photon fluctuation spectroscopy using a two-pulse mode at the Linac Coherent Light Source has great potential for the study of quantum fluctuations in materials as it allows for exploration of low-energy physics. However, the complexity of the data analysis and interpretation still prevent recovering real-time results during an experiment, and can even complicate post-analysis processes. This is particularly true for high-spatial resolution applications using CCDs with small pixels, which can decrease the photon mapping accuracy resulting from the large electron cloud generation at the detector. Droplet algorithms endeavor to restore accurate photon maps, but the results can be altered by their hyper-parameters. We present numerical modeling tools through extensive simulations that mimic previous x-ray photon fluctuation spectroscopy experiments. By modification of a fast droplet algorithm, our results demonstrate how to optimize the precise parameters that lift the intrinsic counting degeneracy impeding accuracy in extracting the speckle contrast. These results allow for an absolute determination of the summed contrast from multi-pulse x-ray speckle diffraction, the modus operandi by which the correlation time for spontaneous fluctuations can be measured.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Interplay of intercalation dynamics and lithium plating in monolithic and architectured graphite anodes during fast charging

Fast charging of high-capacity anodes is challenging due to lithium plating reactions, which lead to poor cycling performance and safety concerns. Thus, accurate predictions of plating onset and an understanding of this electrochemical process are crucial for robust battery design. However, the most commonly used models, based on porous electrode theory (e.g., the pseudo-2D model), are notoriously difficult to calibrate due to their complexity, limiting their predictive power. This work studies the process of lithium plating during fast charging of (small-particle) graphite half-cells by measuring local reaction progression and plating behavior using optical operando techniques. These experiments employ a realistic 1D graphite electrode geometry with commercially-relevant mass loading charged at fast charge rates. It is demonstrated that the local reaction progression and plating onset can not only be predicted accurately with a p2D numerical model, but that these processes follow a simple scaling law. Remarkably, the entire reaction histories of different electrodes charged at different rates (e.g., 160 μm thickness at 0.5C, 111 μm at 1C or 66 μm at 4C) were observed to have self-similar intercalation profiles. It is demonstrated that plating onset is in turn governed by the reaction profile which explains why both processes exhibit the same scaling behavior. Finally, operando measurements of local reaction dynamics are conducted for the first time in electrodes with channeled architectures, quantitatively determining how channels affect reaction uniformity and plating onset. Together, these results reveal underlying simplicity in the complex electrochemical environment of fast charging and lithium plating, improving understanding of this process. These fundamental insights are broadly applicable for design processes, modeling and experimental evaluation of lithium ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage (Task 5 Report Field Scale Stress Modeling)

This report describes research accomplishments achieved with funding provided through Department of Energy (DOE) Contract DE-FE0031686, for the project “A Non-Invasive Approach for Elucidating the Spatial Distribution of In Situ Stress in Deep Subsurface Geologic Formations Considered for CO 2 Storage.” The purpose of this DOE funding is to develop non-invasive methods to obtain important information about subsurface stresses. The overall research goal of the project was to develop and improve methods for determining the spatial distribution of stresses, including magnitude and orientation of the three principal stress components in the subsurface, based on new methods that extract stress information from seismic data combined with well measurements, extended with well tests and logs, and unified in a numerical model that permits computation of the full stress tensor throughout the domains under investigation.

58 GEOSCIENCES↗
DOE OSTI2021

The EGS Collab Project: Status and Accomplishments

The EGS Collab project, supported by the US Department of Energy, is addressing challenges in implementing enhanced geothermal systems (EGS). This includes improving understanding of the stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The project is performing intensively monitored rock stimulation and flow tests at the 10-m scale in an underground research laboratory. Data and observations from the field test are compared to simulations to understand processes and to build confidence in numerical modeling of the processes. In Experiment 1, we examined hydraulic fracturing an underground test bed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km. We drilled eight sub-horizontal boreholes in a well-characterized phyllite. Six of the boreholes were instrumented with many sensor types to allow careful monitoring of stimulation events and flow tests, and the other two boreholes were used for water injection and production. We performed a number of stimulations and flow tests in the testbed. Our monitoring systems allowed detailed observations and collection of numerous data sets of processes occurring during stimulation and during dynamic flow tests. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Data were rapidly analyzed, allowing adaptive control of the tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results. Many simulations were performed in near-real-time in conjunction with the field experiments, with more detailed process study simulations performed on a longer timeframe. Experiment 2 will examine hydraulic shearing in a test bed being built at the SURF at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1. Five sets of fracture orientations were considered in design, and three orientations seem to be consistently observed.

EGS Collab↗
DOE OSTI2023

Energy Model of an Air Source Heat Pump to Explore Performance Improvements under Cold Conditions: A Python Framework: Preprint

Replacing combustion appliances with heat pumps can be part of a strategy to electrify buildings. When their electricity is supplied by renewables, the adoption of heat pumps for space heating may help curb the carbon dioxide emissions inherent to traditional heating systems that burn natural gas or coal. When operating in very cold climates, however, heat pump performance deteriorates, sometimes requiring backup heating supplied by electric resistances or other sources to continue to function. To avoid the use of electric resistances, complementary technologies, like thermal energy storage, solar collectors, or alternative refrigerants, have been explored in recent years. In this study, we present an open-source Python-based numerical model developed to evaluate the annual hourly performance of an air source heat pump (ASHP) operating in the USA. Simulation results include coefficient of performance and heating capacity. The model, validated with experimental results, will be employed in the future to investigate the feasibility of integrating thermal energy storage and other technologies into ASHPs operating in cold climates.

cold climate↗
DOE OSTI2025

Sharp spectroscopic fingerprints of disorder in an incompressible magnetic state

Disorder significantly impacts the electronic properties of conducting quantum materials by inducing electron localization and thus altering the local density of states and electric transport. In insulating quantum magnetic materials, the effects of disorder are less understood and can drastically impact fluctuating spin states like quantum spin liquids. In the absence of transport tools, disorder is typically characterized using chemical methods or by semi-classical modeling of spin dynamics. This requires high magnetic fields that may not always be accessible. Here, we show that magnetization plateaus—incompressible states found in many quantum magnets—provide an exquisite platform to uncover small amounts of disorder, regardless of the origin of the plateau. Using optical magneto-spectroscopy on the Ising-Heisenberg triangular-lattice antiferromagnet K 2 Co(SeO 3 ) 2 exhibiting a 1/3 magnetization plateau, we identify sharp spectroscopic lines, the fine structure of which serves as a hallmark signature of disorder. Through analytical and numerical modeling, we show that these fingerprints not only enable us to quantify minute amounts of disorder but also reveal its nature—as dilute vacancies. Remarkably, this model explains all details of the thermomagnetic response of our system, including the existence of multiple plateaus. Our findings provide a new approach to identifying disorder in quantum magnets.

Infrared spectroscopy↗
DOE OSTI2021

Enhancement of Photovoltaic Current through Dark States in Donor–Acceptor Pairs of Tungsten–Based Transition Metal Di–Chalcogenides

As several photovoltaic materials experimentally approach the Shockley– Queisser limit, there has been a growing interest in unconventional materials and approaches with the potential to cross this efficiency barrier. One such candidate is dark state protection induced by the dipole–dipole interaction between molecular excited states. This phenomenon has been shown to significantly reduce carrier recombination rate and enhance photon-to-current conversion, in elementary models consisting of few interacting chromophore centers. Atomically thin 2D transition metal di-chalcogenides (TMDCs) have shown great potential for use as ultrathin photovoltaic materials in solar cells due to their favorable photon absorption and electronic transport properties. TMDC alloys exhibit tunable direct bandgaps and significant dipole moments. In this work, the dark state protection mechanism has been introduced to a TMDC based photovoltaic system with pure tungsten diselenide (WSe 2 ) as the acceptor material and the TMDC alloy tungsten sulfo-selenide (WSeS) as the donor material. Furthermore, our numerical model demonstrates the first application of the dark state protection mechanism to a photovoltaic material with a photon current enhancement of up to 35% and an ideal photon-to-current efficiency exceeding the Shockley–Queisser limit.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Preliminary Reporting of Thermophysical Property Measurements for the Ghareb Formation

Accurate knowledge of thermophysical properties of rock is vital to develop meaningful models of high level nuclear waste emplacement scenarios. The Israel Atomic Energy Commission is considering storing high level nuclear waste in the Ghareb formation, a porous kerogen bearing chalk. Sandia is supporting this effort with an evolving lab- based geomechanics testing program. We have completed measurements of thermal properties up to 275C and room temperature hydrostatic compaction measurements. We report thermal conductivity, thermal diffusivity, specific heat, and mass loss from our thermal measurements, and we report bulk moduli and porosity loss from our compaction measurements. These values are crucial for the numerical models to simulate heat transfer and formation compressibility around a heat generating repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

IC Project: w19_hossfault “Modelling of stick-slip behavior in sheared granular fault gouge & nonlinear elasticity behavior in cracked solid”

Low frequency earthquakes, non-volcanic tremor, and acoustic emissions are examples of weak seismic signals that may help detect major seismic events, i.e., earthquakes. In a laboratory setting, in which stick-slip events are simulated, acoustic emissions are detected far from the stick-slip events. In both, field and laboratory scale cases, the acoustic emission signals are sourced or detected in a volume remote from the volume that spawns the earthquake. Therefore, it is relevant to establish the causal relationship between signals detected on passive, remote monitors and the dynamics of the elastic structures that launch important seismic events. The work conducted under this IC allocation allowed us to utilize a numerical model that let us follow this causality, i.e., examine and connect the dynamics in a granular system (fault gouge), to signals detected on passive remote monitors. It was demonstrated that stress chains are key in the dynamics of granular systems and that their evolution is the source of the acoustic emission.

58 GEOSCIENCES↗
DOE OSTI2025

Dependence of Convective Cloud Microphysical Properties on Environmental Conditions during the TRACER and ESCAPE Field Campaigns: A Synergistic Approach of Observations, Machine Learning and Parcel Models

The sensitivity of convective clouds to aerosols and their interactions with environment, combined with limited observational constraints in parameterizations, introduces significant uncertainties in atmospheric models. Here, this study investigates the dependence of convective cloud microphysical properties on environmental conditions using a synergistic approach that combines unique observations from the TRACER and ESCAPE field campaigns, machine learning techniques, and parcel model simulations with a super-droplet microphysics scheme. A random forest algorithm identifies in-situ vertical velocity (w), temperature (T), and surface fine-mode aerosol mass concentration as the three most important environmental conditions influencing cloud properties including liquid water content (LWC), number concentration for particles with D max < 50 μm (N c ,<50), 50 μm ≤ D max ≤ 3000 μm (N c,50–3000 ), and droplet effective diameter (D e ). Results show that LWC, N c,<50 , and N c,50–3000 significantly increase with w in updrafts. Across w bins, as T decreases, LWC, D e , and N c,50–3000 increase, while N c,<50 decreases, which are closely linked to the distance above cloud bases. Warmer cloud bases yield higher LWC, greater N c,50–3000 , and smaller N c,<50 , while polluted environments produce greater N c,<50 . Parcel model simulations successfully replicate these observed dependencies. The simulation results indicate that warmer cloud bases enhance condensation generating larger droplets, and differences in droplet sizes are then amplified through collision-coalescence, resulting in a greater N c,50–3000 . Polluted conditions result in a greater N c,<50 primarily due to enhanced cloud condensation nuclei activation despite increased collision-coalescence rates compared to pristine conditions. This study provides observed quantitative patterns characterizing cloud microphysical properties as a function of key environmental parameters, offering valuable constraints for improving physics parameterizations and numerical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Comment on “Review of experimental studies of secondary ice production” by Korolev and Leisner (2020)

Abstract. This is a comment on the review by Korolev and Leisner (2020, hereafter KL2020). The only two laboratory/field studies ever to measure the breakup in ice–ice collisions for in-cloud conditions were negatively criticised by KL2020, as were our subsequent theoretical and modelling studies informed by both studies. First, hypothetically, even without any further laboratory experiments, such theoretical and modelling studies would continue to be possible, based on classical mechanics and statistical physics. They are not sensitive to the accuracy of lab data for typical situations, partly because the nonlinear explosive growth of ice concentrations continues until some maximum concentration is reached. To a degree, the same final concentration is expected regardless of the fragment number per collision. Second, there is no evidence that both lab/field observational studies characterising fragmentation in ice–ice collisions are either mutually conflicting or erroneous such that they cannot be used to represent this breakup in numerical models, contrary to the review. The fact that the ice spheres of one experiment were hail sized (2 cm) is not a problem if a universal theoretical formulation, such as ours, with fundamental dependencies, is informed by it. Although both lab/field studies involved head-on collisions, rotational kinetic energy for all collisions generally is only a small fraction of the initial collision kinetic energy (CKE) anyway. Although both lab/field experiments involved fixed targets, that is not a problem since the fixing of the target is represented via CKE in any energy-based formulation such as ours. Finally, scaling analysis suggests that the breakup of ice during sublimation can make a significant contribution to ice enhancement in clouds, again contrary to the impression given by the review.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A unified thermionic and thermionic-field emission (TE–TFE) model for ideal Schottky reverse-bias leakage current

We present a unified thermionic emission (TE) and thermionic-field emission (TFE) model for the ideal reverse-bias leakage current in Schottky junctions. The unified TE–TFE analytical model advances upon previous analytical TFE models by Murphy–Good and Padovani–Stratton, which are the two most widely adopted models by the community, in two major aspects: (i) the applicability of the TFE expression therein is extended to near-zero surface electric fields by an error-function correction, allowing for the calculation of the total current by a nontrivial sum of TE (over-the-barrier current) and TFE (below-the-barrier current) contributions; therefore, an accurate description of the TE-to-TFE transition region is captured analytically for the first time; (ii) image-force lowering is considered with much-simpler correction terms. Comparisons with the reference numerical model show that the unified TE–TFE model has excellent accuracy, as well as a 10 000× reduction in computation time. The unified model is further tested against experimental data from Schottky barrier diodes based on Si, 4H-SiC, GaN, and β-Ga 2 O 3 , revealing accurate extractions of barrier heights and correct descriptions of the ideal reverse leakage characteristics. With the extended applicable range, improved accuracy, and computational efficiency, the unified TE–TFE model is highly valuable for the design and analysis of devices based on Schottky junctions, as well as for potential integration in technology computer-aided design (TCAD) tools.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Fossil records of early solar irradiation and cosmolocation of the CAI factory: A reappraisal

Calcium-aluminum–rich inclusions (CAIs) in meteorites carry crucial information about the environmental conditions of the nascent Solar System prior to planet formation. Based on models of 50V– 10 Be co-production by in-situ irradiation, CAIs are considered to have formed within ~0.1 AU from the proto-Sun. Here, we present vanadium (V) and strontium (Sr) isotopic co-variations in fine- and coarse-grained CAIs and demonstrate that kinetic isotope effects during partial condensation and evaporation best explain V isotope anomalies previously attributed to solar particle irradiation. We also report initial excesses of 10 Be and argue that CV CAIs possess essentially a homogeneous level of 10 Be, inherited during their formation. Based on numerical modeling of 50V–10Be co-production by irradiation, we show that CAI formation during protoplanetary disk build-up likely occurred at greater heliocentric distances than previously considered, up to planet-forming regions (~1AU), where solar particle fluxes were sufficiently low to avoid substantial in-situ irradiation of CAIs.

58 GEOSCIENCES↗
DOE OSTI2024

Soil Desiccation Treatability Testing at BC Waste Disposal Cribs

During Hanford’s production period, low-level waste products generated from chemical processing of uranium fuel rods were discharged directly to the ground through a system of cribs and trenches located in the 200-BC-1 Operable Unit (OU). The site consists of 6 cribs and 20 trenches that received more than 117,000 m3 of radioactive liquid waste discharged to the soil. These unlined infiltration galleries held volumes of liquid waste while it seeped into the ground, with the understanding that the 100 m (330 ft) thick vadose zone in the area would effectively capture the effluent and prevent groundwater impacts. A conceptual model showing the operation of cribs and trenches is shown in Figure 1. Data show effluent from the 26 cribs and trenches containing about 410 curies of Technetium-99 (Tc-99) is primarily located between 30 m and 70 m (98 ft and 230 ft) depth (Corbin et al., 2005; Ward et al., 2004). Despite no evidence indicating that the contamination has reached the groundwater at BC cribs and trenches, the mobility of Tc-99 had been demonstrated in laboratory tests and was recognized as a threat to groundwater at the site. Using data from numerical models, laboratory analyses, field investigations, and information on historical discharges, the EPA and Ecology identified Tc-99 and U contamination of the vadose zone as a remediation priority. The U.S. DOE was notified by EPA and Ecology regarding risks associated with Tc-99 contamination in a letter requesting development of a strategy for improved methods to understand the nature and extent of vadose zone contamination, specifically Tc-99, and to develop remedial options for addressing such contamination. To develop the appropriate technology for characterizing, remediating, and monitoring the deep vadose zone Tc-99 contamination, the U.S. DOE worked with the EPA and Ecology to create a Treatability Test Plan under a Remedial Investigation/Feasibility Study (RI/FS) for the Hanford 200 Areas. Under this RI/FS, it was determined that a treatability test for soil desiccation should be carried out as it was identified as a promising in-situ treatment technology for mitigating risks posed by Tc-99 contamination to the groundwater table. The BC Cribs and Trenches site was identified as a representative site for Tc-99 and U contamination and selected for the soil desiccation treatability test. In this chapter, we summarize the overlying regulatory framework of RI/FS and treatability tests and illustrate how development and experimentation supported the evaluation of selected remedies. We briefly discuss the RI/FS for the 200 Areas of the Hanford Site and focus on the soil desiccation treatability testing performed at the BC cribs and trenches site under the Deep Vadose Zone Treatability Test Plan for the Hanford Central Plateau (DVZ-TT). The DVZ-TT is one component of the remedial investigation/feasibility study for the Hanford 200 Areas and represents the underlying regulatory framework that drives site operations towards records of decision and site closure.

Mangel, Adam R.↗
DOE OSTI2023

Gallium-doped zinc oxide: nonlinear reflection and transmission measurements and modeling in the ENZ region

Abstract Strong nonlinear materials have been sought after for decades for applications in telecommunications, sensing, and quantum optics. Gallium-doped zinc oxide is a II–VI transparent conducting oxide that shows promising nonlinearities similar to indium tin oxide and aluminum-doped zinc oxide for the telecommunications band. Here we explore its nonlinearities in the epsilon near zero (ENZ) region and show n 2,eff values on the order of 4.5 × 10 −3 cm 2 GW −1 for IR pumping on 200–300 nm thin films. Measuring nonlinear changes in transmission and reflection with a white light source probe in the near-IR while exciting in the near-IR provides data in both time and wavelength. Three films varying in thickness, optical loss, and ENZ crossover wavelength are numerically modeled and compared to experimental data showing agreement for both dispersion and temporal relaxation. In addition, we discuss optimal excitation and probing wavelengths occur around ENZ for thick films but are red-shifted for thin films where our model provides an additional degree of freedom to explore. Obtaining accurate nonlinear measurements is a difficult and time-consuming task where our method in this paper provides experimental and modeled data to the community for an ENZ material of interest.

Ball, Adam (ORCID:0000000318167813)↗
DOE OSTI2022

MODELLING OF ADDITIVE MANUFACTURING – COMPRESSION MOLDING PROCESS USING COMPUTATIONAL FLUID DYNAMICS

A computational fluid dynamics model has been developed to predict the behavior of a printed strands during a novel material extrusion additive manufacturing and compression molding process. While the traditional additive manufacturing process enables control over the fiber orientation within the part, it would also result in high void content. On the other hand, compression molding produces parts with low porosity levels at rapid processing cycle time but lacks control over the microstructure. The novel additive manufacturing – compression molding (AM-CM) integrating both these processes offers control over both the microstructure and porosity to manufacturing high performance composite parts. The numerical model developed here enables to analyze the effect of processing parameters on the behaviour of printed layer that is subsequently compressed and to determine the optimal printing parameters for high-performance composite part design.

Seta, Berin↗