Results for “Numerical Analysi” — Page 7

DOE OSTI2022

Broadband piezoelectric energy harvesting microgyroscopes: Design and nonlinear analysis

Small devices in remote or difficult-to-reach areas can benefit from harvesting energy from mechanical wasted energy, which reduces the requirement for a new power source. A multi-purpose energy harvesting microgyroscope system based on piezoelectric materials is suggested. The necessity of taking spatially varying electrostatic forces is examined. Considering the effects of the system's width and thickness, DC voltage, and angular speed the systems inherent frequencies are found. The partial differential equations describing the system's dynamics are numerically solved by using the differential quadrature method. Further, the numerical analysis enables to identify the optimal system design for broadband energy harvesting. The simulation results reveal that a system with a non-symmetric beam design is adequate for broadband energy harvesting. This is associated with the applied DC voltage, which may be modified to improve the broadband frequency of the system. It is concluded that the nonlinear softening effects create a broadband frequency response with high voltage output. However, DC and AC voltages need to be carefully selected in order to avoid the dynamic pull-in.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Exact WKB and the quantum Seiberg-Witten curve for 4d $N$ = 2 pure SU(3) Yang-Mills. Abelianization

We investigate the exact WKB method for the quantum Seiberg-Witten curve of 4d N = 2 pure SU(3) Yang-Mills in the language of abelianization. The relevant differential equation is a third-order equation on $\mathbb{CP}^1$ with two irregular singularities. We employ the exact WKB method to study the solutions to such a third-order equation and the associated Stokes phenomena. We also investigate the exact quantization condition for a certain spectral problem. Moreover, exact WKB analysis leads us to consider new Darboux coordinates on a moduli space of flat SL(3,$\mathbb{C}$)-connections. In particular, in the weak coupling region we encounter coordinates of the higher length-twist type generalizing Fenchel-Nielsen coordinates. The Darboux coordinates are conjectured to admit asymptotic expansions given by the formal quantum periods series and we perform numerical analysis supporting this conjecture.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Techno-Economic Analysis of Dynamic Building Envelopes Comprising Phase Change Materials and Switchable Insulations

Integrating phase change materials (PCMs) in building envelopes is a recognized technique to reduce the space heating/cooling loads and provide load shedding and shifting capacity. However, PCM benefits have been found to reduce dramatically with increase in the insulation requirements. Dynamic building envelopes that can alter their thermal resistance based on the indoor and outdoor conditions can effectively harness free ambient heating/cooling, thereby greatly enhancing the benefits of the PCM-integrated envelope in managing thermal loads. In this study, we examine various combinations of PCMs and switchable insulations and compare their combined impact on heat flow through the building envelope. Using numerical analysis, we investigate the PCM- and switchable insulation-integrated building walls, calculate the potential energy saving benefits compared to the conventional static walls, and finally perform techno-economic analysis to estimate the acceptable cost of the technology under various payback period scenarios.

dynamic building envelopes↗

DOE OSTI2023

Techno-Economic Analysis of Dynamic Building Envelopes Comprising Phase Change Materials and Switchable Insulations: Preprint

Integrating phase change materials (PCMs) in building envelopes is a recognized technique to reduce the space heating/cooling loads and provide load shedding and shifting capacity. However, PCM benefits have been found to reduce dramatically with increase in the insulation requirements. Dynamic building envelopes that can alter their thermal resistance based on the indoor and outdoor conditions can effectively harness free ambient heating/cooling, thereby greatly enhancing the benefits of the PCM-integrated envelope in managing thermal loads. In this study, we examine various combinations of PCMs and switchable insulations and compare their combined impact on heat flow through the building envelope. Using numerical analysis, we investigate the PCM- and switchable insulation-integrated building walls, calculate the potential energy saving benefits compared to the conventional static walls, and finally perform technoeconomic analysis to estimate the acceptable cost of the technology under various payback period scenarios.

dynamic building envelopes↗

DOE OSTI2020

Thermal Analysis of Large Area Additive Manufacturing Resistance Heating Composites for Out of Oven/Autoclave Applications

Additive Manufacturing (AM) of carbon fiber (CF) reinforced composite has received growing attention because of the design flexibility, superior mechanical properties, improved thermal properties, and weight reduction. Autoclave tooling was proven to be a successful application for large scale AM technology. The capital cost, and cost associated with heating, and cycle time in a conventional autoclave process is relatively high. Thus, an innovative design of AM mold with an efficient heating scheme is essential. This study represents an innovative method of the resistive heating of composite molds which does not require a room size oven for heating during the curing processing. Therefore, it has the potential to reduce the operating cost drastically. For the design validation and feasibility study, we performed a numerical analysis of the wire embedded and AM mold parts. The goal of this study is to determine and optimize the thermal behavior of the printed mold with embedded wire technology. It is anticipated that the larger distance between the embedded wires along the printing direction (z-direction) increase the cold spot, on the other hand, a close distance of the wire can create the unwanted localize heating, thus melting. Constant thermal properties of the 20 wt.% short CF reinforced acrylonitrile butadiene styrene (ABS) was used for the simulation purpose. Thermal characterization was set to 100°C to avoid the thermal deformation or bulging on the part surface.

billah, Kazi↗

DOE OSTI2020

Thermal Analysis of Large Area Additive Manufacturing Resistance Heating Composites for Out of Oven/Autoclave Applications

Additive Manufacturing (AM) of carbon fiber (CF) reinforced composite has received growing attention because of the design flexibility, superior mechanical properties, improved thermal properties, and weight reduction. Autoclave tooling was proven to be a successful application for large scale AM technology. The capital cost, and cost associated with heating, and cycle time in a conventional autoclave process is relatively high. Thus, an innovative design of AM mold with an efficient heating scheme is essential. This study represents an innovative method of the resistive heating of composite molds which does not require a room size oven for heating during the curing processing. Therefore, it has the potential to reduce the operating cost drastically. For the design validation and feasibility study, we performed a numerical analysis of the wire embedded and AM mold parts. The goal of this study is to determine and optimize the thermal behavior of the printed mold with embedded wire technology. It is anticipated that the larger distance between the embedded wires along the printing direction (z-direction) increase the cold spot, on the other hand, a close distance of the wire can create the unwanted localize heating, thus melting. Constant thermal properties of the 20 wt.% short CF reinforced acrylonitrile butadiene styrene (ABS) was used for the simulation purpose. Thermal characterization was set to 100°C to avoid the thermal deformation or bulging on the part surface.

billah, Kazi↗

DOE OSTI2020

Thermal Analysis of Large Area Additive Manufacturing Resistance Heating Composites for Out of Oven/Autoclave Applications

Additive Manufacturing (AM) of carbon fiber (CF) reinforced composite has received growing attention because of the design flexibility, superior mechanical properties, improved thermal properties, and weight reduction. Autoclave tooling was proven to be a successful application for large scale AM technology. The capital cost, and cost associated with heating, and cycle time in a conventional autoclave process is relatively high. Thus, an innovative design of AM mold with an efficient heating scheme is essential. This study represents an innovative method of the resistive heating of composite molds which does not require a room size oven for heating during the curing processing. Therefore, it has the potential to reduce the operating cost drastically. For the design validation and feasibility study, we performed a numerical analysis of the wire embedded and AM mold parts. The goal of this study is to determine and optimize the thermal behavior of the printed mold with embedded wire technology. It is anticipated that the larger distance between the embedded wires along the printing direction (z-direction) increase the cold spot, on the other hand, a close distance of the wire can create the unwanted localize heating, thus melting. Constant thermal properties of the 20 wt.% short CF reinforced acrylonitrile butadiene styrene (ABS) was used for the simulation purpose. Thermal characterization was set to 100°C to avoid the thermal deformation or bulging on the part surface.

billah, Kazi↗

DOE OSTI2020

Non-destructive simulation of node defects in additively manufactured lattice structures

Additive Manufacturing (AM), commonly referred to as 3D printing, offers the ability to not only fabricate geometrically complex lattice structures but parts in which lattice topologies in-fill volumes bounded by complex surface geometries. However, current AM processes produce defects on the strut and node elements which make up the lattice structure. This creates an inherent difference between the as-designed and as-fabricated geometries, which negatively affects predictions (via numerical simulation) of the lattice’s mechanical performance. Although experimental and numerical analysis of an AM lattice’s bulk structure, unit cell and struts have been performed, there exists almost no research data on the mechanical response of the individual as-manufactured lattice node elements. Here we propose a methodology that, for the first time, allows non-destructive quantification of the mechanical response of node elements within an as-manufactured lattice structure. A custom-developed tool is used to extract and classify each individual node geometry from micro-computed tomography scans of an AM fabricated lattice. Voxel-based finite element meshes are generated for numerical simulation and the mechanical response distribution is compared to that of the idealised computer-aided design model. The method demonstrates compatibility with Uncertainty Quantification methods that provide opportunities for efficient prediction of a population of nodal responses from sampled data. Overall, the non-destructive and automated nature of the node extraction and response evaluation is promising for its application in qualification and certification of additively manufactured lattice structures.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗

DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗

DOE OSTI2025

Time correlations from steady-state expectation values

Recovering properties of correlation functions is typically challenging. On one hand, experimentally, it requires measurements with a temporal resolution finer than the system's dynamics. On the other hand, analytical or numerical analysis requires solving the system evolution. Here, we use recent results of quantum metrology with continuous measurements to derive general lower bounds on the relaxation and second-order correlation times that are both easy to calculate and measure. These bounds are based solely on steady-state expectation values and their derivatives with respect to a control parameter, and can be readily extended to the autocorrelation of arbitrary observables. We validate our method on two examples of critical quantum systems: a critical driven-dissipative resonator, where the bound matches analytical results for the dynamics, and the infinite-range Ising model, where only the steady state is solvable and thus the bound provides information beyond the reach of existing analytical approaches. Our results can be applied to experimentally characterize ultrafast systems, and to theoretically analyze many-body models with dynamics that are analytically or numerically hard.

Górecki, Wojciech [INFN, Pavia]↗

DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Exsolution Synthesis of Nanocomposite Perovskites with Tunable Electrical and Magnetic Properties

Nanostructured functional oxides play an important role in enabling clean energy technologies and novel memory and processor devices. Using thin-film La 0.6 Sr 0.4 FeO 3 (LSF) as a model system, the novel utility of exsolution in fabricating self-assembled metal oxide nanocomposites with tunable functionalities is shown. Exsolution triggers the formation of metallic iron (Fe 0 ) nanoparticles, Ruddlesden–Popper domains, and nm-scale percolated Fe-deficient channels in LSF. Combining multimodal characterization with numerical modeling, the chemical, magnetic, and electrical properties of the exsolution-synthesized nanocomposite at different stages of Fe 0 exsolution as well as during redox cycling are assessed. After exsolution, the electronic conductivity of the nanocomposite LSF increased by more than two orders of magnitude. Based on numerical analysis representing all the constituents, it is expected that this increase in conductivity originates mainly from the Fe-deficient percolating channels formed during exsolution. Moreover, the exsolved nanocomposite is redox-active even at moderate temperatures. Such redox capabilities can enable dynamic control of the nanocomposite functionality by tailoring the oxygen non-stoichiometry. This concept is demonstrated with a continuous modulation of magnetization between 0 and 110 emu cm -3 . These findings point out that exsolution may serve as a platform for scalable fabrication of complex metal oxide nanocomposites for electrochemical and electronic applications.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Hybrid Plasmonic Coupled‐Mode Metasurface Design for Shot‐Noise Limited Ppb‐Level Hydrogen Detection

Plasmonic hydrogen sensors have enabled hydrogen detection below parts-per-million (ppm) range by boosting the sensitivity using localized surface plasmonic resonant (LSPR) structures. However, the intrinsic optical losses of Palladium (Pd), the primary plasmonic metal used for hydrogen detection, result in a low quality (Q) factor LSPR, which fundamentally hinders further improvement. In this work, a hybrid plasmonic metasurface is proposed that couples Pd-based LSPR structure with an Au film supporting surface plasmon polariton mode (Au-SPP). The coupled near-perfect absorber resonance yields a spectrally narrow, high Q response that retains strong sensitivity to hydrogen while improving resonance localization. Numerical analysis shows that, under shot-noise-limited conditions, the limit of detection (LoD) can be improved by over threefold compared to the state-of-the-art designs. Furthermore, this hybrid plasmonic coupled-mode metasurface thus presents a promising pathway to achieve parts-per-billion-level (ppb-level) hydrogen detection with enhanced spectral precision and robustness.

coupled modes↗

DOE OSTI2020

Cascade decays of heavy Higgs bosons through vectorlike quarks in two Higgs doublet models

We study cascade decays of heavy neutral Higgs bosons through vectorlike quarks. We focus on scenarios where decay modes into pairs of vectorlike quarks are not kinematically open which extends the sensitivity of the LHC to larger masses. Assuming only mixing with the third family of standard model quarks the new decay modes of heavy Higgs bosons are: H → t 4 t → Wbt, Ztt, htt and H → b 4 b → Wtb, Zbb, hbb, where t 4 (b 4 ) is the new up-type (down-type) quark mass eigenstate. In the numerical analysis we assume the CP even Higgs boson in the two Higgs doublet model type-II but the signatures are relevant for many other scenarios. We identify the region of the parameter space where these decay modes are significant or can even dominate, and thus they provide the best opportunities for the simultaneous discovery of a new Higgs boson and vectorlike quarks. We further explore the reach of the High Luminosity LHC for two representative decay modes,t 4 → Zt → ℓℓt and b 4 → Zb → ℓℓb, and found that cross sections at a 0.1 fb level can be probed with simple cut based analyses. We also find that the rates for Higgs cascade decays can be much larger than the rates for a single production of vectorlike quarks. Furthermore, the reach for vectorlike quarks in Higgs cascade decays and pair production extends to comparable masses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

One-loop charge-breaking minima in the two-Higgs doublet model

We analyze the vacuum structure of the one-loop effective potential in the two Higgs doublet model. We find that electroweak-breaking vacuua can coexist with charge breaking ones, contradicting a theorem valid at tree-level. We perform a numerical analysis of the model and supply explicit parameter values for which charge-breaking vacuua can be the global minimum of the theory, and deeper than charge-preserving ones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗