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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Observation of Road Salt Aerosol Driving Inland Wintertime Atmospheric Chlorine Chemistry

Inland sources of particulate chloride for atmospheric nitryl chloride (ClNO2) formation remain unknown and unquantified, hindering air quality assessments. Globally each winter, tens of millions of tons of road salt are spread on roadways for deicing. Here, we identify road salt aerosol as the primary chloride aerosol source, accounting for 80–100% of ClNO 2 formation, at an inland urban area in the wintertime. This study provides experimental evidence of the connection between road salt and air quality through the production of this important reservoir for nitrogen oxides and chlorine radicals, which significantly impact atmospheric composition and pollutant fates. A numerical model was employed to quantify the contributions of chloride sources to ClNO 2 production. The traditional method for simulating ClNO 2 considers chloride to be homogeneously distributed across the atmospheric particle population; yet, we show that only a fraction of the particulate surface area contains chloride. Our new single-particle parametrization considers this heterogeneity, dramatically lowering overestimations of ClNO 2 levels that have been routinely reported using the prevailing methods. The identification of road salt as a ClNO 2 source links this common deicing practice to atmospheric composition and air quality in the urban wintertime environment.

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

Numerical Study of Mineral Nucleation and Growth on a Substrate

We report nucleation and growth processes of minerals and other crystals can significantly affect one another due to the transport limitations and local depletion of reactive ions in the solution. Most numerical models and experimental measurements have typically focused on either growth or nucleation, but not both. In this work, we incorporate a heterogeneous nucleation process based on classical nucleation theory into a microcontinuum model that implements the Darcy–Brinkman–Stokes approach to study the interplay between nucleation and crystal growth on a substrate in diffusive systems. We demonstrate how the Damköhler number (reaction rate) and nucleation rate prefactor change the effective nucleation rate on a substrate. Higher surface growth rates deplete the solute concentration around the nuclei that appear initially on the substrate, creating islands that screen against further nucleation. The model predicts that measured nucleation rates may be affected by the history of crystal nucleation on the substrate. In the extreme case of high growth rates relative to diffusion, it predicts that the rate of subsequent nucleation is limited by reactant depletion. We introduce a nondimensional number α to represent the relation between surface propagation rate during growth and the heterogeneous nucleation rate. We show that it is important to control Damköhler number and α to achieve similar precipitation regimes at different reaction and nucleation rates. We suggest that the observed universality can guide the interpretation of experimental results on nucleation rates, since matching experiment can be achieved by tuning transport, reaction, and nucleation parameters simultaneously. In addition, we show how the bulk solution concentration affects the structure and topology of precipitation on a substrate.

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

Harnessing Photoelectrochemistry for Wastewater Nitrate Treatment Coupled with Resource Recovery

Wastewater is a misplaced resource well suited to recover nutrients, value-added chemicals, energy, and clean water. A photoelectrochemical device is proposed to transform wastewater nitrates to ammonia and nitrous oxide, coupled with water oxidation. Numerical models were developed to quantify the dependence of process efficiencies and nitrogen-removal rates on light absorber band gaps, electrocatalytic kinetic parameters, competing oxygen reduction and hydrogen evolution reactions, and the reacting nitrate species concentrations that affect the mass-transfer limited current densities. With a single light-absorber and state-of-the-art catalysts, optimal solar-to-chemical efficiencies of 7% and 10% and nitrogen-removal rates of 260 and 395 gN m -2 day -1 are predicted for nitrate reduction to ammonia and nitrous oxide, respectively. The influence of competing reactions on the performance depends on the nitrate concentration and band gap of the light absorber modeled. Oxygen reduction is more dominant than hydrogen evolution to compete with the nitrate reduction reaction, but it is mass-transfer limited. Even with kinetic parameters that enhanced the driving forces for the competing reactions, the performance is only minimally affected by these reactions for optimally selected band gaps and nitrate concentrations larger than 100 mM. Here, theoretically predicted peak nitrogen removal rates and specific energy intensities are competitive with reported estimates for electrochemical and Sharon-Anammox processes for ammonia recovery and nitrogen removal, respectively. This result, together with the added benefit of harnessing sunlight to produce value-added products, indicates promise in the photoelectrochemical approach as a tertiary pathway to recover nutrients and energy from wastewater nitrates.

25 ENERGY STORAGE↗
DOE OSTI2020

Timescales and Processes of Methane Hydrate Formation and Breakdown, With Application to Geologic Systems

Abstract Gas hydrate is an ice‐like form of water and low molecular weight gas stable at temperatures of roughly −10°C to 25°C and pressures of ~3 to 30 MPa in geologic systems. Natural gas hydrates sequester an estimated one sixth of Earth's methane and are found primarily in deepwater marine sediments on continental margins, but also in permafrost areas and under continental ice sheets. When gas hydrate is removed from its stability field, its breakdown has implications for the global carbon cycle, ocean chemistry, marine geohazards, and interactions between the geosphere and the ocean‐atmosphere system. Gas hydrate breakdown can also be artificially driven as a component of studies assessing the resource potential of these deposits. Furthermore, geologic processes and perturbations to the ocean‐atmosphere system (e.g., warming temperatures) can cause not only dissociation, but also more widespread dissolution of hydrate or even formation of new hydrate in reservoirs. Linkages between gas hydrate and disparate aspects of Earth's near‐surface physical, chemical, and biological systems render an assessment of the rates and processes affecting the persistence of gas hydrate an appropriate Centennial Grand Challenge. This paper reviews the thermodynamic controls on methane hydrate stability and then describes the relative importance of kinetic, mass transfer, and heat transfer processes in the formation and breakdown (dissociation and dissolution) of gas hydrate. Results from numerical modeling, laboratory, and some field studies are used to summarize the rates of hydrate formation and breakdown, followed by an extensive treatment of hydrate dynamics in marine and cryospheric gas hydrate systems.

58 GEOSCIENCES↗
DOE OSTI2020

Fate of Internal Waves on a Shallow Shelf

Internal waves strongly influence the physical and chemical environment of coastal ecosystems worldwide. We report novel observations from a distributed temperature sensing (DTS) system that tracked the transformation of internal waves from the shelf break to the surf zone over a narrow shelf slope region in the South China Sea. The spatially continuous view of temperature fields provides a perspective of physical processes commonly available only in laboratory settings or numerical models, including internal wave reflection off a natural slope, shoreward transport of dense fluid within trapped cores, and observations of internal rundown (near-bed, offshore-directed jets of water preceding a breaking internal wave). Analysis shows that the fate of internal waves on this shelf—whether transmitted into shallow waters or reflected back offshore—is mediated by local water column density structure and background currents set by the previous shoaling internal waves, highlighting the importance of wave-wave interactions in nearshore internal wave dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulations of Winter Arctic Clouds and Associated Radiation Fluxes Using Different Cloud Microphysics Schemes in the Polar WRF: Comparisons With CloudSat, CALIPSO, and CERES

Clouds are important for the Arctic climate, but simulating such clouds with numerical models is still challenging. The accuracy of model clouds has not been sufficiently examined due to the harsh Arctic environment obstructing cloud observations, especially during Arctic winters experiencing polar nights. This study compares the Arctic winter clouds simulated by a weather forecast model to cloud observations from active (lidar and radar) satellite instruments. The model successfully produced cloud patterns similar to the satellite observations. However, the choice of the cloud physics module in the model can modify the amount of cloud water significantly enough to affect the simulated surface climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Four Theories of the Madden‐Julian Oscillation

Abstract Studies of the Madden‐Julian Oscillation (MJO) have progressed considerably during the past decades in observations, numerical modeling, and theoretical understanding. Many theoretical attempts have been made to identify the most essential processes responsible for the existence of the MJO. Criteria are proposed to separate a hypothesis from a theory (based on the first principles with quantitative and testable assumptions, able to predict quantitatively the fundamental scales and eastward propagation of the MJO). Four MJO theories are selected to be summarized and compared in this article: the skeleton theory, moisture‐mode theory, gravity‐wave theory, and trio‐interaction theory of the MJO. These four MJO theories are distinct from each other in their key assumptions, parameterized processes, and, particularly, selection mechanisms for the zonal spatial scale, time scale, and eastward propagation of the MJO. The comparison of the four theories and more recent development in MJO dynamical approaches lead to a realization that theoretical thinking of the MJO is diverse and understanding of MJO dynamics needs to be further advanced.

58 GEOSCIENCES↗
DOE OSTI2020

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗
DOE OSTI2021

Methane Hydrate Formation and Evolution During Sedimentation

We explored methane hydrate formation with sedimentation with a newly developed one-dimensional, multiphase flow, multicomponent transport numerical model. Our model couples methane hydrate formation from in situ microbial methane generation within the hydrate stability zone (HSZ), methane recycling, and microbial methane generation below the base of the hydrate stability zone (BHSZ). Both recycled methane and deeply generated methane are transported into the HSZ by buoyancy-driven free gas flow. Free gas flows through the HSZ by both the processes of capillary-dependent pore fillings and by salt exclusion during hydrate formation, with the former being the dominant mechanism. We quantitively illustrated the formation of enriched hydrate in muddy sediments above, and interconnected free gas below, the BHSZ, which are common features along the world's continental margin. In addition, we showed two ways to form concentrated methane hydrate above the BHSZ. The first mechanism is local free gas flow during methane recycling. This happens at sites with sufficient methane generation above the BHSZ. The second mechanism is deep microbial methane generation which is transported into the HSZ by free gas flow. This mechanism plays a more important role at sites with high sedimentation rates. This study provides new insights into methane hydrate formation and distribution below the seafloor. It is important for understanding the carbon cycle and carbon storage below the seafloor and for resource evaluation and exploitation.

58 GEOSCIENCES↗
DOE OSTI2020

A New Crater Near InSight: Implications for Seismic Impact Detectability on Mars

A new 1.5 m diameter impact crater was discovered on Mars only ~40 km from the InSight lander. Context camera images constrained its formation between 21 February and 6 April 2019; follow-up High Resolution Imaging Science Experiment images resolved the crater. During this time period, three seismic events were identified in InSight data. In this paper, we derive expected seismic signal characteristics and use them to evaluate each of the seismic events. However, none of them can definitively be associated with this source. Atmospheric perturbations are generally expected to be generated during impacts; however, in this case, no signal could be identified as related to the known impact. Using scaling relationships based on the terrestrial and lunar analogs and numerical modeling, we predict the amplitude, peak frequency, and duration of the seismic signal that would have emanated from this impact. The predicted amplitude falls near the lowest levels of the measured seismometer noise for the predicted frequency. Hence, it is not surprising this impact event was not positively identified in the seismic data. Finding this crater was a lucky event as its formation this close to InSight has a probability of only ~0.2, and the odds of capturing it in before and after images are extremely low. We revisit impact-seismic discriminators in light of real experience with a seismometer on the Martian surface. Using measured noise of the instrument, we revise our previous prediction of seismic impact detections downward, from ~a few to tens, to just ~2 per Earth year, still with an order of magnitude uncertainty.

58 GEOSCIENCES↗
DOE OSTI2021

Unraveling the Combined Effects of Ice and Permafrost on Arctic Delta Morphodynamics

Arctic delta channel dynamics and coastal deposition patterns are uniquely affected by ice and permafrost. Building on previous experiments that documented the individual effects of ice and permafrost, we conducted new numerical modeling experiments to examine the combined influences of ice thickness and permafrost erodibility. We quantified the effects of ice, permafrost, and their combination on channel migration, avulsion, and persistence; onshore and nearshore deposition; and large-scale delta morphology. Both ice and permafrost limited channel mobility, enhanced overbank deposition, and roughened shorelines. Thick ice increased offshore channel incision and sediment delivery. Permafrost stabilized channels and focused deposition that led to enhanced channel levee development. The increased levee deposition limited widespread channel infilling leading to the preservation of abandoned channels. Levee growth was enhanced by overbank flooding due to ice-driven changes in water surface elevations. Furthermore, this flooding led to more extensive levees on deltas with both ice and permafrost compared to deltas with only permafrost. Ice also preserved small channels, leading to size-dependent channel persistence: the presence of ice in the channels reduced the fraction of time small channels remained active relative to larger channels. Resistant permafrost damped this effect, as smaller and more distal channels were active earlier in the flood on deltas with resistant permafrost likely due to the smaller number of channels on permafrost-laden deltas. Finally, our results suggest that the loss of both ice and permafrost will increase channel mobility, decrease overbank flooding, increase shoreline progradation rates, and decrease offshore transport of riverine fluxes.

58 GEOSCIENCES↗
DOE OSTI2020

Evaluation of Methods for Causal Discovery in Hydrometeorological Systems

Understanding causal relations is of utmost importance in hydrology and climate research for systems identification, prediction, and understanding systems behavior in a changing climate. Traditionally, researchers in hydrometeorology attempted to study causal questions by conducting controlled experiments using numerical models. This approach, however, in most cases of interest provides uncertain results because the models are approximate representation of the natural system. An alternative approach that has recently drawn significant attention in several fields is to infer causal relations from purely observational data. It possesses several traits to its utility particularly in hydrometeorology due to the rapid accumulation of in situ and remotely sensed data records. The first objective of this study is to present a brief description of four causal discovery methods (Granger causality, Transfer Entropy, graph-based algorithms, and Convergent Cross Mapping) with special emphasis on the assumptions on which they are built. Second, using synthetic data generated from a hydrological model, we assess their performance in retrieving causal information taking into account sensitivity to sample size and presence of noise. Last, we use causal analysis to examine and formulate hypotheses on causal drivers of evapotranspiration in a shrubland region during summer and winter seasons. An interpretation of the hypotheses based on canopy seasonal dynamics and evapotranspiration processes is presented. It is hoped that the results presented here can be useful in guiding researchers studying hydrometeorological systems as to which causal method is most appropriate to the characteristics of the system under study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Preferential Solute Transport in Low Permeability Zones During Spontaneous Imbibition in Heterogeneous Porous Media

Multiphase flows in porous media|and the associated solute transport processes|are controlled by a combination of gravity, capillary, and viscous forces. Geologic heterogeneity influences flow and transport processes, yet gaps exist in our understanding of how heterogeneity impacts capillary-driven transport, such as during spontaneous imbibition. Here we use positron emission tomography, combined with a newly developed method for conducting spontaneous imbibition experiments, to observe solute transport into low permeability regions during imbibition. Using an experimentally parameterized 2D numerical model, we demonstrate how capillary-driven ow near the imbibition front can carry solutes into low permeability regions. This process displaces solutes from high permeability zones while the cumulative amount of solute in low permeability zones increases- opposite of what is observed under fully saturated solute transport conditions. These results highlight the complex flow and solute transport behavior that arise during multiphase displacements in heterogeneous geologic porous media.

58 GEOSCIENCES↗
DOE OSTI2021

Tropical Cyclone Frequency

Abstract The frequency with which tropical cyclones (TCs) occur controls all other aspects of tropical cyclone risk since a storm that does not occur can do no harm. Yet this frequency is poorly understood. There is no accepted theory that explains the average number of TCs that occur each year on the Earth, nor how that number will change with global warming. Arguments based on global budgets of heat or moisture do not yet appear helpful, nor does a detailed understanding of the physical processes of TC genesis. Empirical indices that predict TC frequency as a function of large‐scale environmental variables can explain some of its relative variations in space and time, but not its absolute value. Global numerical models with horizontal grid spacings on the order of 25–50 km have allowed much improved simulations of TC activity, however. Many such models project a decrease in frequency with warming, but some project an increase. Idealized simulations, including those at higher resolutions, offer promise by allowing a systematic, deductive investigation of the roles of individual environmental factors. In addition to the larger‐scale environmental modulation of genesis likelihood, precursor disturbances, or “seeds”, may exert an independent influence on TC frequency.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The Global Fingerprint of Modern Ice-Mass Loss on 3-D Crustal Motion

Crustal motion generated by rapid ice-mass loss from Earth's glaciers and ice sheets has previously been considered in Global Navigational Satellite System (GNSS) analyses and numerical models across regions of ice retreat. However, the fingerprint of ice-mass loss is not limited to glaciated areas, but is characterized by a global pattern of 3-D crustal deformation. We compute “far-field” vertical and horizontal deformation rates that occurred in response to early 21st century mass flux from the Greenland and Antarctic Ice Sheets, global glaciers and ice caps, and associated ocean loading. We demonstrate that mass changes in the Greenland Ice Sheet and high latitude glacier systems each generated average crustal motion of 0.1–0.4 mm/yr across much of the Northern Hemisphere, with significant year-to-year variability in magnitude and direction. Horizontal motions associated with ice-mass loss exceed vertical rates in many far-field areas, and both should be considered in future analysis of GNSS measurements.

58 GEOSCIENCES↗
DOE OSTI2021

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗