Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differentiable predictive control”

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 181 records · Page 10

Semitauonic b-hadron decays: A lepton flavor universality laboratory

The study of lepton flavor universality violation (LFUV) in semitauonic b-hadron decays has become increasingly important in light of long-standing anomalies in their measured branching fractions, and the large datasets anticipated from the LHC experiments and Belle II. Here, in this review, a comprehensive survey of the experimental environments and methodologies for semitauonic LFUV measurements at the B factories and LHCb is undertaken, along with an overview of the theoretical foundations and predictions for a wide range of semileptonic decay observables. The future prospects of controlling systematic uncertainties down to the percent level, matching the precision of standard model (SM) predictions, are examined. Furthermore, new perspectives and caveats on combinations of the LFUV data are discussed and the world averages for the R(D(*)) ratios are revisited. Here it is demonstrated that different treatments for the correlations of uncertainties from D** excited states can vary the current 3σ tension with the SM within a 1σ range. Prior experimental overestimates of D**τν contributions may further exacerbate this. The precision of future measurements is also estimated; their power to exploit full differential information, and solutions to the inherent difficulties in self-consistent new physics interpretations of LFUV observables, are explored.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A provisional regulatory gene network for specification of endomesoderm in the sea urchin embryo

We present the current form of a provisional DNA sequence-based regulatory gene network that explains in outline how endomesodermal specification in the sea urchin embryo is controlled. The model of the network is in a continuous process of revision and growth as new genes are added and new experimental results become available; see http://www.its.caltech.edu/~mirsky/endomeso.htm (End-mes Gene Network Update) for the latest version. The network contains over 40 genes at present, many newly uncovered in the course of this work, and most encoding DNA-binding transcriptional regulatory factors. The architecture of the network was approached initially by construction of a logic model that integrated the extensive experimental evidence now available on endomesoderm specification. The internal linkages between genes in the network have been determined functionally, by measurement of the effects of regulatory perturbations on the expression of all relevant genes in the network. Five kinds of perturbation have been applied: (1) use of morpholino antisense oligonucleotides targeted to many of the key regulatory genes in the network; (2) transformation of other regulatory factors into dominant repressors by construction of Engrailed repressor domain fusions; (3) ectopic expression of given regulatory factors, from genetic expression constructs and from injected mRNAs; (4) blockade of the beta-catenin/Tcf pathway by introduction of mRNA encoding the intracellular domain of cadherin; and (5) blockade of the Notch signaling pathway by introduction of mRNA encoding the extracellular domain of the Notch receptor. The network model predicts the cis-regulatory inputs that link each gene into the network. Therefore, its architecture is testable by cis-regulatory analysis. Strongylocentrotus purpuratus and Lytechinus variegatus genomic BAC recombinants that include a large number of the genes in the network have been sequenced and annotated. Tests of the cis-regulatory predictions of the model are greatly facilitated by interspecific computational sequence comparison, which affords a rapid identification of likely cis-regulatory elements in advance of experimental analysis. The network specifies genomically encoded regulatory processes between early cleavage and gastrula stages. These control the specification of the micromere lineage and of the initial veg(2) endomesodermal domain; the blastula-stage separation of the central veg(2) mesodermal domain (i.e., the secondary mesenchyme progenitor field) from the peripheral veg(2) endodermal domain; the stabilization of specification state within these domains; and activation of some downstream differentiation genes. Each of the temporal-spatial phases of specification is represented in a subelement of the network model, that treats regulatory events within the relevant embryonic nuclei at particular stages. (c) 2002 Elsevier Science (USA).

Non-NASA Center↗

Model studies of crosswind landing-gear configurations for STOL aircraft

A dynamic model was used to directly compare four different crosswind landing gear mechanisms. The model was landed as a free body onto a laterally sloping runway used to simulate a crosswind side force. A radio control system was used for steering to oppose the side force as the model rolled to a stop. The configuration in which the landing gears are alined by the pilot and locked in the direction of motion prior to touchdown gave the smoothest runout behavior with the vehicle maintaining its crab angle throughout the landing roll. Nose wheel steering was confirmed to be better than steering with nose and main gears differentially or together. Testing is continuing to obtain quantitative data to establish an experimental data base for validation of an analytical program that will be capable of predicting full scale results.

Stubbs, S. M.↗

RNAseq analysis of Cellvibrio japonicus during starch utilization differentiates between genes encoding carbohydrate active enzymes controlled by substrate detection or growth rate

ABSTRACT Bacterial utilization of starch is increasingly of interest as the importance and contributions of animal gut microbiomes become more defined. Consequently, identifying and characterizing the bacterial enzymes responsible for the degradation, transport, and metabolism of starch will enable developments in pharmaceutical, biotechnological, and culinary industries searching for novel prebiotics, carrier molecules, and low glycemic index sweeteners. The current challenge is that bacteria proficient at starch utilization often have hundreds of carbohydrate active enzymes, and it is unclear which are essential for starch utilization using only homology-based bioinformatics or computational methods. Complementary experimental data are also needed, especially to understand the regulation of bacterial starch utilization. We have completed an RNAseq analysis of the Gram-negative bacterium Cellvibrio japonicus and found that it has sophisticated regulation that includes substrate sensing and growth rate components for genes that encode starch-degrading enzymes. Among the 22 genes predicted to encode starch-active enzymes, C. japonicus has 10 alpha-amylases, 4 alpha-glucosidases, 2 pullulnases, and 2 cyclomaltodextrin glucanotransferases, 15 of which were up-regulated during exponential growth on starch and 8 up-regulated in stationary phase. Growth analyses with an enzyme secretion deficient mutant of C. japonicus suggested that secreted amylases are essential for this bacterium to degrade starch. Our approach of coupling a physiological growth assay with transcriptomic data provides a platform to identify targets for further genetic or biochemical analysis that can be broadly applied to other starch-utilizing bacteria. IMPORTANCE Understanding the bacterial metabolism of starch is important as this polysaccharide is a ubiquitous ingredient in foods, supplements, and medicines, all of which influence gut microbiome composition and health. Our RNAseq and growth data set provides a valuable resource to those who want to better understand the regulation of starch utilization in Gram-negative bacteria. These data are also useful as they provide an example of how to approach studying a starch-utilizing bacterium that has many putative amylases by coupling transcriptomic data with growth assays to overcome the potential challenges of functional redundancy. The RNAseq data can also be used as a part of larger meta-analyses to compare how C. japonicus regulates carbohydrate active enzymes, or how this bacterium compares to gut microbiome constituents in terms of starch utilization potential.

59 BASIC BIOLOGICAL SCIENCES↗

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗

A search for the prewetting line

This paper describes efforts to locate the prewetting line in a binary liquid system (isopropanol-perfluoromethylcyclohexane) at the vapor-liquid interface. Tight upper bounds were placed on the temperature separation (0.2 K) between the prewetting line and the line of bulk liquid phase separation. The prewetting line in systems at equilibrium was not detected. Experimental signatures indicative of the prewetting line occurred only in nonequilibrium situations. Several theories predict that the adsorption of one of the components (the fluorocarbon, in this case) at the liquid-vapor interface should increase abruptly, at a temperature sightly above the temperature at which the mixture separates into two liquid phases. A regular solution calculation indicates that this prewetting line should have been easily detectable with the instruments used in this experiment. Significant features of the experiment are: (1) low-gradient thermostatting, (2) in situ stirring, (3) precision ellipsometry from the vapor-liquid interface, (4) high resolution differential index of refraction measurements using a novel cell design, and (5) computer control.

Schmidt, J. W.↗

Effects of Ionizing Radiation on Murine Gene Expression in Skin and Bone

Long duration spaceflight causes a negative calcium balance and reduces bone density in astronauts. The potential for exposure to space radiation to contribute to lasting decrements in bone mass is not yet understood. Sustained changes to bone mass have a relatively long latency for development, however skin is a radiation sensitive organ and changes in skin gene expression may serve as an early radiation biomarker of exposures and may correlate with adverse effects on skeletal tissue. Previous studies have shown that FGF18 gene expression levels of hair follicles collected from astronauts on the ISS rose over time. In the hair follicle, FGF18 signaling mediates radioresistance in the telogen by arresting the cell cycle, and FGF18 has the potential to function as a radioprotector. In bone, FGF18 appears to regulate cell proliferation and differentiation positively during osteogenesis and negatively during chondrogenesis. Cellular defense responses to radiation are shared by a variety of organs, hence in this study, we examined whether radiation induced gene expression changes in skin may be predictive of the responses of skeletal tissue to radiation exposure. We have examined oxidative stress and growth arrest pathways in mouse skin and long bones by measuring gene expression levels via quantitative polymerase chain reaction (qPCR) after exposure to total body irradiation (TBI). To investigate the effects of irradiation on gene expression, we used skin and femora (cortical shaft) from the following treatment groups: control (normally loaded, sham-irradiated), and TBI (0.5 Gy Fe-56 600 MeV/n and 0.5 Gy H-1 150 MeV/n). Animals were euthanized one and 11 days post-IR. Statistical analysis was performed via a Student's ttest. In skin samples one day after IR, skin expression of FGF18 was significantly greater (3.8X) than sham-irradiated controls (3.8X), but did not differ 11 days post TBI. Expression levels of other radiation related genes (Nfe2l2, Trp53, Cdkn1a, FoxO3, Gadd45g, SOD1), was not different due to TBI at either time point. In bone (femora) TBI significantly increased (3.8X) expression of the pro-bone resorption cytokine, MCP-1, one day after TBI. FGF18 expression in skin and MCP- 1 expression in bone were found to be positively correlated (P less than 0.002, r=0.8779). Further, microcomputed tomography analysis of tibae from these animals showed reduced fractional cancellous bone volume (-21.7%) at 11 days post exposure. These results suggest that early radiation induced changes in FGF18 gene expression in skin may have value for predicting subsequent loss of cancellous bone mass.

space radiation↗

Measurement of 𝑑 2⁢ 𝜎/𝑑⁢|$\vec{q}$|⁢𝑑⁢𝐸 avail in charged current 𝜈 𝜇 -nucleus interactions at ⟨𝐸 𝜈 ⟩=1.86 GeV using the NOvA Near Detector

Double- and single-differential cross sections for inclusive charged-current 𝜈 𝜇 -nucleus scattering are reported for the kinematic domain 0 to 2 GeV/𝑐 in three-momentum transfer and 0 to 2 GeV in available energy, at a mean 𝜈 𝜇 energy of 1.86 GeV. The measurements are based on an estimated 995,760 𝜈 𝜇 charged-current (CC) interactions in the scintillator medium of the NOvA Near Detector. The subdomain populated by 2-particle-2-hole (2p2h) reactions is identified by the cross section excess relative to predictions for 𝜈 𝜇 -nucleus scattering that are constrained by a data control sample. Models for 2-particle-2-hole processes are rated by 𝜒 2 comparisons of the predicted-versus-measured 𝜈 𝜇 CC inclusive cross section over the full phase space and in the restricted subdomain. Shortfalls are observed in neutrino generator predictions obtained using the theory-based València and SuSAv2 2p2h models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Vibrational Power Flow Analysis of Rods and Beams

A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.

Wohlever, James Christopher↗

Smart patch piezoceramic actuator issues

The Phillips Laboratory is undertaking the challenge of finding new and innovative ways to integrate sensing, actuation, and the supporting control and power electronics into a compact self-contained unit to provide vibration suppression for a host structure. This self-contained unit is commonly referred to as a smart patch. The interfaces to the smart patch will be limited to standard spacecraft power and possibly a communications line. The effort to develop a smart patch involves both contractual and inhouse programs which are currently focused on miniaturization of the electronics associated with vibrational control using piezoceramic sensors and actuators. This paper is comprised of two distinct parts. The first part examines issues associated with bonding piezoceramic actuators to a host structure. Experimental data from several specimens with varying flexural stiffness are compared to predictions from two piezoelectric/substructure coupling models, the Blocked Force Model and the Uniform Strain Model with Perfect Bonding. The second part of the paper highlights a demonstration article smart patch created using the insights gained from inhouse efforts at the Phillips Laboratory. This demonstration article has self contained electronics on the same order of size as the actuator powered by a voltage differential of approximately 32 volts. This voltage is provided by four rechargeable 8 volt batteries.

Griffin, Steven F.↗

Scientific Discovery with Physics-Informed System Identification (Abbreviated Report)

My fellowship research focused on making physics-based simulations faster and more useful through machine learning. Many problems in science and engineering are governed by partial differential equations, but high-fidelity simulations are often too expensive to run repeatedly. I worked on improving Latent Space Dynamics Identification (LaSDI), a reduced-order modeling framework that compresses large simulation data sets into a smaller representation and then learns how that representation evolves over time. The motivation was to develop reduced models that remain accurate for more challenging systems, especially when predictions must remain reliable over long time intervals or when the underlying dynamics are more complicated than standard methods can easily handle. I also contributed to related work on Quandary, a high-performance software effort for simulation and control of open quantum systems, before focusing primarily on Latent Space Dynamics Identification methods. The main outcomes of the fellowship were two new algorithms (both of which were published), Rollout-LaSDI and Higher-Order LaSDI, together with supporting work on multi-stage Latent Space Dynamics Identification. Rollout-LaSDI improved long-term prediction by training the model to stay accurate over extended time horizons, and Higher-Order LaSDI broadened the method so it could model systems with higher-order time dynamics. My contributions to multistage Latent Space Dynamics Identification also helped show that its later training stages could be simplified without losing effectiveness, and that this behavior held across different model architectures and training strategies. Taken together, these advances improved the accuracy, flexibility, and practical value of reduced-order modeling tools for computational science.

97 MATHEMATICS AND COMPUTING↗

Four-dimensional phase-space reconstruction of flat and magnetized beams using neural networks and differentiable simulations

Beams with cross-plane coupling or extreme asymmetries between the two transverse phase spaces are often encountered in particle accelerators. Flat beams with large transverse-emittance ratios are critical for future linear colliders. Similarly, magnetized beams with significant cross-plane coupling are expected to enhance the performance of electron cooling in hadron beams. Preparing these beams requires precise control and characterization of the four-dimensional transverse phase space. In this study, we employ generative phase-space reconstruction techniques to rapidly characterize magnetized and flat-beam phase-space distributions using a conventional quadrupole-scan method. The reconstruction technique is experimentally demonstrated on an electron beam produced at the Argonne Wakefield Accelerator and successfully benchmarked against conventional diagnostics techniques. Specifically, we show that predicted beam parameters from the reconstructed phase-space distributions (e.g., as magnetization and flat-beam emittances) are in excellent agreement with those measured from the conventional diagnostic methods. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Blade-mounted trailing edge flap control for BVI noise reduction

Numerical procedures based on the 2-D and 3-D full potential equations and the 2-D Navier-Stokes equations were developed to study the effects of leading and trailing edge flap motions on the aerodynamics of parallel airfoil-vortex interactions and on the aerodynamics and acoustics of the more general self-generated rotor blade vortex interactions (BVI). For subcritical interactions, the 2-D results indicate that the trailing edge flap can be used to alleviate the impulsive loads experienced by the airfoil. For supercritical interactions, the results show the necessity of using a leading edge flap, rather than a trailing edge flap, to alleviate the interaction. Results for various time dependent flap motions and their effect on the predicted temporal sectional loads, differential pressures, and the free vortex trajectories are presented. For the OLS model rotor, contours of a BVI noise metric were used to quantify the effects of the trailing edge flap on the size and directivity of the high/low intensity noise region(s). Average reductions in the BVI noise levels on the order of 5 dB with moderate power penalties on the order of 18 pct. for a four bladed rotor and 58 pct. for a two bladed rotor were obtained.

Hassan, A. A.↗

Commentary: Duckweeds as model organisms for metabolic studies

Duckweeds have many practical applications, for example in human nutrition, as animal feed, in the production of bioplastics or vaccines, and phytoremediation (Acosta et al., 2021). Under most conditions, they reproduce asexually which provides genetically uniform material with predictable patterns of growth that make them ideal as sentinel organisms for phytotoxicity testing (Park et al., 2021). Asexual growth also results in high biomass production which makes duckweeds promising candidates as biofuel feedstocks (Acosta et al., 2021; Liang et al., 2023). In addition, duckweed species like Lemna minor and Spirodela polyrhiza are also reemerging as model organisms in plant biology as high-quality full genome assemblies and other genomic resources become available (Chang et al., 2016; Acosta et al., 2021). We argue that duckweed species are particularly of interest for the study of primary plant metabolism. Primary metabolism concerns the part of metabolism that is directly involved in the growth and development of plants, and which tends to be highly conserved among plant species. What makes duckweeds particularly attractive is that when grown on liquid media more precise control of physiological conditions can be attained relative to growth of plants in soil. Also, due to their relatively simple anatomical structure and asexual reproduction of fronds by budding, precise characterization of the physiological state under study is possible through one simple metric, i.e., the specific growth rate (rate of dry weight increase per existing dry weight), which can be incorporated relatively easily into metabolic models. This is not possible for land plants, such as Arabidopsis, because over the course of their life cycle, they go through multiple growth stages and phases of anatomical differentiation, which are much more complex to quantify. Furthermore, duckweeds can grow on organic substrates under heterotrophic or photomixotrophic conditions that facilitate isotope tracer studies. For example, in a previous study on duckweed by one of the authors, Lemna gibba (L). was grown on glucose with a position-specific 13 C-label that can be detected and resolved by Mass Spectrometry or Nuclear Magnetic Resonance spectrometry. Using this approach, the 13 C-label was traced into biomass compounds formed from glucose, particularly isoprenoid compounds. Some of the resulting labeling patterns were in apparent disagreement with predictions based on known metabolic pathways for the biosynthesis of isopentenyl pyrophosphate, the universal building block for isoprenoids, (Lichtenthaler et al., 1997). From this data it was deduced that isoprenoid compounds such as carotenoids and isoprenoid chains of phytol and plastoquinone, synthesized in the chloroplast, are produced via a previously unreported plant metabolic pathway, now known as the methylerythitol/deoxyxylulose-5-phosphate pathway (Lichtenthaler et al., 1997).

59 BASIC BIOLOGICAL SCIENCES↗

The effects of simulated hypogravity on murine bone marrow cells

Mouse bone marrow cells grown in complete medium at unit gravity were compared with a similar population cultured in conditions that mimic some aspects of microgravity. After the cells adjusted to the conditions that simulated microgravity, they proliferated as fetal or oncogenic populations; their numbers doubled in twelve hour periods. Differentiated subpopulations were depleted from the heterogeneous mixture with time and the undifferentiated hematopoietic stem cells increased in numbers. The cells in the control groups in unit gravity and those in the bioreactors in conditions of microgravity were monitored under a number of parameters. Each were phenotyped as to cell surface antigens using a panel of monoclonal antibodies and flow cytometry. Other parameters compared included: pH, glucose uptake, oxygen consumption and carbon-dioxide production. Nuclear DNA was monitored by flow cytometry. Functional responses were studied by mitogenic stimulation by various lectins. The importance of these findings should have relevance to the space program. Cells should behave predictably in zero gravity; specific populations can be eliminated from diverse populations and other populations isolated. The availability of stem cell populations will enhance both bone marrow and gene transplant programs. Stem cells will permit developmental biologists study the paths of hematopoiesis.

Lawless, Desales↗

Dynamic balance control in elders: gait initiation assessment as a screening tool

OBJECTIVE: To determine whether measurements of center of gravity-center of pressure separation (CG-CP moment arm) during gait initiation can differentiate healthy from disabled subjects with sufficient specificity and sensitivity to be useful as a screening test for dynamic balance in elderly patients. SUBJECTS: Three groups of elderly subjects (age, 74.97+/-6.56 yrs): healthy elders (HE, n = 21), disabled elders (DE, n = 20), and elders with vestibular hypofunction (VH, n = 18). DESIGN: Cross-sectional, intact-groups research design. Peak CG-CP moment arm measures how far the subject will tolerate the whole-body CG to deviate from the ground reaction force's CP; it represents dynamic balance control. Screening test cutoff points at 16 to 18 cm peak CG-CP moment arm predicted group membership. RESULTS: The magnitude of peak CG-CP moment arm was significantly greater in HE than in DE and VH subjects (p<.01) and was not different between the DE and VH groups. The peak CG-CP moment arm occurred at the end of single stance phase in all groups. As a screening test, the peak moment arm has greater than 50% sensitivity and specificity to discriminate the HE group from the DE and VH groups with peak CG-CP moment arm cutoff points between 16 and 18 cm. CONCLUSIONS: Examining dynamic balance through the use of the CG-CP moment arm during single stance in gait initiation discriminates between nondisabled and disabled older persons and warrants further investigation as a potential tool to identify people with balance dysfunction.

NASA Discipline Neuroscience↗

Probabilistic error estimation for non-intrusive reduced models learned from data of systems governed by linear parabolic partial differential equations

This work derives a residual-based a posteriori error estimator for reduced models learned with non-intrusive model reduction from data of high-dimensional systems governed by linear parabolic partial differential equations with control inputs. It is shown that quantities that are necessary for the error estimator can be either obtained exactly as the solutions of least-squares problems in a non-intrusive way from data such as initial conditions, control inputs, and high-dimensional solution trajectories or bounded in a probabilistic sense. Here, the computational procedure follows an offline/online decomposition. In the offline (training) phase, the high-dimensional system is judiciously solved in a black-box fashion to generate data and to set up the error estimator. In the online phase, the estimator is used to bound the error of the reduced-model predictions for new initial conditions and new control inputs without recourse to the high-dimensional system. Numerical results demonstrate the workflow of the proposed approach from data to reduced models to certified predictions.

97 MATHEMATICS AND COMPUTING↗

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗