Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariant”

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,405 records · Page 78

Regression models using shapes of functions as predictors

Functional variables are often used as predictors in regression problems. A commonly used parametric approach, called scalar-on-function regression, uses the $\mathbb L^2$ inner product to map functional predictors into scalar responses. This method can perform poorly when predictor functions contain undesired phase variability, causing phases to have disproportionately large influence on the response variable. One past solution has been to perform phase–amplitude separation (as a pre-processing step) and then use only the amplitudes in the regression model. In this paper, we propose a more integrated approach, termed elastic functional regression model (EFRM), where phase-separation is performed inside the regression model, rather than as a pre-processing step. This approach generalizes the notion of phase in functional data, and is based on the norm-preserving time warping of predictors. Due to its invariance properties, this representation provides robustness to predictor phase variability and results in improved predictions of the response variable over traditional models. We demonstrate this framework using a number of datasets involving gait signals, NMR data, and stock market prices.

97 MATHEMATICS AND COMPUTING↗

Graviton-photon oscillations in an expanding universe

Through the Gertsenshtein effect, the presence of a large external B-field may allow photons and gravitons to mix in a way that resembles neutrino oscillations and is even more similar to axion-photon mixing. Assuming a background B-field (or any field that behaves like one), we examine the Gertsenshtein mechanism in the context of FLRW expanding universe models, as well as de Sitter space. The conformal invariance of Maxwell’s equations and the conformal noninvariance of the Einstein equations preclude the operation of the Gertsenshtein effect at all scales. In general we find for the matter- and radiation-dominated cases, graviton-oscillations are possible only at late conformal times or when the wavelengths are much shorter than the horizon (kn $\gg$ 1), but that the time-dependent oscillations eventually damp out in any case. Further, the presence of charged particles additionally damps out the oscillations. For the de Sitter universe, we find that oscillations are possible only at early conformal times (n $\ll$ H -1 ) and for wavelengths short compared to the Hubble radius, but eventually freeze in when wavelengths become longer than the Hubble radius. In principle a Gertsenshtein-like mechanism might influence the balance of particle species in an inflationary phase before freezing in; however, we find that in all our models the mixing length is larger than the Hubble radius. We discuss several possible remedies to this situation.

79 ASTRONOMY AND ASTROPHYSICS↗

Sequence2Self: Self-supervised image sequence denoising of pixel-level spray breakup morphology

Optical imaging of fast and transient phenomena such as the turbulent breakup of liquid sprays exhibit low signal-to-noise ratios due to the limited illumination intensity relative to the short exposure time. Image denoising is required to facilitate physical studies over these data but is challenging due to the absence of clean ground-truths and the stringency of the denoising task (e.g., strong and complex noise, limited resolution, preserving physical fidelity), preventing supervised and existing un-/self-supervised deep learning methods. To this end, Sequence2Self (Seq2S) is proposed, an extension of Self2Self (S2S) to image sequences that leverages both the signal’s spatial and temporal correlation. Seq2S is demonstrated on time-resolved x-ray phase contrast imaging of liquid jet fuel sprays in a gas turbine combustor, which possesses all of challenges detailed above. Experiments are conducted across four fuels with different breakup morphology using various state-of-the-art methods. Overall, many of the methods failed and Seq2S was most successful: (1) Accurate spray structures were reconstructed with consistent evolution across frames void of artifacts. (2) The performance was robust, invariant to the hyperparameter choice. (3) Computational time is short and can be made eligible for real-time denoising. In particular, the images denoised by Seq2S showed spray droplet diameter distributions with near-zero Kullback–Leibler divergence (0.01 ± 0.01) to a cleaner reference, whereas the second best method yielded 0.06 ± 0.03. In conclusion, this suggests that Seq2S can be reliably used prior to subsequent quantitative spray analyses as it retains (if not, improves) the statistical physical properties of the data.

97 MATHEMATICS AND COMPUTING↗

New data-driven approach to bridging power system protection gaps with deep learning

Protection is a critical function in power systems to avoid equipment damage, maintain personnel safety, and support system reliability. However, current protective relay technology cannot adequately protect equipment and personnel from effects of some events; these deficiencies are termed protection gaps. In this paper, a data-driven approach is proposed to complement traditional protection technology and distinguish fault conditions from transients caused by normal operations. A combined convolutional neural network and long short-term memory (CNN-LSTM) network is implemented to achieve data translation invariance and capture the temporal correlation of the time-series input data. As a result, the data-driven method can accurately detect system faults despite variation and noise in the input data. In addition, using the CNN-LSTM--based method avoids the complicated, manual feature extraction procedure required by many traditional data-driven methods. The effectiveness of the proposed approach is tested on two kinds of protection gaps: high-impedance faults and transformer inter-turn faults. Lastly, a transfer learning method is also proposed to address the common issue of data-driven methods for which real-world training data are scarce. Extensive study results demonstrate that the proposed approach can accurately bridge power system protection gaps.

42 ENGINEERING↗

Hyperelastic nature of the Hoek–Brown criterion

In this article, we propose a nonlinear elasto-plastic model, for which a specific class of hyperbolic elasticity arises as a straight consequence of the yield criterion invariance on the plasticity level. We superimpose this nonlinear elastic (or hyperelastic) behavior with plasticity obeying the associated flow rule. Interestingly, we find that a linear yield criterion on the thermodynamical force associated with plasticity results in a quadratic yield criterion in the stress space. This suggests a specific hyperelastic connection between Mohr–Coulomb and Hoek–Brown (or alternatively between Drucker–Prager and Pan–Hudson) yield criteria. We compare the elasto-plastic responses of standard tests for the Drucker–Prager yield criterion using either linear or the suggested hyperbolic elasticity. Notably, the nonlinear case stands out due to dilatancy saturation observed during cyclic loading in the triaxial compression test. We conclude this study with structural finite element simulations that clearly demonstrate the numerical applicability of the proposed model.

97 MATHEMATICS AND COMPUTING↗

Verification of nonperturbative guiding center theory in symmetric fields

We verify a recently-developed nonperturbative guiding center formalism to charged particle dynamics in fields with two-parameter continuous symmetry groups. This entails finding exact constants of motion, valid in the nonperturbative regime, that agree with Kruskal’s adiabatic invariant series to all orders in the perturbative regime, when the field scale length is large compared with a typical gyroradius. We demonstrate that the nonperturbative guiding center model makes exact predictions in these cases, even though it eliminates the cyclotron timescale, thereby establishing a theoretical baseline for performance of the nonperturbative formalism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Weighted Composition Operators for Learning Nonlinear Dynamics

Operator theoretic methods in dynamical system have been dominated by the use of Koopman operators and their continuous time counterparts, such as Koopman Generators and Liouville Operators. The advantage gained from their use primarily stems from the ability to extract subspaces and eigenfunctions within a space of observables that are invariant with respect to the Koopman operator over that space. When this occurs, a dynamic mode decomposition of the systems state provides a linear model for the dynamical system. Not all Koopman operators have eigenfunctions that may be exploited in this manner. However, the framework can still be leveraged for approximations using other operators. In this setting, we present a different operator for the study of dynamical systems, the weighted composition operator. These operators are compact for a wide range of dynamics and spaces, and through their interactions with occupation kernels and vector valued kernels, they admit an estimation of the underlying dynamics. Here, this manuscript presents a new algorithm for the data driven study of dynamical systems from data, and also provides two numerical experiments where convergence is achieved as a proof of concept.

97 MATHEMATICS AND COMPUTING↗

Physics augmented machine learning discovery of composition-dependent constitutive laws for 3D printed digital materials

Multi-material 3D printing, particularly through polymer jetting, enables the fabrication of digital materials by mixing distinct photopolymers at the micron scale within a single build to create a composite with tunable mechanical properties. Here, this work presents an integrated experimental and computational investigation into the composition-dependent mechanical behavior of 3D printed digital materials. We experimentally characterize five formulations, combining soft and rigid UV-cured polymers under uniaxial tension and torsion across three strain and twist rates. The results reveal nonlinear and rate-dependent responses that strongly depend on composition. To model this behavior, we develop a physics-augmented neural network (PANN) that combines a partially input convex neural network (pICNN) for learning the composition-dependent hyperelastic strain energy function with a quasi-linear viscoelastic (QLV) formulation for time-dependent response. The pICNN ensures convexity with respect to strain invariants while allowing non-convex dependence on composition. To enhance interpretability, we apply $L_0$ sparsification. For the time-dependent response, we introduce a multilayer perceptron (MLP) to predict viscoelastic relaxation parameters from composition. The proposed model accurately captures the nonlinear, rate-dependent behavior of 3D printed digital materials in both uniaxial tension and torsion, achieving high predictive accuracy for interpolated material compositions. This approach provides a scalable framework for automated, composition-aware constitutive model discovery for multi-material 3D printing.

Constitutive modeling↗

A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS

The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.

Nek5000↗

Additional considerations in analytical solution for time-dependent heat conduction in a three-dimensional multilayer sphere

This work presents an analytical method to solve the heat conduction equation in three dimensions for problems consisting of multilayer concentric spheres. The method can be used to treat time-varying heat conduction problems where the heat source that drives the transient is time-invariant. Equally applicable to all Poisson-type problems with concentric spherical geometry, the method consists of representing the solution as a summation of weighted eigenfunctions. The weights for each eigenfunction are computed algebraically. Previous work has already established the core constituents of the methodology. The current work augments the existing methods by including consideration of nonzero interface resistance between layers and explicit discussion on the boundary condition homogenization required to treat inhomogeneous problems. Also, two demonstration problems are presented. One demonstration problem is based on the method of manufactured solutions and therefore allows for comparison with exact expressions for the solution temperature distribution. The second, more complex, demonstration problem relies on the finite element method for comparisons. The expected convergence behavior is observed for both demonstration problems.

97 - MATHEMATICS AND COMPUTING↗

Symmetry structure of a Riccati equation appearing in penetration mechanics

In the design of projectile penetration experiments a matter of considerable interest is scaling: that is, the potential relevance of small-scale experiments to their full-scale counterparts, in a manner analogous to that most often encountered in the context of fluid mechanics. From the theoretical standpoint, phenomena associated with scaling and scalability can be assessed using the well-established tools of dimensional analysis and the Buckingham-Pi Theorem. However, the familiar precepts of dimensional analysis are themselves a specific manifestation of the broader group invariance properties or symmetries of a mathematical model. Here, this work explores these notions – that is, dimensional analysis, the conditions for realizing complete similarity, and any additional symmetry structures – in the context of a Riccati differential equation appearing in the context of penetration mechanics. The aim of the investigation is twofold: 1) to complement existing empirical considerations with a concrete theoretical basis, and 2) to provide a deeper theoretical understanding of the projectile penetration model and its many implications.

97 MATHEMATICS AND COMPUTING↗

Re-investigation of high-temperature phase equilibria in Fe-rich Sm–Fe–Ti alloys

Recently renewed attempts to develop high-performance rare-earth-lean permanent magnets based on the Sm(Fe,Ti) 12 compound have drawn attention to the limited knowledge about the high-temperature phase equilibria in the Sm–Fe–Ti system. Experimental investigation of equilibrated alloys with electron probe microanalysis, X-ray diffraction and thermomagnetic analysis revealed several inaccuracies in the currently accepted phase relations at 1000°C and allowed for a revision of the Fe-rich corner of the Sm–Fe–Ti phase diagram. The Sm(Fe,Ti) 12 and Sm 3 (Fe,Ti) 29 phases were found to have more extended Ti ranges of 5.1–9.7 at% and 2.8–6.9 at%, respectively. With increasing of the Ti content, the Curie temperature of the Sm(Fe,Ti) 12 remains nearly constant at 306–312°C, whereas that of the Sm 3 (Fe,Ti) 29 increases from 188°C to 207°C. The low-titanium Sm 3 (Fe,Ti) 29 phase equilibrates not only with the Sm(Fe,Ti) 12 and Sm2(Fe,Ti) 17 phases, but also with (α-Fe) solid solution. Newly demonstrated equilibrium between Sm 2 (Fe,Ti) 17 and TiFe 2 phases makes impossible the earlier reported equilibrium between the Sm 3 (Fe,Ti) 29 and Sm(Fe,Ti) 11 phases. Because of an invariant reaction at 1000 °C, the revised phase diagram also features a class II four-phase equilibrium Sm 3 (Fe,Ti) 29 + TiFe 2 + Sm(Fe,Ti) 12 + Sm 2 (Fe,Ti) 17 . Peritectic decomposition of the Sm(Fe,Ti) 11 phase, which occurs either at 1075°C or at 1087°C, was found to have among its products the Sm(Fe,Ti) 12 phase. Although no such equilibration was attempted, it must be possible to obtain above 1087°C a two-phase state composed of the Sm(Fe,Ti) 12 phase and a liquid – which is important for manufacturing of the Sm(Fe,Ti) 12 - based permanent magnets via the liquid-phase sintering.

36 MATERIALS SCIENCE↗

Carbidic Mo is the sole kinetically-relevant active site for catalytic methane dehydroaromatization on Mo/H-ZSM-5

Systematic variation of Mo/Al = {0.10, 0.25, 0.35} in Mo/H-ZSM-5 catalysts effects ~20× change of Mo-to-H + ratio during methane dehydroaromatization (DHA) reactions at steady state and on deactivating catalysts. Mo and Brønsted-acid site count, respectively enumerated by ethane hydrogenolysis probe reaction and dimethyl ether chemical titration, evolve disparately during DHA catalysis, indicating the two sites reduce in number by distinct mechanisms. Here, the deactivation of Mo/H-ZSM-5 catalytic activity is uniquely attributed to the loss of active, or accessible, Mo sites, evinced by (i) the non-selective nature of deactivation (i.e. that deactivation occurs without modification of residual active sites), (ii) invariance in product distribution from Mo/H + = 0.10–2.1, and (iii) linear correspondence between Mo-catalyzed ethane hydrogenolysis rate and active site count during DHA on deactivating samples. DHA forward rate – the intrinsic kinetic descriptor of catalyst activity – normalizes by Mo content at steady state and on deactivating catalysts with Mo/Al = 0.10–0.35 for Mo-based contact times in the range τ Mo = 0.79–44 mol Mo s mol –1 C , unequivocally establishing Mo aggregates as the sole kinetically-relevant active site on Mo/H-ZSM-5. Brønsted-acid site content over a ~40× range in τ H+ = 1.4–54 mol H+ s mol –1 C has no discernable correlation with or effect on DHA rate, methane conversion, or product distribution – demonstrating H + does not catalyze any rate- or selectivity-determining steps in the benzene formation pathway and suggesting that the benefits of zeolitic acid sites are only to disperse Mo during catalyst synthesis and, conceivably, to catalyze equilibrated reaction steps.

02 PETROLEUM↗

Rates and reversibility of CO 2 hydrogenation on Cu-based catalysts

Kinetics of reaction pathways involved in the conversion of CO 2 to methanol and CO on Cu/ZnO/Al 2 O 3 are resolved using in situ chemical titration, steady-state kinetic measurements, and mathematical formalisms for reversibility to probe salient species governing methanol selectivity and yield during CO 2 hydrogenation. Across a range of H 2 :CO 2 = 1:1 to H 2 :CO 2 = 80.5:1, active site density determined from in situ chlorine uptake remained invariant; hence, observed trends in rates can be interpreted as only arising from reaction kinetics and not from changing active site density. Kinetic and thermodynamic contributions to rates are decoupled to evaluate forward and reverse rates of methanol synthesis and reverse water-gas shift (RWGS) reactions. These kinetic analyses show that the forward rates of methanol synthesis exhibit persistent first order dependence on hydrogen pressure and are inhibited by water more significantly than the forward rates of RWGS. In contrast, the reverse rates of methanol synthesis and RWGS are both inhibited by H 2 . Consequently, without any modifications to the Cu/ZnO/Al 2 O 3 catalyst formulation, methanol selectivity can be increased to > 80% by increasing inlet H 2 partial pressure and methanol yield can be enhanced by ~20% by adding water adsorbents even under conditions far from equilibrium. Here, the kinetic treatments presented herein demonstrate a dearth of H* species during catalysis, provide thermodynamic constraints precluding sequential RWGS and CO hydrogenation as the pathway for methanol synthesis, reveal P H2 and P H2O as salient in determining methanol selectivity and yield by impacting both the forward and reverse rates of CO 2 hydrogenation on Cu/ZnO/Al 2 O 3 , and explicate the fundamentals of novel sorption-enhanced methanol synthesis, which not only alleviates equilibrium constraints but also alters the intrinsic rate at which the system approaches equilibrium.

42 ENGINEERING↗

Mechanistic and kinetic relevance of hydrogen and water in CO 2 hydrogenation on Cu-based catalysts

Here, we ally steady-state kinetics, kinetic isotope effects, and density functional theory (DFT) calculations to illustrate that Cu-based catalysts remain saturated by H-adatoms (H*) and molecular formic acid (HCOOH**) during CO 2 hydrogenation. High H* coverage under methanol synthesis conditions is evidenced by reverse water-gas shift (RWGS) rates that exhibit positive H 2 reaction orders only at P H2 ≲ 0.5 bar, above which methanol synthesis and RWGS rates exhibit first and zeroth order dependence on P H2 , respectively. HCOOH** also accumulates on the surface with increasing P CO2 as informed by the Langmuir-type dependence on P CO2 (0.25-23 bar) for both methanol synthesis and RWGS. As both HCOOH** and H* have one H-atom per site occupied, the two species share the same P H2 dependence and give rise to CO 2 reaction orders that are independent of P H2 . Surface coverages determined based on kinetic analyses are further corroborated with DFT-derived adsorption energies that show favorable HCOOH** adsorbate-adsorbate interactions as well as repulsive interactions for bidentate formate (HCOO**) on H*-saturated surfaces. Methanol selectivity remains invariant with P CO2 and P CO despite CO inhibiting reaction rates, thereby demonstrating methanol synthesis and RWGS occur on the same active site. In contrast, water preferentially inhibits methanol synthesis rates, increases methanol synthesis H 2 reaction order from 1.0 to 1.5, and alters the methanol synthesis H 2 /D 2 kinetic isotope effect; the inhibitory effect of H 2 O thus cannot be attributed to competitive adsorption alone and instead reflects a change in the rate-determining step for methanol synthesis. The disparate kinetics of methanol synthesis and RWGS evince a branching pathway where methanol is formed from formates and CO is formed from carboxylates. The presented work thus identifies the relevant surface species, underscores the distinct catalytic role of water in branching methanol synthesis and RWGS pathways, and, in doing so, details a mechanistic picture that yields predictable rates and reaction orders for both methanol synthesis and RWGS on Cu-based CO 2 hydrogenation catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy conservative SBP discretizations of the acoustic wave equation in covariant form on staggered curvilinear grids

In this work, we develop a numerical method for solving the acoustic wave equation in covariant form on staggered curvilinear grids in an energy conserving manner. The use of a covariant basis decomposition leads to a rotationally invariant scheme that outperforms a Cartesian basis decomposition on rotated grids. The discretization is based on high order Summation-By-Parts (SBP) operators and preserves both symmetry and positive definiteness of the contravariant metric tensor. To improve accuracy and decrease computational cost, we also derive a modified discretization of the metric tensor that leads to a conditionally stable discretization. Bounds are derived that yield a point-wise condition that can be evaluated to check for stability of the modified discretization. This condition shows that the interpolation operators should be constructed such that their norm is close to one.

97 MATHEMATICS AND COMPUTING↗

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗