SEARCH · Search NASA

Results for “Process settings”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Relaxation time estimation in surface NMR

NMR relaxation time estimation methods and corresponding apparatus generate two or more alternating current transmit pulses with arbitrary amplitudes, time delays, and relative phases; apply a surface NMR acquisition scheme in which initial preparatory pulses, the properties of which may be fixed across a set of multiple acquisition sequence, are transmitted at the start of each acquisition sequence and are followed by one or more depth sensitive pulses, the pulse moments of which are varied across the set of multiple acquisition sequences; and apply processing techniques in which recorded NMR response data are used to estimate NMR properties and the relaxation times T.sub.1 and T.sub.2* as a function of position as well as one-dimensional and two-dimension distributions of T.sub.1 versus T.sub.2* as a function of subsurface position.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Evaluating the Dynamic Behavior of Information Technology Systems in Healthcare using Markov Simulation

In recent years, the safety and reliability of information technology (IT) systems in the healthcare industry are of increasing importance. In this paper, we propose an approach for monitoring and predicting reliability degradation in Health IT (HIT) using Markov chain (MC). A MC model provides an opportunity to represent highly dynamic systems, such as HIT, in a succinct manner to simulate the evolution of the system over time in discrete time steps. The model can also represent system behavior that varies over along duration. Consequently, using electronic health records (EHR) data from systems such as the Veterans Affairs’ Corporate Data Warehouse systems, we defined clinical workflow as a Transaction Process Model (TPM). The TPM represents a set of states in the Consult workflow. It is also an ideal workflow description and has several degrees of freedom. The TPM is then converted into a MC representation and the EHR data is used to compute transition probabilities between the nodes in the MC. The original MC representation is perturbed by changing the transition probabilities to simulate alternative system workflow paths and identifying scenarios that could impact system reliability. We present scenarios that illustrate the proposed approach and discuss some of the insights from the results.

Omitaomu, Olufemi↗

DOE OSTI2023

Configuration of an optical switch fabric using machine learning

An optical switch fabric comprises two or more optical switch elements. The optical switch elements are configured in a topology. A switch control has a plurality of bias control signals. The switch control can address one or more of the optical switch elements and can apply one of the bias control signals to bias of the addressed optical switch element to establish a switch setting. The topology and switch settings determine how each of one of the inputs is connected to each of one of the outputs of the optical switch fabric. The switch settings are determined by a machine learning process which includes a model creation. The model can be made to adapt dynamically during optical switch fabric operation.

Dupuis, Nicolas↗

DOE OSTI2026

Universal Jaynes-Cummings Control of an Oscillator

The Jaynes-Cummings (JC) interaction-the coherent exchange of excitations between a two-level system and a harmonic oscillator-is one of the fundamental interactions of quantum optics, realized across platforms such as cavity quantum electrodynamics, trapped ions, mechanical resonators, and superconducting circuits. Although JC interactions and qubit rotations form a universal gate set for oscillator control, practical implementations have not been demonstrated. Here we develop and experimentally demonstrate universal JC-based oscillator control by compiling arbitrary unitary gates into sequences of JC interactions and qubit rotations. In our experiment, the oscillator is realized using a mode of a high quality factor microwave cavity and the ancilla qubit using a superconducting transmon circuit, with the JC interaction implemented by a sideband interaction enabled by the Josephson nonlinearity. The native gates are constructed to be closed below a chosen cutoff photon number, encoding a qudit with suppressed leakage errors, while ancilla relaxation errors are detectable. We further find that the dispersive shift serves as a compilation resource that reduces circuit depths. We demonstrate universal qudit control and implement a single-qutrit gate set with a mean post-selected process fidelity of 96%, as well as ququart and ququint shift gates. These results establish Jaynes-Cummings control as a practical route to universal oscillator control, enabling programmable bosonic processors across a variety of quantum platforms.

Huang, Jordan [Rutgers U., Piscataway]↗

DOE OSTI2018

Affordable and accurate large-scale hybrid-functional calculations on GPU-accelerated supercomputers

Performing high accuracy hybrid functional calculations for condensed matter systems containing a large number of atoms is at present computationally very demanding or even out of reach if high quality basis sets are used. We present a highly optimized multiple graphics processing unit implementation of the exact exchange operator which allows one to perform fast hybrid functional density-functional theory (DFT) calculations with systematic basis sets without additional approximations for up to a thousand atoms. With this method hybrid DFT calculations of high quality become accessible on state-of-the-art supercomputers within a time-to-solution that is of the same order of magnitude as traditional semilocal-GGA functionals. The method is implemented in a portable open-source library.

density functional theory↗

DOE OSTI2020

Predicting Excitation Energies of Twisted Intramolecular Charge-Transfer States with the Time-Dependent Density Functional Theory: Comparison with Experimental Measurements in the Gas Phase and Solvents Ranging from Hexanes to Acetonitrile

Electronically excited states characterized by intramolecular charge transfer play an essential role in many biological processes and optical devices. The ability to make quantitative ab initio predictions of the relative energetics involved is a challenging yet desirable goal, especially for large molecules in solution. In this work, we present a data set of 61 experimental measurements of absorption and emission processes, both in the gas phase and in solvents representing a broad range of polarities, which involve intramolecular charge transfer mediated by a nonzero, "twisted" dihedral angle between one or more donor and acceptor subunits. Among a variety of density functionals investigated within the framework of linear-response theory, the "optimally tuned" LRC-ωPBE functional, which utilizes a system-specific yet nonempirical procedure to specify the range-separation parameter, emerges as the preferred choice. For the entire set of excitation energies, involving changes in dipole moment ranging from 4 to >20 Debye, the mean signed and absolute errors are 0.02 and 0.18 eV, respectively (compared, e.g., to -0.30 and 0.30 for PBE0, 0.44 and 0.47 for LRC-ωPBEh, 0.83 and 0.83 for ωB97X-V). We analyze the performance of polarizable continuum solvation models available in Q-Chem that partition the solvent response into fast and slow time scales, and clear trends emerge when measurements corresponding to the four small 4-(dimethylamino)benzonitrile (DMABN)-like molecules and a charged species are excluded. We make the case that the large errors found only for small molecules in the gas phase and weak solvents cannot be expected to improve via the optimal tuning procedure, which enforces a condition that is exact only in the well-separated donor-acceptor limit, and present empirical evidence implicating the outsized importance for small donor-acceptor systems of relaxation effects that cannot be accounted for by the linear-response time-dependent density functional theory within the adiabatic approximation. Finally, we demonstrate the utility of the optimally tuned density functional approach by targeting the charge-transfer states of a large biomimetic model system for light-harvesting structures in Photosystem II.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Performance assessment for climate intervention (PACI): preliminary application to a stratospheric aerosol injection scenario

As the prospect of exceeding global temperature targets set forth in the Paris Agreement becomes more likely, methods of climate intervention are increasingly being explored. With this increased interest there is a need for an assessment process to understand the range of impacts across different scenarios against a set of performance goals in order to support policy decisions. The methodology and tools developed for Performance Assessment (PA) for nuclear waste repositories shares many similarities with the needs and requirements for a framework for climate intervention. Using PA, we outline and test an evaluation framework for climate intervention, called Performance Assessment for Climate Intervention (PACI) with a focus on Stratospheric Aerosol Injection (SAI). We define a set of key technical components for the example PACI framework which include identifying performance goals, the extent of the system, and identifying which features, events, and processes are relevant and impactful to calculating model output for the system given the performance goals. Having identified a set of performance goals, the performance of the system, including uncertainty, can then be evaluated against these goals. Using the Geoengineering Large Ensemble (GLENS) scenario, we develop a set of performance goals for monthly temperature, precipitation, drought index, soil water, solar flux, and surface runoff. The assessment assumes that targets may be framed in the context of risk-risk via a risk ratio, or the ratio of the risk of exceeding the performance goal for the SAI scenario against the risk of exceeding the performance goal for the emissions scenario. From regional responses, across multiple climate variables, it is then possible to assess which pathway carries lower risk relative to the goals. The assessment is not comprehensive but rather a demonstration of the evaluation of an SAI scenario. Future work is needed to develop a more complete assessment that would provide additional simulations to cover parametric and aleatory uncertainty and enable a deeper understanding of impacts, informed scenario selection, and allow further refinements to the approach.

climate intervention↗

DOE OSTI2022

Advanced stationary and nonstationary kernel designs for domain-aware Gaussian processes

Gaussian process regression is a widely-applied method for function approximation and uncertainty quantification. The technique has gained popularity recently in the machine learning community due to its robustness and interpretability. The mathematical methods we discuss in this paper are an extension of the Gaussian-process framework. We are proposing advanced kernel designs that only allow for functions with certain desirable characteristics to be elements of the reproducing kernel Hilbert space (RKHS) that underlies all kernel methods and serves as the sample space for Gaussian process regression. These desirable characteristics reflect the underlying physics; two obvious examples are symmetry and periodicity constraints. In addition, non-stationary kernel designs can be defined in the same framework to yield flexible multi-task Gaussian processes. We will show the impact of advanced kernel designs on Gaussian processes using several synthetic and two scientific data sets. The results of our research show that including domain knowledge, communicated through advanced kernel designs, has a significant impact on the accuracy and relevance of the function approximation.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Molten Salt vs. Liquid Sodium Receiver Selection Using the Analytic Hierarchy Process

An international team led by NREL analyzed the favorability of two alternative liquid receiver designs for a 700+ degrees C receiver under the Gen3 CSP Liquid Pathway project. The competing liquid heat transfer fluids were a ternary chloride salt and liquid-metal sodium. The team applied a facilitated analytic hierarchy process (AHP) to arrive at a recommended alternative and set a path forward for the project. The AHP criteria were formulated, weighted, and scored by the project leadership team and technical advisory committee consisting of energy industry and CSP experts. The six- month process culminated with a two-day workshop where the sodium alternative was deemed to have both a significantly higher benefit (19.3%) and a lower LCOE (11.4%), with only a slightly higher risk (~3%) than the salt alternative. Consequently, a sodium-receiver design was selected for the Liquid Pathway project, where it will be used to charge a two- tank chloride salt thermal energy storage system.

analytic hierarchy process↗

DOE OSTI2025

First Search for Neutral Current Coherent Single-Photon Production in MicroBooNE

This article presents the first search for neutrino-induced neutral current coherent single-photon production (NC coherent 1$\gamma$). The search makes use of data from the MicroBooNE 85-tonne active volume liquid argon time projection chamber detector, situated in the Fermilab Booster Neutrino Beam (BNB), with an average neutrino energy of $\langle E_{\nu}\rangle \sim 0.8$ GeV. A targeted selection of candidate neutrino interactions with a single photon-like electromagnetic shower in the final state and no visible vertex activity was developed to search for the NC coherent 1$\gamma$ process, along with two auxiliary selections used to constrain the dominant background from NC$\pi^0$ production. With an integrated exposure of $6.87 \times 10^{20}$ protons on target delivered by the BNB, we set the world's first limit for this rare process, corresponding to an upper limit on the flux-averaged cross section of $\sigma<1.49 \times 10^{-41}\text{cm}^2$ at 90% C.L.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Differentiable, Learnable, Regionalized Process-Based Models With Multiphysical Outputs can Approach State-Of-The-Art Hydrologic Prediction Accuracy

Predictions of hydrologic variables across the entire water cycle have significant value for water resources management as well as downstream applications such as ecosystem and water quality modeling. Recently, purely data-driven deep learning models like long short-term memory (LSTM) showed seemingly insurmountable performance in modeling rainfall runoff and other geoscientific variables, yet they cannot predict untrained physical variables and remain challenging to interpret. Here, we show that differentiable, learnable, process-based models (called δ models here) can approach the performance level of LSTM for the intensively observed variable (streamflow) with regionalized parameterization. We use a simple hydrologic model HBV as the backbone and use embedded neural networks, which can only be trained in a differentiable programming framework, to parameterize, enhance, or replace the process-based model's modules. Without using an ensemble or post-processor, δ models can obtain a median Nash-Sutcliffe efficiency of 0.732 for 671 basins across the USA for the Daymet forcing data set, compared to 0.748 from a state-of-the-art LSTM model with the same setup. For another forcing data set, the difference is even smaller: 0.715 versus 0.722. Meanwhile, the resulting learnable process-based models can output a full set of untrained variables, for example, soil and groundwater storage, snowpack, evapotranspiration, and baseflow, and can later be constrained by their observations. Both simulated evapotranspiration and fraction of discharge from baseflow agreed decently with alternative estimates. The general framework can work with models with various process complexity and opens up the path for learning physics from big data.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Employing a Hardware-in-the-Loop Approach to Realize a Fully Homomorphic Controller for a Small Modular Advanced High Temperature Reactor

This paper addresses the cybersecurity challenges of advanced nuclear reactors by integrating fully homomorphic encryption (FHE) into their control systems, enabling encrypted processing of control signals without compromising functionality. Advanced nuclear reactors, including Small Modular Reactors (SMRs) and microreactors, aim to achieve autonomous and remote operations, reducing costs and enhancing competitiveness. However, these advancements expand the attack surface for cyberattacks, particularly in autonomous and remote operation scenarios. Cyberattacks can exploit vulnerabilities to manipulate physical processes, causing shutdowns, asset damage, or public harm. Such attacks begin with passive reconnaissance, where adversaries intercept communications or observe behaviors to gather information, which is then leveraged to execute cyber-physical attacks by injecting malicious commands. Nuclear power must adopt cybersecurity protection measures to secure the integrity and availability of their digital control systems. This paper demonstrates the application of FHE to secure operations by enabling encrypted processing of sensitive signals and parameters -- ensuring privacy without exposing data. FHE supports secure mathematical operations on encrypted data without requiring decryption. Using a hardware-in-the-loop (HIL) approach, this paper implements an FHE-integrated controller on a BeagleBone Black (BBB) controlling a simulation of the Small Modular Advanced High Temperature Reactor (SmAHTR). By doing so, the encrypted controller protects the integrity of critical set points and control signals during transmission and processing. Thus, FHE-integrated controllers enhance secure operations of advanced nuclear reactors while maintaining functionality.

control systems↗

DOE OSTI2017

Dynamic mesh adaptation for front evolution using discontinuous Galerkin based weighted condition number relaxation

A new mesh smoothing method designed to cluster cells near a dynamically evolving interface is presented. The method is based on weighted condition number mesh relaxation with the weight function computed from a level set representation of the interface. The weight function is expressed as a Taylor series based discontinuous Galerkin projection, which makes the computation of the derivatives of the weight function needed during the condition number optimization process a trivial matter. For cases when a level set is not available, a fast method for generating a low-order level set from discrete cell-centered fields, such as a volume fraction or index function, is provided. Results show that the low-order level set works equally well as the actual level set for mesh smoothing. Meshes generated for a number of interface geometries are presented, including cases with multiple level sets. Lastly, dynamic cases with moving interfaces show the new method is capable of maintaining a desired resolution near the interface with an acceptable number of relaxation iterations per time step, which demonstrates the method's potential to be used as a mesh relaxer for arbitrary Lagrangian Eulerian (ALE) methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

JAX-CanVeg: A Differentiable Land Surface Model

Land surface models consider the exchange of water, energy, and carbon along the soil-canopy-atmosphere continuum, which is challenging to model due to their complex interdependency and associated challenges in representing and parameterizing them. Differentiable modeling provides a new opportunity to capture these complex interactions by seamlessly hybridizing process-based models with deep neural networks (DNNs), benefiting both worlds, that is, the physical interpretation of process-based models and the learning power of DNNs. Here, we developed a differentiable land model, JAX-CanVeg. The new model builds on the legacy CanVeg by incorporating advanced functionalities through JAX in the graphic processing unit support, automatic differentiation, and integration with DNNs. We demonstrated JAX-CanVeg's hybrid modeling capability by applying the model at four flux tower sites with varying aridity. To this end, we developed a hybrid version of the Ball-Berry equation that emulates the water stress impact on stomatal closure to explore the capability of the hybrid model in (a) improving the simulations of latent heat fluxes (LE) and net ecosystem exchange (NEE), (b) improving the optimization trade-off when learning observations of both LE and NEE, and (c) benefiting a multi-layer canopy model setup. Our results show that the proposed hybrid model improved the simulations of LE and NEE at all sites, with an improved optimization trade-off over the process-based model. Additionally, the multi-layer canopy set benefited hybrid modeling at some sites. Anchored in differentiable modeling, our study provides a new avenue for modeling land-atmosphere interactions by leveraging the benefits of both data-driven learning and process-based modeling.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Molten Salt vs. Liquid Sodium Receiver Selection Using the Analytic Hierarchy Process: Preprint

An international team led by NREL analyzed the favorability of two alternative liquid receiver designs for a 700+ °C receiver under the Gen3 CSP Liquid Pathway project. The competing liquid heat transfer fluids were a ternary chloride salt and liquid-metal sodium. The team applied a facilitated analytic hierarchy process (AHP) to arrive at a recommended alternative and set a path forward for the project. The AHP criteria were formulated, weighted, and scored by the project leadership team and technical advisory committee consisting of energy industry and CSP experts. The six-month process culminated with a two-day workshop where the sodium alternative was deemed to have both a significantly higher benefit (19.3%) and a lower LCOE (11.4%), with only a slightly higher risk (~3%) than the salt alternative. Consequently, a sodium-receiver design was selected for the Liquid Pathway project, where it will be used to charge a two-tank chloride salt thermal energy storage system.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2016

Compositional simulation of hydrocarbon recovery from oil shale reservoirs with diverse initial saturations of fluid phases by various thermal processes

We have studied the hydrocarbon production from oil shale reservoirs filled with diverse initial saturations of fluid phases by implementing numerical simulations of various thermal in-situ upgrading processes. We use our in-house fully functional, fully implicit, and non-isothermal simulator, which describes the in-situ upgrading processes and hydrocarbon recovery by multiphase-multicomponent systems. We have conducted two sets of simulation cases—five-spot well pattern problems and Shell In-situ Conversion Process (ICP) problems. In the five-spot well pattern problems, we have analyzed the effects of initial fluid phase that fills the single-phase reservoir and thermal processes by four cases—electrical heating in the single-phase-aqueous reservoir, electrical heating in the single-phase-gaseous reservoir, hot water injection in the single-phase-aqueous reservoir, and hot CO 2 -injection in the single-phase-gaseous reservoir. In the ICP problems, we have analyzed the effects of initial saturations of fluid phases that fill two-phase-aqueous-and-gaseous reservoir by three cases—initial aqueous phase saturations of 0.16, 0.44, and 0.72. Through the simulation cases, system response and production behavior including temperature profile, kerogen fraction profile, evolution of effective porosity and absolute permeability, phase production, and product selectivity are analyzed. In the five-spot well pattern problems, it is found that the hot water injection in the aqueous phase reservoir shows the highest total hydrocarbon production, but also shows the highest water-oil-mass-ratio. Productions of phases and components show very different behavior in the cases of electrical heating in the aqueous phase reservoir and the gaseous phase reservoir. In the ICP problems, it is found that the speed of kerogen decomposition is almost identical in the cases, but the production behavior of phases and components is very different. It is found that more liquid organic phase has been produced in the case with the higher initial saturation of aqueous phase by the less production of gaseous phase.

04 OIL SHALES AND TAR SANDS↗