Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow forming”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Coordinated time variability of multi-phase ultra-fast outflows in J132216.25 + 052446.3

ABSTRACT We present a time variability analysis of broad absorption lines (BAL; spread over the velocity range of 5800–29 000 km s−1) seen in the spectrum of J132216.25 + 052446.3 (zem = 2.04806) at ten different epochs spanning over 19 yr. The strongest absorption component (BAL-A; spread over 5800–9900 km s−1) is made up of several narrow components having velocity separations close to C iv doublet splitting. The C iv, N v, and Si iv absorption from BAL-A show correlated optical depth variability without major changes in the velocity structure. A very broad and shallow absorption (BAL-C; spread over the velocity range 15 000–29 000 km s−1) emerged during our monitoring period coinciding with a dimming episode of J1322 + 0524. All the identified absorption lines show correlated variability with the equivalent widths increasing with decreasing flux. This together with the C iv emission line variability is consistent with ionization being the main driver of the correlated variability. The observed UV-continuum variations are weaker than what is required by the photoionization models. This together with a scatter in the C iv equivalent width at a given continuum flux can be understood if variations of the C iv ionizing photons are much larger than that of the UV continuum, the variations in the ionizing photon and UV fluxes are not correlated and/or the covering factor of the flow varies continuously. We suggest BAL-A is produced by a stable clumpy outflow located beyond the broad emission line region and BAL-C is a newly formed wind component located near the accretion disc and both respond to changes in the ionizing continuum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Progress Report on Performance of A709 and G91 Steels in Sodium

Specimens in six different processing and heat treatment conditions of A709 H58776 were exposed to sodium at 550, 600, and 650°C, respectively for various exposure times. G91 base metal was tested in sodium at 550, 600, and 650°C, and G91 weldment at 550 and 600°C, respectively. Parallel thermal aging experiments were conducted on G91 and A709 to obtain thermal aging data for comparison with sodium exposure results to separate the thermal and sodium effects. The corrosion data obtained on A709 steel continue to show low corrosion rates, which indicate good compatibility of A709 with sodium when oxygen content is controlled. G91 has also shown acceptable corrosion rates over the temperature range of 550-650°C in sodium. Thermal aging or sodium exposures reduced the tensile strength, uniform elongation, and total elongation for A709 ESR, AOD specimens. These specimens consistently showed higher yield stress and tensile strength and lower uniform and total elongations after sodium exposure than after thermal aging under comparable testing conditions. It implies that there was an additional effect of sodium exposure that gave rise to an increase in tensile strength of A709. Dynamic strain aginginduced flow serrations were completely removed by either thermal aging or sodium exposure at 650°C, which imply that carbon or nitrogen in solution was removed from the solution and formed precipitates during aging or sodium exposures. The tensile data ruled out the possibility of decarburization under the sodium exposure conditions. The HOMO specimens behaved somewhat differently from the ESR and AOD specimens in sodium. While thermal aging and sodium exposures at 550 and 600°C decreased the yield stress and the ultimate tensile strength of G91, there was virtually no additional effect resulting from sodium exposures at these two temperatures. In contrast, sodium exposures at 650°C had a drastic effect on the yield stress and the ultimate tensile strength of G91. It is suggested that G91 experienced decarburization in the 650°C sodium environments. Carbon concentrations in sodium in the SMT-1 and SMT-2 loops were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8-1.2 ppm in the SMT-1 loop and 0.3-0.7 ppm in the SMT-2 loop. The carbon activity in sodium in the SMT-2 loop was estimated to be 0.03-0.08 at 600°C, and 0.08-0.2 at 550°C. The carbon activity in sodium at 650°C in the SMT-1 loop was 0.04-0.07. Equilibrium simulation of the carburization – decarburization processes was conducted for A709 and G91 steels exposed in sodium environments at temperatures of 550-700°C. The carbon activity-concentration relationship for G91 was re-evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC carbides dominated the process at low carbon activities. Formation of NbC increases the decarburization resistance of G91 steel. It remains to be understood whether the decarburization resistance provided by MC carbides in G91 can be maintained during long-term operations of SFRs. The carbon concentration-activity relationship for A709 was calculated based on the equilibrium of fcc-austenite and M 23 C 6 carbide phase. The calculations showed that A709 decarburizes at 650°C in the SMT-1 loop environments, which is different from the experimental findings. The carbon concentration-activity relationship in A709 was further evaluated by considering four phases, including fcc-austenite, M 23 C 6 , NbC, and TiC carbide phases. TiC is the most stable carbides among the three carbide phases. The carburization-decarburization process in A709 is dictated by M 23 C 6 carbides at high carbon activities and by TiC at low carbon activities. Formation of TiC increases the decarburization resistance of A709 steel. Because of the complexity of the precipitation process in A709 during sodium exposures, detailed characterization of precipitates in sodium-exposed specimens and improved thermodynamic models are needed to understand the carburization-decarburization behavior of A709. Kinetic analysis of carbon transfer will be investigated for G91 and A709 in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

LGAD-in-CMOS Design Studies

We describe a process flow for a CMOS sensor with intrinsic gain. The design is based on a base wafer which is prepared with closely spaced buried n and p doping layers that form an amplifying high field gain region when the device is biased. A deep n-well forms the collecting electrode. The buried n layer serves to moderate the internal field so that the CMOS transistors within the n-well remain within their operating region.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Filament Formation and Melt Spinning of Coal-Based Mesophase Pitch for Carbon Fiber Production

High-performance carbon fibers excel as high specific strength and modulus materials, and are utilized in composites applications ranging across aerospace, automotive, energy, infrastructure, and sports equipment sectors. Compared to PAN-based carbon fiber, mesophase pitch-based carbon fiber has lower tensile strength but provides a higher modulus, higher thermal conductivity, and a lower cost potential. However, challenges in stable, continuous melt spinning processing are a serious limiting factor. Unlike typical melt spinning of long, linear chain polymers, mesophase pitch is a comprised of relatively shorter polycyclic aromatic hydrocarbons which form a liquid crystal. In its ‘green’ or ‘as-spun’ state, the fibers are very fragile. Moreover, its temperature of processing is quite high, often approaching 400ºC. Complex flow dynamics combined with short length and time scale heat transfer render its melt spinning a formidable processing challenge. Improved understanding of the root causes for nascent filament breakage and spinning instabilities are needed. This work aims to determine fundamental phenomena that govern mesophase pitch filament formation by comparison of required draw force and uninterrupted spinning minutes with applied draw down ratio. The effect of perturbations on spinning stability are analyzed using capillary rheometry, microscopy of elongating filaments, and filtration.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗
DOE OSTI2020

Transcritical diffuse-interface hydrodynamics of propellants in high-pressure combustors of chemical propulsion systems

Rocket engines and high-power new generations of gas-turbine jet engines and diesel engines oftentimes involve the injection of one or more reactants at subcritical temperatures into combustor environments at high pressures, and more particularly at pressures higher than those corresponding to the critical points of the individual components of the mixture, which typically range from 13 to 50 bars for most propellants. This class of trajectories in the thermodynamic space has been traditionally referred to as transcritical. However, the fundamental understanding of fuel atomization, vaporization, mixing, and combustion processes at such high pressures remains elusive. In particular, whereas fuel sprays are relatively well characterized at normal pressures, analyses of dispersion of fuel in high-pressure combustors are hindered by the limited experimental diagnostics and theoretical formulations available. The description of the thermodynamics of hydrocarbon-fueled mixtures employed in chemical propulsion systems is complex and involves mixing-induced phenomena, including an elevation of the critical point whereby the coexistence region of the mixture extends up to pressures much larger than the critical pressures of the individual components. As a result, interfaces subject to surface-tension forces may persist in multicomponent systems despite the high pressures, and may give rise to unexpected spray-like atomization dynamics that are otherwise absent in monocomponent systems above their critical point. In this article, the current understanding of this phenomenon is reviewed within the context of propulsion systems fueled by heavy hydrocarbons. Emphasis is made on analytical descriptions at mesoscopic scales of interest for computational fluid dynamics. In particular, a set of modifications of the constitutive laws in the Navier–Stokes equations for multicomponent flows, supplemented with a high-pressure equation of state and appropriate redefinitions of the thermodynamic potentials, are introduced in this work based on an extended version of the diffuse-interface theory of van der Waals. The resulting formulation involves revisited forms of the stress tensor and transport fluxes of heat and species, and enables a description of the mesoscopic volumetric effects induced by transcritical interfaces consistently with the thermodynamic phase diagram of the mixture at high pressures. Additionally, applications of the theory are illustrated in canonical problems, including dodecane/nitrogen transcritical interfaces in non-isothermal systems. The results indicate that a transcritical interface is formed between the propellant streams that persists downstream of the injection orifice over distances of the same order as the characteristic thermal-entrance length of the fuel stream. The transcritical interface vanishes at an edge that gives rise to a fully supercritical mixing layer.

02 PETROLEUM↗
DOE OSTI2021

Impact of Primary and Secondary ZDDP and Ionic Liquid as Lubricant Oil Additives on the Performance and Physicochemical Properties of Pd-Based Three-Way Catalysts

In the present study, two industry primary and secondary zinc dialkyldithiophosphate standards, ZDDP1 and ZDDP2, respectively, are evaluated for their impact on the performance of Pd-based three-way catalyst and bench-marked against two mixed lubricant additives formed from either ZDDP1 or ZDDP2 with a second-generation oil-miscible phosphoric-containing ionic liquid (IL). The three-way catalysts (TWCs) are exposed to the lubricant additives in an engine bench under four different scenarios: a base case with no additive (NA), ZDDP1, IL+ZDDP1, ZDDP2, and IL+ZDDP2. The engine-aged TWC samples are characterized through a variety of analytical techniques, including evaluation of catalyst reactivity in a bench-flow reactor. With respect to the water–gas shift reaction and the oxygen storage capacity, the ZDDP2- and IL+ZDDP2-aged TWC samples are more degraded than the ZDDP1- and IL+ZDDP1-aged TWC samples. X-ray diffraction (XRD) patterns indicate that phosphorus in the form of CePO4 was found to be present in the washcoat of all TWC samples, with the highest amount found in the ZDDP2-aged TWC sample. The results obtained from XRD are further confirmed by those from inductively coupled plasma-optical emission spectroscopy (ICP-OES), which show that more phosphorus is detected in the washcoat of ZDDP2- and IL+ZDDP2-aged TWC samples than in the ZDDP1- and IL+ZDDP1-aged TWC samples.

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

Hopper flows of mixtures of spherical and rod‐like particles via the multisphere method

Abstract The discharge of granular media from hoppers plays a vital role in many solids handling operations. The dynamic flow of poly‐disperse nonspherical particles differs markedly from mono‐disperse spheres and highlights the need to rigorously account for particle shape effects. The multisphere discrete element method (MS‐DEM) was utilized here to examine the bulk packing fraction, discharge rate, collisional stress tensor, and order parameter of a bi‐disperse sphere‐rod mixture evacuating a hopper. A modified form of the Myers and Sellers correlation with variable shape factor is capable of replicating the MS‐DEM discharge data. The collisional stress tensor is observed to be highly symmetric for mono‐disperse spheres but becomes asymmetric when rod‐like particles are introduced. The order parameter within the hopper increases significantly during discharge and the rod‐like particles undergo shear induced alignment with their long dimensions orthogonal to the outflow plane.

Lattanzi, Aaron M.↗
DOE OSTI2020

Asphaltene-induced spontaneous emulsification: Effects of interfacial co-adsorption and viscoelasticity

Asphaltenes are a class of high molecular weight aromatic compounds found in crude oil. They adsorb onto toluene-water interfaces and induce a spontaneous emulsification phenomenon, whereby stable water-in-oil emulsions form without the need of an external energy input. This work aims to control and understand the factors affecting spontaneous droplet formation in the presence of asphaltene adsorption. This is particularly useful for crude oil refining, where the presence of a stable emulsion hampers the efficiency of downstream processing operations. We explore the effect of the addition of copolymers designed as crude oil flow improvers as a means to control the extent of emulsion formation. We find that the polymers competitively adsorb onto the toluene-water interface and diminish spontaneous emulsification. Additionally, we conduct fluorescence microscopy experiments and measurements of the interfacial energy to determine the mechanism of spontaneous emulsification in asphaltene systems. We conclude that an emulsion forms via the diffusion of molecular water into the oil phase and subsequent binding with asphaltene aggregates, leading to the nucleation of micrometer-sized water droplets. We find that the polymer forms complexes with the dissolved asphaltenes, possibly hampering the ability of diffused water to bind to the asphaltenes and reducing the extent of spontaneous emulsification. Finally, we investigate the role of interfacial shear and dilatational viscoelasticity to better understand which fundamental interfacial properties are important in the emulsification of asphaltene-laden systems. We find that the rate of formation of an interfacial microstructural network is inversely correlated with the extent and rate of spontaneous emulsification.

42 ENGINEERING↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2022

The heat transfer coefficient associated with a moving packed bed of silica particles flowing through parallel plates

Concentrating Solar Power (CSP) with thermal energy storage has the potential to be a renewable energy technology with long duration, inexpensive energy storage. Higher temperature operation increases the conversion efficiency and reduces the cost of energy storage. Several emerging CSP designs utilize particles as the solar receiver due to their high temperature stability and low cost. In some designs, the particles are also used as the thermal energy storage (TES) media. In either configuration, an energy transfer is required between the hot particles and the working fluid in the power cycle in a Particle-to-Fluid Heat Exchanger (PtFHX). Understanding the heat transfer between a moving packed bed and a stationary surface is critical to the successful design of a PtFHX. In this paper, a test facility is described in which a moving packed bed of silica sand with particle size 100–600 μm is introduced into the channel formed by two parallel plates, one of which is heated. The effective static thermal conductivity of the particles used for the test are separately measured over the entire range of test temperatures. The inlet and outlet bulk temperatures of the particle flow are measured as are the surface temperatures at several axial locations along the centerline of the plate. The result is the measurement of heat transfer coefficient as a function of temperature for several velocities. The uncertainty of the measurements is presented and the results are compared to model results found in the literature.

14 SOLAR ENERGY↗
DOE OSTI2021

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Three-Dimensional Free Gas Flow Focuses Basin-Wide Microbial Methane to Concentrated Methane Hydrate Reservoirs in Geological System

Here, we present a systematic model that links the generation, migration, phase partitioning, and accumulation of methane into a closed loop as the sediment is deposited from the seafloor and buried through the base of hydrate stability zone (BHSZ). In our model, methane is generated by biodegradation of organic carbon in muds. Hydrate does not form and methane is not trapped until a coarse-grained layer is deposited, because the small pores prevent hydrate formation in muds. Instead, methane diffuses into sands/silts where methane solidifies into hydrate. As hydrate-bearing sands/silts pass through the BHSZ during sediment burial, methane hydrate dissociates, and releases free gas. The released and the newly generated free gas below the BHSZ concentrates into a vertical/dipping zone with low capillary entry pressure and high permeability and flows upward driven its buoyancy. When free gas reaches the hydrate stability zone (HSZ), capillary forces drive free gas to flow laterally, preferentially enter sands/silts, feed hydrate growth, and elevate hydrate saturation. With three-dimensional focused free gas flow, microbial methane that is generated from a much larger fetch area of the entire basin, both above and below the BHSZ, is concentrated into coarse-grained layers at structural closures for hydrate formation. Our model illustrates how geological evolution, microbial methane generation, and gas flow by buoyancy couple to generate concentrated hydrate deposits in geological system. These insights can be used to explore for high-concentration methane hydrate and are important for understanding the methane budget and carbon cycle under the seafloor.

58 GEOSCIENCES↗
DOE OSTI2022

Shock Formation and Vorticity Creation for 3d Euler

Abstract We analyze the shock formation process for the 3D nonisentropic Euler equations with the ideal gas law, in which sound waves interact with entropy waves to produce vorticity. Building on our theory for isentropic flows in [3, 4], we give a constructive proof of shock formation from smooth initial data. Specifically, we prove that there exist smooth solutions to the nonisentropic Euler equations which form a generic stable shock with explicitly computable blowup time, location, and direction. This is achieved by establishing the asymptotic stability of a generic shock profile in modulated self‐similar variables, controlling the interaction of wave families via: (i) pointwise bounds along Lagrangian trajectories, (ii) geometric vorticity structure, and (iii) high‐order energy estimates in Sobolev spaces. © 2022 Wiley Periodicals LLC.

Mathematics↗
DOE OSTI2020

Shock initiation of the HMX-based explosive PBX 9012: Experiments, uncertainty analysis, and unreacted equation-of-state

In this study, shock initiation experiments have been carried out on the polymer-bonded explosive PBX 9012 [nominally 90% octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazine (HMX), 10% vinylidene–hexafluoropropylene copolymer (Viton A) by weight] in order to provide calibration data for the explosive’s unreacted equation-of-state (EOS) and shock-to-detonation transition for reactive burn rate calibration. The input pressures covered the range of 1.86–4.43 GPa. This provided run-to-detonation depths ranging from > 22.3 mm for the lowest pressure shot to 4.91 mm at the highest pressure. The relative shock sensitivity of PBX 9012 is compared to other HMX-explosives in terms of the Pop-plot, showing that the studied explosive is more sensitive than other similar HMX-based counterparts (with notable exceptions). The uncertainty in the shock velocity determinations from the shock tracker measurements are also investigated, yielding new uncertainty measures in the generated Hugoniot data and run-to-detonation coordinates. Finally, the unreacted equation-of-state is determined using a linear U s – u p Mie–Grüneisen relation and the Davis reactants EOS analytical form, the latter being more suited for reactive burn modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modelling and simulation of brinicle formation

Below the Arctic sea ice, under the right conditions, a flux of icy brine flows down into the sea. The icy brine has a much lower fusion point and is denser than normal seawater. As a result, it sinks while freezing everything around it, forming an ice channel called a brinicle (also known as ice stalactite). In this paper, we develop a mathematical model for this phenomenon, assuming cylindrical symmetry. The fluid is considered to be viscous and quasi-stationary. The heat and salt transport are weakly coupled to the fluid motion and are modelled with the corresponding conservation equations, accounting for diffusive and convective effects. Finite-element discretization is employed to solve the coupled system of partial differential equations. We find that the model can capture the general behaviour of the physical system and generate brinicle-like structures while also recovering dendrite composition, which is a physically expected feature aligned with previous experimental results. This represents, to our knowledge, the first complete model proposed that captures the global structure of the physical phenomenon even though it has some discrepancies, such as brine accumulation.

97 MATHEMATICS AND COMPUTING↗