Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “motion correction”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A Pn magnitude scale m b ( Pn ) for earthquakes along the equatorial Mid-Atlantic Ridge

We developed a short-period Pn magnitude scale m b (Pn) for earthquakes along the equatorial Mid-Atlantic Ridge. Due to low signal-to-noise ratios, teleseismic body wave magnitude and long-period surface wave magnitude cannot be confidently determined for small earthquakes of m b < 4. Local magnitude scales are also not useful for these events because the oceanic environment does not allow the propagation of crustal phases. However, regional high-frequency Pn waves from these small- to moderate-size (m b 3–6) earthquakes are well recorded in the equatorial Atlantic region and can be used to assign magnitudes. We measured over 2041 Pn peak amplitudes on vertical records from about 21 stations in northeastern Brazil and 11 stations in western Africa in the distance range of 700–3700 km. We analysed data from 189 events from the global centroid moment tensor catalogue to tie our m b (Pn) scale to M w so that seismic moments can be readily estimated. Pn arrivals show apparent group velocity between 7.9 km s –1 at short ranges (~1000 km) and up to 9.1 km s –1 at 3500 km. The measured peak amplitudes have a frequency between 0.8 and 3 Hz at 1000–1800 km, but at greater distances, 1800–3700 km, they show a remarkably consistent frequency of about 0.8 Hz. The peak amplitude attenuates at a higher rate at short distances (~0.65 magnitude units between 700 and 2000 km) but attenuates at a lower rate at long distances (~0.35 magnitude units between 2000 and 3700 km). The low rate of amplitude decay with distance and nearly constant frequency content of the peak amplitudes suggest that Pn waves propagate efficiently in the lower part of the upper mantle in the equatorial Atlantic Ocean basins. Furthermore, these are important attributes of oceanic Pn waves that can be used to assign magnitude for small- to moderate-size earthquakes in the equatorial mid-Atlantic region. The estimated station corrections correlate well with upper mantle low-velocity anomalies, especially in Brazil.

58 GEOSCIENCES↗
DOE OSTI2022

Multispectral and Thermal Imager Onboard Aerial Platforms Instrument Handbook

The MicaSense Altum imager is an off-the-shelf, synchronized multispectral and thermal camera, combined with a Global Positioning System (GPS) unit and downwelling light sensor (DLS). The Altum takes images of the land surface within its field of view (FOV) across five visible bands and one longwave infrared thermal band. The sensor records the radiance and converts it to digital numbers. Imagery is radiometrically corrected, taking into account the sensor calibration, lens distortions, vignette effects, sun angle, and atmospheric effects (scattering and absorption). The photogrammetry software Agisoft PhotoScan v 1.4 is used to align and stitch the images into a larger composite image using the technique of structure from motion image capture to construct a dense cloud and 3D model of the surface, which is used to produce a digital elevation model of the terrain surveyed and orthomosaic imagery. Land surface leaf area index, albedo, and surface skin temperature are provided for the user. The user can perform raster calculations on the imagery to produce various other vegetative indices, which are used to indicate plant health and land surface characteristics, including Normalized Difference Vegetation Index (NDVI), Green Normalized Difference Vegetation Index (GNDVI), Normalized Difference Water Index (NDWI), surface temperature, and Enhanced_vegetation_index (EVI). This code is provided with the readme file for each data set. The images have a resolution of 20-60 cm/pixel, which are subsampled to 100 cm/pixel.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Quantifying Site Effects and Their Influence on Earthquake Source Parameter Estimations Using a Dense Array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes ( M L = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures.

58 GEOSCIENCES↗
DOE OSTI2021

Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control

Large-Eddy Simulations (LES) of a gasoline spray, where the mixture was ignited rapidly during or after injection, were performed in comparison to a previous experimental study with quantitative flame motion and soot formation data [SAE 2020-01-0291] and an accompanying Reynolds-Averaged Navier–Stokes (RANS) simulation at the same conditions. The present study reveals major shortcomings in common RANS combustion modeling practices that are significantly improved using LES at the conditions of the study, specifically for the phenomenon of rapid ignition in the highly turbulent, stratified mixture. At different ignition timings, benchmarks for the study include spray mixing and evaporation, flame propagation after ignition, and soot formation in rich mixtures. A comparison of the simulations and the experiments showed that the LES with Dynamic Structure turbulence were able to capture correctly the liquid penetration length, and to some extent, spray collapse demonstrated in the experiments. For early and intermediate ignition timings, the LES showed excellent agreement to the measurements in terms of flame structure, extent of flame penetration, and heat-release rate. However, RANS simulations (employing the common G-equation or well-stirred reactor) showed much too rapid flame spread and heat release, with connections to the predicted turbulent kinetic energy. With confidence in the LES for predicted mixture and flame motion, the predicted soot formation/oxidation was also compared to the experiments. The soot location was well captured in the LES, but the soot mass was largely underestimated using the empirical Hiroyasu model. An analysis of the predicted fuel–air mixture was used to explain different flame propagation speeds and soot production tendencies when varying ignition timing.

42 ENGINEERING↗
DOE OSTI2026

A curated benchmark for cofolding models on kinase conformational states

Abstract Protein kinases are critical drug targets, requiring therapeutics that can modulate their active and inactive conformational states. While cofolding models can generate global folds directly from kinase sequences and ligand SMILES strings, these models have not yet been tested on their ability to recover ligand-induced-fit conformational states of the kinase proteins. Here, we introduce KinConfBench, a curated benchmark of 2225 high-quality human kinase chains to evaluate the ability of four state-of-the-art cofolding models—Boltz-2, Chai-1, Protenix, and RoseTTAFold-All-Atom—to recover both canonical and rare conformational states. We show that geometric success metrics of a ligand pose in the active site do not correlate strongly with the correct kinase conformational state, motivating a new set of dynamical benchmarks for assessing cofolding models. While all four cofolding models achieve ~60–80% prediction accuracy for kinase conformational classification, they exhibit severe mode collapse when performing multiple inferences, show negligible structural diversity in sampling induced-fit motions, and display a prevalent “apo-drift” in which most cofolding models predominantly predict the kinase to be in its ligand-free state. Our results highlight that capturing ligand-induced protein conformational diversity, not just geometric fit, is critical for next-generation structure-based drug discovery.

Sun, Kunyang↗
DOE OSTI2023

A Blueprint for the Milky Way’s Stellar Populations. IV. A String of Pearls—the Galactic Starburst Sequence

We continue our series of papers on phase-space distributions of stars in the Milky Way based on photometrically derived metallicities and Gaia astrometry, with a focus on the halo-disk interface in the local volume. To exploit various photometric databases, we develop a method of empirically calibrating synthetic stellar spectra based on a comparison with observations of stellar sequences and individual stars in the Sloan Digital Sky Survey, the SkyMapper Sky Survey, and the Pan-STARRS1 surveys, overcoming band-specific corrections employed in our previous work. In addition, photometric zero-point corrections are derived to provide an internally consistent photometric system with a spatially uniform metallicity zero-point. Using our phase-space diagrams, we find a remarkably narrow sequence in the rotational velocity (v $\phi$ ) versus metallicity ([Fe/H]) space for a sample of high proper-motion stars (>25 mas yr -1 ), which runs along Gaia Sausage/Enceladus (GSE) and the Splash substructures and is linked to the disk, spanning nearly 2 dex in [Fe/H]. Notably, a rapid increase of v $\phi$ from a nearly zero net rotation to ~180 km s -1 in a narrow metallicity interval (-0.6 ≲ [Fe/H] ≲ -0.4) suggests that some of these stars emerged quickly on a short gas-depletion timescale. Through measurements of a scale height and length, we argue that these stars are distinct from those heated dynamically by mergers. This chain of high proper-motion stars provides additional support for recent discoveries suggesting that a starburst took place when the young Milky Way encountered the gas-rich GSE progenitor, which eventually led to the settling of metal-enriched gas onto the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Review of analytical studies on seismic fluid-structure interaction of base-supported cylindrical tanks

Analytical solutions for seismic fluid-structure interaction (FSI) of tanks are used by analysts and engineers for a number of purposes, including preliminary sizing, risk assessment, and verification of numerical models. Analytical solutions for seismic FSI of tanks were first developed in the 1930s and have been extended and/or modified in the following nine decades to accommodate different seismic inputs and boundary conditions. This report reviews published analytical solutions for rigid and flexible, base-supported, cylindrical tanks subjected to unidirectional horizontal motion of a small amplitude. The studies parsed the fluid-structure (FSI) response into an impulsive and a convective component, and addressed frequencies of lateral motions of the tank (i.e., impulsive frequencies), hydrodynamic pressures, frequencies of waves (i.e., convective frequencies), wave heights, and reactions at the base. The solutions were categorized into exact and approximate solutions based on the methodologies used. Herein, FSI responses are calculated and compared for a range of tank dimensions using the analytical solutions from the different studies. The responses presented here are normalized to be unitless, and can be used for tanks with different dimensions and mechanical properties, and subjected to different input motions. A sample steel water tank is analyzed using the normalized solutions for two earthquake ground motions with different frequency contents. Ignoring tank flexibility may significantly underestimate seismic responses of the tank and its contained fluid. The approximate solutions, which have been widely used for decades, may not be sufficiently accurate in the modern era. Calculation errors in the reviewed studies are identified and corrected, and differences are quantified for the sample tank.

42 ENGINEERING↗
DOE OSTI2021

Lower-hybrid drift waves and their interaction with plasmas in a 3D symmetric reconnection simulation with zero guide field

We investigate lower-hybrid drift waves (LHDW) in symmetric magnetic reconnection with zero guide field using three-dimensional particle-in-cell simulations. The long-wavelength mode with kρiρe∼1 develops in the bifurcated electron current layer around the X-line within the width of the electron meandering motion from the mid-plane, where ρi(e) is the ion (electron) gyroradius. The short-wavelength mode with kρe∼1 develops in the separatrix region downstream of the electron outflow jet, producing electron vortices in the background flow frame. Electrons follow the E × B drift with corrections from the diamagnetic drift and are heated inside the vortices with diverging electric fields. In the vortices, ions have comparable E × B and inertia drifts, which together mostly cancel the diamagnetic drift. Toward the center of diverging field vortices, ions are decelerated, leading to a decrease in the perpendicular temperature, while the loss of low-energy ions results in an increase in the parallel temperature. Parallel electric fields exist as a combination of the LHDW wave field projected to the magnetic field direction and the penetration of whistler waves that are mainly outside of the LHDW layer. The magnetic flux tube is twisted in the vortices. The twist may potentially lead to slippage reconnection, as indicated by the non-uniform parallel potential variation across field lines, while the periodic variations of the twisting directions are a limiting factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mason–Weaver theory: Revised and extended for a semi-infinite domain

Mason and Weaver developed equations to describe small particles settling under the influence of gravity and Brownian motion, including the limiting case for an infinitely deep suspension. We encountered this common convection–diffusion equation and no-flux boundary conditions in a model for dynamics of adsorbed polymers in dead end pores of a depolymerization catalyst. Close examination reveals that the Mason–Weaver solution is not correct for the infinite domain with a non-uniform initial condition. In this paper, we obtain the time dependent Green’s function for the no flux boundary condition and also for a more general reactive boundary condition. We demonstrate how the results provide solutions, via superposition, which provide solutions for several boundary conditions and all initial conditions.

Chen, Ziqiu↗
DOE OSTI2021

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

We report accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume SEVEN-LEAGUE HYDRO code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the α-β and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo–LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Error estimates of finite difference methods for the Dirac equation in the massless and nonrelativistic regime

We present four frequently used finite difference methods and establish the error bounds for the discretization of the Dirac equation in the massless and nonrelativistic regime, involving a small dimensionless parameter 0 < ε &NestedLessLess; 1 inversely proportional to the speed of light. In the massless and nonrelativistic regime, the solution exhibits rapid motion in space and is highly oscillatory in time. Specifically, the wavelength of the propagating waves in time is at O(ε), while in space, it is at O(1) with the wave speed at O(ε -1 ). We adopt one leap-frog, two semi-implicit, and one conservative Crank-Nicolson finite difference methods to numerically discretize the Dirac equation in one dimension and establish rigorously the error estimates which depend explicitly on the time step τ, mesh size h, and the small parameter ε. The error bounds indicate that, to obtain the “correct” numerical solution in the massless and nonrelativistic regime, i.e., 0 < ε &NestedLessLess; 1, all these finite difference methods share the same ε-scalability as time step τ = O(ε 3/2 ) and mesh size h = O(ε 1/2 ). A large number of numerical results are reported to verify the error estimates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code

A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code: Preprint

A passive yaw implementation is developed, validated, and explored for the Wave Energy Converter Simulator (WEC-Sim), an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01 was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than or equal to 1, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2024

Aberration corrected RF flipper for high resolution neutron spectroscopy

Project Summary Company: Adelphi Technology, Inc. Title: Aberration-corrected High Frequency RF Flipper for High-Resolution Neutron Spectroscopy PI: Dr. Jay Theodore Cremer Topic: C55-11 Enhancement of Scattering Instrumentation Technology Used at Pulsed and Continuous Sources Subtopic: d. Other Statement of the problem or situation that is being addressed. The quest to understand heterogeneous and hierarchical materials is gathering momentum, as described in a 2015 report by the Basic Energy Sciences Advisory Committee on Challenges at the Frontiers of Matter and Energy. For the past 40 years a technique called neutron spin echo (NSE) has been used to probe molecular motions in such non-crystalline materials over time scales from 10’s of picoseconds to 100’s of nanoseconds. The method has provided unique information about the dynamics of soft heterogeneous materials, including confirmation of the de Gennes model of polymer reptation and quantitative measurement of bending constants of biologically relevant lipid membranes. However, scientists continue to clamor for even higher resolution than NSE can provide. Biomaterials, polymers, glasses, and artificially nanostructured materials all manifest slow molecular motions because of weak or competing interactions between subunits and are amenable to study with neutrons, provided sufficiently long dynamical correlation times can be achieved. All these materials have important applications to advanced technologies so understanding them is key to technological progress. General statement of how this problem is being addressed. We will address the need for high-resolution neutron spectroscopy by using a technique called Neutron Resonance Spin Echo (NRSE). While similar to NSE in many respects, this method has the potential to exceed the NSE capabilities, if 2 technical hurdles can be overcome. The major impediments to successful high-resolution NRSE are the availability of two technologies: a very high frequency, efficient, radiofrequency (rf) flipper for neutrons and a method to correct certain magnetic aberrations. Based on previous STTR support and follow-on research we have developed a suitable rf flipper and we have invented a method to correct the magnetic aberrations. Both technologies need refinement to make them suitable for implementation at a neutron source such as the Oak Ridge National Laboratory nuclear reactor. In this proposal we seek to perfect the two technologies and to combine them into a single, operationally convenient device. Commercial Applications and Other Benefits In view of the increasing demand for the unique scientific information that high resolution neutron spectroscopy can provide, we expect several major instrumentation upgrades at both U.S. and foreign neutron centers will require make use of the NRSE method over the coming decade, creating a market for the devices we will design. These components will enhance scientists’ abilities to probe the time dependence of density fluctuations in a wide range of hierarchical and heterogeneous materials many of which are vital to existing and future technologies. Key Words – Polarized Neutrons, Neutron Spin Echo, Neutron Scattering, advanced materials. Summary for Members of Congress Neutron beams are a powerful materials-science probe that provide unique information about the structure of matter. The proposed devices will accelerate scientific discoveries required to achieve national goals for new technological materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Yield and PDV Measurements of the Godiva-IV Critical Assembly

A high dynamic range, scinillator based detector, the MHD-240, was used to measure the gamma and neutron flux of several prompt critical bursts of the Godiva-IV assembly. The measurements provided the alpha and FWHM for several different burst temperatures. Calibrations of the MHD240 gamma and neutron sensitivity along with MCNP based simulation of the Godiva-IV assembly allowed us to also determine the total fission yield of the burst integrated over approximately 5 ms. In addition a Photo Doppler Velocimetry measurement was conducted to determine the surface motion vibrational frequency of Godiva-IV assembly during the bursts. The combined data were collected to be used as input to multi-physics simulation of the Godiva-IV assembly. This is a preliminary report, there is an ongoing effort to update the MCNP simulations used in the analysis and develop a correction for the room return component of the measured signal.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Baryogenesis from ultra-slow-roll inflation

The ultra-slow-roll (USR) inflation represents a class of single-field models with sharp deceleration of the rolling dynamics on small scales, leading to a significantly enhanced power spectrum of the curvature perturbations and primordial black hole (PBH) formation. Such a sharp transition of the inflationary background can trigger the coherent motion of scalar condensates with effective potentials governed by the rolling rate of the inflation field. We show that a scalar condensate carrying (a combination of) baryon or lepton number can achieve successful baryogenesis through the Affleck-Dine mechanism from unconventional initial conditions excited by the USR transition. Viable parameter space for creating the correct baryon asymmetry of the Universe naturally incorporates the specific limit for PBHs to contribute significantly to dark matter, shedding light on the cosmic coincidence problem between the baryon and dark matter densities today.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The Pantheon+ analysis: Improving the redshifts and peculiar velocities of Type Ia supernovae used in cosmological analyses

We examine the redshifts of a comprehensive set of published Type Ia supernovae, and provide a combined, improved catalogue with updated redshifts. We improve on the original catalogues by using the most up-to-date heliocentric redshift data available; ensuring all red shifts have uncertainty estimates; using the exact formulae to convert heliocentric redshifts into the Cosmic Microwave Background (CMB) frame; and utilising an improved peculiar velocity model that calculates local motions in redshift-space and more realistically accounts for the external bulk flow at high-redshifts. We review 2607 supernova redshifts; 2285 are from unique supernovae and 322 are from repeat observations of the same supernova. In total, we updated 990 unique heliocentric redshifts, and found 5 cases of missing or incorrect heliocentric corrections, 44 incorrect or missing supernova coordinates, 230 missing heliocentric or CMB frame redshifts, and 1200 missing redshift uncertainties. The absolute corrections range between 10 -8 ≤ Δz ≤ 0.038, and RMS(Δz) ~ 3 × 10 -3 . The sign of the correction was essentially random, so the mean and median corrections are small: 4 × 10 -4 and 4 × 10 -6 respectively. We examine the impact of these improvements for H 0 and the dark energy equation of state w and find that the cosmological results change by ΔH 0 = -0.12 km s -1 Mpc -1 and Δw = 0.003, both significantly smaller than previously reported uncertainties for H 0 of 1.0 km s -1 Mpc -1 and w of 0.04 respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A physics-constrained neural network for multiphase flows

The present study develops a physics-constrained neural network (PCNN) to predict sequential patterns and motions of multiphase flows (MPFs), which includes strong interactions among various fluid phases. To predict the order parameters, which locate individual phases in the future time, a neural network (NN) is applied to quickly infer the dynamics of the phases by encoding observations. The multiphase consistent and conservative boundedness mapping algorithm (MCBOM) is next implemented to correct the predicted order parameters. This enforces the predicted order parameters to strictly satisfy the mass conservation, the summation of the volume fractions of the phases to be unity, the consistency of reduction, and the boundedness of the order parameters. Then, the density of the fluid mixture is updated from the corrected order parameters. Finally, the velocity in the future time is predicted by another NN with the same network structure, but the conservation of momentum is included in the loss function to shrink the parameter space. The proposed PCNN for MPFs sequentially performs (NN)-(MCBOM)-(NN), which avoids nonphysical behaviors of the order parameters, accelerates the convergence, and requires fewer data to make predictions. Numerical experiments demonstrate that the proposed PCNN is capable of predicting MPFs effectively.

Mechanics↗