Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

On the use of adjoints in the inversion of observed quasi-static deformation

An adjoint-based conjugate gradient algorithm provides an efficient means for imaging sources of deformation within the Earth, such as volume stresses associated with fluid flow in aquifers and reservoirs. For time intervals over which the overburden deforms elastically, one can calculate the gradient elements for a single model update using just two numerical simulations. The first is a forward run that is used to compute the residuals associated with the given iteration. The second simulation is to evaluate the application of the adjoint operator to the residuals. In this adjoint calculation, the residual displacements are applied as sources at the measurement locations, driving the deformation in the simulation. The volume stress on the source grid blocks, in response to the residual displacements, provide the gradient components. We apply this technique to satellite-based interferometric synthetic aperture radar (InSAR) line-of-sight displacements that were observed over an oil reservoir in California’s Central Valley. Furthermore, we find that the adjoint-based gradient estimates, requiring 18 CPU seconds, agree with conventional numerical calculations that take over 3700 CPU seconds to compute. Conjugate gradient algorithms utilizing the conventional approach and adjoint-based gradient computations give roughly the same reductions in misfit and similar final estimates of reservoir volume change.

58 GEOSCIENCES↗
DOE OSTI2023

Validation of semi-analytical, semi-empirical covariance matrices for two-point correlation function for early DESI data

ABSTRACT We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel’dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Chemical and structural characterization of EUV photoresists as a function of depth by standing-wave x-ray photoelectron spectroscopy

The success in the miniaturization of the electronic device constituents depends mostly on the photolithographic techniques. Recently, to achieve patterning at the sub-10-nm node, extreme ultraviolet (EUV) lithography has been introduced into high volume production. Continued scaling of EUV via increased numerical aperture to achieve nodes at 3-nm and below requires the development of fundamentally new patterning materials and new characterization methods. Current EUV-resist film thicknesses are in the 20- to 40-nm range, and further thickness reduction is required for the next generation. Therefore, interfaces become exceedingly important, and the properties of the resist film would be dominated by top and bottom interfacial effects. X-ray photoelectron spectroscopy (XPS) combined with standing-wave excitation (SW-XPS), a fairly new method in the EUV lithography field, previously had been largely applied in multilayers and superlattices for characterizing the composition and electronic structure of buried layers and interfaces as a function of depth. We applied the SW-XPS method to organic/inorganic photoresists to provide depth-selective information on their structural and chemical conditions of as a function of temperature, EUV exposure, different underlayers, and other fundamental parameters. As a first attempt, we perform an SW-XPS feasibility study on self-assembled monolayer (SAM) films after exposure to an electron beam. By SW-XPS, we determined that the interface between the Al 2 O 3 underlayer and the SAMs is smooth, with a mean roughness of about 0.2 nm. Moreover, we determined that the SAM chains are, on average, tilted by 1/430 deg off the sample normal. The SW-XPS results also suggest that the SAM is not a perfectly aligned and uniform monolayer, with some areas having thickness higher than a single monolayer. We demonstrated that SW-XPS can provide useful information on ultrathin materials with high potential for being used as a characterization method of organic/inorganic photoresists.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Establishing gold and platinum standards to 1 terapascal using shockless compression

New techniques are advancing the frontier of high-pressure physics beyond 1 terapascal, leading to new discoveries and offering stringent tests for condensed-matter theory and advanced numerical methods. However, the ability to absolutely determine the pressure state remains challenging, and well-calibrated pressure-density reference materials are required. We conducted shockless dynamic compression experiments at the National Ignition Facility and the Z machine to obtain quasi-absolute, high-precision, pressure-density equation-of-state data for gold and platinum. We derived two experimentally constrained pressure standards to terapascal conditions. Establishing accurate experimental determinations of extreme pressure will facilitate better connections between experiments and theory, paving the way toward improving our understanding of material response to these extreme conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Assessments of data centers for provision of frequency regulation

There are numerous opportunities for data centers to participate in demand response programs considering their large energy capacities, flexible working environments and workloads, redundant design and operation, etc. As a type of demand response, frequency regulation requires fast response, and its potential is not fully explored by data centers yet. This paper proposes a synergistic control strategy for data center frequency regulation which uses both IT and cooling systems. It combines power management techniques at the server level with control of the chilled water supply temperature to track the regulation signal from the electrical market. A frequency regulation flexibility factor is also proposed to increase the IT capacity for frequency regulation. The performance of the control strategy is studied through numerical simulations using an equation based object-oriented Modelica platform designed for data centers. Simulation results show that with well-tuned control parameters, data centers can provide frequency regulation service in both regulation up and down. The performance of data centers in providing frequency regulation service is largely influenced by the regulation capacity bid, frequency regulation flexibility factor, workload condition, and cooling mode of the cooling system, and not significantly influenced by the time constant of chillers. In addition, compared with a server-only control strategy, the proposed synergistic control strategy can provide an extra regulation capacity of 3% of the design power when chillers are activated. Here, when chillers are deactivated, both strategies have a similar regulation capacity.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Generation of selective single-mode guided waves by d 36 type piezoelectric wafer

In general, mixed modes and dispersion of guided wave occur at any excitation frequency, while single modes are needed for effective nondestructive evaluation of structures. In this work, we present an approach to generate a selective single-mode guided wave in plate-like structures by exploiting the unique directionality of the d 36 type piezoelectric wafer and the symmetry of fundamental guided wave modes. Specifically, we devise a unified fundamental shear horizontal (SH 0 ) wave/fundamental antisymmetric mode (A 0 ) wave directional transducer pair by attaching the d 36 type piezoelectric wafers back to back on the opposite sides of the plate to provide selective, directional single-mode guided waves for actuation and sensing. Originally, the d 36 -type piezoelectric wafer, poled and cut from the lead magnesium niobate-lead titanate crystal, produces a mixture of the in-plane SH 0 (symmetric) mode and out-of-plane A 0 (asymmetric) and S 0 (symmetric) modes. By applying the in-phase or out-of-phase applied electrical field to the coupled d 36 -type piezoelectric wafer pair, our devised approach generates selective single modes, SH 0 or A 0 , respectively, to propagate in the structure. We describe the theoretical development of our approach and conduct both numerical simulations and laboratory experiments for validations.

42 ENGINEERING↗
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimization of flow in additively manufactured porous columns with graded permeability

Chemical engineering systems often involve a functional porous medium, such as in catalyzed reactive flows, fluid purifiers, and chromatographic separations. Ideally, the flow rates throughout the porous medium are uniform, and all portions of the medium contribute efficiently to its function. The permeability is a property of a porous medium that depends on pore geometry and relates flow rate to pressure drop. Additive manufacturing techniques raise the possibilities that permeability can be arbitrarily specified in three dimensions, and that a broader range of permeabilities can be achieved than by traditional manufacturing methods. Using numerical optimization methods, we show that designs with spatially varying permeability can achieve greater flow uniformity than designs with uniform permeability. We consider geometries involving hemispherical regions that distribute flow, as in many glass chromatography columns. By several measures, significant improvements in flow uniformity can be obtained by modifying permeability only near the inlet and outlet.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A survey on degradation modeling, prognosis, and prognostics-driven maintenance in wind energy systems

Wind energy generation proliferated over the past decades, introducing unique challenges and opportunities for failure prediction, operation and maintenance. Decision-makers are continuously looking into new methods to infer failure mechanisms and behaviors of wind turbine components to detect and intervene in the failures before they happen. Evidently, degradation modeling and prognosis become engaging topics for researchers and practitioners to prevent catastrophic failures. Prognostics-driven approaches predict the time of failure for the components (e.g., predicting remaining useful life), which provides significant insights for scheduling of operations and maintenance activities. Integrating these prognostics-driven insights into wind farm operations and maintenance presents a substantial challenge, demanding careful consideration of numerous factors such as accessibility, crew routing, and spare part logistics. This study provides state-of-the-art review for degradation modeling, prognosis, and prognostics-driven maintenance techniques for wind energy systems. The discussed techniques align with the United Nations' sustainable development goals, in particular Goal 7 (Affordable and Clean Energy), by enhancing effectiveness and sustainability of wind energy operations. This work also showcases open research questions related to degradation modeling, prognosis, and prognostics-driven maintenance.

Altinpulluk, Nur Banu↗
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
DOE OSTI2020

CRNT4SBML: a Python package for the detection of bistability in biochemical reaction networks

Motivation: Signaling pathways capable of switching between two states are ubiquitous within living organisms. They provide the cells with the means to produce reversible or irreversible decisions. Switchlike behavior of biological systems is realized through biochemical reaction networks capable of having two or more distinct steady states which are dependent on initial conditions. Investigation of whether a certain signaling pathway can confer bistability involves a substantial amount of hypothesis testing. The cost of direct experimental testing can be prohibitive. Therefore, constraining the hypothesis space is highly bene?cial. One such methodology is based on Chemical Reaction Network Theory which uses computational techniques to rule out pathways that are not capable of bistability regardless of kinetic constant values and molecule concentrations. Although useful, these methods are complicated from both pureandcomputationalmathematicsperspectives.Thus,theiradoptionisverylimitedamongstbiologists. Results: We brought Chemical Reaction Network Theory approaches closer to experimental biologists by automating all the necessary steps in CRNT4SMBL. The input is based on SBML format which is the community standard for biological pathway communication. The tool parses SBML and derives C-graph representations of the biological pathway with mass action kinetics. Next steps involve an ef?cient search for potential saddle-node bifurcation points using an optimization technique. This type of bifurcation is important as it has the potential of acting as a switching point between two steady states. Finally, if any bifurcation points are present, numerical continuation analysis extends the equilibria branches for generating the diagram. Presence of an S-shaped bifurcation diagram indicates that the pathway acts as a bistable switch for the given optimization parameters. Availability: CRNT4SBML is available via the Python Package Index. The documentation can be found at https://crnt4sbml.readthedocs.io. CRNT4SBML is licensed under the Apache Software License 2.0. Contact: vladislav.petyuk@pnnl.gov

Reyes, Brandon C.↗
DOE OSTI2022

Spectral density reconstruction with Chebyshev polynomials

Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semiexclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. Here, in the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to obtain the local density of states for graphene in a magnetic field as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Charge per Micro-Pulse Calculation for Ideal Lujan and WNR beams

A complete redesign of the LANSCE front-end is currently under development for the Los Alamos Modernization Project (LAMP). This includes the replacement of the traditional Cockroft-Walton injection system to the newer radiofrequency quadrupole (RFQ) standard. LANSCE accelerates both H- and H+ beams and therefore requires an unconventional application for RFQ systems. A new facility could consider using an independent RFQ for each beam species. However, with the tight space of the LANSCE injector facility, it would prove difficult to instrument. Our initial studies have shown that a single RFQ could be used. The negatively- and positively-charged direct current (DC) beams are shaped, bunched, and accelerated out of phase to each other, with very little interaction. This would work for all possible H+ beam gates and most H- Beam gates. However, a challenge arises with the Weapons Neutron Research (WNR) beam structure that requires the acceleration of bunched beams at higher charge density per pulse. The Long Bunch Enable Gate (LBEG) used for the Lujan center is the best example of a standard injector beam gate. The beam gate, or Macro-Pulse (MP), is 625-us long, separated by at least 8.3 ms. A chopper is used to dice the MP into 1750 mini-Pulses (mP), which are 290-ns wide. After chopping, the MP contains 100,000 micro-Pulses (uP), which are generated by the linac acceleration structure modulated at 201.25 MHz; that is, the micro-Pulses are separated 5 ns apart from each other. The H+ beam structure will be very similar once a chopper is added to its transport. However, the Multi-Pulse Enable Gate (MPEG) used for WNR is significantly different in structure. In this case, the mP is chopped to be only 35-ns wide. This width is then rotated into a single 201.25 MHz RF bucket using the Low Frequency Buncher (LFB). This is proving difficult for the RFQ design in the LAMP upgrade. We began investigating the amount of current sent to WNR for two reasons. Firstly, we wanted to find a solution to avoid using the Low Frequency Buncher to simplify the LAMP RFQ design. It is simple to inject and accelerate a 35 mA DC beam in the LAMP RFQ (this is done in numerous accelerator facilities). However, injecting a short bunch with a larger current (once compressed to 5 ns to meet the RFQ frequency) may be difficult. Second, we have a technique called PSR2WNR in which we do not employ the Low Frequency Buncher, but rather accumulate 5 LBEG uP (5 ns) in the PSR and deliver the accumulated beam to WNR every 1.8 us. As a result, we could provide to WNR with five accumulated LBEG uPs every shot. The goal here is to create a rapid kicker with a frequency response of 555 kHz. DARHT-II and other accelerator facilities have constructed kickers of this type for other uses, and our engineers consider this technology as feasible.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Scattering Observables from Few-Body Densities and Compton Scattering on $^6$Li

The dynamics of scattering on light nuclei is numerically expensive using standard methods. Fortunately, recent developments allow one to factor the relevant quantities for a given probe into a convolution of an n -body Transition Density Amplitude (TDA) and the interaction kernel for a given probe. These TDAs depend only on the target, and not the probe; they are calculated once for each set of kinematics and can be used for different interactions.\ % in the same kinematics. The kernels depend only on the probe, and not on the target; they can be reused for different targets and different kinematics. The calculation of TDAs becomes numerically difficult for more than four nucleons, but we discuss a new solution through the use of a Similarity Renormalization Group transformation, and a subsequent back-transformation. This technique allows for extending the TDA method to heavier nuclei such as 6 Li. We present preliminary results for Compton scattering on 6 Li and compare with available data, anticipating an upcoming, more thorough study. We also discuss ongoing extensions to pion-photoproduction and other reactions on light nuclei.

Long, Alexander [George Washington University, Was↗
DOE OSTI2022

Parametric dynamic mode decomposition for reduced order modeling

Dynamic Mode Decomposition (DMD) is a model-order reduction approach, whereby spatial modes of fixed temporal frequencies are extracted from numerical or experimental data sets. The DMD low-rank or reduced operator is typically obtained by singular value decomposition of the temporal data sets. For parameter-dependent models, as found in many multi-query applications such as uncertainty quantification or design optimization, the only parametric DMD technique developed was a stacked approach, with data sets at multiple parameter values were aggregated together, increasing the computational work needed to devise low-rank dynamical reduced-order models. Here in this paper, we present two novel approach to carry out parametric DMD: one based on the interpolation of the reduced-order DMD eigen-pair and the other based on the interpolation of the reduced DMD (Koopman) operator. Numerical results are presented for diffusion-dominated nonlinear dynamical problems, including a multiphysics radiative transfer example. All three parametric DMD approaches are compared.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2022

Machined phononic crystals to block high-order Lamb waves and crosstalk in through-metal ultrasonic communication systems

For systems that require complete metallic enclosures (e.g., containment buildings for nuclear reactors), it is impossible to access interior sensors and equipment using standard electromagnetic techniques. A viable way to communicate and supply power through metallic barriers is the use of elastic waves and ultrasonic transducers, introducing several design challenges that must be addressed. Specifically, the use of multiple communication channels on the same enclosure introduces an additional mechanism for signal crosstalk between channels: guided waves propagating in the barrier between channels. This work numerically and experimentally investigates a machined phononic crystal to block MHz Lamb wave propagation between ultrasonic communication channels, greatly reducing wave propagation and the resulting crosstalk voltage. Blind grooves are machined into one or both sides of a metallic barrier to introduce a periodic unit cell, greatly altering the guided wave dispersion in the barrier. Numerical simulations are used to determine a set of groove geometries for testing, and experiments were performed to characterize the wave-blocking performance of each design. The best-performing design was tested using piezoelectric transducers bonded to the barrier, showing a 14.4 dB reduction in crosstalk voltage. Overall, the proposed periodic grooving method is a promising technique for completely isolating ultrasonic power/data transfer systems operating in a narrow frequency range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗