SEARCH · Engineering Papers

Results for “multiscale problem”

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

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Hydrodynamics of countercurrent flows in a structured packed column: effects of initial wetting and dynamic contact angle

The structured packed column has gained considerable attention as a prospective technology to mitigate carbon dioxide emissions from thermal power plants. The computational flow investigation of such systems is a multiscale problem, and local hydrodynamics plays a key role in overall column efficiency. Accordingly, multiphase flow studies in the representative elementary unit (REU) of a packed column can offer insights into the local flow behavior such as flow pattern, interfacial area, etc. The interfacial area plays a key role in the overall interphase mass transfer and thereby in separation efficiency. The predicted effective area matches well with corresponding ones in experiments for Mellapak 250.Y packings. The contact angle (i.e., solid substrate behavior) is one of the critical factors that dictates wettings, and thereby the interfacial area. The dynamic contact angle (DCA) was considered in order to explore the effects of contact angle hysteresis on the interfacial area. DCA has more pronounced impact on the interfacial area for a solvent possessing a higher value of surface tension and equilibrium contact angle. The interfacial area shows temporal undulation and does not achieve a pseudo steady state due to the capillary force. In contrary, the interfacial area gets a net value for a solvent having low surface tension value. Further, specification of the static contact angle at the solid substrate shows the least value of the interfacial area. We also extensively studied the effect of the initial sheet condition (dry vs wet) on the interfacial area. The initially wetted sheets show slightly higher value of the interfacial area as compared to the initially dry sheets at a fixed flow rate for a solvent having high surface tension value. Wetting hysteresis decreases with increasing liquid loads. On the other hand, wetting hysteresis does not occur for a solvent having value of low surface tension and equilibrium contact angle.

structured packings, dynamic contact angle, interf↗

DOE OSTI2018

Residence time distribution in a structured packing unit for monitoring aerosol emissions

Solvent absorption in a structured packed column has been demonstrated as a viable option for post-combustion carbon capture. Amine emissions in the form of aerosol particles pose a major challenge that severely impacts the separation efficiency of the packed column. The aerosol particles quickly achieve the motion of the gas flow. In this context, residence time distribution (RTD) of the gaseous phase will provide insights for managing aerosol particles, thereby controlling amine loss in the packed column. Computational fluid dynamics (CFD) modeling of such columns is a multiscale problem because of a wide range of disparity in length scales. Accordingly, three-dimensional CFD simulations in the representative elementary unit of a structured packing are conducted using OpenFOAM. The RTD is analyzed in terms of residence time distribution function, mean residence time, and variance. Effects of the diffusivity and flow rates on RTD are extensively studied, and the dispersion coefficient is calculated using a one-dimensional axial dispersion model. Furthermore, RTDs computed in a smaller domain are extrapolated to a larger domain using a convolution integral. The computed and corresponding CFD models derived from the convolution integral match well at different flow rates and domain sizes. A comparison of RTD among different packing designs also is presented at a fixed number of repeating units and different flow rates.

structured packings, residence time distribution↗

DOE OSTI2015

A coupled ordinates method for solution acceleration of rarefied gas dynamics simulations

Non-equilibrium rarefied flows are frequently encountered in a wide range of applications, including atmospheric re-entry vehicles, vacuum technology, and microscale devices. Rarefied flows at the microscale can be effectively modeled using the ellipsoidal statistical Bhatnagar–Gross–Krook (ESBGK) form of the Boltzmann kinetic equation. Numerical solutions of these equations are often based on the finite volume method (FVM) in physical space and the discrete ordinates method in velocity space. However, existing solvers use a sequential solution procedure wherein the velocity distribution functions are implicitly coupled in physical space, but are solved sequentially in velocity space. This leads to explicit coupling of the distribution function values in velocity space and slows down convergence in systems with low Knudsen numbers. Furthermore, this also makes it difficult to solve multiscale problems or problems in which there is a large range of Knudsen numbers. In this paper, we extend the coupled ordinates method (COMET), previously developed to study participating radiative heat transfer, to solve the ESBGK equations. In this method, at each cell in the physical domain, distribution function values for all velocity ordinates are solved simultaneously. This coupled solution is used as a relaxation sweep in a geometric multigrid method in the spatial domain. Enhancements to COMET to account for the non-linearity of the ESBGK equations, as well as the coupled implementation of boundary conditions, are presented. The methodology works well with arbitrary convex polyhedral meshes, and is shown to give significantly faster solutions than the conventional sequential solution procedure. Acceleration factors of 5–9 are obtained for low to moderate Knudsen numbers on single processor platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling

Neoclassical tearing modes are macroscopic ($L$ ~ 1 m) instabilities in magnetic fusion experiments; if unchecked, these modes degrade plasma performance and may catastrophically destroy plasma confinement by inducing a disruption. Fortunately, the use of properly tuned and directed radiofrequency waves ($λ$ ~ 1 mm) can eliminate these modes. Numerical modeling of this difficult multiscale problem requires the integration of separate mathematical models for each length and time scale [T. G. Jenkins, S. E. Kruger, Phys. Plasmas 19 (2012) 122508]; the extended MHD model captures macroscopic plasma evolution while the RF model tracks the flow and deposition of injected RF power through the evolving plasma profiles. The scale separation enables use of the eikonal (ray-tracing) approximation to model the RF wave propagation. In this work we demonstrate a technique, based on methods of computational geometry, for mapping the ensuing RF data (associated with discrete ray trajectories) onto the finite-element/pseudospectral grid that is used to model the extended MHD physics. In the new representation, the RF data can then be used to construct source terms in the equations of the extended MHD model, enabling quantitative modeling of RF-induced tearing mode stabilization. Though our specific implementation uses the NIMROD extended MHD [C. R. Sovinec et al., J. Comp. Phys. 195 (2004) 355] and GENRAY RF [A. P. Smirnov, R. W. Harvey, K. Kupfer, Bull. Amer. Phys. Soc. 39 (1994) 1626] codes, the approach presented can be applied more generally to any code coupling requiring the mapping of ray tracing data onto Eulerian grids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Modeling Nanoconfinement Effects Using Active Learning

Predicting the spatial configuration of gas in nanopores of is relevant in applications such as fluid flow forecasting and hydrocarbon reserves estimation. For example, shale reservoirs have suffered from computationally intractable multiscale problems, since fluid properties such as viscosity, density, and adsorption must be calculated by using expensive molecular dynamics (MD) simulations within each nanopore, whereas flow through these connected nanopores must be simulated at the micrometer scale. We utilize machine learning techniques to quickly and accurately model nanoscale confinement effects as an important step toward bridging the nano and micro scales. Our workflow is based on building and training physics-based deep-neural-networks models by learning from a database of MD calculations. The model accounts for the adsorption phenomenon by predicting the statistical distribution of gas inside nanopores. Because large databases of MD calculations are expensive to create, we investigate active learning (AL) as a data set construction strategy. In this workflow, new data are selected based on the model uncertainty via the query-by-committee approach. We show that our workflow obtains accurate models that generalize to real scanning electron microscopy geometries with 1/10th of the number of MD calculations required vs random data set generation. Our method enables the possibility of modeling nanoconfinement effects at the mesoscale, where complex connected sets of nanopores affect flow.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Structure and dynamics of H{sub 2}O vis-á-vis phenylalanine recognition at a DPPC lipid membrane via interfacial H-bond types: Insights from polarized FT-IRRAS and ADMP simulations

Preferential and enantioselective interactions of L-/D-Phenylalanine (L-Phe and D-Phe) and butoxycarbonyl-protected L-/D-Phenylalanine (LPA and DPA) as guest with 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (L-DPPC) as host were tapped by using real time Fourier transform infrared reflection absorption spectroscopy (FT-IRRAS). Polarization-modulated FT-IRRAS of DPPC monolayers above the phenylalanine modified subphases depicted fine structure/conformation differences under considerations of controlled 2D surface pressure. Selective molecular recognition of D-enantiomer over L-enantiomer driven by the DPPC head group via H-bonding and electrostatic interactions was evident spectroscopically. Accordingly, binding constants (K) of 145, 346, 28, and 56 M{sup −1} for LPA, DPA, L-Phe, and D-Phe, respectively, were estimated. The real time FT-IRRAS water bands were strictly conformation sensitive. The effect of micro-solvation on the structure and stability of the 1:1 diastereomeric L-lipid⋯, LPA/DPA and L-lipid⋯, (L/D)-Phe adducts was investigated with the aid of Atom-centered Density Matrix Propagation (ADMP), a first principle quantum mechanical molecular dynamics approach. The phosphodiester fragment was the primary site of hydration where specific solvent interactions were simulated through single- and triple- “water-phosphate” interactions, as water cluster’s “tetrahedral dice” to a “trimeric motif” transformation as a partial de-clusterization was evident. Under all the hydration patterns considered in both static and dynamic descriptions of density functional theory, L-lipid/D-amino acid enantiomer adducts continued to be stable structures while in dynamic systems, water rearranged without getting “squeezed-out” in the process of recognition. In spite of the challenging computational realm of this multiscale problem, the ADMP simulated molecular interactions complying with polarized vibrational spectroscopy unraveled a novel route to chiral recognition and interfacial water structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

A high-fidelity approach towards simulation of pool boiling

A novel numerical approach is developed to simulate the multiscale problem of pool-boiling phase change. The particular focus is to develop a simulation technique that is capable of predicting the heat transfer and hydrodynamic characteristics of nucleate boiling and the transition to critical heat flux on surfaces of arbitrary shape and roughness distribution addressing a critical need to design enhanced boiling heat transfer surfaces. The macro-scale of the phase change and bubble dynamics is addressed through employing off-the-shelf Computational Fluid Dynamics (CFD) methods for interface tracking and interphase mass and energy transfer. The micro-scale of the microlayer, which forms at early stage of bubble nucleation near the wall, is resolved through asymptotic approximation of the thin-film theory which provides a closed-form solution for the distribution of the micro-layer and its influence on the evaporation process. In addition, the sub-grid surface roughness is represented stochastically through probabilistic density functions and its role in bubble nucleation and growth is then represented based on the thermodynamics of nucleation process. This combination of deterministic CFD, local approximation, and stochastic representation allows the simulation of pool boiling on any surface with known roughness and enhancement characteristics. The numerical model is validated for dynamics and hydrothermal characteristics of a single nucleated bubble on a flat surface against available literature data. In addition, the prediction of pool-boiling heat transfer coefficient is verified against experimental measurements as well as reputable correlations for various roughness distributions and different surface orientations. Finally, the model is employed to demonstrate pool-boiling phenomenon on enhanced structures with reentrance cavities and to explore the effect of enhancement feature design on thermal and hydrodynamic characteristics of these surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Turbulent mixing and transition criteria of flows induced by hydrodynamic instabilities

In diverse areas of science and technology, including inertial confinement fusion (ICF), astrophysics, geophysics, and engineering processes, turbulent mixing induced by hydrodynamic instabilities is of scientific interest as well as practical significance. Because of the fundamental roles they often play in ICF and other applications, three classes of hydrodynamic instability-induced turbulent flows—those arising from the Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities—have attracted much attention. ICF implosions, supernova explosions, and other applications illustrate that these phases of instability growth do not occur in isolation, but instead are connected so that growth in one phase feeds through to initiate growth in a later phase. Essentially, a description of these flows must encompass both the temporal and spatial evolution of the flows from their inception. Hydrodynamic instability will usually start from potentially infinitesimal spatial perturbations, will eventually transition to a turbulent flow, and then will reach a final state of a true multiscale problem. Indeed, this change in the spatial scales can be vast, with hydrodynamic instability evolving from just a few microns to thousands of kilometers in geophysical or astrophysical problems. These instabilities will evolve through different stages before transitioning to turbulence, experiencing linear, weakly, and highly nonlinear states. Here, the challenges confronted by researchers are enormous. The inherent difficulties include characterizing the initial conditions of such flows and accurately predicting the transitional flows. Of course, fully developed turbulence, a focus of many studies because of its major impact on the mixing process, is a notoriously difficult problem in its own right. In this pedagogical review, we will survey challenges and progress, and also discuss outstanding issues and future directions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Adjoint-based deviational Monte Carlo methods for phonon transport calculations

In the field of linear transport, adjoint formulations exploit linearity to derive powerful reciprocity relations between a variety of quantities of interest. We develop an adjoint formulation of the linearized Boltzmann transport equation for phonon transport. We use this formulation for accelerating deviational Monte Carlo simulations of complex, multiscale problems. Benefits include significant computational savings via direct variance reduction, or by enabling formulations which allow more efficient use of computational resources, such as formulations which provide high resolution in a particular phase-space dimension (e.g., spectral). We show that the proposed adjoint-based methods are particularly well suited to problems involving a wide range of length scales (e.g., nanometers to hundreds of microns) and lead to computational methods that can calculate quantities of interest with a cost that is independent of the system characteristic length scale, thus removing the traditional stiffness of kinetic descriptions. Applications to problems of current interest, such as simulation of transient thermoreflectance experiments or spectrally resolved calculation of the effective thermal conductivity of nanostructured materials, are presented and discussed in detail.

42 ENGINEERING↗

DOE OSTI2023

Multiscale evolution of charmed particles in a nuclear medium

Parton energy-momentum exchange with the quark gluon plasma (QGP) is a multiscale problem. In this work, we calculate the interaction of charm quarks with the QGP within the higher twist formalism at high virtuality and high energy using the Modular All Twist Transverse-scattering Elastic-drag and Radiation (MATTER) model, while the low-virtuality and high-energy portion is treated via a linearized Boltzmann transport formalism. Coherence effect that reduces the medium-induced emission rate in the MATTER model is also taken into account through a virtuality-dependent qˆ, leaving the simultaneous dependence of qˆ on heavy quark mass and virtuality for future studies. The interplay between these two formalisms is studied phenomenologically and used to produce a first description of the D-meson and charged hadron nuclear modification factor R AA across multiple centralities. As a result, all calculations were carried out utilizing the Jet Energy-loss Tomography with a Statistically and Computationally Advanced Program Envelope framework.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1994

Dynamic modeling of orographically induced precipitation

Local orography governs the triggering of cloud formation and the enhancement of processes such as condensation and hydrometeor nucleation and growth in mountainous regions. Intense, lengthy precipitation events are typical upwind of the topographic divide, with sharply decreasing magnitude and duration on the lee side. Differences in mean annual precipitation of several hundred percent between windward slopes of orographic barriers and adjacent valleys or lee side slopes are not unusual. Because much of the streamflow in areas such as the western United States is derived from mountainous areas that are remote and often poorly instrumented, modeling of orographic precipitation has important implications for water resources management. Models of orographically induced precipitation differ by their treatment of atmospheric dynamics and by the extent to which they rely on bulk parameterization of cloud and precipitation physics. Adiabatic ascent and a direct proportionality between efficiency and orographically magnified updrafts are the most frequent assumptions in orographic precipitation modeling. Space-time discretization (i.e., resolution) is a major issue because of the high spatial variability of orographic precipitation. For a specific storm, relative errors as large as 50 to 100% are common in the forecast/hindcast of precipitation intensity and can be even larger in the case of catastrophic storms. When monthly or seasonal timescales are used to evaluate model performance, the magnitude of such errors decreases dramatically, reaching values as low as 10 to 15%. Current research is focusing on the development of data assimilation techniques to incorporate radar and satellite observations, and on the development of aggregation and disaggregation methodologies to address the implications of modeling a multiscale problem at restricted spatial and temporal resolutions.

Barros, Ana Paula↗

NASA NTRSDate not supplied

Application of a Rapid Design Tool to a 3D Woven Structural Joint

The optimization of composite structural joints is an iterative process and a multiscale problem. High fidelity finite element modeling of joints with 3D woven and laminated materials can become computationally expensive. The aim of this paper is to establish a reliable analysis process for the optimization of a composite Y-joint (curved Pi-joint), to be used in an aircraft fuselage, using commercial rapid joint design, analysis, and optimization software. The rapid joint design tool was investigated as a substitute and/or complement to a finite element analysis software. A composite Pi-joint and a composite Y-joint were evaluated using the rapid joint design tool to determine the applicability and limits of the software. Furthermore, three main preliminary parametric studies were performed to better understand the capability of the tool in predicting the stress distributions and trends in the Y-joint. The parameters investigated were the joint curvature, the laminated skin thickness, the adhesive systems, and the ply composition. Lastly, trends in predicted failure load were produced as a function of skin thickness (16ply, 24ply and 32ply). Failure loads were found with varying joint curvature and skin thickness using stress-based adherend failure criterion. The rapid joint design tool was also validated against existing experimental results.

Bonded Joints↗

DOE OSTI2024

A simplified all-frequency stable formulation with an implicit Coulomb gauge

A potential-based finite element formulation has been developed in the past to circumvent the low-frequency breakdown issues commonly encountered in electromagnetic (EM) simulations of low-frequency and multiscale problems. In this formulation, the magnetic vector and electric scalar potentials have been employed to express the electric field and magnetic flux, leading to a set of two equations that represents all four Maxwell’s equations and the current continuity equation. To enforce the Coulomb gauge, an auxiliary potential has been introduced, which results in a total of three equations to be solved simultaneously. To reduce the number of equations and unknowns needed in a simulation, a simplified formulation is proposed in this paper to enforce the Coulomb gauge implicitly without the need for the auxiliary potential. A numerical example is given to demonstrate the accuracy of the proposed formulation.

Mekonnen, Minyechil↗

DOE OSTI2024

PhILMs: Collaboratory on Mathematics and Physics-Informed Learning Machines for Multiscale and Multiphysics Problems

The landscape of computational science and engineering is continually evolving, with the challenge of high-dimensional regression problems standing as a significant hurdle in numerous scientific endeavors. Addressing this challenge, our research, funded by this award, has led to the development of an innovative computational framework known as Probabilistic Partition of Unity Networks (PPOU-Nets). This initiative represents a collaborative effort to harness the potential of mathematics and physics-informed machine learning in tackling multiscale and multiphysics problems prevalent in high-dimensional spaces. Through this work, we have proposed a novel methodology that seamlessly integrates adaptive dimensionality reduction and a mixture of experts model, thereby facilitating a more efficient and accurate approximation of complex functions. This research effort has not only advanced the state of computational science but also opened new avenues for exploration in quantum computing and beyond. This report outlines the motivation, methodology, key findings, and implications of our work, underscoring our contributions to the broader scientific community and the potential pathways for future research.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Physics-Informed Learning Machines for Multiscale and Multiphysics Problems (PHILMS) (Technical Report)

The research work at University of California Santa Barbara (UCSB) resulted in several new developments in the areas of scientific machine learning, numerical analysis, and practical methods for data-driven modeling, prediction, reductions, and simulation. Many of the projects were carried out in collaboration with members of the national laboratories at Sandia National Laboratories (SNL), Pacific Northwestern National Laboratories (PNNL), and other institutions. Results included developing new scientific machine learning methods, related theory and mathematical frameworks for analysis and training, data-driven numerical solvers, and related tools and software for scientific computation. During the support period, over 16+ papers were submitted for publication, and 4 open-source software packages were developed and released (available at http://atzberger.org/). In addition, 7+ students and 2 post-docs were mentored in collaboration with the laboratory staff for future careers in academia, government labs, and industry.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Multiscale Mesh Adaptation for Transonic Aeroelastic Flutter Problems

This work applies multiscale mesh adaptation with refine to reduce spatial discretization error of aeroelastic computational fluid dynamics (CFD) simulations. Benchmark flutter models, such as the pitch and plunge NACA64A-010 airfoil and the benchmark supercritical wing, are studied with both a linearized frequency-domain solver and time-marching CFD coupled to a modal structural solver in FUN3D. The undeformed NASA Common Research Model (CRM), an aeroelastic jig shape variant of the CRM, is also studied with the linearized frequency-domain approach. For these cases, the adaptation process converges to comparable flutter predictions to hand-generated meshes but with smaller node counts. However the additional disciplines of the linearized frequency-domain analysis, the mesh deformation, and the unsteady finite-volume solver create robustness challenges that need to be addressed before it can be applied as a fully automated process for complex transonic aeroelastic problems. In particular, negative volumes are observed to be an issue for FUN3D’s linear elasticity mesh deformation solver when moving the adapted meshes.

Aeroelasticity↗

DOE OSTI2023

An integrated coupling model for solving multiscale fluid-fluid coupling problems in SAM code

In this study, an integrated coupling method has been developed for solving multiscale fluid-fluid coupling problems in plant-scale safety analysis models in SAM (System Analysis Module) code. In this method, a higher-fidelity multi-dimensional (3D) flow module is used for reactor components of complex flow features (e.g., reactor core) and a lumped parameter one-dimensional (1D) flow module for plant-scale flow loops (e.g., primary loop pipe network), respectively. In this method, the 3D fluid equation/domain and 1D fluid equation/domain are tightly coupled at the residual level and solved simultaneously using the Newton’s method to overcome the convergence issues typically seen in existing approaches like separate domain approach, where the 3D fluid equation and 1D fluid equation are solved separately. Extensive and successful code verifications and demonstrations have been performed for this newly developed method. This new modeling approach significantly simplify the work flow in developing high-fidelity plant-scale safety analysis model, e.g. for pool-type reactors and pebble-bed reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗