Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix Models”

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 1,099 records · Page 61

Impact of Concurrent Solubilization and Fines Migration on Fracture Aperture Growth in Shales during Acidized Brine Injection

Acidic hydraulic fracturing fluid chemically and physically alters shale rock fabric during injection and shut-in, creating a “reaction-altered zone” along the fracture faces. To better characterize the variable thickness and composition of this reaction-altered zone under advective flow, we take a coupled experimental and modeling approach. A fluidic cell, with six fiducial markers, is first fabricated to keep the rock sample in place during the core floods and to allow image alignment of acquired images. Then, we conduct a series of reactive core floods in a clay-rich siliceous Wolfcamp shale sample with 10 wt % carbonate, using a synthetic fracturing fluid under no confining stress and at room temperature. High-resolution computed tomography (CT) scans are periodically conducted to observe the spatial alteration of the fracture network. We then perform scanning electron microscopy (SEM-EDS) on the two orthogonal surfaces (fracture surface and freshly cut profile face) to generate high-resolution elemental maps that show the change in mineralogy, both with distance along a given flow path along a fracture surface and with depth from the fracture surface into the shale matrix. These results are contrasted against a two-dimensional (2D) advection-diffusion reaction model developed previously for batch reactions between shale and synthetic fracturing fluids. The model simulates the geochemical interaction occurring at the fracture/matrix interface and penetrating into the shale matrix during the reactive core flood. Both model and experimental results show that the acidic brine is neutralized during the core flood, corresponding to an increase in fracture aperture as a function of fluid volume injected with the greatest change near the inlet. SEM-EDS scans reveal significant dissolution of carbonates on the fracture surface without pyrite oxidation. The reactive transport model indicates that carbonate depletion into the shale interior should be observable, yet SEM-EDS shows no discernible loss of carbonate in the orthogonal profile face. The combination of these observations suggests an additional fracture evolution mechanism in the reactive system, i.e., fines migration. We show that fines migration enhances the access of fracturing fluid to the matrix resulting in a more pronounced fracture widening. Finally, we conclude that coupled mineral dissolution and fines migration govern fracture aperture growth during acidized brine injection. In this work, we effectively show the underlying risk of relying solely on models that do not include an important (transport) process that can alter the system significantly and propose a combined chemomechanical mechanism for fracture evolution appropriate for this shale mineralogy.

58 GEOSCIENCES↗

Secondary production of massive quarks in thrust

We present a factorization framework that takes into account the production of heavy quarks through gluon splitting in the thrust distribution for e{sup +}e{sup −} → hadrons. The explicit factorization theorems and some numerical results are displayed in the dijet region where the kinematic scales are widely separated, which can be extended systematically to the whole spectrum. We account for the necessary two-loop matrix elements, threshold corrections, and include resummation up to N{sup 3}LL order. We include nonperturbative power corrections through a field theoretical shape function, and remove the O(Λ{sub QCD}) renormalon in the partonic soft function by appropriate mass-dependent subtractions. Our results hold for any value of the quark mass, from an infinitesimally small (merging to the known massless result) to an infinitely large one (achieving the decoupling limit). This is the first example of an application of a variable flavor number scheme to final state jets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Temperature dependent dielectric and conductivity studies of polyvinyl alcohol-ZnO nanocomposite films by impedance spectroscopy

Dielectric and conductivity behaviors of nano ZnO doped polyvinyl alcohol (PVA) composites for various concentrations of dopant were investigated using impedance spectroscopy for a wide range of temperatures (303 K–423 K) and frequencies (5 Hz–30 MHZ). The dielectric properties of host polymer matrix have been improved by the addition of nano ZnO and are found to be highly temperature dependent. Anomalous dielectric behavior was observed in the frequency range of 2.5 MHz–5 MHz. Increase in dielectric permittivity and dielectric loss was observed with respect to temperature. The Cole-Cole plot could be modeled by low resistance regions in a high resistance matrix and the lowest resistance was observed for the 10 mol. % films. The imaginary part of the electric modulus showed asymmetric peaks with the relaxation following Debye nature below and non-Debye nature above the peaks. The ac conductivity is found to obey Jonscher's power law, whereas the variation of dc conductivity with temperature was found to follow Arrhenius behavior. Two different activation energy values were obtained from Arrhenius plot indicating that two conduction mechanisms are involved in the composite films. Fitting the ac conductivity data to Jonscher's law indicates that large polaron assisted tunneling is the most likely conduction mechanism in the composites. Maximum conductivity is observed at 423 K for all the samples and it is optimum for 10 mol. % ZnO doped PVA composite film. Significant increase in dc and ac conductivities in these composite films makes them a potential candidate for application in electronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A multidimensional approach to quantum state tomography of photoelectron wavepackets

There is a growing interest in reconstructing the density matrix of photoelectron wavepackets, in particular in complex systems where decoherence can be introduced either by a partial measurement of the system or through coupling with a stochastic environment. To this end, several methods to reconstruct the density matrix, quantum state tomography protocols, have been developed and tested on photoelectrons ejected from noble gases following absorption of extreme ultraviolet (XUV) photons from attosecond pulses. It remains a challenge to obtain model-free, single scan protocols that can reconstruct the density matrix with high fidelities. Current methods require extensive measurements or involve complex fitting of the signal. Efficient single-scan reconstructions would be of great help to increase the number of systems that can be studied. We propose a new and more efficient protocol that is able to reconstruct the continuous variable density matrix of a photoelectron in a single time delay scan. It is based on measuring the coherences of a photoelectron created by absorption of an XUV pulse using a broadband infrared (IR) probe that is scanned in time and a narrowband IR reference that is temporally fixed to the XUV pulse. We illustrate its performance for a Fano resonance in He as well as mixed states in Ar arising from spin-orbit splitting. We show that the protocol results in excellent fidelities and near-perfect estimation of the purity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

MicroRNA-130a and -130b enhance activation of hepatic stellate cells by suppressing PPARγ expression: A rat fibrosis model study

Hepatic stellate cells (HSCs) are the primary sources of extracellular matrix (ECM) in normal and fibrotic liver. Peroxisome proliferator-activated receptor gamma (PPARγ) maintains HSCs in a quiescent state, and its downregulation induces HSC activation. MicroRNAs (miRNAs) can induce PPARγ mRNA degradation, but the mechanism by which miRNAs regulate PPARγ in rat HSCs is unclear. This study aimed to investigate some miRNAs which putatively bind to the 3′-untranslated region (3′-UTR) of PPARγ mRNA, and increase expression of ECM genes in rat HSCs. In carbon tetrachloride injection (CCl{sub 4}) and common bile duct ligation (CBDL) liver fibrosis models, miRNAs miR-130a, miR-130b, miR-301a, miR-27b and miR-340 levels were found to be increased and PPARγ expression decreased. Overexpression of miR-130a and miR-130b enhanced cell proliferation by involving Runx3. MiR-130a and miR-130b decreased PPARγ expression by targeting the 3′-UTR of PPARγ mRNA in rat HSC-T6 cells. Transforming growth factor-β1 (TGF-β1) may mediate miR-130a and miR-130b overexpression, PPARγ downregulation, and ECM genes overexpression in cell culture. These findings suggest that miR-130a and miR-130b are involved in downregulation of PPARγ in liver fibrosis. - Highlights: • MiR-130a and miR-130b are increased and PPARγ is decreased in liver fibrosis models. • MiR-130a and miR-130b decreased PPARγ by targeting the 3′-UTR of PPARγ mRNA. • MiR-130a and miR-130b enhanced HSC cell proliferation by involving Runx3. • TGF-β1 may mediate miR-130a and miR-130b overexpression.

60 APPLIED LIFE SCIENCES↗

Constraining astrophysical neutrino flavor composition from leptonic unitarity

The recent IceCube observation of ultra-high-energy astrophysical neutrinos has begun the era of neutrino astronomy. In this work, using the unitarity of leptonic mixing matrix, we derive nontrivial unitarity constraints on the flavor composition of astrophysical neutrinos detected by IceCube. Applying leptonic unitarity triangles, we deduce these unitarity bounds from geometrical conditions, such as triangular inequalities. These new bounds generally hold for three flavor neutrinos, and are independent of any experimental input or the pattern of lepton mixing. We apply our unitarity bounds to derive general constraints on the flavor compositions for three types of astrophysical neutrino sources (and their general mixture), and compare them with the IceCube measurements. Furthermore, we prove that for any sources without ν{sub τ} neutrinos, a detected ν{sub μ} flux ratio < 1/4 will require the initial flavor composition with more ν{sub e} neutrinos than ν{sub μ} neutrinos.

79 ASTRONOMY AND ASTROPHYSICS↗

A numerical manifold method model for analyzing fully coupled hydro-mechanical processes in porous rock masses with discrete fractures

A numerical manifold method (NMM) model was developed for fully coupled analysis of hydro-mechanical (HM) processes in porous rock masses with discrete fractures. Using an NMM two-cover-mesh system of mathematical and physical covers, fractures are conveniently discretized by dividing the mathematical cover along fracture traces to physical cover, resulting in a discontinuous model on a non-conforming mesh. In this model, discrete fracture deformation (e.g. open and slip) and fracture fluid flow within a permeable and deformable porous rock matrix are rigorously considered. For porous rock, direct pore-volume coupling was modeled based on an energy-work scheme. For mechanical analysis of fractures, a fracture constitutive model for mechanically open states was introduced. For fluid flow in fractures, both along-fracture and normal-to-fracture fluid flow are modeled without introducing additional degrees of freedom. When the mechanical aperture of a fracture is changing, its hydraulic aperture and hydraulic conductivity is updated. At the same time, under the effect of coupled deformation and fluid flow, the contact state may dynamically change, and the corresponding contact constraint is updated each time step. Therefore, indirect coupling is realized under stringent considerations of coupled HM effects and fracture constitutive behavior transfer dynamically. To verify the new model, examples involving deformable porous media containing a single and two sets of fractures were designed, showing good accuracy. Finally, the model was applied to analyze coupled HM behavior of fractured porous rock domains with complex fracture networks under effects of loading and injection.

58 GEOSCIENCES↗

Quantum spin liquid with emergent chiral order in the triangular-lattice Hubbard model

The interplay between spin frustration and charge fluctuation gives rise to an exotic quantum state in the intermediate-interaction regime of the half-filled triangular-lattice Hubbard model, while the nature of the state is under debate. Using the density matrix renormalization group with SU(2) spin ⓍU(1) charge symmetries implemented, we study the triangular-lattice Hubbard model defined on the long cylinder geometry up to circumference W=6. A gapped quantum spin liquid, with on-site interaction 9≲U/t≲10.75, is identified between the metallic and the antiferromagnetic Mott insulating phases. In particular, we find that this spin liquid develops a robust long-range spin scalar-chiral correlation as the system length L increases, which unambiguously unveils the spontaneous time-reversal symmetry breaking. In addition, the degeneracy of the entanglement spectrum supports symmetry fractionalization and spinon edge modes in the obtained ground state. The possible origin of chiral order in this intermediate spin liquid and its relation to the rotonlike excitations have also been discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Where Does Water Go During Hydraulic Fracturing?

Typically, during hydraulic fracturing millions of gallons of water are injected at high pressure into deep shale formations. This water can be housed in fractures, within the shale matrix, and can potentially migrate beyond the shale formation via fractures and/or faults raising environmental concerns. Here, we describe a generic framework for producing estimates of the volume available in fractures and undamaged shale matrix where water injected into a representative shale site could reside during hydraulic fracturing, and apply it to a representative site that incorporates available field data. Moreover, the amount of water that can be stored in the fractures is estimated by calculating the volume of all the fractures associated with a discrete fracture network (DFN) based on real data and using probability theory to estimate the volume of smaller fractures that are below the lower cutoff for the fracture radius in the DFN. The amount of water stored in the matrix is estimated utilizing two distinct methods—one using a two-phase model at the pore-scale and the other using a single-phase model at the continuum scale. Based on these calculations, it appears that most of the water resides in the matrix with a lesser amount in the fractures.

58 GEOSCIENCES↗

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Matrix-assisted relaxation in Fe(phen){sub 2}(NCS){sub 2} spin-crossover microparticles, experimental and theoretical investigations

In this study, we present the influence of the embedding matrix on the relaxation of Fe(phen){sub 2}(NCS){sub 2} (phen = 1,10-phenanthroline) spin-transition microparticles as revealed by experiments and provide an explanation within the framework of an elastic model based on a Monte-Carlo method. Experiments show that the shape of the high-spin → low-spin relaxation curves is drastically changed when the particles are dispersed in glycerol. This effect was considered in the model by means of interactions between the microparticles and the matrix. A faster start of the relaxation for microparticles embedded in glycerol is due to an initial positive local pressure acting on the edge spin-crossover molecules from the matrix side. This local pressure diminishes and eventually becomes negative during relaxation, as an effect of the decrease of the volume of spin-crossover microparticles from high-spin to low-spin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Magnetic properties of alternating Hubbard ladders

We explore the Hubbard Hamiltonian on ladders where the number of sites per rung alternates between two and three. These geometries are bipartite with nonequal or equal number of sites on the two sublattices. Thus they share a key feature of the Hubbard model in a class of lattices which Lieb has shown analytically to exhibit long-range ferrimagnetic order while being amenable to powerful numeric approaches developed for quasi-one-dimensional geometries. The density matrix renormalization group (DMRG) method is used to obtain the groundstate properties, e.g., excitation gaps, charge and spin densities as well as their correlation functions at half filling. We show the existence of long-range ferrimagnetic order in the one-dimensional ladder geometries. Our work provides detailed quantitative results which complement the general theorem of Lieb for generalized bipartite lattices. It also addresses the issue of how the alternation between quasi-long-range order and spin liquid behavior for uniform ladders with odd and even numbers of legs might be affected by a regular alternation pattern.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kaon-nucleon scattering in three-dimensional technique

Kaon-nucleon (KN) scattering is formulated in the three-dimensional (3D) momentum space, in which the basis state is not expanded into partial waves. Based on this basis the Lippmann-Schwinger equation for the T-matrix is evaluated. We obtain as final equation for the T-matrix elements a set of two coupled integral equations in two variables, which are the momentum’s magnitude and the scattering angle. Calculations for the differential cross section and some spin observables are shown, for which we employ a hadrons exchange model with the second order contributions only.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An in vivo dose verification method for SBRT–VMAT delivery using the EPID

Purpose: Radiation treatments have become increasingly more complex with the development of volumetric modulated arc therapy (VMAT) and the use of stereotactic body radiation therapy (SBRT). SBRT involves the delivery of substantially larger doses over fewer fractions than conventional therapy. SBRT–VMAT treatments will strongly benefit from in vivo patient dose verification, as any errors in delivery can be more detrimental to the radiobiology of the patient as compared to conventional therapy. Electronic portal imaging devices (EPIDs) are available on most commercial linear accelerators (Linacs) and their documented use for dosimetry makes them valuable tools for patient dose verification. In this work, the authors customize and validate a physics-based model which utilizes on-treatment EPID images to reconstruct the 3D dose delivered to the patient during SBRT–VMAT delivery. Methods: The SBRT Linac head, including jaws, multileaf collimators, and flattening filter, were modeled using Monte Carlo methods and verified with measured data. The simulation provides energy spectrum data that are used by their “forward” model to then accurately predict fluence generated by a SBRT beam at a plane above the patient. This fluence is then transported through the patient and then the dose to the phosphor layer in the EPID is calculated. Their “inverse” model back-projects the EPID measured focal fluence to a plane upstream of the patient and recombines it with the extra-focal fluence predicted by the forward model. This estimate of total delivered fluence is then forward projected onto the patient’s density matrix and a collapsed cone convolution algorithm calculates the dose delivered to the patient. The model was tested by reconstructing the dose for two prostate, three lung, and two spine SBRT–VMAT treatment fractions delivered to an anthropomorphic phantom. It was further validated against actual patient data for a lung and spine SBRT–VMAT plan. The results were verified with the treatment planning system (TPS) (ECLIPSE AAA) dose calculation. Results: The SBRT–VMAT reconstruction model performed very well when compared to the TPS. A stringent 2%/2 mm χ-comparison calculation gave pass rates better than 91% for the prostate plans, 88% for the lung plans, and 86% for the spine plans for voxels containing 80% or more of the prescribed dose. Patient data were 86% for the lung and 95% for the spine. A 3%/3 mm χ-comparison was also performed and gave pass rates better than 93% for all plan types. Conclusions: The authors have customized and validated a robust, physics-based model that calculates the delivered dose to a patient for SBRT–VMAT delivery using on-treatment EPID images. The accuracy of the results indicates that this approach is suitable for clinical implementation. Future work will incorporate this model into both offline and real-time clinical adaptive radiotherapy.

60 APPLIED LIFE SCIENCES↗

Atomic data for beam-stimulated plasma spectroscopy in fusion plasmas

Injection of high energy atoms into a confined plasma volume is an established diagnostic technique in fusion research. This method strongly depends on the quality of atomic data for charge-exchange recombination spectroscopy (CXRS), motional Stark effect (MSE) and beam-emission spectroscopy (BES). We present some examples of atomic data for CXRS and review the current status of collisional data for parabolic states of hydrogen atoms that are used for accurate MSE modeling. It is shown that the collisional data require knowledge of the excitation density matrix including the off-diagonal matrix elements. The new datasets for transitions between parabolic states are used in an extended collisional-radiative model. The ratios between the {sigma}- and {pi}-components and the beam-emission rate coefficients are calculated in a quasi-steady state approximation. Good agreement with the experimental data from JET is found which points out to strong deviations from the statistical distribution for magnetic sublevels.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Latent-Variable Formulation of the Poisson Canonical Polyadic Tensor Model: Maximum Likelihood Estimation and Fisher Information

We establish parameter inference for the Poisson canonical polyadic (PCP) tensor model through a latent-variable formulation. Our approach exploits the observation that any random PCP tensor can be derived by marginalizing an unobservable random tensor of one dimension larger. The loglikelihood of this larger dimensional tensor, referred to as the “complete” loglikelihood, is comprised of multiple rank one PCP loglikelihoods. Using this methodology, we first derive maximum likelihood estimators for the PCP model and demonstrate that several existing algorithms for fitting non-negative matrix and tensor factorizations are Expectation-Maximization algorithms. Next, we derive the observed and expected Fisher information matrices for the PCP model. The Fisher information provides us crucial insights into the well-posedness of the tensor model, such as the role that tensor rank plays in identifiability and indeterminacy. For the special case of rank one PCP models, we demonstrate that these results are greatly simplified.

97 MATHEMATICS AND COMPUTING↗

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗

Precipitate growth kinetics under inhomogeneous concentration fields using a phase-field model

We investigate precipitation dynamics in the presence of a local solute gradient using phase-field simulations. During the homogenization heat treatment of the solidified Inconel 718 alloy, high Nb concentration within the Laves phases or at the core of the secondary arms results in Nb diffusion into the γ matrix. The volume fraction and spatial distribution of precipitation during subsequent annealing can be controlled by tailoring the Nb concentration gradient in the matrix during homogenization. We use a surrogate Ni-Fe-Nb alloy for Inconel 718 to explore the growth dynamics of δ precipitates related to the local Nb concentration levels. The simulations indicate that in the presence of a Nb concentration gradient the growth rate of δ precipitates is higher than in a matrix of uniform average Nb concentration. The higher growth rate is a result of the higher local thermodynamic driving force at the interface between the solute-rich matrix and the δ interface. Here, we propose a phenomenological model to describe the diffusion-controlled growth kinetics of the δ phase under a solute concentration gradient.

36 MATERIALS SCIENCE↗