Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LID”

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 289 records · Page 16

A Semi-Implicit Lattice Method for Simulating Flow

We propose a new semi-implicit lattice numerical method for modeling fluid flow that depends only on local primitive variable information (density, pressure, velocity) and not on relaxed upstream distribution function values. This method has the potential for reducing parallel communication and permitting larger time steps compared to the lattice Boltzmann method. The lid-driven cavity is modeled to demonstrate the accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A spectral-element discontinuous Galerkin lattice Boltzmann method for incompressible flows.

We present a spectral-element discontinuous Galerkin lattice Boltzmann method for solving nearly incompressible flows. Decoupling the collision step from the streaming step offers numerical stability at high Reynolds numbers. In the streaming step, we employ high-order spectral-element discontinuous Galerkin discretizations using a tensor product basis of one-dimensional Lagrange interpolation polynomials based on Gauss-Lobatto-Legendre grids. Our scheme is cost-effective with a fully diagonal mass matrix, advancing time integration with the fourth-order Runge-Kutta method. We present a consistent treatment for imposing boundary conditions with a numerical flux in the discontinuous Galerkin approach. We show convergence studies for Couette flows and demonstrate two benchmark cases with lid-driven cavity flows for Re = 400-5000 and flows around an impulsively started cylinder for Re = 550-9500. Computational results are compared with those of other theoretical and computational work that used a multigrid method, a vortex method, and a spectral element model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical schemes for dynamically orthogonal equations of stochastic fluid and ocean flows

The quantification of uncertainties is critical when systems are nonlinear and have uncertain terms in their governing equations or are constrained by limited knowledge of initial and boundary conditions. Such situations are common in multiscale, intermittent and non-homogeneous fluid and ocean flows. The dynamically orthogonal (DO) field equations provide an adaptive methodology to predict the probability density functions of such flows. The present work derives efficient computational schemes for the DO methodology applied to unsteady stochastic Navier-Stokes and Boussinesq equations, and illustrates and studies the numerical aspects of these schemes. Semi-implicit projection methods are developed for the mean and for the DO modes, and time-marching schemes of first to fourth order are used for the stochastic coefficients. Conservative second-order finite-volumes are employed in physical space with new advection schemes based on total variation diminishing methods. Other results include: (i) the definition of pseudo-stochastic pressures to obtain a number of pressure equations that is linear in the subspace size instead of quadratic; (ii) symmetric advection schemes for the stochastic velocities; (iii) the use of generalized inversion to deal with singular subspace covariances or deterministic modes; and (iv) schemes to maintain orthonormal modes at the numerical level. To verify our implementation and study the properties of our schemes and their variations, a set of stochastic flow benchmarks are defined including asymmetric Dirac and symmetric lock-exchange flows, lid-driven cavity flows, and flows past objects in a confined channel. Different Reynolds number and Grashof number regimes are employed to illustrate robustness. Optimal convergence under both time and space refinements is shown as well as the convergence of the probability density functions with the number of stochastic realizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An efficient particle Fokker–Planck algorithm for rarefied gas flows

This paper is devoted to the algorithmic improvement and careful analysis of the Fokker–Planck kinetic model derived by Jenny et al. [1] and Gorji et al. [2]. The motivation behind the Fokker–Planck based particle methods is to gain efficiency in low Knudsen rarefied gas flow simulations, where conventional direct simulation Monte Carlo (DSMC) becomes expensive. This can be achieved due to the fact that the resulting model equations are continuous stochastic differential equations in velocity space. Accordingly, the computational particles evolve along independent stochastic paths and thus no collision needs to be calculated. Therefore the computational cost of the solution algorithm becomes independent of the Knudsen number. In the present study, different computational improvements were persuaded in order to augment the method, including an accurate time integration scheme, local time stepping and noise reduction. For assessment of the performance, gas flow around a cylinder and lid driven cavity flow were studied. Convergence rates, accuracy and computational costs were compared with respect to DSMC for a range of Knudsen numbers (from hydrodynamic regime up to above one). In all the considered cases, the model together with the proposed scheme give rise to very efficient yet accurate solution algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Variance reduction for Fokker–Planck based particle Monte Carlo schemes

Recently, Fokker–Planck based particle Monte Carlo schemes have been proposed and evaluated for simulations of rarefied gas flows [1–3]. In this paper, the variance reduction for particle Monte Carlo simulations based on the Fokker–Planck model is considered. First, deviational based schemes were derived and reviewed, and it is shown that these deviational methods are not appropriate for practical Fokker–Planck based rarefied gas flow simulations. This is due to the fact that the deviational schemes considered in this study lead either to instabilities in the case of two-weight methods or to large statistical errors if the direct sampling method is applied. Motivated by this conclusion, we developed a novel scheme based on correlated stochastic processes. The main idea here is to synthesize an additional stochastic process with a known solution, which is simultaneously solved together with the main one. By correlating the two processes, the statistical errors can dramatically be reduced; especially for low Mach numbers. To assess the methods, homogeneous relaxation, planar Couette and lid-driven cavity flows were considered. For these test cases, it could be demonstrated that variance reduction based on parallel processes is very robust and effective.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Discontinuous Petrov-Galerkin Methodology for Adaptive Solutions to the Incompressible Navier-Stokes Equations

The discontinuous Petrov-Galerkin methodology with optimal test functions (DPG) of Demkowicz and Gopalakrishnan [18, 20] guarantees the optimality of the solution in an energy norm, and provides several features facilitating adaptive schemes. Whereas Bubnov-Galerkin methods use identical trial and test spaces, Petrov-Galerkin methods allow these function spaces to differ. In DPG, test functions are computed on the fly and are chosen to realize the supremum in the inf-sup condition; the method is equivalent to a minimum residual method. For well-posed problems with sufficiently regular solutions, DPG can be shown to converge at optimal rates—the inf-sup constants governing the convergence are mesh-independent, and of the same order as those governing the continuous problem [48]. DPG also provides an accurate mechanism for measuring the error, and this can be used to drive adaptive mesh refinements. We employ DPG to solve the steady incompressible Navier-Stokes equations in two dimensions, building on previous work on the Stokes equations, and focusing particularly on the usefulness of the approach for automatic adaptivity starting from a coarse mesh. We apply our approach to a manufactured solution due to Kovasznay as well as the lid-driven cavity flow, backward-facing step, and flow past a cylinder problems.

adaptive finite elements↗

Importance sampling variance reduction for the Fokker–Planck rarefied gas particle method

The Fokker–Planck approximation to the Boltzmann equation, solved numerically by stochastic particle schemes, is used to provide estimates for rarefied gas flows. This paper presents a variance reduction technique for a stochastic particle method that is able to greatly reduce the uncertainty of the estimated flow fields when the characteristic speed of the flow is small in comparison to the thermal velocity of the gas. The method relies on importance sampling, requiring minimal changes to the basic stochastic particle scheme. We test the importance sampling scheme on a homogeneous relaxation, planar Couette flow and a lid-driven-cavity flow, and find that our method is able to greatly reduce the noise of estimated quantities. Significantly, we find that as the characteristic speed of the flow decreases, the variance of the noisy estimators becomes independent of the characteristic speed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A general CFD framework for fault-resilient simulations based on multi-resolution information fusion

We develop a general CFD framework for multi-resolution simulations to target multiscale problems but also resilience in exascale simulations, where faulty processors may lead to gappy, in space-time, simulated fields. We combine approximation theory and domain decomposition together with statistical learning techniques, e.g. coKriging, to estimate boundary conditions and minimize communications by performing independent parallel runs. To demonstrate this new simulation approach, we consider two benchmark problems. First, we solve the heat equation (a) on a small number of spatial “patches” distributed across the domain, simulated by finite differences at fine resolution and (b) on the entire domain simulated at very low resolution, thus fusing multi-resolution models to obtain the final answer. Second, we simulate the flow in a lid-driven cavity in an analogous fashion, by fusing finite difference solutions obtained with fine and low resolution assuming gappy data sets. We investigate the influence of various parameters for this framework, including the correlation kernel, the size of a buffer employed in estimating boundary conditions, the coarseness of the resolution of auxiliary data, and the communication frequency across different patches in fusing the information at different resolution levels. In addition to its robustness and resilience, the new framework can be employed to generalize previous multiscale approaches involving heterogeneous discretizations or even fundamentally different flow descriptions, e.g. in continuum-atomistic simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preconditioning a Newton-Krylov solver for all-speed melt pool flow physics

n this paper, we introduce a multigrid block-based preconditioner for solving linear systems arising from a Discontinuous Galerkin discretization of the all-speed Navier-Stokes equations with phase change. The equations are discretized in conservative form with a reconstructed Discontinuous Galerkin (rDG) method and integrated with fully-implicit time discretization schemes. To robustly converge the numerically stiff systems, we use the Newton-Krylov framework with a primitive-variable formulation (pressure, velocity, and temperature), which is better conditioned than the conservative-variable form at low-Mach number. In the limit of large acoustic CFL number and viscous Fourier number, there is a strong coupling between the velocity-pressure system and the linear systems become non-diagonally dominant. To effectively solve these ill-conditioned systems, an approximate block factorization preconditioner is developed, which uses the Schur complement to reduce a 3 x 3 block system into a sequence of two 2 x 2 block systems: velocity-pressure,v P , and velocity-temperature, v T . We compare the performance of the v P -v T Schur complement preconditioner to classic preconditioning strategies: monolithic algebraic multigrid (AMG), element-block SOR, and primitive variable block Gauss-Seidel. The performance of the preconditioned solver is investigated in the limit of large CFL and Fourier numbers for low-Mach lid-driven cavity flow, Rayleigh-Bénard melt convection, compressible internally heated convection, and 3D laser-induced melt pool flow. Here, numerical results demonstrate that the v P -v T Schur complement preconditioned solver scales well both algorithmically and in parallel, and is robust for highly ill-conditioned systems, for all tested rDG discretization schemes (up to 4th-order).

97 MATHEMATICS AND COMPUTING↗

A fully-integrated lattice Boltzmann method for fluid–structure interaction

Here we present a fully-integrated lattice Boltzmann (LB) method for fluid–structure interaction (FSI) simulations that efficiently models deformable solids in complex suspensions and active systems. Our Eulerian method (LBRMT) couples finite-strain solids to the LB fluid on the same fixed computational grid with the reference map technique (RMT). An integral part of the LBRMT is a new LB boundary condition for moving deformable interfaces across different densities. With this fully Eulerian solid–fluid coupling, the LBRMT is well-suited for parallelization and simulating multi-body contact without remeshing or extra meshes. We validate its accuracy via a benchmark of a deformable solid in a lid-driven cavity, then showcase its versatility through examples of soft solids rotating and settling. The LBRMT achieves a spatial convergence rate between first-order and second-order for FSI simulations and is designed for low to intermediate Reynolds number flows with finite inertia at small Mach numbers. With simulations of complex suspensions mixing, we highlight the potential of the LBRMT for studying collective behavior in soft matter and biofluid dynamics.

97 MATHEMATICS AND COMPUTING↗

The Mechanism and Function of Group II Chaperonins

We report protein folding in the cell requires the assistance of enzymes collectively called chaperones. Among these, the chaperonins are 1-MDa ring-shaped oligomeric complexes that bind unfolded polypeptides and promote their folding within an isolated chamber in an ATP-dependent manner. Group II chaperonins, found in archaea and eukaryotes, contain a built-in lid that opens and closes over the central chamber. In eukaryotes, the chaperonin TRiC/CCT is hetero-oligomeric, consisting of two stacked rings of eight paralogous subunits each. TRiC facilitates folding of approximately 10% of the eukaryotic proteome, including many cytoskeletal components and cell cycle regulators. Folding of many cellular substrates of TRiC cannot be assisted by any other chaperone. A complete structural and mechanistic understanding of this highly conserved and essential chaperonin remains elusive. However, recent work is beginning to shed light on key aspects of chaperonin function and how their unique properties underlie their contribution to maintaining cellular proteostasis.

59 BASIC BIOLOGICAL SCIENCES↗

Surface chemistry and physics of deuterium retention in lithiated graphite

Lithium wall conditioning in TFTR, CDX-U, T-11M, TJ-II and NSTX is found to yield enhanced plasma performance manifest, in part, through improved deuterium particle control. X-ray photoelectron spectroscopy (XPS) experiments examine the affect of D irradiation on lithiated graphite and show that the surface chemistry of lithiated graphite after D ion bombardment (500 eV/amu) is fundamentally different from that of non-Li conditioned graphite. Instead of simple LiD bonding seen in pure liquid Li, graphite introduces additional complexities. XPS spectra show that Li-O-D (533.0 {+-} 0.6 eV) and Li-C-D (291.4 {+-} 0.6 eV) bonds, for a nominal Li dose of 2 {micro}m, become 'saturated' with D at fluences between 3.8 and 5.2 x 10{sup 17} cm{sup -2}. Atomistic modeling indicate that Li-O-D-C interactions may be a result of multibody effects as opposed to molecular bonding.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microstructural features and deuterium diffusion in lithium penta-aluminate pellets under He + and D + ion irradiation

Lithium (Li) penta-aluminate (LiAl 5 O 8 ) is investigated as a potential tritium (T) breeding material, with a focus on microstructural response to ion irradiation and deuterium (D) diffusion behavior. Under high-fluence ion irradiation (2 x 10 17 (He + +D + )/cm 2 ) at 773 K, LiAl 5 O 8 exhibits significant disorder on the Li sublattice, as revealed by atomic-resolution scanning transmission electron microscopy, while the Al and O sublattices remain stable, demonstrating strong resistance to structural amorphization. Irradiation induces the formation of platelet-shaped antiphase boundaries (APBs), which may serve as effective D trapping sites. Atom probe tomography suggests the presence of 6 LiD clusters in the mass spectra, though definite conclusions regarding APB composition are hindered by signal overlap and limited data statistics. Time-of-flight secondary ion mass spectrometry reveals that D retention approaches to saturation at 3 x 10 17 (He + +D + )/cm 2 . Isothermal and isochronal annealing studies determine an average diffusivity of 1.6 x 10 -13 at 773 K and an effective activation energy of 0.8 ± 0.1 eV for D migration. Compared to γ-LiAlO 2 , LiAl 5 O 8 demonstrates superior irradiation resistance, minimal Li loss, and enhanced D retention, underscoring its potential as a durable breeder material for T production. In conclusion, these findings provide key insights into the microstructural evolution, defect dynamics, and D retention mechanisms in LiAl 5 O 8 under reactor-relevant conditions.

42 ENGINEERING↗

Specific corrosion product on interior surface of a bronze wine vessel with loop-handle and its growth mechanism, Shang Dynasty, China

In this paper, a kind of specific stalactitic product was found on the interior surface of a covered bronze wine vessel with loop-handle (Chinese name is you), which was fabricated in Shang Dynasty (1700 B.C.-1100 B.C.) and now is collected in Xiaogan Museum, Hubei province of China. The microstructures of the product were characterized systematically by using optical microscopy, scanning electron microscope, transmission electron microscope, X-ray diffraction, and Raman microscopy. The experimental results revealed that the product belonged to a kind of malachite with high purity and high crystallinity. The growth of the product was considered to be a possible reason that the vessel was overly airtight within a museum display cabinet besides a lid of the vessel, which made the excess of H{sub 2}O and CO{sub 2} gas concentrations inside the vessel during long-term storage. This corrosion product is very harmful to bronze cultural relics, because of a large amount of copper consumption from the matrix which will reduce its life. The growth mechanism of the specific stalactitic product and the suggestions for preservation of the similar bronze relics in museum were proposed. - Highlights: Black-Right-Pointing-Pointer The stalactitic product was the high purity and good crystallinity malachite. Black-Right-Pointing-Pointer Its growth was related to the excess of H{sub 2}O and CO{sub 2} gas concentrations in museum. Black-Right-Pointing-Pointer It is harmful to the bronzes, because copper will be consumed from the matrix. Black-Right-Pointing-Pointer The suggestions for preservation of the similar bronzes in museum were proposed.

36 MATERIALS SCIENCE↗

Innovative laser-based methods for monitoring fuel retention in ITER

This paper addresses the challenge of tritium inventory management in ITER and future fusion reactors, highlighting the importance of accurate tritium measurement and its spatial distribution within the vacuum vessel. Given ITER’s operational constraints, especially the limit on tritium retention, precise measurement is essential for both safety and regulatory compliance. To tackle these questions, the paper presents the T-monitor diagnostic system developed by Forschungszentrum Jülich, which uses Laser-Induced Desorption (LID) in combination with Diagnostic Residual Gas Analysis (DRGA) to measure hydrogen isotope concentrations on the surface of divertor tiles. The system integrates a high-power laser, advanced optical components, and a Fast Scanning Mirror Unit (FSMU) for accurate laser spot positioning with rapid response. Designed to measure in situ tritium retention, the diagnostic provides high-resolution spatial mapping, vital for evaluating detritiation strategies. The laser heating process increases the divertor surface temperature to 1600 K within the laser spot, promoting hydrogen isotope desorption. Accurate measurements require the precise control of laser parameters, including pulse duration and spot size, with a target relative accuracy of 20%. The optical design includes both in-vessel and ex-vessel components, such as durable high-reflectivity mirrors made of gold and copper, selected not only for their infrared performance but also for their transmission of visible wavelengths for observation purposes. To protect optical components from contamination, a pneumatic shutter is used.

Hydrogen isotopes↗

Three-dimensional flow model development for thermal mixing and stratification modeling in reactor system transients analyses

Mixing, thermal-stratification, and mass transport phenomena in large pools or enclosures play major roles for the safety of reactor systems. Depending on the fidelity requirement and computational resources, various modeling methods, from the 0-D perfect mixing model to 3-D Computational Fluid Dynamics (CFD) models, are available. Each is associated with its own advantages and shortcomings. It is very desirable to develop an advanced and efficient thermal mixing and stratification modeling capability embedded in a modern system analysis code to improve the accuracy of reactor safety analyses and to reduce modeling uncertainties. An advanced system analysis tool, SAM, is being developed at Argonne National Laboratory for advanced non-LWR reactor safety analysis. While SAM is being developed as a system-level modeling and simulation tool, a reduced-order three-dimensional module is under development to model the mull-dimensional flow and thermal mixing and stratification in large enclosures of reactor systems. This paper provides an overview of the three-dimensional finite element flow model in SAM, including the governing equations, stabilization scheme, and solution methods. Additionally, several verification and validation tests are presented, including lid-driven cavity flow, natural convection inside a cavity, laminar flow in a channel of parallel plates. Here, based on the comparisons with the analytical solutions and experimental results, it is demonstrated that the developed 3-D fluid model can perform very well for a wide range of flow problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Structural Basis of Lysosomal Phospholipase A 2 Inhibition by Zn 2+

Lysosomal phospholipase A 2 (LPLA 2 /PLA2G15) is a key enzyme involved in lipid homeostasis and is characterized by both phospholipase A2 and transacylase activity and by an acidic pH optimum. Divalent cations such as Ca 2+ and Mg 2+ have previously been shown to have little effect on the activity of LPLA 2 , but the discovery of a novel crystal form of LPLA 2 with Zn 2+ bound in the active site suggested a role for this divalent cation in regulating enzyme activity. In this complex, the cation directly coordinates the serine and histidine of the α/β-hydrolase triad and stabilizes a closed conformation. This closed conformation is characterized by an inward shift of the lid loop, which extends over the active site and effectively blocks access to one of its lipid acyl chain binding tracks. Therefore, we hypothesized that Zn 2+ would inhibit LPLA 2 activity at a neutral but not acidic pH because histidine would be positively charged at lower pH. Indeed, Zn 2+ was found to inhibit the esterase activity of LPLA 2 in a noncompetitive manner exclusively at a neutral pH (between 6.5 and 8.0). Because lysosomes are reservoirs of Zn 2+ in cells, the pH optimum of LPLA 2 might allow it to catalyze acyl transfer unimpeded within the organelle. Here, we conjecture that Zn 2+ inhibition of LPLA 2 at higher pH maintains a lower activity of the esterase in environments where its activity is not typically required.

59 BASIC BIOLOGICAL SCIENCES↗

Crystal Structure of d -Ornithine/ d -Lysine Decarboxylase, a Stereoinverting Decarboxylase: Implications for Substrate Specificity and Stereospecificity of Fold III Decarboxylases

A newly discovered Fold III pyridoxal 5'-phosphate (PLP)-dependent decarboxylase, d-ornithine/lysine decarboxylase (DOKDC), catalyzes decarboxylation of d-lysine and d-ornithine with inversion of stereochemistry. The X-ray crystal structure of DOKDC has been determined to 1.72 Å. DOKDC has a low level of sequence identity (<30%) with meso-diaminopimelate decarboxylase (DAPDC) and l-lysine/ornithine decarboxylase (LODC), but its three-dimensional structure is very similar. The distal binding site of DAPDC contains a conserved arginine that forms an ion pair with the l-carboxylate end of DAP. In both LODC and DOKDC, this distal site is modified by replacement of the arginine with aspartate, changing the substrate specificity. l-Ornithine decarboxylase (ODC) and LODC have a conserved phenylalanine on the re-face of the PLP complex that has been found to play a key role in the decarboxylation mechanism. We have found that both DAPDC and DOKDC have tyrosine instead of phenylalanine at this position, which precludes the binding of l-amino acids. Because the PLP-binding lysine in ODC, LODC, DAPDC, and DOKDC is located on the re-face of the PLP, we propose that this is the acid group responsible for protonation of the product, thus resulting in the observed retention of configuration for decarboxylation of l-amino acids and inversion for decarboxylation of d-amino acids. Here, the reactions of DAPDC and DOKDC are likely accelerated by positive electrostatics on the re-face by the lysine ε-ammonium ion and on the si-face by closure of the lid over the active site, resulting in desolvation and destabilization of the d-amino acid carboxylate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗